Plattenbiegung

Prof. Dr. K. Weinberg

Universitt Siegen, Lehrstuhl fr Festkrpermechanik

Beispiele zur Plattendurchbiegung

Vergleich Plattentheorie-Balkentheorie
Gegeben sei eine rechteckige Platte der Dicke h, die in y-Richtung eine sehr groe Ausdehnung
hat ( y ), bei x = 0 fest eingespannt und bei x = a durch eine Randquerkraft
qRand und ein Randmoment mRand belastet ist. Man berechne fr diesen Plattenstreifen die
Verformungen und Schnittgren.
Da die Durchbiegung fr dieses Beispiel eine reine Funktion der Variable x ist, reduziert sich
die Differentialgleichung der Plattendurchbiegung w (x) zu:
w,xxxx = 0
Fr die Schnittgren erhalten wir
mx = Kw,xx ,

my = Kw,xx ,

mxy = 0,
11111
00000
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111

Die obige, gewhnliche Differentialgleichung lsst sich durch elementares Integrieren lsen. Es folgt:
2

w (x) = C0 + C1 x + C2 x + C3 x

11111
00000
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111
00000
11111

qx = Kw,xxx ,

qy = 0
qRand

mRand
a

x
y

Die unbekannten Konstanten C0 bis C3 erhalten wir aus den folgenden Randbedingungen:
w (0) = 0 = C0
w,x (0) = 0 = C1
mx (a) = mRand = K (2C2 + 6C3 a)
qx (a) = qRand = K6C3
und C2 =
Mit C3 = qRand
6K

mRand
2K

qRand a
2K

erhalten wir insgesamt

qRand 3
1
(mRand + qRand a) x2
x
2K
6K
mx (x) = mRand qRand (a x) my (x) = mx (x)
qx (x) = qRand
w (x) =

Die grte Durchbiegung tritt bei x = a auf

wmax = w (a) =

mRand a2 qRand a3
+
2K
3K
1

23. Juni 2010

Plattenbiegung
Prof. Dr. K. Weinberg
Universitt Siegen, Lehrstuhl fr Festkrpermechanik

Bei dem Vergleich zur Balkentheorie geht man bei der Berechnung der zugehrigen Biegelinie
ebenfalls von einer Dgl. vierter Ordnung aus. Fr einen Balken der Lnge a mit der Biegesteifigkeit EI, der einseitig eingespannt und am anderen Ende durch eine Kraft F und ein Moment
M belastet wird, ergibt sich mittels des Superpositionsprinzips:
wmax =

M a 2 F a3
+
2EI
3EI

An dieser Stelle muss man also die Biegesteifigkeit des Balkens EI [N mm2 ] mit der Plattensteifigkeit K [N mm] vergleichen. Die Plattensteifigkeit K ist definiert als
Eh3
.
K=
12 (1 2 )
3

b
Ein Vergleich mit einem Rechteckbalken der Biegesteifigkeit Eh
fhrt zu dem Ergebnis, dass
12
wb
nEI
sich eine Platte der Breite nb um den Faktor wp = nbK = 1 2 weniger durchbiegt, als n
nebeneinanderliegende Balken der Breite b. Allgemein gilt zwischen der Biegeflche des Plattenstreifens und der Biegelinie des Rechteckbalkens folgender Zusammenhang:

wPlattenstreifen (x) = 1 2 wBalken (x)

Recheckplatte unter einer konstanten rechteckigen Flchenlast
u

Gegeben sei eine Rechteckplatte der Dicke h, die bei x = 0, x = a, y = 0 und

y = b gelenkig gelagert ist. Diese wird mit
einer rechteckflchigen (2c 2d), konstanten Flchenlast p mit Mittelpunkt x = u
und y = v belastet.

Allgemeine Lsung
Im Wesentlichen nutzt man aus, dass eine Fourierreihe durch Anwendung des Bipotentialoperators ebenfalls in eine Fourierreihe berfhrt wird. Man bedenke,
dass beispielsweise die vierte Ableitung
der Funktion sin (x) erneut die Funktion
sin (x) liefert. Allgemein gilt z.B.:

111111111
000000000
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111
000000000
111111111

a0
(ak cos (kx) + bk sin (ky)) +
2
k=1

d
d
b

p (x, y)

(k cos (kx) + k sin (ky))

k=1

Aus diesem Grund drcken wir die Belastung p (x, y) als eine ungerade Funktion von x und y
in einer Fourierreihe mit den Perioden 2a und 2b aus:
p (x, y) =

X
X

m=1 n=1

pmn sin

mx
ny
sin
a
b

XX
m

