FEM细部分析方法.ppt

1、有限元分析方法,北京迈达斯技术有限公司 (http:/www.midasIT.com),桂 满 树,顺 序,板单元/实体单元的特点及正确使用方法 通过例题说明各种建模方法 建立板单元网格的方法 使用扩展功能建立实体单元的方法 实际模型例题 实际工程中细部精密分析的方法,板单元的特点 (1),平面应力zz = xz = yz = 0 薄板 (Kirchhoff Plate) 忽略剪切变形的影响 1-D: Euler-Bernoulli Beam 厚板 (Mindlin Plate) 考虑剪切变形的影响 1-D: Timoshenko Beam 大部分情况可选用厚板(误差不到 2%), 非常薄的板

2、应使用薄板 防止Shear Locking,板单元的特点 (2),板单元没有面内旋转(In-plane Rotation)自由度-五个自由度/每个节点 由厚度来体现面内和面外的刚度 挠度比板单元厚度薄的时候, 可忽略面内变形,MIDAS的板单元,平面内特性- 三角形: LST (Linear Strain Triangle)- 四角形: Plane Stress Formulation with Incompatible Modes 平面外特性 薄板- DKT/DKQ (Discrete Kirchhoff Tria./Quad.)- DKQ: Taylor Kinematic Couplin

3、g)的功能是在不太重要的位置上将结构连接起来(相对运动)，并传递荷载。 使用刚性连接时, 在连接位置在某一方向上位移不是连续的(相同)，应力分布也不是很圆滑 从属节点本应该依靠外部荷载而产生位移, 但因为被设置为从属于主节点，所以不能产生与相邻节点的正常位移 位移不连续 应力的不连续发生在距连接位置单位特性长度(一般为厚度或高度尺寸)的局部范围内, 该范围内的应力不可信 刚性连接应尽量使用于距重要位置23倍长度范围以外 在受扭(Torsion)位置最好不要使用刚性连接。因为刚性连接约束了截面的翘曲(Warping)，所以会夸大结构的抗扭刚度 管型,对称条件,对称结构最好利用结构的对称性进行分析

4、 建模简便，结果对称 在MIDAS/Civil中可以将简化后的模型按对称条件输出整体模型 对称条件 几何形状、材料、荷载、边界条件均应对称 边界条件应设置为不能让结构的变形越过对称面 特征值分析/屈曲分析中不能使用对称条件 因为模态不是对称的 荷载的大小也应满足对称条件,荷载的处理 (1),在节点处作用集中荷载时，在节点处容易发生应力奇异(Stress Singularity)现象 平面弹性问题、节点支承 集中荷载作用下的应力奇异性随网格密度的增加而增加直至。 板单元/实体单元网格中的集中荷载的处理 垂直于面的荷载: 压力荷载 板单元端部的竖向荷载 压力荷载中的边压力荷载 虚拟梁和梁单元荷载

5、其他(如: 平面内荷载)可将相应节点刚性连接后， 在主节点处加集中荷载,荷载的处理 (2),当不可避免地需要加集中荷载时, 较理想的处理方法如下(实际使用起来同样有难度) 不使用相应位置的分析结果 只使用St. Venant原理适用的范围的结果 在周边建立较细的三角形单元网格, 忽略相应位置的分析结果 在非常小的范围内用均布荷载替代集中荷载 在集中荷载位置删掉非常小的单元，用均布荷载代替。,单元应力和节点应力,以位移连续为基础的有限元分析中，节点应力值是不精确的(最大误差 30%). 中心点的应力是较为准确的 单元应力是不连续的, 节点应力(绕节点平均值)是单元应力的平均值 使用节点应力应该比

6、使用单元应力更合理一些 最大应力值: 节点应力 单元应力,轴力(q)作用下的桁架单元的节点位移,板单元的应力特点,计算板单元的节点应力时，板顶应力与板底应力分别取平均值 相邻板单元的法线方向(z轴)不同时，绕节点平均法计算的节点平均应力可能会计算有误。,在MIDAS/Civil中查看板单元应力结果时，在坐标系选项中选择整体坐标系时，因为输出的是整体坐标系方向的应力，所以可防止因法线方向的不同引起的错误。 (默认选项为整体坐标).,根据材料特性的应力评价 (1),根据材料特性的应力评价(2),von Mises Stress (Maximum Distorsion Energy Theory)

7、组合应力作用下的变形能与抗拉试验的变形能相同时，材料发生屈服 与试验结果最接近 Tresca Stress (Maximum Shear Stress Theory) 最大剪切应力与抗拉试验的最大剪切应力(屈服应力的1/2)相同时，材料发生屈服 接近试验结果，提供相对保守的结果，所以一般用于实际设计(特别是危险部位设计)以及抵抗剪切的部位(锚拴焊接部位等) Principal Stress (Maximum Principal Stress Theory) 最大主应力超过极限强度(u)时，材料发生屈服 与脆性材料试验结果最接近,根据材料特性，使用下列应力评价结构的性能,实用的分析结果 (1),

8、应力向量输出三个主应力方向的应力向量可以查看应力的向量图，可以检查出应力集中的位置可以很方便地查出超过容许应力的位置，以便于加强该位置的截面或采取相应措施。 剖断面可以给出任意剖断面上的实体单元的应力 特别适用于实体单元，与动画功能相结合可方便地查看实体单元内部的应力。 等值面可以输出试题单元内任意应力值的等值面。,实用的分析结果(2),剖断面图可以输出任意剖断线或剖断面上的板单元的应力图和内力图，对于计算配筋和制作计算书非常实用特别是可依输出任意方向的内力，对查看斜桥主应力方向的内力非常实用 局部方向内力的合力利用某平面上板单元或实体单元上各节点的内力，可求出作用于选择单元形心位置上的各方向

