有限元第二章-平面问题有限元方法.pptVIP

二、常应变三角形单元分析 5.单元载荷及等效节点力 同单元刚度矩阵的组集一样，将单元载荷矩阵按扩阶的方法进行按点叠加，以获得结构的载荷列阵 ●结构载荷列阵的形成 此时所不同的是：载荷应含三项内容：集中力、面力、体力。对于平面问题，结构载荷为2n阶列阵，即 二、常应变三角形单元分析 6.结构平衡方程 外载荷的势能 结构的总变形能为 于是，弹性体的总位能为 二、常应变三角形单元分析 6.结构平衡方程 根据最小势能原理 可以推导出 它是结构的静力控制方程，亦称刚度方程、结构平衡方程，求解上述节点位移，即可获得问题的解。 二、常应变三角形单元分析 7.边界条件引入 由于总体刚度矩阵结构的总变形能为奇异阵，不存在逆阵，从数学上讲，结构平衡方程不能直接求解，必须引入一定的补充条件，消除弹性体由于外力作用而产生的刚体位移，使刚度阵从半正定成为正定矩阵，从而获得唯一解。 边界条件的种类分3类： 力边界 位移边界 混合边界 有限元位移法中仅涉及位移边界。位移边界又分成零位移和非零位移，通常将位移边界条件简称为边界条件。 边界条件的引入最常用的有主元改“1”法、主元乘大数法、方程降阶法3种。 二、常应变三角形单元分析 7.边界条件引入 边界条件引入最常用3种方法 有主元改“1”法 主元乘大数法 方程降阶法 二、常应变三角形单元分析 7.边界条件引入 二、常应变三角形单元分析 7.边界条件引入 ●主元改“1”法 设第i个节点x方向的位移为ui=0,将结构平衡方程作如下处理，对总刚阵 的第2i-1行和第2i-1列的对角线元素改为“1”，其它各元素改为“0”，节点载荷向量中第2i-1分量也改为“0”。 二、常应变三角形单元分析 经过上述修改后的方程为 于是原来的第2i-1个方程变为 7.边界条件引入 二、常应变三角形单元分析 问题与主元改“1”法相同，将第2i-1行对角线元素乘上一个大数，例如1010以上，这时修改后的方程为 7.边界条件引入 ●主元乘大数法 二、常应变三角形单元分析 7.边界条件引入 于是原来的第2i-1个方程变为 ●主元乘大数法 用K2i-1，2i-1 ×1010除以第2i-1行和第2i-1列的元素后，近似有 二、常应变三角形单元分析 7.边界条件引入 ●非零位移的引入 对于非零指定位移边界条件的引入，完全可以用上述两种方法，处理时稍微复杂些。 对于主元改“1”法，其余步骤相同，只是将节点载荷向量中第2i-1分量改为“ub”。 对于主元乘大数法，其余步骤相同，只是将节点载荷向量中第2i-1分量改为“K2i-1，2i-1 ×1010 × ub”。 二、常应变三角形单元分析 7.边界条件引入 将所有位移为零的节点自由度对应的所在所有行和列的元素划去，使总体刚阵降阶，从而实现边界条件的引入。 ●方程降阶法 二、常应变三角形单元分析 8.平衡方程求解及后处理 求的平衡方程的解后，还需要进行后处理，包括计算单元应变、单元应力以及强度和刚度分析等，也包括应力、应变的三维可视化。 三、有限元分析应该注意的问题 1.节点选择与单元划分 1.集中力作用点、分布力突变点、支撑点应选做节点 2.不同厚度、不同材料的部分不应划在同一单元 3.应力变化大处单元应密集一些，节点的多少与疏密要考虑计算机的容量和计算精度 4.单元边长的边长之比应尽可能接近1 5.相邻单元的尺寸尽可能接近 6.节点所连接的单元个数尽可能一致 宜 不宜 三、有限元分析应该注意的问题 1.节点选择与单元划分 三、有限元分析应该注意的问题 2.节点编号 尽可能使相关节点的节点编码差值最小 总刚阵半带宽=（节点编码最大差值+1）＊节点自由度数 总刚阵半带宽=（7+1）＊2=16 总刚阵半带宽=（2+1）＊2=6 总刚阵需占用的存储空间为 总刚阵需占用的存储空间为 16＊14＊2=448 6＊14＊2=168 三、有限元分析应该注意的问题 3.结构对称性 三、有限元分析应该注意的问题 4.应力结果的处理 位移的计算结果一般比应力、内力结果精度高。位移达到满意结果，由几何方程求应变，再由物理方程球应力，结果的精度较差。上述三角形单元为常应变和常应力，矩形单元应力线性变化，而工程问题的应力是比较复杂的。为更好地反映实际应力情况，需要对计算结果进行整理。常用处理方法有两种：