ANSYS 整体刚度矩阵提取子程序——

合集下载

Ansys用于提取刚度矩阵的命令

矩阵类型用3个字符表示，第1个字符可取：R---实数矩阵、C---复数矩阵、P---仅矩阵结构（无元素数值）；第2个字符可取：S---对称矩阵、U---不对称矩阵、H---Hermitian矩阵、Z---病态对称矩阵；R---带状矩阵；第3个字符可取：A---组装的矩阵、E---单元矩阵（未组装）。对称矩阵只存储下三角元素，如结构刚度矩阵为对称矩阵，Harwell-Boeing格式则仅记录下三角元素。
第1行：格式（A72），为文件头的字符型解释，如刚度矩阵或质量矩阵等标题。
第2行：格式（5I14），分别表示该文件的总行数（不包括文件头）、矩阵列指针的总行数、矩阵行索引的总行数、矩阵元素数值的总行数、右边项总行数。
Ansys用于提取刚度矩阵的命令
命令：HBMAT,fname,ext,--,form,matrx,rhs
其中：
Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。
=MASS：输出质量矩阵。可用于特征值屈曲、子结构分析、模态分析。
=DAMP：输出阻尼矩阵。仅用于有阻尼的模态分析。
ቤተ መጻሕፍቲ ባይዱ
第5行：格式（A3,11X,2I14），A3各列分别表示右边项格式、应用高斯起始矢量、应用eXact求解矢量；两个整数分别表示右边项列数、行索引数。三个字符中的第1个字符可取：F---全部存贮（如节点荷载向量的全部元素）、M---与系数矩阵相同方法。
第6行后：矩阵元素值（单列）。
第3行：格式（A3,11X,4I14），分别为矩阵类型、矩阵行数、矩阵列数、矩阵行索引数（对组装后的矩阵，该值等于矩阵行索引数）、单元元素数（对组装后的矩阵此值为0）。
第4行：格式（2A16,2A20），分别表示列指针格式、行索引格式、系数矩阵数值格式、右边项数值格式。

ansys提出刚度和质量矩阵

一、连续钢梁的刚度和质最矩阵弹性模量：EX=2.0x 10n Pa 泊松比：PRXY=0.3 密度：DENS=7811Kg/m3 截面特性如右图。