9、的内力特别是对实体单元，一般只输出应力结果，用于实际配筋计算很不方便。使用该功能可以使设计人员非常方便地计算出内力。,二维单元网格的类型 (1),结构化网格 (Structured Mesh) 内部节点的Valence均为4的网格(内部节点的位相相同) 目标的限制(1) 必须为由四条曲线组成的封闭区域 (2) 相对的两条曲线上有相同的节点数量 FEmodeler: Map-mesher,Mapping,目标曲线不是四个的例子,二维单元网格的类型 (2),非结构化网格 (Unstructured Mesh) 内部节点的Valence维持在35较好 对目标没有限制 FEmodeler: Auto-

10、mesher (Grid Mesher, Flex Mesher) 内部可以包含点和线,一般使用Grid Mesher, 当尺寸控制超过1:2时可使用Flex Mesher,建模时的参考事项,分割 尽量细分模型 尽量先建立最小的标准单元或对称单元，然后利用复制和镜像功能 查看是否能由低次单元扩展为高次单元 (扩展后可删除低次单元) 使用最小的标准单元建模 事先建立建模计划 越是复杂的部分，越要努力寻找复合规律的方案 尽量使用结构化网格 注意倾斜角(Skew Angle) 要控制网格密度，不要建立过多的单元 只细分重要的部位 检查自由边和自由面，确认单元之间的连接,二维单元网格的划分 (1),W

11、ireframe 模型 (平面),(1),(2),Auto-mesh Planar Domain,(3),Extrude Planar Linear Array Curve,Auto-mesh + Extrude,二维单元网格的划分 (2),二维单元网格的划分 (3),Wireframe 模型,Curve Project on Plane,(1),(2),(3),Map-mesh 4-curve Surface,(4),Map-mesh 4-curve Surface,Mapped-mesh,三维单元网格划分时的参考事项,必须决定对复杂位置的处理方法后再开始建模 检查可使用的扩展方法 首先检查

12、能否使用Stack Method 在一个平面内画出所有截面(也可使用投影命令), 然后生成二维单元网格, 然后扩展成为三维网格 选择各种扩展方法 复制、旋转、投影(平面/实体)、滑行 沿一个方向扩展较难时，查看是否能向多个方向扩展 检查是否能组合使用各种扩展方法 变截面时，检查是否可先扩展为等截面，然后调整(投影)节点位置 检查分区域扩展的可能性，特别是分区域后使用方法3的可能性,三维单元网格的扩展方法,Planar Mesh,Translate Planar Mesh,Project SolidMesh,Translate Planar Mesh + Mirror,建立三维网格 (1),St

13、ack Method,首先在平面上画出所有截面,Auto-mesh Planar Domain,Extrude Solid Linear Array Planar Element,(1),(2),(3),第一次扩展,第二次扩展,建立三维网格 (2),建立三维网格 (3-1), 对称模型,Wireframe 模型,Auto-mesh Planar Domain,Extrude Solid Glide Planar Element,Extrude Solid Glide Planar Element,(1),(2),(3),(4),建立三维网格 (3-2),Extrude Solid Linear

14、 Array Planar Element,(5),(6),Curve/Element Translate (Copy) Auto-mesh Planar Element,(7),Extrude Solid Linear Array Planar Element,建立三维网格 (4-1), 对称模型,Wireframe 模型,(1),(2),Auto-mesh Planar Domain,建立三维网格 (4-2),(3),Extrude Solid Project PlanarElement,(4),Node Translate Node Project on Plane,(5),Extrud

15、e Solid Linear Array ProjectPlanar Element,(6),Extrude Solid Project Solid Element,建立三维网格 (5),建立三维网格 (6), 对称模型,(1),Wireframe 模型,(2),Map-mesh 12-curve Volume,(3),Auto-mesh Planar Domain,(4),Extrude Solid Glide Planar Element,建立三维网格 (7),建立细部分析模型时的注意事项,利用梁单元模型分析结果建立细部分析(使用强制位移)模型时，细部模型的各端部位置(1) 要与整体模型的

16、节点位置一致,(2) 要距离重要位置相当远的距离 参照前面刚性连接说明 要注意使细部模型的坐标系与整体模型的坐标系一致 强制位移的各位移分量方向一致 细部分析不正确最常见的原因就是细部模型与整体模型的坐标系不一致所致。 输入强制位移时，要输入六个位移分量，并尽可能加大有效数值位数 特别是对于旋转自由度分量,将细部分析模型的各端部的中心位置建立一个主节点，然后将各端部分别设置刚性连接，最后将整体模型分析中得到的位移值(6个位移分量)作为强制位移加到各端部的主节点上。,标准方法,细部分析方法 (1),细部分析方法 (2),当细部分析模型的端部位置不在节点上时，可在主节点处建立一段梁单元，梁单元长度为从主节点到整体模型的节点位置。,强制位移加载位置,细部分析方法 (3),直接将细部分析模型插入到杆系模型中的方法,细部分析方法 (4),插入,Top View,梁+板单元 (B)模型,细部分析板单元 (A)模型,板单元(A)、 (B)部分靠螺栓连接为一体。两个模型分别在螺栓位置建立节点，然后用刚性连接在相应螺栓节点位置将两个模型连接起来。将对应的螺栓节点的一点设置为主节点，另一点设置为从属节点。,使用了U型加劲肋的敞开型桥梁 的横向支撑连接和横膈板细部分析模型,细部分析实例 (1),拱肋和横向支