FEChFORTRAN程序设计.ppt

Chap. 7 典型有限单元程序设计 §7.7 程序结构和特点 教学程序结构 3．Case 循环 f90 中的新加入 integer day read (*,*) day selectcase(day) case(1) write(*,*) mon’ case(2) write(*,*) tues case(3) write(*,*) wend …… case default write(*,*) wrong! endselect 用FORTRAN语言编制,该程序的主要特点是采用结构化设计方法，即采用模块化方法，把各个主要有限元模拟步骤分别用独立的子程序来实现，而后主要的子程序又依靠辅助子程序来完成次级运算。一个辅助子程序可以为一个以上的主要子程序所调用，主要子程序最后由一个主程序控制。结构化程序条理清楚，具有易读易维护的优点。主要分为三大模块： 1.前处理，即有限元模型的建立。主要包括输入初始数据，形成有限元网格、节点约束的确定及边界节点的形成等。 2.求解。采用三角分解方法进行刚度矩阵组装和求解，3.后处理。根据得到的节点位移,计算单元应变、应力并输出和变形网格。 4.子模块（子程序） subroutine a.??? 特点：独立性强，只要输入输出接口，象一个黑匣子，与外界无关。 b.?? 作用：完成一个独立的功能（求应力，矩阵分解，投放等） 格式： subroutine或function(ip1,ip2-,rp1,rp2,io1,io2,ro1,ro2).(其中ip1,ip2,rp1,rp2,是输入形参，io1,io2.rp1,rp2是输出形参) INITIAL 数组初始化 INPUT 输入几何尺寸、边界条件和材料性能常数 开 始 SOLVER 求解线性方程组 STRESS 求解应力 SUBROUTINE SFR2(ETASP,EXISP) INCLUDE COEFEM C--------------------------------------------------------------------------C EVALUATES SHAPE FUNCTIONS AND THEIR DERIVATIVES FOR LINEAR,QUADRATIC LAGRANGIAN AND SERENDIPITY ISOPARAMETRIC 2-D ELEMENTS. C-------------------------------------------------------------------------C S=EXISP T=ETASP IF (NNODE.GT.4) GOTO 10 ST=S*T SHAPE(1)=(1.0-T-S+ST)*0.25 SHAPE(2)=(1.0-T+S-ST)*0.25 SHAPE(3)=(1.0+T+S+ST)*0.25 SHAPE(4)=(1.0+T-S-ST)*0.25 C C*** SHAPE FUNCTION DERIVATIVES. C DERIV(1,1)=(-1.0+T)*0.25 DERIV(1,2)=( 1.0-T)*0.25 DERIV(1,3)=( 1.0+T)*0.25 DERIV(1,4)=(-1.0-T)*0.25 DERIV(2,1)=(-1.0+S)*0.25 DERIV(2,2)=(-1.0-S)*0.25 DERIV(2,3)=( 1.0+S)*0.25 DERIV(2,4)=( 1.0-S)*0.25 RETURN 10 IF (NNODE.GT.8) GOTO 30 S2=S*2.0 T2=T*2.0 SS=S*S TT=T*T ST=S*T SST=S*S*T STT=S*T*T ST2=S*T*2.0 C** SHAPE FUNCTIONS FOR 8 NODED ELEMENT. SHAPE(1)=(-1.0+ST+SS+TT-SST-STT)*0.25 SHAPE(2)=( 1.0-T-SS+SST)*0.5 SHAPE(3)=(-1.0-ST+SS+TT-SST+STT)*0.25 SHAPE(4)=( 1.0+S-TT-STT)*0.5 SHAPE(5)=(-1.0+ST+SS+TT+SST+STT)*0.25 SHAPE(6)=( 1.0+T-SS-SST)*0.5