第十章 薄板弯曲

合集下载

薄板弯曲问题

薄板弯曲问题

物理方程
应变
位移函数
薄板在弯曲变形后,薄板的法线没有伸缩;
w z 0 z
w wx, y
位移函数
薄板的法线,在薄板弯扭以后,保持为薄 板弹性曲面的法线;
xz yz 0
w u 0 x z
w v 0 y z
位移函数
u w z x
利用12个结点位移条件,由广义坐标法可 建立形函数,显然十分麻烦。
位移函数
w( x, y ) 1 2 x 12 xy
3
f x, y
w f x, y x y y
w f x, y y x x
D Dz
薄板弯曲问题的有限元法
结点 位移函数 位移 用插值方法求 内部各点位移
几何方程
结点力
平衡方程
应力
物理方程
应变
内力与应力的关系
薄板内力微元体如图所示。

h/2
- h/2
yx zdxdz
h/2 - h/2
y
h/2
- h/2
x zdydz

h/2
- h/2
x xy zdydz
该转角的确定包含了单元全部结点位移参数,由于非公共 边上结点位移的协调关系不能保证,因此一般
综上所述,本节构造的位移场不能完全满足收敛的协调性 准则,具体为挠度及切向转角跨单元协调,法向转角跨单 元不协调,因此该单元不是完全协调元。
弹性薄板矩形(R12)单元
4) 非完全协调元的收敛性
4 i 1
w N i d i N d
已知支座位移问题时
薄板弯曲问题的有限元法

弹性力学:平板弯曲问题 (2) 薄板弯曲经典解法

弹性力学:平板弯曲问题 (2) 薄板弯曲经典解法

16q0
6
Dmn
m2 a2
n2 b2
2
(m 1,3,5, ; n 1,3,5, )
代入式(10.22),即得挠度的表达式 (受均布载荷)
m x n y
w 16q0
sin sin
a
b
D 6 m1,3,5,n1,3,5,
mn
m2 a2
n2 b2
2
(10.24)
由此可以用公式(10.11)求得内力的表达式。
y
2
w
t2 4
z 2
(10.5)
其中,D称为板的抗弯刚度,其表达式为
D Et3
12(1 2 )
(10.6)
最后,次要应力分量σZ,可根据z方向的平衡方程求得。
z xz yz
z
x y
将式(10.5)代入上式得
x
z
6D t3
4
w
t2 4
z 2
积分上式得
z
6D t3
4
w
t2 4
在边界上
w n 0
D 4 w q
将式(10.18)代入式(10.8)得
D
24 m a4
16 m a2b2
24m
b4
q
解得m并代入式(10.18)得
w
q
x a
2 2
y2 b2
2 1
8D
3 a4
2 a2b2
3 b4
这就是夹支边椭圆薄板在均布载荷作用下的挠度 表达式。
有了挠度表达式,就可以求的内力。
y2 b2
2
1
(10.18)
o
a
y 图10.6 椭圆板

薄板弯曲

薄板弯曲
e
k e e {B}T D{B}dxdy
S
(3)节点荷载 当单元上作用有分布载荷 p( x, y ) 单元等效结点力
Q e N pdxdy
e S
T
见书5.17
z, w n (wn,xn,yn) 2b 2a k (wk,xk,yk) y() 3m ) (wm,xm,ym) x() l (wl,xl, yl)
2 2w z 2 2 x x x 2w 2 y z 2 z w 2 y y xy 2 2w 2 z 2 xy xy
二、几何关系
法线转角和挠度的关系
w w z x y x u z y w v z z w x y w w x w x y
M x x zdz
h 2 h 2
M y
M xy M yx xy zdz
Mx M My M xy
1 Mx h 3 3 h 1 h 2 M y h zσdz D p 2 12 12(1 ) 2 M xy 0
2w 2 x 0 2w 1 0 2 y 1 0 2w 2 2 xy
h3 1 1 M Dp D 12
其中 称为板的弯曲刚度,D为弹性薄板的 弹性系数矩阵。
由于是平面的小变形,称

1 为曲率。
三、物理关系
x y z xy yz zx T

薄板的小挠度弯曲问题及经典解法

薄板的小挠度弯曲问题及经典解法

(z2

d2
4
)
y
2 w

(9-5)
(4)用w表示应力分量z
由平衡方程(7-1)式的第三式有(取 fz=0):
z zx yz
z
x y
(c)
若体力不为零,可把薄板单位面积内的体力及面力归入薄板上面的
面力,并用 q表示。
d
q ( f )z zd
FRB

2D(1
)
2w xy
B
(9-18)
集中剪力或集中反力的正负号决定于角点处的扭矩的正负号, 而不能另行规定。据此,A点和C点处的剪力以沿z轴的正方向为正, 而O点和B点处的剪力以沿 z轴的负向时为正。
如果点B是自由边AB和自由边BC的交点,而点B并没有任何支 柱对薄板施以此向集中力,则应有FRB=0 ,亦即:
z
w
w(x, y)即在垂直于中面的任一法线
上,薄板全厚度内各点的挠度相同。
2)由几何方程, zy

w v y z
0
, zx

u z
w x
0
,得
v w , u w z y z x (2) z 引起的形变可以不计。
(9-1)
由物理方程(7-12),有:
(3)应用时可查相关手册,若是双向配筋时,扭矩的影响 也可不考虑。
§9-4 边界条件 扭矩的等效剪力
矩形薄板,OC边简支;OA边固支;AB和BC边自由。
1. 固支边,OA边(x = 0)
(w) x ( w )
x
0 x0
0 0
(9-13)
2. 简支边,OC边 (y = 0)
x y

薄板弯曲问题

薄板弯曲问题

