FE_Ch05.1_3有限元应用中的实际考虑.pptVIP

Chap 5 有限单元法应用中的若干实际考虑; Wilson非协调元的特点和分片试验的意义 及实施方法 子结构方法的特点、使用条件和实施步骤 有限元建模中有效利用结构对称性和周期 性的方法和实施步骤。 ; 分析过程的有效性 计算结果的可靠性;一.单元类型和形状的选择;单元形状;例题：如图是一悬臂梁， E=10MPa ，ν=0.3，L/B=100，端部有一集中力P=10KN作用。计算精度与单元网格划分及网格畸变间的关系。;此悬臂梁在端部的垂直位移需要考虑横向剪切的影响，可按弹性力学解出：;单元（T：三角形 Q：矩形）;网格划分：M1至M5， ——三角形单元数是矩形单元数的2倍 ——网格M5，T3、Q4的自由度为16，T6的自由度为42， Q8的自由度为36。;;;NDLT：总自由度数;网格划分与计算结果;二、网格的划分的考虑;采用自适应分析方法;三种方法：;结点;结 点 线;三、疏密网格的过渡;(b) 采用三角形单元过渡;(c) 采用多点约束方程过渡;区域III;方案1：采用三角形与四边形等参单元过渡;方案2：采用形状不规则的四边形等参单元过渡;方案3：采用四边形等参元附加多点约束方程过渡。其中：@表示4边形5节点Serendipity过渡单元;四、无限域问题;§5.2 应力近似解的性质与处理 ;;故，得到：;显然求 的极小，即：求位移变分 所引起的应变能为极小值的问题；;§5.３ 子结构 (Substructures);例如，对一个实际物体按照3维问题进行离散处理;当实际物体太大，可以考虑先离散成子结构(Substructures)，用超级单元处理;;;一.子结构法的理论根据：;转化矩阵：;2.子结构集成方程的讨论;ii) 尚未集成完毕;二.内部自由度的凝聚;1. 对总体方程的贡献仅外部自由度有关项;2. 由(*)式可以计算 ，并计算各单元应力。;b) 、 正是集成总体方程所需的项。;§5.４ 非协调元与分片试验; Wilson平面非协调单元的思路：引入内部位移自由度，使位移模式具有完全二次式的形式 。然而，内部位移自由度并不和实际的节点相联系，因此可以理解为增加了虚拟的内部节点．这种节点不能保证两个单元的界面位移连续。因此，称为“非协调单元”; 在矩形双线性单元基础上，增加一位移修正项：; 根据这一位移模式，可得应变矩阵为;双线性; 由于[α]e是单元内部位移参数，与单元集合体（整体）的其他单元无关，因此在集成之前可做如下消去(称:静态凝聚); 由此可见做了上述消去后，单元刚度方程形式上和双线性单元完全一样。;二. 三维非协调元;其位移模式包含：;此处：;2.Wilson非协调元的有限元一般格式：;代入? 可得：;可以得到：;二.Wilson非协调元收敛性与分片试验：;则分片试验通过;证明：;被积函数中包含的项为;小结：;§ 5.4 应力计算结果的性质和处理 ;所以，位移元得到的位移解总体不大于真解。;二、应力精度的改进;三、应力抹平的概念;位移元的有限单元法;应力解的近似性：;2、 应力在单元内分布;定义：;4、等参元的最佳应力点;;可视单元内应力平均值，或形心处应力。 所以，平均应力 =;（b）加权抹平 ;(i)构造一个改进的应力解， ;;因此;得到各个结点改进的应力值，如果 是协调， 单元交界面上是连续。;应力在单元内抹平，减少工作量。;或 只算一个方向上的应力；;对于等参元：Gauss积分点的应力值求得后，