有限元程序设计.pptVIP

有限元程序设计 输出结果（场变量和场变量导数等），打印结果。 场变量包括：位移、温度、流场的速度势等； 场变量导数包括：应力、应变、热流、流场速度势等 §1-2 有限元求解器 带宽法：只存储总刚矩阵的半带宽内的元素。 等带宽法：每行的带宽相等，常采用二维数组存储 变带宽法：每行的带宽不等，常采用一维数组存储 波前法是另一种分块储存的变带宽法，其消元次序是按单元编号进行，边组集边消元。 调入内存的单元所保留的节点称为波前节点，所消去的节点称为消元节点，消元节点是与以后调入的单元没有联系的节点，即该节点的方程已组集完。 波前法的回代按消元节点的反序进行。对内存的需求取决与最大波前宽。 V-Fortran 第二章 FEP2 程序的总体设计与输入数据 FEP2 的设计任务 FEP2 的结构设计 FEP2 的主程序 FEP2 的主控程序 FEP2 的数据输入格式与程序实现 第三章 总刚度矩阵的变带宽存贮与求解 总刚度矩阵的变带宽存贮与对角线地址 LDLT 三角分解 自由度指示矩阵ID与对角线地址矩阵 §3-1 总刚度矩阵的变带宽存贮与对角线地址 §3-2 LDLT 三角分解 §3-3 自由度指示矩阵ID与对角线地址矩阵 第四章 单刚与载荷的组集，指定约束的处理，单元分析的预处理 子程序 PFORM 的功能 子程序 MODIFY 完成将指定位移转化成等价内力 将单刚组集成总刚子程序 ADDSTF 子程序 BASBLY 组集载荷与节点反力向量 PRTREAC 输出节点反力 PRTDIS 输出节点位移 ELEMLIB 单元库管理程序 §4-1 子程序 PFORM 的功能 §4-2 子程序 MODIFY 完成将指定位移转化成等价内力 §4-3 将单刚组集成总刚子程序 ADDSTF §4-4子程序 BASBLY 组集载荷与节点反力向量 §4-5 PRTREAC 输出节点反力 §4-6 PRTDIS 输出节点位移 §4-7 ELEMLIB 单元库管理程序 第五章 平面单元设计 平面四节点单元的基本公式 数值积分 形函数及其导数 四节点单刚的构造 3 – 9 节点平面单元 平面单元 PLANE 子程序 §5-1平面四节点单元的基本公式 §5-2 数值积分 §5-3 形函数及其导数 §5-4 四节点单刚的构造 §5-5 3–9 节点平面单元 §5-6 平面单元 PLANE 子程序 节点位移按节点次序重新排列 注意： 1）解得得位移向量为活化自由度，编号小于NEQ 2）约束 (包括指定位移)得值 存放在二维数组 F（载荷）中 SUBROUTINE PRTDIS(ID,U,F) COMMON/CDATA/NDF,NDM,NEN,NJ,NE,NEQ DIMENSION ID(NDF,*),U(*),F(NDF,*),UL(6) DO 102 II=1,NJ,50 WRITE(6,2000)(K,K=1,NDF) JJ=MIN0(NJ,II+49) DO 102 N=II,JJ DO 100 I=1,NDF K=ID(I,N) IF(K.LT.0) U(NEQ-K)=F(I,N) IF(K.LT.0) UL(I)=U(NEQ-K) 100 IF(K.GT.0) UL(I)=U(K) WRITE(6,2001) N,(UL(I),I=1,NDF) 102 CONTINUE 2000 FORMAT(//1X,NODAL DISPLACEMENTS/2X,NODE,6(I8,DISP)) 2001 FORMAT(I6,6E12.4) END FEP2 中已有两种单元： BEAM 、 PLANE 若要增加单元，可修改第七句。例如，已编好一个板单元，名为：PLATE（形参），则程序修改如下： … … GOTO (1, 2, 3) IEL … … CALL PLATE（D , U, X, IX, S, P, I, J, K, ISW) RETURN … … SUBROUTINE ELEMLIB(D,U,X,IX,S,P,I,J,K,ISW) REAL P(K),S(K,K),U(1) DIMENSION IX(1),D(1),X(1) COMMON/CDATA/NDF,NDM,NEN,NJ,NE,NEQ COMMON /ELDATA/