第五章薄板弯曲问题机场学院2011/11/21CAUCCAUC两个平行面和垂直于这两个平行面的柱面或棱柱面所围成的物体,称为平板,简称为板。

bhyxzCAUCCAUC垂直于板面——平板弯曲问题byxzCAUCCAUC1、小变形假设:虽然板很薄,但它的挠度远小于板的厚度。

byxz)(0==z u 0)(0==z v 因为:2、板中面各点都没有平行于中面的位移,只发生弯曲变形。

x u x ∂∂=εy v y ∂∂=εyu x v xy ∂∂+∂∂=γ所以:0)(0==z x ε0)(0==z y ε0)(0==z y x γCAUC CAUC3、沿板的厚度方向挤压变形忽略不计。

byxz=∂∂=zw z ε所以:),(y x w w =在薄板中面的任一根法线上,薄板全厚度内的所有各点都具有相同的挠度。

CAUCCAUC保持在挠曲面法线上。

byxz应力分量:zx τzy τzσ远小于其余三个应力分量,其引起的形变忽略不计。

0=zx γ0=zx γ0=∂∂+∂∂xw z u 0=∂∂+∂∂yw z v 即:等价于:这样=∂∂=z w z ε0=zx γ0=zx γ中面法线不伸缩,仍为变形后曲面的法线CAUC CAUCxyxy x y y y x x EEE τµγµσσεµσσε)1(2)(1)(1+=−=−=薄板弯曲与平面应力问题有相同的物理方程。

CAUCCAUC1、几何方程byxz0=∂∂+∂∂x w z u 0=∂∂+∂∂y w z v xw z u ∂∂−=∂∂y w z v ∂∂−=∂∂),(2y x f z yw v +∂∂−=),(1y x f z xwu +∂∂−=0)(0==z u 0)(0==z v 因为:),(),(21==y x f y x fCAUCCAUCzxu ∂−=zyv ∂−=zxwx u x 22∂∂−=∂∂=εzyw y v y 22∂∂−=∂∂=εz yx w y u x v xy∂∂∂−=∂∂+∂∂=22γ221xw x ∂∂−=ρ221ywy ∂∂−=ρyx wxy ∂∂∂−=221ρ令:xx zρε=yy z ρε=xyxyz ργ=得:CAUCCAUCw y x y x xy y x ⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂−=⎪⎪⎪⎭⎪⎪⎪⎬⎪⎪⎪⎩⎪⎪⎪⎨=⎭⎬⎫⎩⎨⎧222221111ρρρρ{}w y x y x z xy y x ⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂−=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=222222γεεε写成列阵形式:应变列阵:CAUCCAUCxyxy x y y y x x EEE τµγµσσεµσσε)1(2)(1)(1+=−=−=xyxy x y y y x x EEE γµτµεεµσµεεµσ)1(2)(1)(122+=+−=+−={}w y x y x z xy y x ⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂−=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=222222γεεεyx w Ez x w y w Ez y wx w Ez xy y x ∂∂∂+−=∂∂+∂∂−−=∂∂+∂∂−−=222222222221)(1)(1µτµµσµµσCAUCCAUCyx w Ez xw y w Ez yx xyy x ∂∂∂+−=∂∂+∂∂−−=∂+∂−−=2222222221)(1)(1µτµµσµµσ其它几项应力:w yh z E w xh z E zy zx22222222)4()1(2)4()1(2∇∂∂−−=∇∂∂−−=µτµτw hz h z Eh z 4223)1()21()1(6∇+−−−=µσCAUCCAUC在薄板的上表面有:qh z z −==2)(σ得:q w Eh =∇−423)1(12µ令:)1(1223µ−=Eh D qw D =∇42、微分方程CAUCCAUC xyab边界条件:0)(,0)(0)(,0)(0)(,0)(0)(,0)(220220220220=∂∂==∂∂==∂∂==∂∂=========b y b y y y a x a x x x xww x ww x ww x w w qw D =∇4微分方程:四边简支矩形薄板的重三角级数解答——纳维叶解法CAUCCAUC设重三角级数解为:b yn a x m A w m n mn ππsinsin 11∑∑∞=∞==代入微分方程:qb yn a x m A b n am D m n mn =+∑∑∞=∞=πππsin sin )(1122224b yn a x m C q m n mn ππsinsin 11∑∑∞=∞==将),(y x q q =也展成重三角级数:CAUCCAUC222226)(16bn a m Dmn q A mn +=π(m=1,3,5, m=1,3,5, ………… n=1,3,5, n=1,3,5, …………)∑∑∞=∞=+=...5,3,1,...5,3,12222260)(sin sin 16m n bn a m mn b yn a x m D q w πππ得挠度的表达式:CAUC CAUC荷代替q ,得:dxdyP q =b n a m bn a m abD P dxdy b n a m dxdy P b n a m abD A mn ηπξππηπξππsin sin )(4sin sin )(4222224222224+=+=CAUC CAUC集中载荷作用下的简支矩形板挠度表达式:b y n a x m bn a m b m a m abD P w m n ππηπξππsin sin )(sin sin 411222224∑∑∞=∞=+=M x yxzM y{}[]zDxyyx⎭⎬⎫⎩⎨⎧=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=ρτσσσ1zdzMhhxx∫−=22σ1、弯曲应力zdzMhhyy∫−=22σzdzMhhxyxy∫−=22τCAUC CAUCCAUC CAUC{}zdzM M M M h xy y x ∫−=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=22}{σ完成积分:⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧=ρρ1][1][12}{3D D hM ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡−−=21000101)1(12][23µµµµEh DCAUCCAUC2b2ayxzlmn kw θ yθ x(1)节点位移单元任一节点有三个位移分量:{}⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂−∂∂=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧=i i i yi xi i i x w y w w w )()(θθδ{}{}Tyk xk k ynxn n ymxm m yl xl li w w w w θθθθθθθθδ={}{}T T kT nT mTli δδδδδ=CAUCCAUC31231131029283726524321xya y x a y a xy a y x a xa y a xy a x a y a x a a w +++++++++++=写成矩阵形式:{}a xy yx yxyyx xy xy xy xw ]1[33322322=或:{}a y x M w )],([=CAUCCAUC{}a xy yx yxy yx xy xy xy xw ]1[33322322={}a xy xyxy xy x yw x ]332020100[2322=∂∂=θ{}a y y x y xy xy x xw y ]302302010[3222−=∂∂−=θCAUC CAUC⎪⎪⎪⎪⎭⎪⎪⎬⎪⎪⎪⎪⎩⎪⎪⎨⎥⎥⎥⎥⎥⎥⎦⎢⎢⎢⎢⎢⎢⎣=⎪⎪⎪⎪⎭⎪⎪⎬⎪⎪⎪⎪⎩⎪⎪⎨654310000110000001a a a a y x y x y x y x v u v u n nn n m m m m n n m m {}[]{}a A e=δ[]{}[][]{}a A A A e 11−−=δ{}[]{}eA a δ1−=[]{}[][]{}{}eey x N A y x M a y x M w δδ)],([),(),(1===−A[][]k nm lN N N N y x N =),(形函数CAUCCAUC⎥⎥⎦⎤⎟⎠⎞⎜⎝⎛−⎟⎠⎞⎜⎝⎛+⎟⎠⎞⎜⎝⎛+−⎟⎠⎞⎜⎝⎛−⎟⎠⎞⎜⎝⎛+⎟⎠⎞⎜⎝⎛+⎢⎣⎡⎟⎟⎠⎞⎜⎜⎝⎛−−++⎟⎠⎞⎜⎝⎛+⎟⎠⎞⎜⎝⎛+=111,111,21181][2222222222222222a x x b y y a x x x b y y b y y a x x y b y a x b y y a x x b y y a x x N i i i i i i i i ii i i i (i =l ,m ,n ,k )单元刚度阵:ee xy y x B N y x y x w y x y x }]{[}]{[2211112222222222δδρρρρ=⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂−=⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂−=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧=⎭⎬⎫⎩⎨⎧CAUCCAUC][][k n m l B B B B B =单元内力:eB D M }]{][[}{δ=[][][][]dxdy B D B k Ts ee∫=单元刚度阵:[]{}{}Q K =δ整体方程:。

薄板的弯曲破坏分析与预测

薄板的弯曲破坏分析与预测

薄板的弯曲破坏分析与预测薄板是一种常见的结构材料,广泛应用于建筑、航空航天、汽车等领域。

然而,在使用过程中,薄板可能会遭受弯曲破坏,导致结构的失效。

因此,对薄板的弯曲破坏进行分析与预测,对于设计和使用薄板结构具有重要意义。

首先,我们来探讨薄板弯曲破坏的原因。

薄板在受到外力作用时,会发生弯曲变形。

当外力超过薄板的承载能力时,薄板可能会发生破坏。

薄板的弯曲破坏主要包括弯曲变形和局部破坏两个方面。

在弯曲变形方面,薄板在受到外力作用时,会发生曲率变化,即薄板的中部会凸起或凹陷。

这种变形会导致薄板的强度和刚度下降,进而影响结构的稳定性和安全性。

因此,对于薄板的弯曲变形进行分析与预测,可以帮助我们更好地评估薄板结构的承载能力。

而在局部破坏方面,薄板在受到外力作用时,可能会出现局部的破坏现象,如薄板的边缘开裂、孔洞扩展等。

这种局部破坏会导致薄板的强度降低,进而引发整体结构的失效。

因此,对于薄板的局部破坏进行分析与预测,可以帮助我们更好地评估薄板结构的寿命和可靠性。

接下来,我们来探讨薄板弯曲破坏的分析与预测方法。

薄板的弯曲破坏是一个复杂的力学问题,需要运用弹性力学、塑性力学、断裂力学等多个学科的知识进行分析。

其中,有限元分析是一种常用的分析方法,可以通过建立薄板的数值模型,计算薄板的应力和变形,进而评估薄板的弯曲破坏情况。

此外,实验方法也是分析薄板弯曲破坏的重要手段。

通过设计合适的试验装置和加载方式,可以模拟薄板在实际使用中的受力情况,从而观察薄板的弯曲变形和破坏过程。

通过实验数据的分析,可以得到薄板的弯曲破坏特征和破坏机制,为薄板结构的设计和使用提供参考依据。

此外,还可以借助计算机模拟和人工智能等新技术手段,对薄板的弯曲破坏进行预测。

通过建立合适的模型和算法,可以预测薄板在不同工况下的弯曲破坏情况,从而指导薄板结构的设计和使用。

这种方法不仅可以提高分析和预测的准确性,还可以节省时间和成本,提高工作效率。

综上所述,薄板的弯曲破坏分析与预测对于设计和使用薄板结构具有重要意义。

薄板弯曲问题有限元法

薄板弯曲问题有限元法

T
wl xl yl
Fzl M zl M yl T
j
xj
yj
wj
7
第8页/共24页
薄板弯曲时,只有w(x,y)是薄板变形的未知基本函数,而其它量,如u,v 等都是w(x,y)的函数,故薄板矩形单元的位移函数的选择实际就是w(x,y) 的选取。注意单元有12个自由度,则
w(x, y) 1 2x 3 y 4x2 5xy 6 y2
1 2
(w,
Ljj
w, Ljm
),
a5
1 2
(w,Lii
w, Lim
),
6
1 2
(w,Lii
w, Lij
w, Lji
w,Ljj
),
7
wj
wm
1 2 (w,Ljj
w, Ljm
)
8
wi
wm
1 2
(w,Lii
w, Lim
)
w,Lij 表示w对Li的 偏导数在j点的值。
9
wi
wj
1 2
(w,Lii
角形和矩形。为了使相邻单元间同时可传递力和力矩,节点当作刚性节点
,即节点处同时有节点力和节点力矩作用。每个节点有三个自由度,即一
个扰度和分别绕x,y轴的转角。 1.设位移函数
l
xl
yl wl
m
xm ym wm
节点位移分量和节点力分量
i
xi
yi
wi
q e wi xi yi F e Fzi M xi M yi
w(x, y) c1 c2 x c3x2 c4 x3
四个系数刚好通过i,j两个端点的扰度值和绕y轴的两个转角值唯一确定 ;同时,相邻单元在此边界上也能通过i,j的值唯一确定,故连续。

薄板弯曲有限元法

薄板弯曲有限元法

3 边界条件(续):
(3) 自由边
a
x
b
在x = a 处自由:
M x (a, y) 0,
2w 2w 2 0 x 2 y x a
y
M xy (a, y) 0, Qx (a, y) 0
M xy 0 Qy x x a
在y = b 处自由:
3w 3w (2 ) 0 2 x3 xy x a
M y (a, y) 0, M xy ( x, b) 0, Qy ( x, b) 0
3w 3w (2 ) 2 0 y3 x y y b
o
x
h b
y
a
z
u( x, y,0) v( x, y,0) 0, w( x, y, z) w( x, y,0)
基本位移函数为中面挠度: w w( x, y)
o
x
一、薄板弯曲理论基础
y
h
b
a
2、基本方程 位移:
w( x, y, z ) w( x, y ) w u ( x, y , z ) z x w v ( x, y , z ) z y
T
实外力: q( x, y), V ( s), M n ( s)
实应力: { } { x , y , xy }T
虚应变:
实内力:{M } {M x , M y , M xy }
虚功方程:
* * * { } { M }d x d y w q d x d y w V ds A A S2 T
M1 y 1 y
单元插值函数:
w a1 a2 x a3 y a4 x 2 a5 xy a6 y 2 a7 x 3 a8 x 2 y a9 xy 2 a10 y 3 a11 x 3 y a12 xy 3

弹性力学_第10章_薄板弯曲

弹性力学_第10章_薄板弯曲
t 2 t 2
dy
dx
xy在板厚上的总和为零,只能分别合 成为弯矩 M x和扭矩 M xy ;而 xz只能合
显然,在垂直于x 轴的横截面 上,每单位宽度之值如下:
成横向剪力Qx 。
M M Q
x


t 2

t 2 t 2
x zdz
t 2
xy

t 2


xz
xy
zdz
x



t 2

《弹塑性力学》课件
内容提要
§10-1 有关概念和基本假设 基本方程 薄板横截面上的内力 薄板的边界条件 薄板弯曲的直角坐标求解 圆形薄板的轴对称弯曲
2 48
薄 板 弯 曲 问 题
§10-2 §10-3 §10-4 §10-5 §10-6
2012/5/31
§10-1 有关概念和基本假设
薄板区别于厚板。通常情况下,板的厚度t与板面的 最小尺寸b的比值满足如下条件。
2
将(1)式对z 积分,得:
Et z 6 1 2
3


1 z z 4 1 2 t t
2
设在薄板顶面上每单位面积作用的载荷q(包括横向面 力和横向体力),板上面的边界条件为:
z z t
q
2
将 z的表达式代入该边界条件,得薄板挠曲微分方程:
zx z t 0, zy z t 0
2 2
将上列二式对z 进行积分,得:
zx
E 2 1 2


zy
E 2 1 2


2 t2 2 z 4 x 2 t2 2 z 4 y

机械工程用有限元法学习笔记(四)

机械工程用有限元法学习笔记(四)

薄板弯曲问题的有限元法一、 薄板弯曲问题的基本方程什么是薄板?薄板就是指厚度t 远小于其长度、宽度的板。

1. 三个基本假设(克希霍夫假设): (1) 法线假设,εz =0,γyz =γzx =0 (2) 正应力假设,σz <<σx ,σy ,τxy (3) 小挠度假设,w<t/4根据假设,可以得到位移分量()()()()()(),,,,,,,,,,, x y z u x y z z x x y z v x y z z y x y z x y ωωωω∂⎧=-⎪∂⎪∂⎪=-⎨∂⎪⎪=⎪⎩式4-1图 1 薄板弯曲后某点B 的位移2. 应变分量{}222222x y z x z y x y ωεωεεεω⎧⎫∂-⎪⎪∂⎪⎪⎧⎫⎪⎪∂⎪⎪⎪⎪==-⎨⎬⎨⎬∂⎪⎪⎪⎪⎩⎭⎪⎪∂-⎪⎪∂∂⎪⎪⎩⎭式4-23. 曲率{}222222x y z x y x y ωχωχχχω⎧⎫∂-⎪⎪∂⎪⎪⎧⎫⎪⎪∂⎪⎪⎪⎪==-⎨⎬⎨⎬∂⎪⎪⎪⎪⎩⎭⎪⎪∂-⎪⎪∂∂⎪⎪⎩⎭式4-3 22=x x ωχ∂-∂——薄板弹性曲面在x 方向的曲率22=y yωχ∂-∂——薄板弹性曲面在y 方向的曲率2=z x yωχ∂-∂∂——薄板弹性曲面在x 方向和y 方向的扭率4. 应力分量与应变分量间的关系:{}[]{}2222222222221 11D Ez xy Ez x y Ez x y σεωωμμωωμμωμ=⎧⎫⎛⎫∂∂-+⎪⎪ ⎪-∂∂⎝⎭⎪⎪⎪⎪⎛⎫∂∂⎪⎪=-+⎨⎬ ⎪-∂∂⎝⎭⎪⎪⎪⎪∂⎪⎪--∂∂⎪⎪⎩⎭式4-4 5. 线力矩{}()2222222101012110022x y z x M Et M M y M x y ωμωμμμω⎧⎫∂-⎪⎪⎡⎤∂⎪⎪⎢⎥⎧⎫⎪⎪⎢⎥∂⎪⎪⎪⎪==-⎨⎬⎨⎬⎢⎥∂-⎪⎪⎪⎪⎢⎥-⎩⎭⎪⎪⎢⎥∂⎣⎦-⎪⎪∂∂⎪⎪⎩⎭式4-5a广义应力与广义应变之间的关系式{}[]{}D M χ= 式4-5b式中:[D]—薄板弯曲问题的弹性矩阵6. 薄板弯曲问题的基本方程(双调和方程)()32222222121Et p xx y y ωωωμ⎛⎫∂∂∂++= ⎪∂∂∂∂-⎝⎭ 式4-6()32121Et μ-——薄板弯曲刚度 二、 矩形薄板单元分析 1、矩形薄板单元图 2 矩形薄板单元2、位移函数22123456322333789101112 a a x a y a x a xy a y a x a x y a xy a y a x y a xy ω=+++++++++++ 式4-73、形状函数[]{}k i i xi xi yi yi j j xj xj yj yj k kxk xk yk y l l xl xl yl yl N N N N N N N N N N N N N q ωωθθωθθωθθωθθ=+++++++++++= 式4-8式中:i,j,k,l ——节点号N i ,N xi ,N yi ,……,N yl ——形状函数()()()()()()()()()()2211128N 111 ,,,8111 8y i i i i i xi i i iyi i i i b N i i j h l N a x a b ξξηηξξηηξηηξξηηηξξξηηξξη⎧⎫++++--⎪⎪⎧⎫⎪⎪⎪⎪⎪⎪=-++-=⎨⎬⎨⎬⎪⎪⎪⎪⎩⎭⎪⎪++-⎪⎪⎩⎭==, 式4-94、单元刚阵[][][][]S K TB D B dxdy =⎰ 式4-10式中:[]22222222222222222222 2222yi yl i xiyi yl ixi yi yl i xi N N N N x x x x N N NN B y y y y N N N N x yx yx yx y ⎡⎤∂∂∂∂⎢⎥∂∂∂∂⎢⎥⎢⎥∂∂∂∂=⎢⎥∂∂∂∂⎢⎥⎢⎥∂∂∂∂⎢⎥⎢⎥∂∂∂∂∂∂∂∂⎣⎦式4-11 5、节点力与节点位移的关系式{}[]{}F K q = 式4-12三、 三角形薄板单元分析1、三角形薄板单元当薄板具有斜交边界或曲线边界时,可采用三角形单元较好地反映边界形状。

弹性薄板的小挠度弯曲课件

弹性薄板的小挠度弯曲课件
践指导。
06
参考文献
参考文献
总结词:详细描述了弹性力学的基本 原理,包括应力和应变的关系,以及 弹性薄板在受到外力作用时的弯曲变 形规律。
详细描述:在弹性力学中,薄板的小 挠度弯曲是指薄板在受到外力作用时 发生的弯曲变形,其弯曲变形程度较 小,可以忽略不计薄板的剪切变形和 转动惯性。这种变形情况下,薄板的 弯曲变形可以通过挠度(即变形量) 来描述。在弹性力学中,应力和应变 之间的关系由胡克定律(Hooke's Law)描述,即应力与应变成正比, 比例系数为材料的弹性模量。
详细描述
圆形薄板在受到垂直于其平面的力时,会在力的方向上发生弯曲,形成弧形。与矩形薄板类似,这种弯曲程度较 小,也称为小挠度弯曲。在圆形薄板中,各个方向的弯曲程度基本相同,因此圆心位置的应力最大。
实例三:不规则形状薄板的弯曲
总结词
不规则形状薄板在受到垂直于其平面的力时,会发生小挠度弯曲。
详细描述
不规则形状薄板在受到垂直于其平面的力时,会在力的方向上发生弯曲,形成弧形。与矩形和圆形薄 板类似,这种弯曲程度较小,也称为小挠度弯曲。不规则形状薄板的弯曲情况较为复杂,需要考虑各 个方向的弯曲程度以及应力分布。
05
结论与展望
研究结论
结论一
弹性薄板在受到小挠度弯 曲时,其弯曲行为与材料 属性、几何尺寸等因素密 切相关。
结论二
通过理论分析和数值模拟, 我们得到了弹性薄板在小 挠度弯曲下的变形规律和 应力分布。
结论三
实验结果与理论预测和数 值模拟结果基本一致,验 证了理论的正确性和数值 方法的可靠性。小的单元,对每 个单元进行弯曲分析,通过求解每个 单元的平衡方程得到整体的挠度分布。
对于某些特定形状和载荷条件的薄板, 可以通过解析方法直接求解弯曲微分 方程,得到挠度分布的精确解。

薄板弯曲问题弹性理论分析及数值计算资料

薄板弯曲问题弹性理论分析及数值计算资料

薄板弯曲问题弹性理论分析及数值计算课程设计指导教师:孙秦学院:航空学院姓名:程云鹤学号: 2011300092班级: 01011105薄板弯曲问题弹性理论分析及数值计算一、一般三维体弹性系统求解微分方程体系总结1、弹性力学中的基本假定(1)连续性,即假定整个物体的体积都被组成这个物体的介质所填满。

(2)完全弹性,物体在引起形变的外力被除去后可完全恢复原形 (3)均匀性,即假定物体是由同一材料组成的。

(4)各向同性,物体的弹性在所有各个方向都相同。

(5)和小变形假定,即假定位移和形变是微小的。

2、平衡微分方程在一般空间问题中,包含15个未知函数,即6个应力分量、6个形变分量和3个位移分量,它们都是x,y,z 坐标变量的函数。

对于空间问题,在弹性体区域内部,考虑静力学、几何学和物理学三方面条件,分别建立平衡微分方程、几何方程和物理方程;并在给定约束面或面力的边界上,建立位移边界条件或应力边界条件。

然后在边界条件下根据所建立的三套方程求解应力分量、形变分量和位移分量。

在物体内的任一点P ,割取一个微小的平行六面体,如图1-1所示。

根据平衡条件即可建立方程。

(1)分别以连接六面体三对相对面中心的直线为矩轴,列出力矩的平衡方程0=∑M ,可证明切应力的互等性:yx xy xz zx zy yz ττττττ===,,(2)分别以轴轴、轴、z y x 为投影轴,列出投影的平衡方程0=∑x F ,0=∑y F ,0=∑z F ,对方程进行约简和整理后,得到空间问题的平衡微分方程如下⎪⎪⎪⎭⎪⎪⎪⎬⎫=+∂∂+∂∂+∂∂=+∂∂+∂∂+∂∂=+∂∂+∂∂+∂∂000z yzxz z y xyzy y x zx yx x f y x z f x z y f z y x ττσττσττσ (1-1)3、物体内任一点的应力状态现在,假定物体在任一点P 的6个直角坐标面上的应力分量 ,,z y x ,σσσyx xy xz zx zy yz ττττττ===,,为已知,试求经过P 点的任一斜面上的应力。

薄板弯曲问题-浙江大学共36页PPT资料

薄板弯曲问题-浙江大学共36页PPT资料

从薄板内取出一个平行六面体,
它的三边长度分别为d x , d y和板的厚度
图(9-2)
在x为常量的横截面上,作用着 x , y , xz 在该截面的每单位宽度上,应力分量 x
a
对中面合成为弯矩 M x
2 a
z
x
dz
2
将式(9-4)中的第一式代入并对z进行积分,得
M x 1 E 2 2 x w 2 2 y w 2 a 2 a 2 z 2 d z 1 2 ( 1 E 32 ) 2 x w 2 2 y w 2
0 取 z
由几何方程的第三式得 w0wwx,y
z
结论:中面的任一根法线上的各点都有相同的横向位移,也就等于挠度
2)应力分量 xz , yz , z 远小于其余的3个应力分量
所引起的形变可以忽略不计
z 0,zx 0,yz 0
从而有 u w,v w z x z y
可见:中面的法线在薄板弯曲时保持不伸缩,并且成为弹性曲面的法线
x
12M 3
x
z,
y
12M 3
y
z,
xy
yx
12M 3
xy
z,
zx
6 FSx 3
2
4
z2
,
yz
6 FSy 3
2 4
z2
,
z
2
q
1 2
z
2
1
z

(9-11)
• 由内力表示的平衡微分方程
Qx Qy q0 x y
Mx x
M yyxQx
0
MxyMy x y
Qy
0
2M x2x22x M yxy 2M y2yq0
薄板的物

薄板弯曲问题

薄板弯曲问题

z
y
x
将x、y、xy的表达式代入得:
zx Ez ( 3w 3w ) Ez 2w z 1 2 x3 xy2 1 2 x
zy
z

1
Ez

2
(
3w y3

3w yx2
)

1
Ez

2
2w y
2020年3月10日星期二
2020年3月10日星期二
专题:薄板弯曲问题
12
最后得
z


Et3
6(1 2 )
(1 2

z )2 (1 t
z )4w t
边界条件: z t q 板上面 2
将z的表达式代入此边界条件,得:
Et 3
12(1
2
)
4w

q
or或写成: D4w q
其中:D

Et 3
2
二、薄板弯曲的基本假定 (1)薄板弯曲时,中面为曲面,称为弹性曲面或中曲面, 中面内各点在垂直方向的位移称为挠度.
(2)小挠度弯曲:挠度t,(本节讨论) 大挠度弯曲:挠度=t 薄膜问题 :挠度t
(3)薄板弯曲的基本假定:(Kirchhoff-Love假定)
a、假定应变分量z=0,xz=0,yz=0
19
六、边界条件
求薄板的小挠度弯曲问题,就是在满足板边的边界条件
下,由方程
D4w q
求出挠度w
下面以矩形板为例:
O
如图所示矩形板,OA边固定
a
x C
,OC边简支,AB、BC边是自由 b

OA边 w 0 x0
w 0 A x x0

【弹塑性力学】6 薄板弯曲

【弹塑性力学】6 薄板弯曲
2 Qx D w x
Qy D
2 w y
• 应力与内力的关系
12 z x 3 M x t
6 t2 yz 3 ( z 2 )Q y t 4
12 z y 3 M y t
xy
12 z 3 M xy t
6 t2 zx 3 ( z 2 )Qx t 4
Et 3 4 wq 2 12(1 ) Et 3 D= 12(1 v 2 )
D是板的弯曲刚度,板厚的三次方成正比,与弹模成正比,与 梁的弯曲刚度类似
D4 w q
• 薄板横截面上的内力
x
M x z x dz
M xy
t 2 t 2
t 2 t 2
z xy dz
3w 3w x 3 (2 ) xy 2 =0 x a
3w 3w y 3 (2 v) x 2y =0 y b
2w 2w y 2 x 2 =0 y b
x
z
My Mxy Myx Qy Qx Mx
y
• 内力由挠度表示
将应力的表达式代入积分得到
2w 2w M x D( 2 v 2 ) D( K x vK y ) x y 2w 2w M y D( 2 v 2 ) D( K y vK x ) y x 2w M xy M yx D1 xy
3w 3w V y D 3 (2 v) 2 x y y
O
2w RB 2 D(1 v)( ) xy
C
x
b
A y
a
B
(1)自由边
弯矩和合成剪力为零,因此,
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

取薄板的中面为xy面,如图10-1薄板的计算假定如下 (1)垂直于中面方向的线应变,即 ε z 可不计,取 ε z = 0, ∂ω = 0 从而得 ω = ω ( x, y ) 则由几何方程得
∂z
这就是说,横向位移w只是x,y的函数,不随z而变。因 此,在中面的任一根法线上各点都具有相同的横向位移 ,也就等于挠度。
2 (1 − μ
)⎣ 4
3⎝
8 ⎠⎦
Eδ 3 ⎛ 1 z ⎞ =− − ⎟ 2 ⎜ 6 (1 − μ ) ⎝ 2 δ ⎠
2
z⎞ ⎛ 1 + ⎟ ∇ 4ω ⎜ ⎝ δ⎠
(10-6)
现在来导出w的微分方程。由薄板的上板面的边界条件
(σ z ) z =−δ
= −q
(f)
2
其中q是薄板每单位面积内的横向载荷,包括横向面力 及横向体力,如式(d)所示。
§10-2 弹性曲面的微分方程
注意 τ xz = τ zx ,τ yz = τ zy ,将这两个应力分量的表达式(10-5 )代入式(c)得 ⎛δ2 ∂σ z E 2⎞ 4 = − z ⎟∇ ω 2 ⎜ ∂z 2 (1 − μ ) ⎝ 4 ⎠
积分
⎡δ 2 E z3 ⎤ 4 σz = z − ⎥ ∇ ω + F3 ( x, y ) 2 ⎢ 3⎦ 2 (1 − μ ) ⎣ 4
第十章 薄板弯曲
目 录
§10-1 有关概念及计算假定 §10-2 弹性曲面的微分方程 §10-3 薄板横截面上的内力 §10-4 边界条件 扭矩的等效剪力 §10-5 四边简支矩形薄板的重三角级数解 §10-6 矩形薄板的单三角级数解 §10-7 薄板弯曲的直角坐标求解 §10-8 圆形薄板的弯曲 §10-9 圆形薄板的轴对称弯曲 §10-10 变分法求薄板的位移
其中的待定函数 F1 ( x, y ) , F2 ( x, y ) ,可以根据薄板的上, 下板面的边界条件求出,即
(τ zx ) z =± δ
= 0,
2
(τ )
zy z =± δ 2
=0
§10-2 弹性曲面的微分方程
应用这两个边界条件求出 F1 ( x, y ) , F2 ( x, y ) 以后,即得 τ zx ,τ zy 的表达式
1 ⎡σ x − μσ y ⎤ ⎦ E⎣ 1 ε y = ⎡σ y − μσ x ⎤ ⎦ E⎣ 2 (1 + μ ) γ xy = τ xy E
εx =
(2)
§10-1 有关概念及计算假定
(3)薄板中面内的各点都没有平行于中面的位移即
( u ) z =0 = 0,
εx =
∂u ∂v ∂u ∂v , ε x = , γ xy = + ∂x ∂y ∂y ∂x
§10-1 有关概念及计算假定
对于薄板的弯曲问题,可以用来计算工程上的问题。对于 厚板,但还不便应用于工程实际问题。当薄板受有一般载 荷时,总可以把每个载荷分解为两个分荷载,一个是平行 于中面的所谓纵向荷载,另一个是垂直于中面的所谓横向 荷载。对于纵向荷载,可以认为它们沿薄板厚度均匀分布 ,因而它们所引起的应力,形变和位移,可以按平面应力 问题进行计算,横向载荷将使薄板弯曲,它们所引起的应 力,形变和位移,可以按薄板弯曲问题进行计算。
§10-1 有关概念及计算假定
两个平行面和垂直于这两个平行面的柱面或棱柱面 所围成的物体,称为平板,或简称为板,如图10-1。这 两个平行面称为板面,而这个柱面或棱柱面称为侧面或 板边。两个板面之间的距离 δ 称为板的厚度,而平分厚 度 δ 的平面称为板的中间平面,或简称为中面。ቤተ መጻሕፍቲ ባይዱ果板 的厚度 δ 远小于中面的最小尺寸b,这个板就称为薄板, 否则就称为厚板。
∂ 2ω ⎞ E ⎛ ∂ 2ω +μ 2 ⎟z σx = − 2 ⎜ 2 1 − μ ⎝ ∂x ∂y ⎠ ∂ 2ω ⎞ E ⎛ ∂ 2ω +μ 2 ⎟z σy = − 2 ⎜ 2 1 − μ ⎝ ∂y ∂x ⎠ E ∂ 2ω τ xy = − z 1 + μ ∂x∂y
(4)
由于w不随z而变可见这三个主应力分量都和z成正比。
∂w u=− z, ∂x
∂w v=− z ∂y
§10-2 弹性曲面的微分方程
(2)将主要应变分量ε x ,ε y , γ xy 用w表示。把上式的u,v 代入几何方程就得到
∂u ∂2w εx = =− 2 z ∂x ∂x ∂v ∂2w εy = = − 2 z ∂y ∂y ∂u ∂v ∂2w = + = −2 z ∂y ∂x ∂x∂y
§10-2 弹性曲面的微分方程
(4)将次要应力分量τ xz ,τ yz用w表示。由于次要应力分 量 τ xz ,τ yz 所引起的形变滤去不计,相应的物理方程也 已放弃。为了求出 τ xz ,τ yz ,可以应用平衡微分方程,并 由于不存在纵向荷载,体力分量 f x = 0, f y = 0,由此得到
( v ) z =0 = 0
(3)
(ε x ) z =0 = 0,
(ε )
y z =0
= 0,
(γ )
xy z = 0
=0
这就是说,中面的任一部分,虽然弯曲成为弹性曲面 的一部分,但它在 xy 面上的投影形状保持不变。
§10-2 弹性曲面的微分方程
薄板的小挠度弯曲问题是按位移求解的,只取挠度 ω = ω ( x, y ) 作为基本未知函数。将其他未知函数——纵向位移u,v 主要应变分量 ε x ,ε y , γ xy ,主要的应力分量σ x ,σ y ,τ xy , 次要应力分量 τ xz ,τ yz 及更次要的应力分量 σ z ,分别都 用挠度 w 来表示,并导出求解挠度的方程。
(a)
γ xy
(3)将主要应力分量用w表示。由薄板的物理方程求解 应力分量,得
E ⎡ε + με y ⎤ 2 ⎣ x ⎦ 1− μ E ⎡ε y + με x ⎤ σy = ⎦ 1− μ2 ⎣ E τ xy = γ xy 2 (1 + μ )
σx =
(b)
§10-2 弹性曲面的微分方程
再把式(a)的应变分量代入上式,就得出
(e)
其中待定函数 F3 ( x, y ) 可由薄板的下板面的边界条件来 确定,即
(σ z ) z = δ
=0
2
§10-2 弹性曲面的微分方程
将式(e)代入,求出 F3 ( x, y ) ,再代回到式(e),即 得 σ z 的表达式 ⎡δ 2 E 1 ⎛ 3 δ 3 ⎞⎤ 4 σz = ( z − δ ) − ⎜ z − ⎟⎥ ∇ ω 2 ⎢
§10-1 有关概念及计算假定
不计 τ xz ,τ yz 所引起的形变
∂u ∂w + = 0, ∂z ∂x ∂v ∂w + =0 ∂z ∂y
γ zx = 0, γ yz = 0
几何方程
∂u ∂w =− , ∂z ∂x
∂v ∂w =− ∂z ∂y
(1)
由于 ε z = 0, γ zx = 0, γ yz = 0 ,可见中面的法线在薄板弯 曲时保持不伸缩,并且成为弹性曲面的法线。
§10-1 有关概念及计算假定
在上述计算假定中虽然采用了 ε z = 0, γ zx = 0, γ yz = 0 ,但 在以后的考虑平衡条件时,仍然必须计入3个次要的应力分 量τ xz ,τ yz , σ z。因此,在薄板的小挠度弯曲理论中,放弃了 关于 ε z , γ zx , γ yz 的物理方程。 因为不计σ z 所引起的形变,所以薄板的物理方程为
⎛ 2 δ2 ⎞ ∂ 2 E τ zx = z − ⎟ ∇ω 2 ⎜ 4 ⎠ ∂x 2 (1 − μ ) ⎝ ⎛ 2 δ2 ⎞ ∂ 2 E z − ⎟ ∇ω τ zy = 2 ⎜ 4 ⎠ ∂y 2 (1 − μ ) ⎝
(10-5)
这两个切应力沿横向为抛物线分布。
§10-2 弹性曲面的微分方程
(5)最后,将更次要的应力分量σ z 用w表示。用平衡微 分方程,并取体力分量 f z = 0
σ x ,σ y ,τ xy
∂2 ∂2 ∇ = 2+ 2 ∂x ∂y
2

§10-2 弹性曲面的微分方程
将上两式对z积分,得
Ez 2 ∂ 2 τ zx = ∇ ω + F1 ( x, y ) 2 2 (1 − μ ) ∂x Ez 2 ∂ 2 τ zy = ∇ ω + F2 ( x, y ) 2 2 (1 − μ ) ∂y
§10-1 有关概念及计算假定
如果薄板的弯曲刚度很小,以致挠度远大于厚度,则薄 板成为薄膜。薄板问题属于空间问题。为了建立薄板的 小挠度弯曲理论,除了引用弹性力学的5个基本假设外, 还补充提出了3个计算假定,用来简化空间问题的基本方 程(这些计算假定,已被大量的实验证实是合理的)。
§10-1 有关概念及计算假定
称为薄板的弹性曲面微分方程,或挠曲微分方程。
§10-3 薄板横截面上的内力
薄板横截面上的内力,称为薄板内力,是指薄板横截面 的每单位宽度上,由应力合成的主矢量和主矩。 为了求出薄板横截面上的内力,从薄板内取出一个 平行六面体,它的三边的长度分别为dx, dy 和板的厚 度δ
§10-3 薄板横截面上的内力
∂τ zx ∂σ x ∂τ yx =− − ∂z ∂x ∂y ∂τ zy ∂σ y ∂τ xy =− − ∂z ∂y ∂x
∂τ zx Ez ⎛ ∂ 3 w ∂ 3 w ⎞ Ez ∂ 2 = ∇ω ⎜ 3 + ⎟= 2 2 2 ∂z 1 − μ ⎝ ∂x ∂x∂y ⎠ 1 − μ ∂x ∂τ zy Ez ⎛ ∂ 3 w ∂ 3 w ⎞ Ez ∂ 2 = + = ∇ω ⎜ ⎟ ∂z 1 − μ 2 ⎝ ∂y 3 ∂y∂x 2 ⎠ 1 − μ 2 ∂y
§10-2 弹性曲面的微分方程
(1)将纵向位移u,v用挠度w表示。把此式对z积分,并注 意w是x,y的函数,即得
相关文档
最新文档