fortran调用ansys做桁架结构优化例子.pdfVIP

fortran 调用ansys 做桁架结构优化例子 ANSYS 命令流部分： /clear B=1000 !基本尺寸 !进入PREP7 并建模 /PREP7 *dim,A,,3,1 ！定义一个3*1 的数组 *VREAD,A,E: \project\trussdata_1,txt,,IJK,3,1 !IJK 表示读入的顺序 (F7.2) A1=A(1,1) A2=A(2,1) A3=A(3,1) *cfopen,E: \project\trussdata_3,txt *vwrite,A1,A2,A3 (F7.2,/F7.2,/F7.2) *cfclos ET,1,LINK1 !二维杆单元 R,1,A1 !以参数形式的实参 R,2,A2 R,3,A3 MP,EX,1,2.1E6 !杨氏模量 N,1,-B,0,0 !定义结点 N,2,0,0,0 N,3,B,0,0 N,4,0,-B,0 E,1,4 !定义单元 REAL,2 E,2,4 REAL,3 E,3,4 FINISH ! !进入求解器，定义载荷和求解 /SOLU D,1,ALL,0,,3 !结点UX=UY=0 F,4,FX,200000 !结点4 上的X 方向载荷分量(工况1) F,4,FY,-200000 !结点4 上的Y 方向载荷分量(工况2) SOLVE FINISH !进入POST1 并读出状态变量数值 /POST1 SET,LAST ETABLE,EVOL,VOLU !将每个单元的体积放入ETABLE SSUM !将单元表格内数据求和 *GET,VTOT,SSUM,,ITEM,EVOL !VTOT=总体积 RHO=2.85E-4 WT=RHO*VTOT !计算总体积 ETABLE,SIG,LS,1 !将轴向应力放入ETABLE *GET,SIG1,ELEM,1,ETAB,SIG !SIG1=第一个单元的轴向应力 *GET,SIG2,ELEM,2,ETAB,SIG !SIG2=第二个单元的轴向应力 SIG1=ABS(SIG1) !计算轴向应力的绝对值 SIG2=ABS(SIG2) !/ESHAPE,2 !以实体单元模式显示壳单元 !/VIEW,1,1,1,1 !轴测视图 !EPLOT !画单元 FINISH *cfopen,E:\project\trussdata_2,txt *vwrite,SIG1,SIG2,WT (F7.2,/F7.2,/F7.2) *cfclos FORTRAN 程序部分： PROGRAM TRUSS USE DFLIB IMPLICIT NONE character(len=20) :: filename1=trussdata_2.txt character(len=20) :: filename2=trussdata_1.txt REAL,parameter :: e=0.01 !e 为计算精度 REAL,parameter :: xu=400 !xu 为许用应力 INTEGER :: fileid,k INTEGER i LOGICAL(4) result REAL u1(50),u2(50),ww(50),umax !u1,u2 为应力比,umax 为最大应力 REAL sig1,sig2,W,ap !sig1 为杆1 的最大应力,sig2 为杆2 的最大应力 REAL x1(50),x2(50),xx1,xx2 !x 为杆截面面积 i=1 k=1 x1(1)=707 !定义杆的初始值 x2(1)=707 ap=10 DO WHILE( ap=e ) !调用ANSYS 做结构分析 result=SYSTEMQQ( C:\Program Files\Ansys Inc\