AN SYS命令流：/FILNAM.BEAM MASS AND KNIFF ANALYSIS/TITLE, BEAM MASS AND KNIFF ANALYSIS/PREP7EI1BEAM3MBEX470E11MRNUXYX0.3MRDENS 丄7800SECTYPE,1,BEAM」,,3SECO 阡SET.CENTSECDATAO150・：15,030・02O02,0・(HO0O0 K10O0K210O0一、框架的刚度和质呈:矩阵的提取模型采用右图的集中质屋模型，Ml=2762Kg ,M2=2760Kg, M3=2300Kg,层间刚度分别为：248500、192100. 152200命令流如下:finish/clear/prep7 et4z combinl4keyopt,1,2,1 !—维弹簧单元（Ux平动）et2mass21 keyopt,2,3/2 !3・D mass没有转动惯最0,1041-1）*3,0 *enddo!质最常数r z2,2762r,3,276054,2300type,2*do,i,2,4,lrealje,i*enddo!层间刚度^12,248500G1*******G14J52200 type」•doj,134 realjl+i eJJ+1 *enddo 哟束nselsnode,24JL d，all,uz，0d,all,uy,0 allsel d,l,all,0 /solu antype z7 seoptsubmatB m,all z all solve selistsubmat.B使用该命令流町以得到结构的刚度矩阵和质最矩阵结果: 刚度矩阵：440600. 00 -192100. 000. 0000000 质量矩阵：2762. 00000. 00000000. 0000000•192100.00344300. 00-152200.000. 0000000-152200. 00152200. 000.00000002760. 00000.00000000.00000000.00000002300. 0000imIT12mik-ki附件1：（运行ANSYS命令的输出结果，最后面是刚度和质最矩阵）PRINT CONTENTS OF SUPERELEMENT submatPRINT OPTION = 3HEADER =8 3 2 3 34 0 4 3 01 1 0 0 10 0 1 12 11550 1101 166 103 submat 1078 0922 928 946 953 976982 1027 1072 0 934940 106004406 0 0 30 0419 0 0 0 00 0 0 0HEADER SUMMARY:NUMBER OF ROWS = 3NUMBER OF MATRICES = 2NUMBER OF EDGE PLOT LINES = 3NUMBER OF DEGREES OF FREEDOM PER NODE = 3NUMBER OE DEGREE OF FREEDOM INDICES = 4NUMBER OF NODES = 3NUMBER OF LOAD VECTORS = 1NUMBER OE TRANSFORMATIONS = 0BASE FILE NAME二submatMAXIMUM STIFFNESS二0. 44060E+06DEGREES OF FREEDOM PER NODE =DEGREES OF FREEDOM PER NODE LABELS =UX UY UZDEGREES OF FREEDOM AS GENERATED =4 7 10REORDERED DOF POSITIONSORIGINAL DOF ORDER =DEGREE OF FREEDOM INDICES =1 2 3 4 TITLE =NODES =23 4NODE X Y THXZ2 0. 0000 3. 0000 0. 00003 0. 0000 6. 0000 0. 00004 0. 0000 9. 0000 0. 0000EDGE PLOT DATA 二EDGE XI Y1Z21 0. 0000 0. 0000 0. 00002 0. 0000 3. 0000 0. 00003 0. 0000 6. 0000 0. 0000GLOBAL DOF SET 二33 65 97GLOBAL DOF SET NODES, LABELS2 UX3 UX4 7 10MASS INFORMATION:Z THXY THYZ 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000 0. 0000Z1 X2 Y20. 0000 0. 0000 3. 0000 0. 0000 0. 0000 6. 0000 0. 0000 0. 0000 9. 00004 UXTOTAL MASS = 7822.0CENTROID (X,Y,Z) = 0. 0000 5. 8228 MOMENT OE INTERTIA ABOUT ORIGIN:IXX = 0.31052E+06 IYY = 0. 0000IXY = 0.0000 IYZ = 0.0000 0. 00001ZZ = 0.31052E+06 IZX = 0. 0000RON 1 MATRIX440600. 00 -192100. 00 0. 0000000ROW 1 MATRIX 22762. 0000 0. 0000000 0. 0000000ROW 2 MATRIX 1-192100. 00 344300. 00 -152200. 00ROW 2 MATRIX 20. 0000000 2760. 0000 0. 0000000ROW 3 MATRIX 10. 0000000 -152200. 00 152200. 00ROW 3 MATRIX 20. 0000000 0. 0000000 2300. 0000LOAD VECTOR 10. 0000000 0. 0000000 0. 0000000。

ansys提取刚度矩阵案例(命令流及矩阵文件)

ansys提取刚度矩阵的三种方法注：本案例借鉴了王新敏老师编著的《ansys工程结构数值分析》部分内容结构刚度矩阵K(1)单元刚阵的提取-----详见page356/PREP7ET,1,BEAM3MP,EX,1,2E5R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLU/OUTPUT,ELEMSTIFF,TXT/DEBUG,-1,,,1SOLVE/OUTPUTFINISH用任一文本编辑器打开ELEMSTIFF.TXT文件可得到单元刚度矩阵ELEMSTIFF.TXT(2)用HBMAT提取原始刚度矩阵和节点载荷FINISH/CLEAR/FILNAME，HBFILE/PREP7ET,1,BEAM3MP,EX,1,2E5R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3/SOLUWRFULL,1SOLVEFINISH/AUX2FILE,HBFILE,FULLHBMAT,HBFILE,TXT,ASCII,STIFF,YESFINISH用文本编辑器打开HBFILE.TXT可看到用Harwell-Boeing格式记录的文件（仅含非零项的上三角阵，刚度矩阵为对称矩阵）。

矩阵有12行12列，33个非零元素HBFILE.TXT（3）用HBMAT提取结构刚度矩阵与（2）相同，但施加约束条件即可。

所生成的HBMAT.TXT的前5行为/PREP7ET,1,BEAM3MP,EX,1,2E5R,1,1E-2,32E-5,0.5N,1N,2,0,4N,3,4,4N,4,4,0EN,1,2,3EN,2,1,2EN,3,4,3F,2,FX,5SFBEAM,1,1,PRES,10,,,,2,-1SFBEAM,2,1,PRES,3D,1,ALLD,4,ALL/SOLUWRFULL,1SOLVEFINISH/AUX2FILE,文件名缺省,FULLHBMAT,文件名缺省,TXT,ASCII,STIFF,YESFINISH其意义同上，但数值有变化。

ansys提出刚度和质量矩阵

2762.0000 0.0000000 0.0000000
附件 1：（运行 ANSYS 命令的输出结果，最后面是刚度和质量矩阵） PRINT CONTENTS OF SUPERELEMENT submat PRINT OPTION = 3 HEADER = 8 4 1 0 0 at 922 982 940 0 419 0 3 0 1 0 1101 1078 928 1027 106004406 0 0 2 4 0 1 166 0 946 1072 0 0 0 3 3 0 12 103 953 0 0 0 0 3 0 1 1155 subm 976 934 3 0 0
一、连续钢梁的刚度和质量矩阵弹性模量：EX=2.0× 1011 Pa 泊松比：PRXY=0.3 密度：DENS=7811Kg/m3 截面特性如右图。 ANSYS 命令流： /FILNAM,BEAM MASS AND KNIFF ANALYSIS /TITLE, BEAM MASS AND KNIFF ANALYSIS /PREP7 ET,1,BEAM3 MP,EX,1,2.0E11 MP,NUXY,1,0.3 MP,DENS,1,7800 SECTYPE,1,BEAM,I,,3 SECOFFSET,CENT SECDATA,0.15,0.15,0.3,0.02,0.02,0.01,0,0,0,0 K,1,0,0,0 K,2,10,0,0
GLOBAL DOF SET = 33 65 97 GLOBAL DOF SET NODES, LABELS = 2 UX 3 UX MASS INFORMATION: TOTAL MASS = 7822.0 CENTROID (X,Y,Z) = 0.0000 5.8228 MOMENT OF INTERTIA ABOUT ORIGIN: IXX = 0.31052E+06 IYY = 0.0000 IXY = 0.0000 IYZ = 0.0000 ROW 1 MATRIX 1

ansys质量矩阵刚度矩阵提取说课材料

a n s y s质量矩阵刚度矩阵提取ansys质量矩阵刚度矩阵提取看了这么久了都没人回，查了一些质料终于找到答案了，，下面提供三种方法：方便与其他程序进行接口编程1.Which matrix you would like? element stiffness matrix or full stiffness matrix?element stiffness is within file.emat. full stiffness matrix is within file.fullA simple way to dump the matrix is as follow:-------------------/aux2fileaux2,file,ematform,longdump,all-------------------2.可以使用/DEBUG命令来得到。

ansys质量矩阵刚度矩阵提取教学提纲

ansys质量矩阵刚度矩阵提取说课材料

an s y s 质量矩阵刚度矩阵提取ansys质量矩阵刚度矩阵提取看了这么久了都没人回，查了一些质料终于找到答案了，，下面提供三种方法：方便与其他程序进行接口编程 1.Which matrix you would like? element stiffness matrix or full stiffness matrix?element stiffness is within file.emat. full stiffness matrix is within file.fullA simple way to dump the matrix is as follow: /aux2 fileaux2,file,emat form,l ong dump,all2.可以使用/DEBUG命令来得到。

详细步骤参见下面的宏文件fin ish/clearPl=3.1415926w1=3w2=10w3=6w4=1.2r=.8t=0.08/PREP7!*ET,1,SHELL63R,1,tET,2,MASS21R,2,500,500,500,2000,2000,2000,!*UIMP,1,EX, , ,2e11UIMP,1,NUXY, , ,0.3,UIMP,1,DAMP, , ,0.2,UIMP,1,DENS, , ,7800,BLC4,0,0,w2,w1ESIZE,1.5,0,AMESH,allNSEL,S,LOC,X,0.0D,all, , , , , ,ALL,,,,,allsel,allSFA,all,1,PRES,12FINISH/OUTPUT,cp,out,, !将输出信息送到cp.out文件/debug,-1,,,1 !指定输出单元矩阵/SOLUSOLVEfin ish/OUTPUT, TERM !将输出信息送到output windows 中!这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵3.其原理很简单，即使用an sys的超单元即可解决问题。

提取刚度矩阵

1在ANSYS建立有限元模型2建立好模型后，求解时，采用下面的命令流/soluallselantype,staticeqslv,sparsensubst,1wrfull,1ematwrite,yessolvefini求解后，得到的文件包括.emat .full 文件将这两个文件拷贝到mor4ansys.exe 所在文件夹下mor4ansys.exe的路径为：本机D:\mor4ansys\mor4ansys\mor4ansys-1.8-windows\mor4ansys.exe在DOS系统中进入D:\mor4ansys\mor4ansys\mor4ansys-1.8-windows 文件夹在文件夹1．进入D盘D：回车2 进入D:\mor4ansys\mor4ansys\mor4ansys-1.8-windows 文件夹在文件夹键入cd D:\mor4ansys\mor4ansys\mor4ansys-1.8-windows并回车然后键入：mor4ansys -f 文件名.full 文件名.emat –w例如我用ANSYS提出的文件分别为即文件名为sphere然后键入：mor4ansys -f sphere.full sphere.emat –w 回车即可得出文件然后将mminfo mmread mmwrite MATLAB_READ TTFORM 等文件拷贝到该文件夹打开matlab执行MATLAB_READ命令得出K.mat M.mat文件在打开TTFORM.m文件编辑for k=1:2316TT(fullmodel(k,1),fullmodel(k,2))=fullmodel(k,3);End其中2316为结构的自由度总数，可根据实际情况修改修改后保存，并运行，得出TT变量导入K.mat M.mat文件进行运算K1=T*M*T’M1=T*M*T’得出的K1，M1即为刚度，质量矩阵。

ansys 提取质量刚度阻尼矩阵

ansys 提取质量刚度阻尼矩阵
在ANSYS中，可以使用以下步骤来提取质量、刚度、阻尼矩阵：
打开ANSYS软件，并导入需要分析的模型。

在主界面中，选择“工具”菜单中的“质量特性”，打开“质量特性”对话框。

在“质量特性”对话框中，选择“全局质量矩阵”，然后单击“计算”按钮。

在计算完成后，将显示全局质量矩阵、全局刚度矩阵和全局阻尼矩阵。

可以使用类似的方法提取其他相关矩阵，例如局部质量矩阵、局部刚度矩阵和局部阻尼矩阵。

需要注意的是，在提取矩阵之前，需要对模型进行正确的材料属性设置和边界条件设置。

此外，对于复杂模型，可能需要使用其他工具或方法来提取矩阵。

NO2 如何提取Ansys整体刚度阵及转换成Matlab数据

3. 利用 matlab 程序“hb_to_msm.m”即可得到 matlab 稀疏矩阵形式的 MCK 阵，程序见附录2
function [a, rhsval] = hb_to_msm (input_file) M、C、K：稀疏矩阵； rhsval：载荷右端项
input_file：步骤2导出的数据文件
6. 利用 matlab 程序“readAnsysSort2.m”即可读取 EMAT 文件，程序见附录 3
function NodeDofs = readAnsysSort2(file1,file2) file1：通过5导出的EMAT文件，即“DUMP.lis” file2：模型约束文件
节点号约束自由度
4. 该方法适用于大型结构，复杂模型。 5. 该方法生成的整体刚度阵，是按优化后节点序列排序的，需要找到各节点在刚度阵中的
对应的位置，此处通过 EMAT 文件获取，其中 EMAT 文件的格式说明见附录。
FINISH /AUX2 FORM,LONG FILEAUX2,'file','emat',' ' DUMP,5, FINISH
FINISH /AUX2 FILE,'file','full',' ' HBMAT,'k','txt',' ',ASCII,STIFF,YES FINISH
HBMAT,'m','txt',' ',ASCII,MASS,YES HBMAT,'c','txt',' ',ASCII,DAMP,NO YES-列出载荷右端项；NO-不列出载荷右端项；

2016新编ANSYS中导出刚度矩阵

[转贴]ansys中导出刚度矩阵ansys中导出刚度矩阵可以使用/DEBUG命令来得到。

详细步骤参见下面的宏文件finish/clearPI=3.1415926w1=3w2=10w3=6w4=1.2r=.8t=0.08/PREP7!*ET,1,SHELL63R,1,tET,2,MASS21R,2,500,500,500,2000,2000,2000,!*UIMP,1,EX, , ,2e11UIMP,1,NUXY, , ,0.3,UIMP,1,DAMP, , ,0.2,UIMP,1,DENS, , ,7800,BLC4,0,0,w2,w1ESIZE,1.5,0,AMESH,allNSEL,S,LOC,X,0.0D,all, , , , , ,ALL, , , , ,allsel,allSFA,all,1,PRES,12FINISH/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中! 这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵干部教育培训工作总结[干部教育培训工作总结] 年干部教育培训工作，在县委的正确领导下，根据市委组织部提出的任务和要求，结合我县实际，以兴起学习贯彻“三个代表”重要思想新高潮为重点，全面启动“大教育、大培训”工作，取得了一定的成效，干部教育培训工作总结。

现总结报告如下：一、基本情况全县共有干部**人，其中中共党员**人，大学本科以上学历**人，大专学历**人，中专学历**人，高中及以下学历**人。

**年，以县委党校、县行政学校为主阵地，举办各类培训**期，培训在职干部**人，占在职干部总数的**.*%，培训农村党员、干部**人，其中：举办科级领导干部轮训班*期，培训**人；举办科级领导干部“三个代表”重要思想专题学习班*期，培训**人；举办科级以下公务员培训班*期，培训**人；举办企业经营管理者培训班*期，培训**人；举办专业技术人员培训班*期，培训**人；举办非中共党员干部培训班*期，培训**人；举办理论骨干培训班*期，培训**人；举办妇女干部培训班*期，培训**人；举办基层团干培训班*期，培训**人；举办农村党支部书记、村主任培训班各*期，培训**人，达到了每年培训在职干部五分之一的要求，超额完成了培训任务。

ANSYS中整体、单元刚度和质量矩阵的提取

ANSYS中整体、单元刚度和质量矩阵的提取1、整体刚度和质量矩阵的提取。

1.1、用户程序法，需要二次开发该功能需要进行二次开发，由ansys形成的二进制文件.full提取整体刚度和质量矩阵。

基于ansys的一个用户开发程序例子编了一个程序（附件中）。

开发环境：compaqfortran6.5运行环境：win2000。

一个主文件self.for；另一个文件matrixout.f90用于矩阵输出；binlib.lib 为ansys提供的库文件，将其引入项目中（也可直接扔进debug目录）；.full文件由子空间迭代模态分析获得。

运行编译后的可执行文件.exe其他文件见/f/EE24A2F87F524606.html1.2、超单元法其原理很简单，即使用ansys的超单元即可解决问题。

定义超单元，然后列出超单元的刚度矩阵即可。

下面是一个小例题，自可明白。

/prep7k,1k,2,3000l,1,2et,1,beam3mp,ex,1,2e5mp,prxy,1,0.3r,1,5000,2e7,200lesize,all,,,10lmesh,allfinish!----以上正常建立模型，不必施加约束和荷载/soluantype,7 !substructuring分析类型seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵1.3、HBMAT命令法提取整体矩阵13.1、命令说明命令：HBMAT,fname,ext,--,form,matrx,rhs其中：Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。

ext---输出矩阵文件的扩展名，缺省为.matrix。

form---定义输出矩阵文件的格式，其值可取：=ASCII：ASCII码格式；=BIN：二进制格式。

Ansys中单元刚度矩阵总刚矩阵的提取

Ansys中单元刚度矩阵、总刚矩阵的提取一、单元刚度矩阵的提取/DEBUG命令详细说明：finish/clearPI=3.1415926w1=3w2=10w3=6w4=1.2r=.8t=0.08/PREP7!*ET,1,SHELL63R,1,tET,2,MASS21R,2,500,500,500,2000,2000,2000,!*UIMP,1,EX, , ,2e11UIMP,1,NUXY, , ,0.3,UIMP,1,DAMP, , ,0.2,UIMP,1,DENS, , ,7800,BLC4,0,0,w2,w1ESIZE,1.5,0,AMESH,allNSEL,S,LOC,X,0.0D,all, , , , , ,ALL, , , , ,allsel,allSFA,all,1,PRES,12FINISH/OUTPUT,cp,out,, ! 将输出信息送到cp.out文件/debug,-1,,,1 ! 指定输出单元矩阵/SOLUSOLVEfinish/OUTPUT, TERM ! 将输出信息送到output windows中这时用编辑器打开cp.out文件，可以看到按单元写出的质量、刚度等矩阵二、整体刚度矩阵的提取（有三种方法：用户程序法、超单元法、HBMAT命令法）1、用户程序法：需要二次开发（略）2、超单元法/soluantype,7 !substructuring分析类型seopt,matname,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)nsel,all !选择所有节点m,all,all !定义所有节点自由度为主自由度solve !求解selist,matname,3 !列出整体刚度矩阵3、HBMAT命令法提取整体矩阵命令：HBMAT,fname,ext,--,form,matrx,rhs其中：Fname---输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。

