下载文档原格式
第五章 薄板的弯曲薄板的概念：厚度t<<Min(B,L)()L B Min t 81~51<中厚板 ()L B Min t 81~51> 厚板()()L B Min t L B Min 81~511001~801<< 薄板()L B Min t 1001~801< 薄膜作用在其上的载荷分解为平行于板面和垂直于板面，当仅有平行于板面的力时，就是我们前面讲到的平面应力问题。
现在我们要解决的就是当有垂直于板面的载荷时（板受弯曲作用时），应该如何计算。
两者都有时，又应该如何考虑。
§5.1 薄板弯曲的基本方程一，基本概念1，中面：变形前平分板厚的平面。
2，挠度：中面上各点在垂直于中面上的位移w 。
3小挠度：通常w/t<1/5。
二，基本假定1，变形前垂直于中面上的直线，变形后仍为直线，且仍垂直于弯曲的中面。
该假定类似与材料力学中梁的平面假定。
它确保与中面平行的的各面之间不存在剪应变。
0==zy zx γγ 2，变形前后，板的厚度不变，即0=z ε。
板内各点的挠度值仅为x 、y 的函数，而与z 轴无关。
()y x w w ,=。
3，薄板中面内的各点没有平行于板面的位移()00==z u 、()00==z v ，只有z 方向的位移。
4，平行于中面的各层之间互不挤压。
0=z σ三，基本方程利用空间的三大方程和以上4个假定，我们可以推求出适用薄板的基本方程。
1，几何方程由假定○1，0=∂∂+∂∂=x w z u zx γ，0=∂∂+∂∂=ywz v zy γ，就有： x w z u ∂∂-=∂∂，ywz v ∂∂-=∂∂，积分可得： ()y x f xwzu ,1+∂∂-= ()y x f ywzv ,2+∂∂-=再由假定○3，()00==z u 、()00==z v ，就是中面上各点没有板面的位移，代入上式，可得()()0,,21==y x f y x f 所以x w zu ∂∂-=，ywz v ∂∂-=。
薄板弯曲问题的有限元求解1.问题描述如图所示，已知悬臂矩形薄板，其几何尺寸为20m×10m×1m,左边固定，右上角节点上作用有向下垂直于板中面的集中载100N。
材料的弹性模量为Ex=300GPa,泊松比μ=0.3，求薄板的位移、应力及固定端反力。
2.分析步骤（1）进入Ansys（设定工作目录和工作文件）；（2）设置计算类型为Structural；（3）选择单元类型shell63，选择与厚度有关，在Real constants中定义厚度参数为1；（4）定义材料参数弹性模量为EX:3e11;泊松比PRXY：0.3；（5）建立几何模型生成节点和单元。
此题结构简单，受力也简单，因此可用4个单元来分析。
首先创造节点，节点的坐标是：1（0，0，1）2（0，5，1）3(0,10,1) 4(10,0，1) 5(10,5,1) 6(10,10,1) 7(20,0,1) 8(20,5,1) 9(20,10,1)，操作如下：GUI：Preprocessor>Modeling>Create>Nodes>In Active CSGUI：Preprocessor>Modeling>Create>Elements>AutoNumbered>ThruNodes，逆时针方向依次连接这几个点形成4个4节点四边形单元（6）施加载荷与约束加载与施加边界条件板的左边完全被固定，其自由度为0；右边第9节点施加了一个垂直方向的集中力（7）求解 （8）查看结果 1）变形结果可得最大变形为51011.0-⨯m 2）查看节点位移3）查看等效应力可以看出3节点受最大应力1248Pa，节点7所受应力最小。
4）查看节点力及力矩可以看出节点1、2、3既受到Z轴的集中力又受到X、Y的弯矩。
节点9只受外载作用。
3.如果将例题中的受力作如下图的改变，则此时单元的计算应为薄壳问题。
按照前面的计算方法可得出节点的线位移、角位移及力和力矩。
第五章 弹性薄板小挠度弯曲问题的变分原理平分板厚度的平面称为板的中面，一般地，当板的厚度t 不大于板中面最小尺寸的5/1时的板称为薄板，薄板的中面是一个平面。
薄板在垂直于中面的载荷作用下发生弯曲时，中面变形所形成的曲面称为弹性曲面或挠度面，中面内各点在未变形中面垂直方向的位移称为板的挠度。
薄板弯曲的精确理论应是满足弹性力学的全部基本方程，但这在数学上将会遇到很大的困难。
1850年，G .R.Kirchhoff 除采用弹性力学的基本假设外，还提出了一些补充的假设，从而建立起了薄板小挠度弯曲的近似理论。
这些假设是：第一，变形前垂直于板中面的直线，在板变形后仍为直线，并垂直于变形后的中面，而且不经受伸缩；第二，与中面平行的各面上的正应力z σ与应力x σ,y σ和xy τ相比属于小量；第三，在横向载荷作用下板发生弯曲时，板的中面并不伸长，这也就是说，薄板中面内各点都没有平行于中面的位移分量。
用变分法可以导出薄板弯曲问题的平衡微分方程和边界条件。
当板的形状和边界条件较复杂时，直接求解偏微分方程时比较困难的，以变分法为基础的各种近似解是求解这类问题的一个重要途径。
本章讨论了用于薄板小挠度弯曲问题的一些基础变分原理，这包括虚功原理、最小位能原理、最小余能原理、两类自变量广义变分原理并推广到三类自变量广义变分原理。
§5.1 基本方程与边界条件回顾取坐标平面oxy 与中面重合，z 轴垂直于中面，x ,y 和z 轴构成一个右手直角笛卡儿坐标系。
变形后的板内各点沿x ,y 和z 轴方向的位移分别用u ,v 和w 表示。
由Kirchhoff 假设，可以得到x w zz y x u ∂∂-=),,(，yw z z y x v ∂∂-=),,(，),(),,(y x w z y x w = （5-1） 并利用弹性力学中位移与应变之间的关系式，可以得到薄板中任意点的应变分量为22x w z x ∂∂-=ε，22ywz y ∂∂-=ε，y x w z xy ∂∂∂-=γ22 （5-2）其余3个应变分量z ε,xz γ和yz γ根据假设都等于零，即0=εz ，0=γxz ，0=γyz （5-3）由薄板的平衡关系，可以确定板的横向分布载荷),(y x q 与剪力x Q ,y Q 以及弯矩x M ,y M 和扭矩xy M （x M ,y M ,xy M 统称为内力矩）与x Q ,y Q 之间的关系式。
§4-2 空间等参数单元的数学分析在进行空间等参数单元的力学分析时，需要用到（1） 各个形函数对整体坐标的导数（求应变）；（2） 局部坐标系中微分体的体积及微分面的面积（载荷移置、刚度矩阵）； （3） 局部坐标面的法向余弦（载荷移置）。
现在来分别导出这些参数的表达式。
一、形函数对整体坐标的导数● 由复合函数的求导法则，有xN x N x N x N i i i i ∂∂∙∂∂+∂∂∙∂∂+∂∂∙∂∂=∂∂ζζηηξξ (x,y,z)不过，由上节(4-17)和(4-18)式知，等参数单元的形函数N i 及x,y,z 均只是局部坐标ξ，η，ζ的显函数，所以，利用上式无法求出形函数N i 对整体坐标x,y,z 的导数。
● 我们如将Ni 理解为局部坐标ξ，η，ζ的复合函数，则有ξξξξ∂∂∙∂∂+∂∂∙∂∂+∂∂∙∂∂=∂∂yy N y y N x x N N i i i i （ξ，η，ζ） 等等，所以有⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂∂∂∂∂=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂∂∂∂∂⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂∂∂∂∂z N y N x N J z N y N x N z y x z y x z yx N N N i i i i i i i i i ][ζζζηηηξξξζηξ (4-19) 从而有⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂∂∂∂∂=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∂∂∂∂∂∂-ζηξi i i i i i N N N J z N y N x N 1][ (4-20) 这里⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂∂=ζζζηηηξξξz y x z y x z yxJ ][ (4-21)称为雅可比矩阵。
将(4-18)式代入上式。
于是得⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂∂∂∂∂∂∂=⎢⎢⎢⎢⎢⎢⎢⎣⎡∂∂∂∂∂∂=∑∑∑======ni n i i i ni ni ii ni n i ii N N N N N N x N x Nx N J 212121111111][ζζηηξξζηξ● 能够保证式(4-20)● 0°〈θ〈180°，一般应尽量控制在● 在单元内的任一点P ，沿局部坐标ξ、ζ的方向分别作微分矢量a 、b 、c 对应局部坐标的d ξ，d η，d ζ)，间形成一平行六面体(图4-9)。
第五章薄板弯曲问题有限元法
第一节薄板弯曲问题的有关概念
一、基本概念
1．薄板的定义：薄板是由上下两个平行的表面所构成的片状结构，其间距称为板厚。
同时，定义等分板厚的面为中面，当中面为平面时，称为平板，当中面为曲面时则称为壳体。
2．挠度; 板结构在承受横向载荷（弯矩、扭矩和横向剪力）作用下，发生弯扭而使薄板中面上各个点沿垂直中面方向发生的横向变形称为挠度，记为w。
3．薄板的两类问题：
（1）平面应力板问题，载荷作用于板面内—（薄膜单元）；在拉、压力和面内切力作用下，板内将产生薄膜内力，从而使板产生面内变形。
（2）薄板弯曲问题：其特点为：
a) 几何尺寸：板的厚度远较长与宽的几何尺寸为小(一般厚度与板面最小尺寸之比小于1/5-1/10)；（否则称为厚板）
b) 载荷条件：结构仅承受垂直于板中面的横向载荷作用。
c) 小挠度条件；即挠度与板厚之比值较小，一般为w/t ≤1/5。
研究薄板弯曲问题时，通常以未变形的板的中面为xoy平面，厚度方向为z轴方向，3．板的一般问题：
一般情况下，板既可承受横向载荷作用，也可同时承受平行于板中面的膜载荷作用。
(1) 薄板：在小挠度情况下，当两种载荷同时作用时，可认为两种变形互不影响，因此膜载荷的作用可按平面应力问题进行处理，而横向载荷的作用则按薄板弯曲问题来分析，两种问题引起的薄膜内力和弯曲内力的叠加便是一般载荷综合作用的结果。
（2）厚板：当1<w／t<5时为大挠度板,w／t≥5时为特大挠度板。
在大挠度情况下，薄板面内变形和弯扭变形之间将相互影响，即横向载荷也可能产生膜内力和面内变形，而膜载荷也可能产生弯曲内力和弯曲变形。
这时描述薄板变形的数学方程是非线性的，应采用更复杂的理论分析方法。
二．薄板弯曲问题求解的假设：（克希霍夫假设）
1．法线假设
垂直板中面的法线在板变形后仍垂直于弯曲的挠曲面，且法线线段没有伸缩，板的厚度无变化。
这样，垂直于中面的正应变便可忽略，即εz=0
根据几何方程，可得
因此挠度只是x,y的函数，表示为w=w(x,y)，也即薄板中面上法线的各点都有相同位移。
2．正应力假设
在平行于中面的截面上，应力分量σz、τzx及τyz远小于其他三个应力分量，可忽略不计。
3．小挠度假设
板中面只发生弯曲变形而没有面内变形，即中面内各点没有平行于中面的位移，表示为：
在这些假设前提下，薄板的位移、应变和应力都可用挠度w表示。
由于w仅为x、y的函数，因此实际薄板的复杂三维问题就可简化为二维问题来进行分析和计算。
其中位移d 与挠度w的关系为
应变ε与w的关系为：
应力σ与w的关系为：
[D]称为薄板弯曲问题的弹性矩阵，该矩阵和平面应力问题的弹性矩阵完全相同。
三、工程中的薄板弯曲问题
车辆工程中的车体地板、高速车辆的顶板及墙板，发动机缸体、齿轮箱箱体，建筑结构的楼板、桥梁桥面等都属于薄板弯曲结构，这类结构的几何模型是由薄板中面构成的平面模型，可用相应的弯曲板单元划分网格。
现在，一些专业有限元软件在弯曲板单元和平面应力单元的基础上发展了6自由度的板单元，这种单元可用于离散受一般载荷（包括面内与垂直板面载荷）的作用。
第二节 弹性薄板弯曲的能量泛函和微分方程式 一、位移向量
根据克希霍夫假设中的小挠度假设，薄板中任意一点没有平行于中面(即xoy 平面)的位移，即0;00
====z z v
u
,而在垂直于板面的载荷作用下也没有绕轴的转动，即θz=0。
在变
形前垂直于x-y 面的任一直线m-n ，在弯曲变形后在垂向产生了一个位移w, 同时也分别绕y 轴产生了一个转角
x
w
∂∂和绕x 轴产生了一个转角y w ∂∂(图中完整没有表示出来)。
这时直线m
—n 上任意一点既有垂向位移，又有沿x 轴和y 轴的转角位移：所以在进行有限元分析时，
薄板任一点有三个位移分量：挠度w 、法线绕x 轴的转动θx 和绕y 轴的转动θy ，即
对薄板进行网格划分时，通常采用矩形或三角形单元。
矩形板单元的节点位移向量为
二、广义应变分量和曲率
根据克希霍夫的法线假设，可知，所以此时薄板的应
变只有εx, εy,νxy 三个分量，表达式分别为：
令：
式中，称为薄板弯曲的广义应变，和分别称为板中面在x、y方向的变形曲率和在xy面上的扭率，并且有
三、应力—应变关系
由克希霍夫法线假设和正应力假设，有，由弹性力学物理方程得由上述三式即得到薄板弯曲的应力—应变关系为
四、广义应力
从上式中可以看到，当z=0时，即在板的中面上，所有的应力均为零，而且在板的截面上，应力沿板厚按线性关系分布(见图左)，这些应力在各截面上形成了相应的力偶，这些力偶就是各截面上所作用的弯矩和扭矩(图右)。
若用Mx，My和Mxy．，分别表示单位长度上的弯矩和扭矩，根据图可以得到:
同理有：
将前述薄板内应力表达式代入上面表达式并对z积分，可得：
因而，薄板弯曲的广义应力为：
式中，Mx、My分别是绕x轴和y轴的截面上单位长度的弯矩；Mxy是垂直于z轴截面上单位面积内的扭矩；[D]是薄板弯曲问题的弹性矩阵，其表达式为
五、能量泛函和微分方程式
弹性体的应变能为：
对于薄板弯曲问题，上式中的应变{ε}与应力{σ}应分别使用广义应变{}χ和广义应力{M}来代替。
由于求{M}时，已经对微分体的厚度积分，所以此时应只在板面范围内对dxdy 进行积分，即
设平板受有强度为Ps的均布载荷作用，此时外力位能为
将应变能与位能相加，得到板弯曲时的总位能为
将上式进行变分，并令δⅡ=0，从而可得出薄板弯曲问题的基本微分方程为
上式是双调和方程，在理论上可以通过该微分方程用解析法求解薄板问题。
但实际上由于边界条件及板形状的复杂性，对着以微分方程直接求解析解是十分困难的。
这类方程所描述的薄板弯曲问题采用有限元法求解则是较为方便和收敛的。
有限元法求解时，一般将薄板分割为矩形板单元和三角形板单元。
第三节薄板弯曲问题有限元法
一、结构离散
由于矩形板单元列式简单，精度也较好，所以求解薄板弯曲问题时，常用四节点矩形单元离散结构。
单元形状如图所示。
二、单元分析
1．位移函数
由第一节可知，薄板的位移、应变和应力都可用挠度w(x,y) 表示，因此将w(x,y)作为薄板弯曲的位移函数。
对于单元中的每个节点只考虑w,θx和θy三个自由度，则四节点矩形板单元共有12个自由度，所以位移函数应取成12项多项式，设为
把矩形板单元四节点的局部坐标值代人上式及其导数表达式（转角按右手螺旋法则定正负）
弯曲板矩形单元的节点自由度和相应的节点力列阵为：
运用边界条件联立求解可以得到α1•••α12的值，代入式(6—14)，经整理得
式中，[N]为形函数矩阵，
其中Ni，Nj, Nm, Np为形函数，其表达式分别为
应该说明，这样设定的。
在单元内及边界上是连续的，但θx和θy在边界上不一定连续，所以称采用这种位移函数的矩形板单元为非协调单元。
但实践证明，利用这种矩形板
单元计算薄板弯曲问题能很好地收敛于精确解。
2．单元刚度矩阵
将式(6—24)代人式(6—12)得弯曲板单元广义应变与位移之间的关系：
式中
则单元刚度矩阵有：
3．单元载荷向量
薄板弯曲问题中最常见的载荷形式是垂直于板中面的分布载荷Ps ，此时单元载荷向量可按下式计算
当Ps为常数，即薄板上作用有均布载荷时，有
三、总刚矩阵集成
按第二章介绍的方法可集成得到结构的总刚矩阵[K]。
(详见书本)
四、载荷移置
将各种载荷移置为等效节点载荷，形成节点载荷列阵}。
载荷移置的方法同书本第二章。
五、边界条件处理
总体刚度矩阵[K]和载荷列阵的建立过程中，还没有考虑结构的边界约束，因此
[K]是奇异阵。
为了求出惟一的节点位移解，必须对边界条件进行处理，排除结构的刚体
运动，消除[K]的奇异性。
六、求解线性方程组
对矩阵[K]{q}进行约束处理之后，就可求解结构的刚度方程，即线性方程
[K](q}={R}
式中，[K]、(q}分别为经约束处理后的总刚矩阵和节点载荷列阵。
求解单元任意点的位移和其他物理量的具体方法同前几章。
