弹性力学及有限元基础 上机实践报告 指导老师： 班级： 机械（茅）101 学号： 1004010510 姓名： 一．程序功能 （1）本程序只用于三节点三角形单元，它可计算平面应力问题，亦可计算平面应变问题。这两类问题用类型码IND 来区别： IND 输入0 的数据——为解平面应力问题。 IND 输入1 的数据——为解平面应变问题。 （2）程序中的载荷包括结点载荷和自重载荷两种，如有其他非节点载荷，则应事先换算成等效结点载荷。 （3）支撑方式，可以是在任一结点的水平和数值支撑，支撑数量不限。但平面弹性体至少应具有保证弹性体几何不变性的三个独立的支撑。 （4）弹性体由单一的材料组成，即只有一组 E,G,u 弹性常数。 （5）弹性体计算简图的厚度是常量。 （6）直要计算机内存允许，本程序对对单元总数，结点总数等没有限制。 （7）本程序输出全部结点位移，对每个但愿输出应力分量，两个主应力及其主平面角。 二．程序中的记号说明 nj 结点个数 ne 单元个数 nz 约束个数 ndd 半带宽 ind 问题类型码 ind=0 平面应力问题 ind=1 nj2 位移分量个数 eo 弹性模量 un 泊松比 gama 材料容重比γ te 单元厚度 ae 单元面积 jm(ne,3)单元结点码数组 nzc(nz)支承数组 cjz(nj,2)结点坐标数组 pj(npj,2)结点载荷数组 b(3,6)几何矩阵 d(3,3)弹性矩阵 s(3,6)应力矩阵 tkz(nj2,ndd)半带储存的整刚矩阵 eke(6,6)单元刚度矩阵 p(nj2)载荷向量，后来存放结点位移 npj1 结点载荷个数加 1。是载荷数组的最大行代码 meo 任一单元码 iask 子程序elest 中的形参，它是计算信息码 iask=1 求单元面积 ae iask=2 求应力矩阵[s] iask=3 求单元刚度[k] ie,je,me 单元meo 的三个角点码 cm,bm,cj,bj 计算形函数中的常cm，bm，cj，bj I,j,k,ii,jj 循环参数 lh 单元行码 ldh 半带存储的整体刚度矩阵的行码 l 单元列码 lz 整体刚度矩阵的列码 ld 半带存储的整体刚度矩阵的列码 i1 工作单元 pe 自重的等效结点荷载 mz 支杆相应的位移分量码 jo 最大列码 j1 工作单元 im 最大行码 c 系数比值 ld1 m 工作单元 wy(6)单元结点位移向量 yl(3)应力向量 sigx sigy toxy 应力分量 “d”x dy 套 xy pyl 平均应力 ryl 应力圆半径 sig 工作单元 sig1 sig2 最大、最小主应力 geta1 工作单元 geta 主平面角 三．子程序名 data 输入 elest（meo，iask） iask=1 ...
