07弹性平面问题程序设计研讨.pptx

弹性平面问题程序设计 有限元法课件 石油大学力学系 周 博 一、弹性平面问题的程序结构 [enter] 二、主控变量和数据描述 [enter] 三、弹性平面问题的功能函数 [enter] 四、弹性平面问题的计算程序 [enter] 五、弹性平面问题的数值算例 [enter] 弹性平面问题程序设计 一、弹性平面问题的程序结构 与平面桁架的程序设计相似，在进行弹性平面问题有限元法程序设计时，应首先设计出反映计算内容和顺序的程序结构图，以便于程序的顺利编制和调试。 根据弹性平面问题的特点，其程序结构图如右图所示。 二、主控变量和数据描述 根据弹性平面问题的几何和力学特性，本章编制的弹性平面问题的有限元分析计算程序的主控变量包括6个，它们分别为： N —— 结点总数； LR —— 已知结点位移分量总数； MT —— 三角形单元总数； LTB —— 单元的材料特征类型数； LPQ —— 外载荷总数； XEV —— 应力状态信息（0为平面应力，1为平面应变）。 弹性平面问题的主控变量 以图示弹性平面应力问题的有限元模型为例，假设板厚均匀、上半部分和下半部分材料不同，其主控变量： 结点总数 N=9； 已知结点位移分量总数 LR=4； 三角形单元总数 MT=8； 单元的材料特征类型数 LTB=2； 外载荷总数 LPQ=4； 应力状态信息 XEV=0。 二、主控变量和数据描述 弹性平面问题的数据描述 二、主控变量和数据描述 1. 结点坐标数组 在计算单元刚度矩阵、单元应力及结点载荷时须用到结点坐标。将结点坐标数组命名为XY，它是一个N行3列的二维数组，其中N为主控变量-结点总数。 如图所示弹性平面问题有限元模型的结点数N=9，结点坐标数组 2. 已知结点位移数组 结点受约束，即结点位移是给定的，属于结点具有已知位移的情况。将已知结点位移数组命名为SU，它是一个LR行、2列的二维数组，其中LR为主控变量-已知结点位移分量总数，它应大于或等于3。 对于图示的弹性平面应力问题有限元模型，已知结点位移分量总数LR=4，已知结点位移数组 二、主控变量和数据描述 3. 单元信息数组 三结点三角形单元的结点的编号应按逆时针方向进行，否则会导致非正定的单元刚度矩阵，得到错误的计算结果。在单元定义中，单元的任一结点都可以作为第一个结点。 将单元定义数组命名为ELT，它是 一个MT行、5列的二维数组，其中MT为主控变量-单元总数。 对于图示弹性平面应力问题有限元模型，单元总数MT=8，材料特征类型总数LTB=2，单元信息数组 二、主控变量和数据描述 4. 材料特征数组 材料特征数组存放的内容包括弹性摸量，泊松比、单元厚度和比重的x、y方向分量Wx和Wy。 材料特征数组命名为EMUTW，它是一个LTB行、5列的二维数组，其中LTB为主控变量-材料特征类型总数。 对于图示的弹性平面应力问题有限元模型，材料特征类型总数LTB=2，若不计体力，材料特征数组 二、主控变量和数据描述 5. 载荷数组 将载荷数组命名为ELPQ，它是一个LPQ行、6列的二维数组，其中LPQ为主控变量-外载荷总数。 第1列：受载单元边界的结点号i； 第2列：受载单元边界的结点号j ； 第3列：集中力的大小或分布力的载荷集度qi ； 第4列：集中力到结点i的距离li或分布力的载荷集度qj ； 第5列：载荷方向与x方向的夹角； 第6列：载荷信息类型，集中力为0、分布力为1。 二、主控变量和数据描述 对于图示的弹性平面应力问题有限元模型，主控变量 LPQ=4，其载荷数组 三、弹性平面问题的功能函数 单元刚度矩阵函数 用MATLAB编写如下一个用单元子方阵函数krst，生成三结点三角形单元刚度矩阵的子方阵。 function krs=krst(xye,e,t,mu,r,s) %*******变量说明**************** %xye(3,2)----单元结点坐标数组 %e----------弹性模量 %t----------单元厚度 %mu-------- 泊松比 %r,s --- 1,2,3 %******************************** 根据上图所示单元刚度矩阵函数的程序结构图，并利用单元子方阵函数krst，用MATLAB语言编制平面3结点三角形单元刚度矩阵函数dan_gang。 function ke=dan_gang(xye,e,t,mu) %*******变量说明************* %ke(6,6)--单元刚度矩阵 %xye(3,2)---单元结点坐标数组 %e----------弹性模量 %t----------单元厚度 %mu--------泊松比 %*************************** 三、弹性平面问题的功能函数 %---