有限元分析及工程应用216第四章.ppt

4.4 典型实例分析及详解 ， ， ， 、 、 E,1,2,3 ! 由3个节点生成一个单元，共生成4个单元 E,2,4,5 E,2,5,3 E,3,5,6 FINISH ! 退出当前处理器 /SOL ! 进入求解器 F,1,FY,-1000 ! 在节点1上施加集中力 NSEL,S,LOC,X,0 ! 选取X=0的所有节点 D,ALL, , , , , ,UX, , , , , ! 在选取节点上施加UX约束 NSEL,S,LOC,Y,0 ! 选取Y=0的所有节点 D,ALL, , , , , ,UY, , , , , ! 在选取节点上施加UY约束 ALLSEL,ALL ! 选取所有实体 SOLVE ! 求解计算 4.4 典型实例分析及详解 ， ， ， 、 、 FINISH ! 退出当前处理器 /POST1 ! 进入后处理器 PRNSOL,U,COMP ! 列表节点位移 PRESOL,S,COMP ! 列表单元应力 FINISH ! 退出当前处理器 谢 谢! 4.2 等参数单元与数值积分 （2）单元矩阵的变换 2) 平面等参单元的刚度矩阵 2.几何矩阵 将位移表达式代入到几何方程: 4.2 等参数单元与数值积分 （2）单元矩阵的变换 2) 平面等参单元的刚度矩阵 3.单元刚度矩阵 4.2 等参数单元与数值积分 （2）单元矩阵的变换 2) 平面等参单元的刚度矩阵 3.单元刚度矩阵 单元的刚度矩阵为： 单元厚度 4.等效节点载荷 单元体积力的等效节点载荷计算公式为： 单元边界力的等效节点载荷只在该单元的受载边界上存在，其积分变量变换式为 4.2 等参数单元与数值积分 （2）单元矩阵的变换 2) 平面等参单元的刚度矩阵 4.等效节点载荷 （3）高斯积分公式 等参数单元刚度矩阵和节点载荷矩阵列阵中，需要进行如下形式的积分运算： 和 一元的高斯积分公式可表述为： 当被积函数比较复杂时，必须要采用数值积分来代替函数积分。 4.2 等参数单元与数值积分 （3）高斯积分公式 一元高斯积分公式中: 对于二维(矩形)区域，通过先对一个坐标进行积分后再对另一个坐标进行积分，就可以得到高斯积分公式为: 4.2 等参数单元与数值积分 （3）高斯积分公式 对于矩形区域四个高斯积分点位于 权系数 4.3 平面问题的较精密单元 由于3节点三角形单元采用的位移模式是线性的，从而其形函数也是线性。 ， ， （1）6节点三角形单元 取位移函数 在单元的边界上，位移分量是按抛物线变化的，每条公共边界上有3个公共节点，正好可以保证相邻两单元位移的连续性。 用面积坐标表示的二次形函数： 角点： 4.3 平面问题的较精密单元 ， ， （1）6节点三角形单元 边中点： 利用上述形函数，可将位移函数表示如下： 或 4.3 平面问题的较精密单元 ， ， （2） 10节点三角形单元 单元位移函数 每边4个节点，可以保证相邻单元之间的位移连续性。 以面积坐标表示的三次形函数 角点： 边三分点： 内点： ， 4.3 平面问题的较精密单元 ， ， （3） 4节点矩形单元 ， 位移函数: 将相应的节点坐标代入上式，即可求出这8个未知量,代入上式 矩阵形式 4.3 平面问题的较精密单元 ， ， （3） 4节点矩形单元 ， 双线性函数 代入弹性力学的几何方程式，可求得单元的几何矩阵为： 4.3 平面问题的较精密单元 ， ， （4） 8节点等参数单元 ， 位移模式 形函数分别取为 4.4 典型实例分析及详解 ， ， ， 解：由于两个对角线均为该板的对称轴，所以只需取板的四分之一作为计算对象即可。根据板的受力情况，应按平面应力问题处理。将讨论对象划分为4个三角形单元。根据计算简图。约束条件为： 选择对称轴作为坐标轴，则各节点的坐标值如表所示: 4.4 典型实例分析及详解 ， ，