ansys提取刚度阵质量阵的快速方法.docxVIP

HBMAT命令：命令：HBMAT,fname,ext,--,form,matrx,rhs其中：? ? ? ? Fname输出矩阵的路径和文件名，缺省为当前工作路径和当前工作文件名。? ? ? ? ext输出矩阵文件的扩展名，缺省为.matrix。? ? ? ? form定义输出矩阵文件的格式，其值可取：? ? ? ? ? ? ? ? =ASCII：ASCII码格式；? ? ? ? ? ? ? ? =BIN：二进制格式。? ? ? ? matrix定义输出矩阵的类型，其值可取：? ? ? ? ? ? ? ? =STIFF：输出刚度矩阵。可用于写入了.FULL文件的任何类型的分析。? ? ? ? ? ? ? ? =MASS：输出质量矩阵。可用于特征值屈曲、子结构分析、模态分析。? ? ? ? ? ? ? ? =DAMP：输出阻尼矩阵。仅用于有阻尼的模态分析。? ? ? ? rhs右边项输出控制（右边项指用矩阵所表示方程的等号右端矢量，这里可为节点荷载向量），如rhs=YES则输出，如rhs=NO则不输出。? ? ? ? 模态分析时，因仅LANB和QR法可生成完整的质量矩阵，因此也仅采用这两种方法时才可使用HBMAT命令得到质量矩阵文件。Harwell-Boeing文件格式：? ? ? 用HBMAT命令可输出结构刚度矩阵、质量矩阵和阻尼矩阵，其文件记录格式为大型稀疏矩阵的标准交换格式，采用索引存储方法仅记录矩阵的非零元素。文件基本格式是前面有4或5行描述数据，其后为单列矩阵元素值，说明如下：? ? ? ? 第1行：格式（A72），为文件头的字符型解释，如刚度矩阵或质量矩阵等标题。? ? ? ? 第2行：格式（5I14），分别表示该文件的总行数（不包括文件头）、矩阵列指针的总行数、矩阵行索引的总行数、矩阵元素数值的总行数、右边项总行数。? ? ? ? 第3行：格式（A3,11X,4I14），分别为矩阵类型、矩阵行数、矩阵列数、矩阵行索引数（对组装后的矩阵，该值等于矩阵行索引数）、单元元素数（对组装后的矩阵此值为0）。? ? ? ? 第4行：格式（2A16,2A20），分别表示列指针格式、行索引格式、系数矩阵数值格式、右边项数值格式。? ? ? ? 第5行：格式（A3,11X,2I14），A3各列分别表示右边项格式、应用高斯起始矢量、应用eXact求解矢量；两个整数分别表示右边项列数、行索引数。三个字符中的第1个字符可取：F全部存贮（如节点荷载向量的全部元素）、M与系数矩阵相同方法。? ? ? ? 第6行后：矩阵元素值（单列）。? ? ? ? 矩阵类型用3个字符表示，第1个字符可取：R实数矩阵、C复数矩阵、P仅矩阵结构（无元素数值）；第2个字符可取：S对称矩阵、U不对称矩阵、HHermitian矩阵、Z病态对称矩阵；R带状矩阵；第3个字符可取：A组装的矩阵、E单元矩阵（未组装）。对称矩阵只存储下三角元素，如结构刚度矩阵为对称矩阵，Harwell-Boeing格式则仅记录下三角元素。? ? ? ? 根据Harwell-Boeing文件格式，可读取矩阵的任意行列元素的数值，也可编程还原为满矩阵存储，以便它用，很显然这种提取方式比较方便。如当生成.FULL文件后，可采用命令：/AUX2FILE,mywork,fullHBMAT,mystiff,txt,ASCII,STIFF,YESFINISH将二进制mywork.full文件输出为ASCII码文件mystiff.txt，并输出右边项。Menu PathsMain MenuGeneral PostprocData File Opts Main MenuTimeHist PostproSettingsFile Utility MenuFileListBinary Files Utility MenuListFilesBinary Files 这一步一定要先读取数据，Main MenuGeneral PostprocData File Opts不然导不出来。然后在计算的工作目录下找到mystiff.txt，OK了，只不过这是Harwell-Boeing文件格式的，还要利用程序把它转换了满矩阵存储，希望用过的高手把程序贴出来，让大家共同进步。急求把Harwell-Boeing转化为满矩阵存储的命令或程序，请大家帮帮忙。一个存储例子：Stiffness matrix from ANSYS FULL file dumped into Harwell-Boeing format 43 7 15 15 6RSA 6