我用ansys建了模型并做了有限元分析，現在需要把它用到matlab里，遇到了matlab與ansys互相調用的問題，我先假設2和7兩個單元的損傷，損傷程度分別為0.3和0.2，我把它們用matlab寫入了txt文件中，然后再用ansys調用計算這種損傷情況下的頻率，這是我的matlab程序

a=[2;7;0.3;0.2];

fid = fopen('b.txt','wt');

fprintf(fid,'%6.1f\n',a);

fclose(fid); 

下面是ansys命令流

finish

/clear

/prep7

L=0.2

H=0.1

B=0.1

blc4,,,L,H,B

ET,1,solid65 !定義單元類型 

!KEYOPT,1,1,1   !對指定單元類型設置所需的特性，與具體的單元類型有關,設置單元類型1的第1個選項的值為1     ？                              

!KEYOPT,1,5,2  

mp,ex,1,210e9

mp,prxy,1,0.3

mp,dens,1,7850

mp,ex,2,147e9

mp,prxy,2,0.3

mp,dens,2,7850

mp,ex,3,168e9

mp,prxy,3,0.3

mp,dens,3,7850

el=0.1/3

ESIZE,el  !設置單元邊的長度為e1

MSHKEY,1  !采用映射網格劃分方式

MSHAPE,0  !2D：四邊形，3D：六面體單元

VSEL,all  !選擇一組體的子集

VMESH,all!在體上生成節點和單元

nsel,s,loc,x,0$nsel,r,loc,y,0$nsel,r,loc,z,B  !選擇一組節點子集

d,all,all                                     !在節點上施加約束

nsel,a,loc,x,0$nsel,r,loc,y,H$nsel,r,loc,z,B

d,all,all 

nsel,a,loc,x,0$nsel,r,loc,y,H$nsel,r,loc,z,0

d,all,all  

nsel,a,loc,x,0$nsel,r,loc,y,0$nsel,r,loc,z,0 

d,all,all

allsel,all   !選擇所有實體

*dim,eleinfo,,4,1    !定義一個4*1的數組

*vread,eleinfo(1,1),G:\ANSYSPROGRAM\b,txt,,JIK,1,4

(f6.1)

a1=eleinfo(1,1)

a2=(eleinfo(1,1)+1)

*cfopen,E:\ANSYSPROGRAM\unrelated1\eleinfofile,txt

*Vwrite,eleinfo(1,1)

(f6.1)

*cfclos

*if,a1,eq,1,then

*do,i,1,a1

esel,s,elem,,eleinfo(i+1,1)   !選擇一組單元子集

emodif,all,mat,2

!esel,all

*enddo

*else

*do,i,1,a1

esel,s,elem,,eleinfo(i+1,1)

emodif,all,mat,2

esel,all

*enddo

*do,i,a2,4

esel,s,elem,,eleinfo(i,1)

emodif,all,mat,3

esel,all

*enddo

*endif

esel,all  

問題就出現在上面這部分，我需要讀入matlab寫出的txt文件，并對2和7號單元做單元屬性修改，急，請各位大神指教！