ext---输出矩阵文件的扩展名，缺省为.matrix。

提取刚度矩阵

============为什么折腾这个文档========我有一个计算线性动力学方程组的瞬态、谐响应和静力学的python程序，现希望开发一个将ANSYS组集好的总体矩阵导入该PYTHON程序中的接口。

该问题可分解为：[STEP1] [ANSYS]->[包含矩阵信息的文件][STEP2] [包含矩阵信息的文件]->[python通用数据对象][STEP3] [python通用数据对象]->[程序特定数据对象]->[进行计算]因此检索了一些帖子，基本上完成了这项工作，本文是对[STEP1]和[STEP2]的整理，并且利用[STEP3]对结果进行了验证============主要内容==================1，了解从ANSYS中提取总体矩阵和载荷向量的方法；2，了解提取出来的矩阵是怎样表示的；3，说明在Python中，如何读取这样的矩阵；4，构造一个简单的算例，说明整个【建模】-【提取】-【读取】过程及其正确性；=========站内检索综述====================检索词：提取矩阵得到21个结果，代表性的帖子有下面这9个：编号[1]标题：ansys中怎样提取质量，刚度，阻尼矩阵？地址：/forum-vi ... fromuid-159019.html 要点：pengweicai给出了一段网上最常见的提取代码,该程序以fortran 写成，可以利用.full文件以及一些列约定将ANSYS中的总体矩阵读入FORTRAN中。