pmn sin (m x) sin (n y) ,

wobei m = m
und n = n
.
a
b
Anmerkung: Die Periode der Glieder der Fourierreihe wird so gewhlt, dass die homogenen
2

23. Juni 2010

Plattenbiegung
Prof. Dr. K. Weinberg
Universitt Siegen, Lehrstuhl fr Festkrpermechanik

Randbedingungen des Problems erfllt werden.

Ein hnlicher Ansatz kann nun fr die Durchbiegung w (x, y) gewhlt werden:

X
X

mx
ny
wmn sin
w (x, y) =
sin
a
b
m=1 n=1

XX
m

wmn sin (m x) sin (n y)

Die Koeffizienten wmn erhlt man nun aus der Differentialgleichung fr die Plattendurchbiegung
p (x, y)
K


XX
X X pmn
2
2 2
wmn m + n sin (m x) sin (n y) =
sin (m x) sin (n y)
m n
m n K
2 w (x, y) =

Aufgrund der linearen Unabhngigkeit der trigonometrischen Polynome folgt mittels Koeffizientenvergleich:
pmn
wmn =
2 + 2 )2
K (m
n
Somit lautet die exakte Lsung dieses Problems:
XX
pmn
w (x, y) =
sin (m x) sin (n y)
2
2 2
m n K (m + n )

Lastentwicklung
Wir entwickeln zunchst die Rechtecklast der Aufgabenstellung als ber die Breite a verteilte
Last in y-Richtung und schreiben allgemein

ny
a0 X
ny
an cos
+ bn sin
p (y) =
+
2
b
b
n=1

Bei der Entwicklung einer ungeraden Funktion gilt:

1Z b
ny
dy
p (y) sin
a0 = 0, an = 0 und bn =
b b
b
Wir erhalten unter Anwendung der folgenden trigonometrischen Eigenschaften
cos (x) = cos (x) sin (x) = sin (x)
cos (a + b) = cos (a) cos (b) sin (a) sin (b)
fr den Fourierkoeffizient bn :
1Z b
ny
bn =
p (y) sin
b b
b

"

p
ny
b
=
cos
b
n
b
"

"

Z v+d
p
ny
dy =

sin
b
b
vd

b
ny
+
cos
n
b
vd

#v+d

#v+d

"

vd

dy +

v+d

vd

ny
sin
b

dy

n(v + d)
n(v d)
n(v + d)
n(v d)
cos
cos
cos
+ cos
b
b
b
b
"
!
!#
n(v d)
n(v + d)
2 cos
2 cos
b
b
"
!
!
!






nd
nd
nd
nv
nv
nv
cos
sin
cos
2 cos
+ 2 sin
2 cos
b
b
b
b
b
b
!#
!




nd
nd
4p
nv
nv
sin
sin
+2 sin
=
sin
b
b
n
b
b

p
=
n
p
=
n
p
=
n

!#

23. Juni 2010

Plattenbiegung
Prof. Dr. K. Weinberg
Universitt Siegen, Lehrstuhl fr Festkrpermechanik

Damit ist
!

nd
4p X
1
nv
ny
sin
sin
sin
p (y) =
n=1 n
b
b
b

p (y)

11111111111111111
00000000000000000
00000000000000000
11111111111111111
00000000000000000
11111111111111111
00000000000000000
11111111111111111

1111111111111111
0000000000000000
0000000000000000
1111111111111111
0000000000000000
1111111111111111
0000000000000000
1111111111111111
0000000000000000
1111111111111111

y
v

Eine Entwicklung in x-Richtung erbigt dann ganz analog

mc
mx
1
mu
4p (y) X
sin
sin
sin
m=1 m
a
a
a
X

16p X
1
= 2
sin (m u) sin (m c) sin (n v) sin (n d) sin (m x) sin (n y)
m=1 n=1 mn

p (x, y) =

Es gilt somit:

16p
sin (m u) sin (m c) sin (n v) sin (n d)
2 mn
Fr den Sonderfall der vollstndig belasteten Rechteckplatte erhlt man mit
a
b
c = u = und d = v =
2
2

0
fr m oder n gerade
16p
n
m
pmn = 2
sin2
= 16p
sin2
mn
2
2
fr m und n ungerade
2 mn
pmn =

Die Durchbiegung lautet dann

w (x, y) =

16p
K 2

sin (m x) sin (n y)
2 + 2 )2
mn (m
m,n ungerade
n
X

Fr die Schnittmomente erhalten wir daraus

2
X
16p
(m
+ n2 ) sin (m x) sin (n y)
mx (x, y) = K (w,xx + w,yy ) = 2
2 + 2 )2
m,n ungerade
mn (m
n
my (x, y) = K (w,yy + w,xx ) =

2
X
(n2 + m
) sin (m x) sin (n y)
16p
2
2 + 2 )2
m,n ungerade
mn (m
n

mxy (x, y) = (1 ) Kw,xy = (1 )

X
m n cos (m x) cos (n y)
16p
2
2 + 2 )2
m,n ungerade
mn (m
n

23. Juni 2010

Plattenbiegung
Prof. Dr. K. Weinberg
Universitt Siegen, Lehrstuhl fr Festkrpermechanik

Die Doppelreihe fr die Durchbiegung weit ein gutes Konvergenzverhalten auf. Zur Berechnung
der grten Durchbiegung einer quadratischen Platte (a = b) bei drei gltigen Ziffern bentigt
man trotzdem etwa 8 Reihenglieder.
wmax = w

sin

16pa
a a
=
,
2 2
K 6

m,n ungerade

m
2

sin

n
2
2 2

mn (m2 + n )

0, 00406

pa4
K

Fr den Sonderfall einer Punkt-/Einzellast F = p 2c 2d im Punkt (x = u, y = v) erhlt man

F
ber einen Grenzwertprozess:
die bentigten Fourierkoeffizienten pmn mit p = 4cd
16p
sin (m u) sin (m c) sin (n v) sin (n d)
2 mn
4F
sin (m u) sin (m c) sin (n v) sin (n d)
=
mcnd




mc
nd




4F
nv sin a sin b
mu
=
sin
sin
mc
nd
ab
a
b
a
b

pmn =

Mit lim

sin

= 1 folgt fr pmn mit c 0 und d 0:

pmn

Fr u =

a
2

und v =

b
2

nv
mu
4F
sin
sin
=
ab
a
b

erhlt man

pmn = 4F
ab

und

w (x, y) =

fr m oder n gerade
fr m und n ungerade

X
sin (m x) sin (n y)
4F
2 + 2 )2
Kab m,n ungerade
(m
n

Im Vergleich zur vollbelasteten Rechteckscheibe ergibt sich mit F = pab, dass die Durchbie2
gung im Falle einer betragsgleichen Einzellast etwa um den Faktor 4 grer ist als bei einer
gleichmigen Lastverteilung.
Durchbiegung von Kreisplatten
Als ein Modellbeispiel betrachten wir eine am Rand (r = R) eingespannte Kreisplatte konstanter Dicke, die mit einer Flchenlast p (r, ) belastet ist. Die Belastung sei gegeben durch

p (r, ) = p0 1 +

r
cos ()
R

Gesucht ist nun die entsprechende Biegeflche w (r, ). Ausgangspunkt unserer Berechnungen
ist die Plattengleichung in ebenen Polarkoordinaten:
p (r, )
K
!
!


2
2
2
2
r
()
1 () 1 ()
1 w 1 w
p0
w
1 + cos ()
+ 2
+
+ 2 2 +
=
r2
r 2
r r
r2
r
r r
K
R
2 w (r, ) =

23. Juni 2010

Plattenbiegung
Prof. Dr. K. Weinberg
Universitt Siegen, Lehrstuhl fr Festkrpermechanik

Die durch die feste Einspannung erzeugten Randbedingungen lauten:

1111111111111111111111111111111111111
0000000000000000000000000000000000000
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
r
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111

0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
p (r, )
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111
0000000000000000000000000000000000000
1111111111111111111111111111111111111

w (R) = 0 und w,r (R) = 0

Es ist zweckmig, die Gesamtlsung aus den Teilproblemen
p0
(rotationssymmetrisches Teilproblem)
K
p0 r
2 w 1 =
cos () (restliches Problem)
KR

2 w 0 =

zusammenzusetzen.

111
000
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111

111
000
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111
000
111

2R
Die homogene Lsung des rotationssymmetrischen Teilproblems ergibt sich anhand der Lsungstabelle fr den Bipotentialoperator in ebenen Polarkoordinaten zu:
w0h = A0 + B0 r2 + C0 ln (r) + D0 r2 ln (r)
Bei der vorliegenden Vollplatte mssen die Konstanten C0 und D0 Null sein, da sonst unendlich
groe Schnittgren fr r = 0 auftreten wrden.
Berechnung der Partikulrlsung fr das rotationssymmetrische Teilproblem
Bei einem rotationssymmetrischen Problem hngen die Belastung p und die Durchbiegung w
nur noch von dem Radialparameter r ab. Dadurch vereinfacht sich der Laplaceoperator zu
d2 () 1 d ()
d ()
1 d
() =
r
. ()
+
=
2
dr
r dr
r dr
dr
!

Die Plattengleichung lautet somit

d2 () 1 d ()
+
dr2
r dr

d2 w 1 dw
+
dr2
r dr

p0
p (r)
=
K
K

Die Partikulrlsung kann nun unmittelbar durch Integration ermittelt werden. Zu diesem
Zwecke spalten wir die Plattengleichung in zwei Differentialgleichungen
M = p0

w =

M
K

auf. Mit Hilfe der Gleichung () erhlt man durch zweifache unbestimmte Integration nach r
und entsprechendes Umformen:
!

1 d
dM
d
dM
r
= p0
r
= p0 r
r dr
dr
dr
dr
R
R
()dr
p0 r 2
p0 r 2
dM
p0 r ()dr
dM
M=
=

=
r
dr
2
dr
2
4
6

23. Juni 2010

Plattenbiegung
Prof. Dr. K. Weinberg
Universitt Siegen, Lehrstuhl fr Festkrpermechanik

Analog ergibt sich mit diesem M die Partikulrlsung zu:

1 d
dw0p
p0 r 2
dw0p
p0 r 3
d
w0p =
r
=
r
=

r dr
dr
4K
dr
dr
4K
R
4
4
3
()dr
p0 r
p0 r
dw0p
p0 r
dw0p
w0p =
=

=
r
dr
16K
dr
16K
64K
!

()dr

Somit lautet die Lsung des rotationssymmetrischen Teilproblems

w0 (r) = w0p + w0h =

p0 r 4
+ A0 + B0 r2
64K

Die homogene Lsung des restlichen Teilproblems ergibt sich anhand der Lsungstabelle fr
den Bipotentialoperator in ebenen Polarkoordinaten zu:

w1h = A1 r + B1 r3 + C1 r1 + D1 r ln (r) cos ()

Auch hier gilt, wegen der Beschrnktheit der Lsung C1 = D1 = 0.
Berechnung der Partikulrlsung fr das restliche Teilproblem
Fr die Partikulrlsung w1p des restlichen Teilproblems gilt:
2 w1p =

p0 r
cos ()
KR

Da der Term cos () unter dem Bipotentialoperator unverndert bleibt muss der Ansatz fr die
Partikulrlsung von der Art
w1p = C1p rn cos () , n N
sein. Einsetzen dieses Ansatzes in den Laplaceoperator liefert:

w1p = C1p n (n 1) rn2 cos () rn2 cos () + nrn2 cos ()

= C1p n2 1 rn2 cos ()

h

C1p n2 1 rn2 cos () = C1p n2 1 ((n 2) (n 3) 1 + (n 2)) rn4 cos ()

= C1p n2 1

n2 4n + 3 rn4 cos ()

! p0 r
cos ()
n2 4n + 3 rn4 cos () =
KR
p0
=
192KR

2 w1p = C1p n2 1
n=5

C1p

Somit folgt fr die Gesamtlsung des restlichen Teilproblems

w1 = w1p + w1h =

p0 r 5
+ A1 r + B1 r3 cos ()
192KR
!

Die Konstanten A0 , B0 , A1 , B1 erhlt man nun aus den Randbedingungen des Problems:
w (R) = 0
p0 R 4
w0 (R) =
+ A0 + B0 R2 = 0
64K

w1 (R) =
7

p0 R 5
+ A1 R + B1 R3 cos () = 0
192KR
!

23. Juni 2010

Plattenbiegung
Prof. Dr. K. Weinberg
Universitt Siegen, Lehrstuhl fr Festkrpermechanik

w,r (R) = 0
5p0 R4
+ A1 + 3B1 R2 cos () = 0
w1,r (R) =
192KR
5p0 R3
A1 =
3B1 R2
192K
!

p0 R 3
w0,r (R) =
+ 2B0 R = 0
16K
p0 R 2
B0 =
32K
Insgesamt folgt somit
A0 =

p0 R 4
,
64K

B0 =

p0 R 2
,
32K

A1 =

p0 R 3
,
192K

B1 =

p0 R
96K

Zusammenfassend erhlt man fr die Teillsungen

p0 R 4
r
w0 =
1
64K
R
"

2 #2

p0 R 4 r
r
w1 =
1
192K R
R
"

2 # 2

cos ()

und endgltig
r
p0 R 4
1
w (r, ) =
192K
R
"

2 # 2

r
3 + cos () .
R

Fr die Absenkung des Mittelpunktes der Kreisplatte gilt

wmax = w (r = 0) =

p0 R 4
64K

23. Juni 2010