ANSYS 提供了VRead 命令.docVIP

ANSYS 提供了*VRead 命令，可以从文本文件读取数据到数组中。 *Vread 命令的格式如下： *VREAD, ParR, Fname, Ext, Dir, Label, n1, n2, n3, NSKIP-- Reads data and produces an array parameter vector or matrix. ParR: 事先定义好的用于保存数据的数组，数组可以是 1维或 2 维或者 3 维的； Fname: 要读取的数据文件的文件名； Ext:数据文件名的后缀； Dir:数据文件所在的目录；（目录中最好别包含中文和空格以及太长的目录名） Label:指定读取方式,一般用 IJK,JIK 这两种方式； n1,n2,,n3：Read as (((ParR(i,j,k), k = 1,n1), i = 1, n2), j = 1, n3) for Label = KIJ. n2 and n3 default to 1.(具体什么意思，没 想明白) NSKIP :指定从文件的哪一行开始读取数据，默认 NSKIP=0,从第一行读取； 根据命令中的 Label 参数常用的 2 种读取方式来讨论。 1. Label=JIK 这是最常用的读取方式。如果希望数组中的数据的排列方式和数据在文件中的 排列方式保持一致，就应该用这种读取方式。 比如数据文件 data.txt 中的数据的排列方式如下： 1.1,1.2,1.3,1.4 2.1,2.2,2.3,2.4 3.1,3.2,3.3,3.4 4.1,4.2,4.3,4.4 如果希望读取数据后，数组 A 中的数据按照相同的方式排列即： A(1,1)=1.1 A(1,2)=1.2 A(1,3)=1.3 A(1,4)=1.4 A(2,1)=2.1 A(2,2)=2.2 A(2,3)=2.3 A(2,4)=2.4 A(3,1)=3.1 A(3,2)=3.2 A(3,3)=3.3 A(3,4)=3.4 就应该用 Label=JIK 的方式读取。 一般而言，假定数组是A(m,n)，则命令流应该为： …….. *Vread,ParR,Fname,Ext,Dir,JIK,n,m (nFw.d) ……….. 需要注意的地方： 1．矩阵是A(m,n)，命令行中是”n,m”,而不是”m,n” 2．格式说明必须为 nFw.d,不能为 mFw.d, n,w,d 在具体使用时，用具体的数字 代替； 现以上面的数据为例，下面是对应的命令流：（在 ANSYS5.7 中验证通过） finish /clear /prep7 *DIM, A,array,3,4 *vread,A(1,1),data,txt,g:\ansysuser,JIK,4,3 !注意顺序是 4,3，而不是 3,4,和数组的定义相反； (4F6.2) !注意是 4Fw.d,而不能用 3Fw.d 应该根据自己的实际情况修改上面的命令流中的文件名和存放路径； 执行命令流后，可以得到上面所预期的结果。 2. Label=IJK 这种情况，读取数据的规律非常复杂，读取数据的方式会受到到说明格式 中的 xFw.d 中的 x 的值的影响，还会收到文件中数据的排列方式的影响； 这里只说明一种常用的情况： 如果希望数据在数组中的排列方式和数据在文件中的排列方式正好是矩阵 转置关系的时候，即数组为 mXn, 数据为 nXm,应该用下面固定的格式： 假定数组为 A(m,n); 则命令格式为： *Vread,ParR,Fname,Ext,Dir,IJK,m,n !要注意是”m,n”，和数组的定义一 致； (mFw.d) !格式说明是 mFw.d，不能是 nFw.d; 注意：文件中数据的排列必须是 nXm，即数据共 n 行，每行 m 个数值，正好是矩 阵的转置； 例子：数据文件中数据的排列方式为4X3, 定义数组为A(3,4) 1.1,1.2,1.3 2.1,2.2,2.3 3.1,3.2,3.3 4.1,4.2,4.3 执行下面的命令流之后，矩阵 A 的值刚好是上面的矩阵的转置，即： A(1,1