编号[2]标题：如何得知HBMAT命令提取的质量、刚度矩阵对应的自由度？地址：/forum-vi ... fromuid-159019.html 要点：提出了使用HBMAT命令提取稀疏矩阵时常见的问题：我们如何知道提取出来的信息是怎么储存的呢？编号[3]标题：[分享]ANSYS中整体、单元刚度和质量矩阵的提取地址：/forum-vi ... fromuid-159019.html 要点：在该帖子的7楼，其实已经给出了帖子[2]中问题的解答，即HBMAT 中提取出来的矩阵是Harwell-Boeing格式的，并且给出了该格式的细节，可惜是英文的，没引起多少关注。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

*deck,rdfull USERDISTRIB ansys
program rdfull
c *** primary function: demonstrates use of binary access routines
c
c ********** define the file names **********定义文件名
pname = 'file.full' !ANSYS的结果文件，存储了整体刚度矩阵和质量刚度矩阵
mname = 'MASS.MATRIX' !用来存放质量刚度矩阵信息的文件的名字
external binset, binini, binrd8, binclo, biniqr, bintfo
!打开二进制文件，初始化输入输出设备，读取二进制文件中的记录
!关闭打开的二进制文件，二进制文件或系统的查询，定义二进制文件头的必要信息
integer binset, biniqr
c ********** define the i/0 buffer page length (integer*4 words)
reclng = biniqr (0,1)
!返回值为16384，表示标准输入输出缓冲区域页的长度为16384个字符=16kb=系统块的大小
c most analyses. This program demonstrates how to read 本程序展示了怎样去对.FULL文件进行读取和改变格式
c and reformat the .FULL file that was created using 后写入到新的文件，其中这里的.FULL文件可由symbolic
units = 0 !新定义的文件头记录中保存的单位制
c --- the value for code below should be changed to appropriate 3rd party code
code = 200 !定义第三方供应商的代码
c *** secondary function: Read and reformat full file
c
c *** copyright(c) 2009 SAS IP, Inc. All rights reserved.
c *** ansys, inc.
c
c NOTICE - A new assembly process, termed 'symbolic assembly', has 一个被称为symbolic assembly的新的组装过程取代了
c (nmode is not used - it can be any value) nmode没有被使用，可以是任何值
c If the free-free stiffness and mass matrices are desired, 如果想获得无约束的刚度矩阵或质量矩阵，则
title(1) = 'New title as given by ROM' !新定义的文件头记录中保存的主标题
title(2) = 'New subtitle as given by ROM' !新定义的文件头记录中保存的副标题
jobnam = 'file ' !新定义的文件头记录中保存的工作名
c replaced the old assembly process, termed 'frontal 原来被称为 frontal assembly的旧的组装过程，这个
c assembly', and is now the default assembly process for 新的组装过程是现在大多数分析默认的组装过程。
c sequence is used. ANSYS will also write the .FULL file 写出.FULL文件。
c when the sparse, ICCG, or JCG solver is used, as well as
integer iarray(ROWLENG)
double precision darray(ROWLENG)
equivalence (darray(1),iarray(1))
c ********** define the unit numbers **********定义设备编号
c
c ********** initialize the bin routines **********
call binini (iout) !初始化输入输出设备
c ********** initialize the header common **********
kname = 'STIFFNESS.MATRIX' !用来存放整体刚度矩阵信息的文件的名字
write (iout,2000)
2000 format (/' ***** WRITE OUT ANSYS MATRICES FROM file.full',
x ' *****'//' MASS MATRIX ON FILE = mass.matrix'/
x irow,node,kdof,lll(ROWLENG),
x mr,nterms,indx(ROWLENG),l(WAVEFRONT)
integer buffer(IOLENG)
double precision krow(ROWLENG), mrow(MAXEQN)
x ' STIFFNESS MATRIX ON FILE = stiffness.matrix'//
x ' Only the symmetric part of the matrices is written')
!仅仅写出矩阵的对称部分
c make sure there are no constraints on the model. 确保没有对模型施加约束。
#include "impcom.inc"
external ihsort !将数据表中的数据按第某一行中的数据从小到大排列
iout = 6 !表示输入输出设备为显示器
kunit = 2 !新建文件unit=kunit=2，用来存放整体刚度矩阵信息
munit = 3 !新建文件unit=munit=3，用来存放质量刚度矩阵信息
c ********** define the number of integers per double precision
c frontal assembly or symbolic assembly. ANSYS writes the assembly组装过程或frontal assembly组装过程产生。
c .FULL file if the PSOLVE,ELFORM ANSYS将在PSOLVE,ELFORM;PSOLVE,ELPREP;PSOLVE,TRIANG三组命令被依次执行
integer IOLENG, ROWLENG, WAVEFRONT, MAXNODE, MAXEQN
parameter (IOLENG=16384, ROWLENG=1000, WAVEFRONT=10000,
x MAXNODE=100000, MAXEQN=500000)
LONGINT jloc
integer kunit,munit,units,code,nunit,nbuf,lbuf,npage,keyrw,j,
x kext,i,n,ivect(100),kbf,assemb,numdof,lenbac,nontp,
c when most mode extraction methods are used.
c Be sure to set up for modal ANTYPE,2 ？？？？？？
c and Block Lanczos MODOPT,LANB,nmode,0,0, ,OFF ？？？？？？
call bintfo (title,jobnam,units,code)!为子程序binhed8写一个新的文件头记录提供必要的信息
c subroutine bintfo (title,jobnam,units,code)
c *** primary function: set information necessary for binhed
intpdp = biniqr (0,2) !返回值=2，表示每个双精度数所占的空间等于2个整数所占的空间
c ********** define the number of characters in the file name
lenfnm = biniqr (0,3) !返回值=260，表示所打开的二进制文件的名字最多占用260个字节
x nmatrx,lumpm,baclst(MAXNODE),nmax,sortlist(2,MAXNODE),
x nmass,nstif,numNodes,nrow,nDofEachNode(MAXEQN),
x dof(MAXEQN),idof,prevDof,currDof,dofMap(MAXEQN),
c 为子程序binhed8写一个新的文件头记录提供必要的信息
c --- This routine is intended to be used in standalone programs.
character*80 title(2)
character*106 pname,mname,kname