薄板弯曲问题弹性理论分析及数值计算课程设计论文 毕业设计.doc

薄板弯曲问题弹性理论分析及数值计算 课程设计 指导教师： 学 院： 航空学院 姓 名： 学 号： 班 级： 薄板弯曲问题弹性理论分析及数值计算 一般三维体弹性系统求解微分方程体系总结 弹性力学中的基本假定 连续性，即假定整个物体的体积都被组成这个物体的介质所填满。 完全弹性，物体在引起形变的外力被除去后可完全恢复原形 均匀性，即假定物体是由同一材料组成的。 各向同性，物体的弹性在所有各个方向都相同。 和小变形假定，即假定位移和形变是微小的。 平衡微分方程 在一般空间问题中，包含15个未知函数，即6个应力分量、6个形变分量和3个位移分量，它们都是x,y,z坐标变量的函数。对于空间问题，在弹性体区域内部，考虑静力学、几何学和物理学三方面条件，分别建立平衡微分方程、几何方程和物理方程；并在给定约束面或面力的边界上，建立位移边界条件或应力边界条件。然后在边界条件下根据所建立的三套方程求解应力分量、形变分量和位移分量。 在物体内的任一点P，割取一个微小的平行六面体，如图1-1所示。根据平衡条件即可建立方程。 分别以连接六面体三对相对面中心的直线为矩轴，列出力矩的平衡方程，可证明切应力的互等性： (2)分别以为投影轴，列出投影的平衡方程，，，对方程进行约简和整理后，得到空间问题的平衡微分方程如下 物体内任一点的应力状态 现在，假定物体在任一点P的6个直角坐标面上的应力分量 为已知，试求经过P点的任一斜面上的应力。为此，在P点附近取一个平面ABC，平行于这一斜面，并与经过P点而平行于坐标面的三个平面形成一个微小的四面体PABC，如图1-2所示。当四面体PABC无限减小而趋于P点时，平面ABC上的应力就成为该斜面上的应力。 命平面ABC的外法线为，则其方向余弦为 三角形ABC上的全应力在坐标轴上的投影用代表.根据四面体的平衡条件进行推到，可以得出 (1-2) 设三角形ABC上的正应力为，则，将式1-2代入，并分别用代替，即得 (1-3) 设三角形ABC上的切应力为，则由于，得 (1-4） 由式1-3和1-4可见，6个应力分量完全决定了一点的应力状态。在特殊情况下，如果ABC是物体上受面力作用的边界面，则成为面力分量 ，于是由式1-2得空间问题的应力边界条件 (1-5) 应力状态有三种表示方式如下: (1)如图1-2,在图中表示 (2)应力状态矩阵 该矩阵为一对称阵。 应力向量 物体内任一点的应变状态 过空间一点P所有方向上的线应变和角应变的集合称为P点的应变状态，通过该点作三个相互垂直的线元。该三线元长度改变（线应变）和线元间夹角改变（角应变）的集合就完整地代表了P点的应变状态。三个线应变为，三个角应变为：. 应变状态的表示方式如下: 向量形式 (2)矩阵形式 几何方程和物理方程 空间问题的几何方程 （1-6） 几何方程的矩阵形式为（在V内），其中L为微分算子 空间问题的物理方程，在材料力学中根据胡克定律导出如下 ，，， ，， 根据关系，其中为体应变，为体积应力，与间的比例常数称为体积模量，可推得物理方程的另一种形式 物理方程的矩阵形式为或，其中D为弹性矩阵，C为柔度矩阵，两矩阵为互逆关系。 边界条件 根据物体内任一点的应力状态可得空间问题的应力边界条件，即式1-5 空间问题的位移边界条件为 (1-9) 按位移求解空间问题 按位移求解问题，是取位移分量为基本未知函数，并要通过消元法，导出弹性体区域内求解位移的基本微分方程和相应的边界条件。将几何方程代入物理方程，得出用位移分量表示应力分量的弹性方程，再将该弹性方程代入平衡微分方程得按位移求解时所需用的基本微分方程。 按应力求解空间问题 按应力求解空间问题，是取应力分量为基本未知函数。对空间问题来说就是，就是要从15个基本方程中消去位移分量和形变分量，得出只包含6个应力分量方程，进行求解。 板弯问题基本概念及微分方程 有关概念 (1)在弹性力学里，两个平行面和垂直于这两个平行面的柱面或棱柱面所围成的物体，称为平板，或简称板，如下图所示。这两个平行面称为板面，而这个柱面或棱柱面称为侧面或板边。两个板边之间的厚度称为板的厚度，而平分厚度的平面称为板的中间平面，或简称为中面。如果板的厚度远小于中面的最小尺寸b，这个板就称为薄板，否则就称为厚板。 (2)当薄板受有一般载荷时，总可以把每个载荷分解为两个分载荷，一个是平行于中面的所谓纵向载荷，另一个是垂直于中面的所谓横向载荷。对于纵向载荷，可以认为它们沿薄板厚度均匀分布，因而它