有限元薄板弯曲问题分析.doc.docxVIP

《有限元》讲义 《有限元》讲义 PAGE PAGE # 第4章 弹性薄板弯曲问题的有限元法 薄板弯曲问题在理论上和应用上都具有重要意义 ，并有专门著作加以论述（如 杨耀乾《平板理论》）。象其它弹性力学问题一样，用微分方程、差分法等经典方 法所能求解的薄板问题很有限，一般只能解决等厚、小孔口、支承情况较简单的单 跨板。故工程设计中以往多采用简化、近似、图表等方法来解决板的设计问题。 在板的分析中，常取板的中面为 xoy平面（如图）。平板结构按其厚度 t与短边a TOC \o 1-5 \h \z 的比值大小而分为 ： 厚板（Thick plate ）和 属平面,具有梁的受力特征，即常说的弯曲问W 属平面 ,具有梁的受力特征，即常说的弯曲问 W 1时称为薄板 平板上所承受的荷载通常有两种 ： 面内拉压荷载。由面内拉压刚度承担 应力问题。 垂直于板的法向荷载 ，弯扭变形为主 题。平板在垂直于板面的荷载作用下产生挠度 当最大挠度 w远小于t时，称为小挠度问题 （or刚性板stiffness plate） 当最大挠度 w与 t相差不大时，称为大挠度问题 （or柔性板）（flexure plate） （工程定义w （工程定义 ~ 5为柔性板; w L —5 为绝对柔性板。 ） 4.1 基本理论 一、基本假定 1、 略去垂直于中面的法向应力。 （z 0 ）,即以中面上沿Z方向的挠度W代表板的挠度） 2、 变形前垂直中面的任意直线，变形后仍保持为垂直中面的直线。（一法向假定 zx 0，zy 0） 3、 板弯曲时，中面不产生应力。 （一中面中性层假定） 上述假定常称为薄板小挠度问题假定（ or柯克霍夫假定）。符合上述假定的平板即 为刚性板。 二、基本方法 以上述假定为基础，板分析中常用挠度 W作为基本未知量，下面介绍以 W为基本未知 量所导出的有关方程。 1、几何方程（应变一挠度关系） ①弹性曲面沿 x, y 方向的倾角 从中面取出一微小矩形 ABCD,如图所示，设其边长 为dx, dy ，变形后弯曲成曲面 ABCD 设A点挠度W,则沿x方向倾角（绕y轴） W y匚（B 点绕度 — dx)x 沿y方向倾角（绕x轴） W x (D 点绕度 y —dy)y ②沿x, y 方向位移 作平行于xoz平面，设中面上点A到A1的距离为 Z,变形后,A点有挠度 W,同时发生弯曲, 曲面沿x方向的倾角为 w ,根据法线假定 x ,则A点沿 7S x方向的位移 同理取yoz平面得: （负号为方向与 x相反） (4-1-1) ③Z平面的应变分量和曲、扭率 基本假定，由于 z zx zy 故板内任意点的应变与平面问题相同 xy 将U .V代入 xy 2 W z r x 2 W Z 亍(4-1-2) 2 W 2z— x y 2 2 2 w w w 此为Z平面的应变一挠度度几何方程。上式中的2 , —r , 此为Z平面的应变一挠度度几何方程。上式中的 x y x y (4-1-3)X,Y方向的曲、扭率，记为： (4-1-3) x y xy 所以， 2、物理方程（应力一挠度关系） 由于忽略b z对变形的影响，因此Z平面的应力一应变关系具有与平面问 题相同的形式： 将（4-1-2）代入得Ez2w2 将（4-1-2）代入得 Ez 2 w 2 w 1 2 2 2 x y x Ez 2 2 w w y 1 2 2 2 y x xy Ez 2 w 1 x y xy 2 1 xy 或简写为: Do z Do x (4-1-4) 式中弹性矩阵D 式中弹性矩阵 Do 0 0 1 2 从板内取微元体 tdxdy ,由其上正应力x , y和剪应力 xy 从板内取微元体 tdxdy ,由其上正应力 x , y和剪应力 xy 可在截面上合成合力矩： Mx( yOz面上由 x产生的绕丫轴弯矩) M y( x0z面上由 y产生的绕x轴弯矩) 扭矩：Mxy (由剪应力产生，如图) 假定Mx,M y,M xy分别表示单位宽度上的内力矩。如是，内力矩阵: MxM yM xyt2tz Mx M y M xy t 2tz 2 dz t 2t z2 Do 2 t3 dz 12 Do 2：；2 简写成 比较(4-1-4) t3F - Do 12 和(4-1-5)可得用内力矩表示的平板应力 (4-1-5) 由此可见，平板上、下表面处的应力最大: 以上是薄板弯曲问题中的基本公式 ，从中可见其挠度V是弯曲问题中的基本未知函数。 且由于忽略了 z方向的变化,因此它只是x，y的函数：w=w(x, y )。若w已知，则位移，内 力、应力均可按上述相应公式求出。在经典解析法中， W(x, y)常设为三角级数形式。例 如，四边简支矩形板的 W(x, y)设为：(纳维尔解) w x, ym x n y Amn s