计算固体力学第四章.ppt

记为 一般情况下，采用高斯积分求解。 内力矩计算公式则为： 三、平面壳体单元 壳也是应用十分广泛的一种基本构件。它的中面是一个曲面，在变形时，弯曲变形和中面变形是耦合在一起的，因此与板有很大的不同。当壳体的厚度跟其他尺寸相比是小量时，称为薄壳。 分析时，可以用薄板单元组成的折板系统去代替原来的薄壳，由平面应力状态和板弯曲应力状态加以简单组合而得到薄壳的应力状态。 通常对于任意形状的壳体，采用三角形单元比较方便。如果在壳体上容易找到同一平面上的四个点，也可以采用四边形单元。 （1） 平面壳体单元 壳体单元的特点：包括两部分，其一为弯曲变形，其二为中面变形，因而壳体中内力包括弯曲内力和中面内力。 （2） 曲面厚壳单元 位移函数： 8. 修正的余能和杂交应力元 * * * * * 系数ai可以用节点坐标和未知的节点位移表示. 上式中形函数由下式表示 单元刚度矩阵: 考虑均匀压力作用下, 确定本题的四面体单元的刚度矩阵,设E=30×106psi, v=0.3 采用刚刚介绍的方法来计算本例题的单元刚度矩阵. 建议采用计算机计算. 习题1: 7. 其他等参单元 板壳单元 板和壳是指厚度比其他尺寸要小得多的平面或曲面构件，在工程中应用广泛。 但是采用前面给出的三维单元来近似，需要采用很细密的网格，这样将导致有限元模型的自由度太大，需要大量的计算时间和前后处理时间。因此需要开发适合于板壳结构的专用单元。 仿照梁理论建立梁单元的思路，采用板理论建立板单元。这里讨论两种板理论： 一是薄板理论，也称为 Kirchhoff 板理论，它忽略了板的横向剪切变形； 另一种是Mindlin 板理论，它考虑了板的横向剪切变形的影响，适合半于板的厚跨比较大的情形，这个也称为 Reissner板理论或中厚板理论。 对于线弹性材料，板内的应力应变关系为： 板的弹性矩阵 D跟平面应力状态是一致的。 在板理论中，一般假设板的中面是一中性面，也就是在没有面内力时，中面上的三个应变εx=εy=γxy=0。另一个基本假设即为所谓的直法线假定：变形前垂直于中面的法线变形后仍然保持直线，但是不一定仍然垂直于变形后的中面。 如图b，这条直线有绕 y和x轴的转角分别为ψx和ψy。则距离中面距离为z的任意点的位移和应变分别为： 弹性板弯曲理论： 这里w是板的横向挠度，假设它沿板厚度方向不变，即w,z=0。上式是Mindlin 板理论的基本假定。 如果假定变形后的法线仍然是变形后中面的法线，即w,x=ψx和 w,y=ψy，则前面方程中的两个横向剪切应变γyz和γxz 为零，这就是退化的Kirchhoff板理论（适用于薄板）。 由于忽略横向剪切变形，即w,x=ψx和 w,y=ψy，因此板内所有的力学量都能用挠度w 表示，比如，应变可以表示为挠度的曲率，即： 根据Mindlin板理论的假设，中面法线在变形后不再垂直于中面，因此必须要用三个位移分量 w, ψx和 ψy来描述板内的变形，即 另外为了修正横向剪应力沿板厚均匀分布导致的误差，引入了所谓剪切修正因子 k 来修正剪力， 由于考虑了横向剪切变形，Mindlin板理论包含了三个互相耦合的中面挠度和中面法线的两个转角位移。这样，Mindlin板理论看上去比Kirchhoff板理论更为复杂，但是在构造有限单元时，有时候比Kirchhoff理论更容易作到位移协调。当然，这是以克服当板的厚跨比很小时出现的“剪切自锁”现象为代价的。 应用有限元法求解板弯曲问题时，用离散的板单元代替原来连续的结构。每个节点有三个广义位移分量，即挠度w,绕x 轴的转角 θx 和绕y 轴的转角 θy。 矩形12自由度板单元 4节点矩形板单元有12个节点位移分量，采用局部坐标，可得 式中弹性矩阵D取平面应变的弹性矩阵，由于矩阵 Bi 是局部坐标的函数，因此在单元内内力矩不是常量。而且由不同单元计算得到的节点内力矩一般也不一致，这跟前面的平面或三维单元一样，即单元间的应力和应变是不连续的。可以通过绕节点平均法或者其他磨平方法来改善内力矩的精度。 9自由度三角形薄板弯曲单元 三个节点，每个节点三个自由度，故位移函数可以使用9个参数。 根据单元刚度矩阵的普遍计算公式，将要对面积坐标的多项式在三角形上进行积分，从而得到显式的刚度矩阵，但是很复杂，一般可以采用数值积分。 设单元受横向的分布载荷 p 的作用，于是等效节点力是 内力矩的计算公式在形式跟矩形板单元的类似，即 前面讨论的矩形单元和三角形单元都是Kirchhoff薄板理论的，它忽略了剪切变形的影响。由于Kirchhoff板理论要求挠度的导数连续，给构造协调单元带来不少的麻烦，为此，采用考虑剪切变形的Mindlin板理论来克服。这种方法比较简单，精度较好，并且能利用等参变换，得到任意四边形甚至曲边四边形单元，实用价值较高。 面积