有限元法课程设计题目.docVIP

《计算力学》课程设计任务书 1．用三角元程序计算求如下各单元的应力。 （图中几何尺寸、力学量自定） ` 2．用三角元程序计算求如下各单元的应力。 （图中几何尺寸、力学量自定） 3．用三角元程序计算求如下各单元的应力。 （图中几何尺寸、力学量自定） 4．用平面桁架程序计算求如下各单元的应力。 （图中几何尺寸、力学量自定） 5．用平面刚架程序计算求如下各单元的应力。 6．利用有限元软件进行： 带孔平板的有限元建模与受力分析 7．利用任意有限元软件进行： 受内外压筒体的有限元建模与应力变形分析 8、直杆受集中拉伸外载情形下圣维南原理的数值验证及分析  9、自选题目 各组可自行选择题目, 但自选题目需在经老师认可。 报告格式： 封皮（“《计算力学》课程设计”，班级，完成人员，指导教师） 内容：（设计题目，内容；设计用到的有限元法基本内容，包括原理，公式等；程序，结果；如用软件做设计，写出软件操作过程，截图），A4输出，左侧面订。 以班级为单位，6月21日前交到建业楼409。 三角元程序 4.3.1本节利用平面三角元，对结构体进行离散处理，通过程序运算，可得出单元应力与节点处位移。程序适用于平面应力问题，如果对形成弹性矩阵子程序MODPS作适当修改，就可用于计算平面应变问题。程序适用于各种边界条件及各种载荷。 4.3.2平面三角元程序框图 4.3.3子程序说明 本程序包含四个子程序，其主要功能分别是： 1.MODPS（DMATX，YOUNG，POISS） 该子程序是形成弹性矩阵[D]的，[D]阵的平面应力问题表达式为： (4-56) 如果要计算平面应变问题，可将上式中弹性模量、泊松比值分别用、值代替。 2、BMATPS（INELE，COORD，LNODE，BMATX，AREA，NPOIN，NELEM） 该子程序用来计算应变矩阵[B]，[B]的表达式为： (4-57) 3、DBE（DMATX，BMATX，DMATX） 该子程序用来形成应力矩阵[S]，[S]阵的表达式为 (4-58) 4、GUASS（NZERO，LZERO，ASTIF，NHBW，ALOAD，NP2，NPOIT） 该子程序是用来进行支座处理和解方程的，最后输出节点处位移分量。支座处理的基本原理是如果位移受到约束，则将其对应刚度对角元素变为1，该对角元素所对应的行和列的其余元素变为0，并且将右端项约束位移对应元素充为0。解方程应用带消去法，解出位移后，将节点位移分量存放在原来存放右端项的标识符ALOAD中，最后将ALOAD（节点位移）输出。 4.3.4数组及变量说明 数组说明 COORD 节点坐标数组 LNODE 单元节点数组 ALOAD 荷载（右端项）数组，解出位移后，将节点处的位移充于其中 ESTIF 单元刚度矩阵，为6×6阶方阵 ASTIF 总刚度矩阵，带状存贮 DBMAT 应力矩阵[S]阵，为3×6阶 DMATX 弹性矩阵[D]阵，为3×3阶 BMATX 应变矩阵[S]阵，为3×6阶 STRES 单元应力矩阵，为3×1阶 JJS 受载节点矩阵 ZX Ｘ方向已知载荷向量 ZY Y方向已知载荷向量 DISP 单元节点位移矩阵 LZERO 约束位移矩阵 变量说明 NPOIT 最大节点数 NELEM 最大单元数 NNODE 单元节点数（三角元中为3） NLOAD 最大受载节点数 NZERO 最大约束位移个数 YOUNG 弹性模量 POISS 泊松比 TE 板厚（假设t=1） NHBW 半带宽 NP2 位移总数 AREA 单元面积 已知一对角受压的正方形薄板，厚度 为1cm，荷载沿厚度均匀分布，为2Kg/cm2, 泊松比，求板内点的应力与位场。 题中由于XZ 面及YZ面均为该板的对称面，所以只选取1/4部分作为计算模型 如图4-12 计算模型的输入数据有以下各量: NPOIN NELEM NNODE NLOAD NZERO YOUNG POISS TE 6 4 3 1 6 210000 0.3 1.0 JJS 1 ZX 0 ZY -1.0 I COORD 1 2 3 4 5 6 COORD(I,1) 0.0 0.0