一个ansys中用fortran77编的routine的范例.docVIP

一个ansys中用fortran77编的routine的范例 发信人: abaya (abaya), 信区: FEA 标 题: Re: 谁能贴一个ansys中用fortran77编的routine的范例 发信站: BBS 水木清华站 (Thu Jun 6 15:53:47 2002) 我为了能够将自己编写的有限元分析计算结果与ANSYS比较，曾经改编过ANSYS的一个用户 接口文件，该文件用于读取ANSYS模态分析（子空间迭代方法）形成的.full文件，提取整 体质量和刚度矩阵并输出至文件。介绍如下： 共有两个文件： Extract.for (fortran 77,主文件) matrixout.f90 (fortran 90,自己编的子程序) 开发环境：COMPAQ FORTRAN 6.5 操作系统：WIN2000 Extract.for: ****************************************************************************** program Extract c *** primary function: Extract global mass and stiffness matrix from full file. c *** secondary function: Read and reformat full file c c *** zx(r) copyright(c) 2001 c *** Xin Zhao c c NOTICE - The .full file as output of ANSYS is needed. c Be sure to set up for modal ANTYPE,2 c and full subspace MODOPT,SUBSP,nmode,0,0, ,OFF c (nmode is not used - it can be any value c If the free-free stiffness and mass matrices are desired, c make sure there are no constraints on the model. c #include impcom.inc c**ZHAOXIN *Use dflib *USE NUMERICAL_LIBRARIES *USE MATRIXOUT_ C*END OF ZHAOXIN external ihsort external binset, binini, binrd, binclo, biniqr, bintfo integer binset, biniqr integer IOLENG, ROWLENG, WAVEFRONT, MAXNODE parameter (IOLENG=16384, ROWLENG=1000, WAVEFRONT=10000, x MAXNODE=100000) integer ivect(100), npage, keyrw, units, code, x nbuf,i, j, nunit, jloc, kbf, kext, n, lbuf, nmatrx, x lll(ROWLENG), indx(ROWLENG),lenbac, nontp, numdof, irow, x l(WAVEFRONT), nrow, mr, nterms, kunit, munit, kdof, x baclst(MAXNODE), node, idof, sortlist(2,MAXNODE), x nmass, nstif, nmax, lumpm integer buffer(IOLENG) double precision krow(ROWLENG), mrow(ROWLENG) character*80 title(2) character*106 pname,mname,kname *character*106 zmname,zkname !zhao xin character*32 jobnam integer