众所周知，MATLAB是matrix&laboratory的组合，翻译成中文也就是矩阵实验室。它的基本数据单位是矩阵。在模型转化这个项目中，我们主要运用了它的文本读入、文本输出以及矩阵变换的功能。

目前业内运用较多的设计软件有PKPM、YJK、MIDAS、ETABS、SAP2000等，其中MIDAS/GEN、ETABS和SAP2000都支持形成直接可读的模型信息文本。在这些文本中，我们可以找到模型的所有信息。我们指的模型信息包括：点信息、单元信息、截面信息、单元与截面的映射信息、单元与材料的映射信息、单元方向信息、材料信息、荷载信息等。下面例举在MIDAS/GEN中部分信息。

点信息

单元关于截面、材料映射、方向的信息

材料信息

板厚信息(属截面信息)

截面信息

支座信息

首先我们需要了解这些数据的意思，在文本的台头均有数据说明。接下来我们需要做的就是把这些数据装入MATLAB的矩阵之中。由于受到个人MATLAB技术限制，我把各个信息做成单个文本，便于后面的操作。通过函数fopen()读入文本，利用函数fgetl()、fscanf()以及一些循环语句完成数据信息的矩阵化。然后按照需要对矩阵进行控制输出。输出的时候主要用到函数strcmp()、正则表达式regexpi()、fprintf()等。

下面附一截面信息读取例子：

文本名为 beamsection.txt。以下为文本内容

56 B-0.80x0.80 B 0.8 0.8 0.04 0.04

57 B-0.40x0.40 B 0.4 0.4 0.016 0.016

501 幕墙转换桁架_方800x600x40x40 B 0.6 0.8 0.04 0.04

502 幕墙转换桁架_方400x700x30x30 B 0.4 0.7 0.03 0.03

544 屋盖格构柱_方350x12 B 0.35 0.35

0.012 0.012

601 屋盖封边_方250*8 B 0.25 0.25

0.008 0.008

608 屋盖桁架弦杆200x8(格构柱) B 0.2 0.2 0.008 0.008

651 屋面桁架腹杆_方120*4 B 0.12 0.12

0.004 0.004

653 屋盖桁架腹杆_方180x6 B 0.18 0.18

0.006 0.006

654 屋盖桁架腹杆_方80x4 B 0.08 0.08

0.004 0.004

656 屋盖桁架腹杆_方180x6 B 0.18 0.18

0.006 0.006

1005 封边次梁口600X200X10X18 B 0.6 0.2 0.01 0.018

1008 封边次梁口600X300X10X20 B 0.6 0.3 0.01 0.02

1012 斜梁口600X600X20X40 B 0.6 0.6 0.02 0.04

1007 檩条C200X70X20X3 C型钢 0.2 0.07

0.003 0 0.02

41 H-0.01x0.65 H 0.65 0.3 0.012 0.016

42 H-0.01x0.50 H 0.5 0.25 0.01 0.014

43 H-0.01x0.30 H 0.3 0.2 0.008 0.012

1002 柱顶纵梁H700X300X12X20 H 0.7 0.3 0.012 0.02

1006 二级次梁H450X200X8X18 H 0.45 0.2 0.008 0.018

1011 挑梁H600X550X14X32 H 0.6 0.55

0.014 0.032

606 屋盖斜拉杆_圆89x3

P 0.089 0.003

652 屋盖桁架腹杆_圆200x8 P 0.2 0.008

655 屋盖桁架腹杆_圆168x5 P 0.168 0.005

1 SB-0.30x0.80 SB 0.8 0.3

2 SB-0.30x0.90 SB

0.9 0.3

3 SB-0.40x1.20 SB

1.2 0.4

44 SB-0.10x0.10 SB

0.1 0.1

45 SB-1.50x1.50 SB

1.5 1.5

47 SB-0.80x0.80 SB

0.8 0.8

54 SB-1.30x1.30 SB

1.3 1.3

50 SR-1.00 SR

1

46 RH2T-1.50x1.50 RH2T 1.5 0.9 0.5 0.04

48 RH2T-1.30x1.30 RH2T 1.3 0.9 0.5 0.04

53 RH2T-1.50x1.50 RH2T 1.5 1.2 0.5 0.03

55 RH2T-1.30x1.30 RH2T 1.3 0.9 0.5 0.04

543 屋盖柱_圆800x25 EPC 0.8 0.025

END

MATLAB文本如下所示：

fid1=fopen('beamsection.txt','r');

stemp1='END';

stemp='START';

itter=0;

%定义元组beaminf{},第一列为梁截面编号，第二列为梁截面名称，第三列为截面参数信息

while ~strcmp(stemp,stemp1)

itter=itter+1;

beaminf{itter,4}=[];

stemp=fgetl(fid1);

pat='S*?(?=s)|(?<=s)S*+';

a=regexpi(stemp,pat,'match');

numba=size(a,2);

for i=1:numba

if i==1

beaminf{itter,1}=a(i);

elseif i==2

beaminf{itter,2}=a(i);

elseif i==3

beaminf{itter,3}=a(i);

else

beaminf{itter,4}=[beaminf{itter,4},a(i)];

end

end

end

fclose(fid1);

%开始读取壳截面参数信息

fid1=fopen('shellsection.txt','r');

shellinf=fscanf(fid1,'%g' ,[2,inf]);

shellinf=shellinf';