应用有限元工程实例-3解答.pptVIP

* * * * * * 三角形壳单元刚度方程 * 是转换矩阵 * * 93个节点 * * * * * * * * * * 单元节点力： * * * * * * 节点4 * * 将式位移u,v两式合写在一起，以矩阵形式表示，则有 * * * * 剪切变形影响 * * 局部坐标系下的单元刚度矩阵 * * * * 单元节点力： * * 最大应力在板上下表面 * * 3. 单元刚度矩阵 单元刚度方程： 令单元刚度矩阵 单元刚度矩阵显式表达式： 式中， A是单元的横截面面积， I 为横截面惯性矩，E 为材料的弹性模量。 第1行、第1列对应于ui,….。 其中， 为剪切影响系数， A为梁的横截面面积，E和G是弹性摸量和剪切摸量， 为截面剪切因子。矩形和圆形截面，文献[3]建议 分别取5/6和10/9。 对于高为h 宽为b 的矩形截面，剪切影响系数 ，所以，当梁的高度 相对于跨度 很小时，可以忽略剪切变形的影响。 4. 考虑剪切变形时的单元刚度矩阵 当梁的高度相对于跨度不是很小时，需要考虑剪切变形的影响，为此，需要对 式（6.63）的刚度矩阵做如下修正 例 按考虑和不考虑剪切变形影响两种情况，计算图示矩形截面悬臂梁在端部横向力 F作用下的端点挠度。截面高为 h，宽为b 。忽略轴向变形。 用一个单元，只有两个节点。 忽略轴向变形，每个节点只有2个自由度， 即y 方向的移动和绕z 轴的转动。 节点1 固定， 总体平衡方程只有与节点2对应的2个方程。 按不考虑剪切变形影响， 根据式（6.63）的单元刚度矩阵，可得总体平衡方程： 使用一个单元，就可以得到与材料力学的理论解相同 的位移结果，这是因为这种梁单元的挠度位移函数是 三次的，悬臂梁的挠度函数也是三次的。 考虑剪切变形影响 总体平衡方程： 矩形截面： 当 ，括号中第二项等于0.54%。 平面梁单元和三维梁单元 平面梁单元：可以处在平面内的任意 方位，与坐标轴平行，或不平行。 取两个直角坐标系： 夹角为 。 （1）节点力与节点位移 局部系 总体系 （2）节点力与节点位移转换关系 令 局部系下的节点位移用总体系下的表示 同理，局部系下的节点力用总体系下的表示 由式（6.62） 局部系下刚度方程 将式（6.69） （6.70）代入（6.65），得到 用 左乘式（6.71）两端，得总体坐标系下、任意方位平面梁单元的节点力与节点位移关系： （1）局部坐标系下的三维梁单元 三维梁单元能够承受： 轴力： 剪力： 弯矩： 扭矩： 每个节点6个自由度（广义）： 三个方向的位移： 绕三个轴的转角： 刚度方程 式（6.62）给出了局部系下节点力 与节点位移 的关系 同理，可以给出局部系下节点力 与节点位移 的关系 局部系下空间梁单元刚度方程（节点力与节点位移关系）： 节点力： 节点位移： 刚度矩阵： 空间梁单元 （2）总体坐标系下的三维梁单元 总体坐标系下的梁单元刚度矩阵 可通过如下转换得到 [T]是三维转换矩阵，可以根据单元的总体节点坐标求出，具体计算方法可见文献[1]。 5.3 板壳单元 * 对于板壳结构,原则上说，可以采用二维或三维实 体单元对其进行结构分析，但是，可能会出现： （1）单元的三个方向的尺寸相差很大，从而导致形成的有限元总体方程为病态或奇异，最后使解丧失精度，或无法求解。 （2）为避免上述问题，若使单元在各个方向尺寸相近，将导致单元总数和节点总数过分庞大，也会使计算无法进行。 * 与板理论和壳理论对应，在有限元法中，有板单元和壳单元。 * 对板壳结构应该采用板壳单元。 板与壳的区别： 板只考虑弯曲变形，壳考虑弯曲和中曲面内的变形。 本章只介绍比较简单的板单元和壳单元。 5.3.1 薄板的弯曲变形 工程中平板构件的长度或宽度要小得多，如果厚度与板面内的其他尺寸之比小于1/5，可以认为是薄板。 位于板厚度中点的平面称为板的中面， 薄板所受的载荷可以分为两类： 作用在中面内的纵向载荷，认为沿厚度均匀分布，可以按平面应力进行计算， 垂直于中面的所谓横向载荷。横向载荷使板弯曲，可以按薄板弯曲问题进行计算。 这里的薄板弯曲是指弹性薄板小挠度弯曲。 在忽略横向剪切变形的前提下， 克希霍夫 Kirchhoff 薄板弯曲理论假设： * 直法线假设，即变形前垂直于中面的法线， 变形后仍垂直于中面。 * 板的中面只发生弯曲变形，没有面内的 伸缩变形，即中面无变形。 * 忽略沿厚度方向的挤压变形， 即 ： 由上述假设，可以写出薄板弯曲问题的位移、应变、应力和内力矩： （1）位移 ， z方向挠度只是 x,y的函数， 板内各点的位移为 2 应变 由于 w仅是坐标x,y 的函数，虑 三个所以只需要考应变分量 根据几何方程和式，应变可以