利用Fortran讀取并修改ANSYS的APDL
2016-08-30 by:CAE仿真在線 來源:互聯(lián)網(wǎng)
利用Fortran讀取并修改ANSYS的APDL
APDL是ANSYS的命令文件,可以執(zhí)行ANSYS分析。利用Fortran可以實現(xiàn)APDL的參數(shù)化修改并且調(diào)用ANSYS的批處理bat進行有限元分析。本次采用Fortran對APDL進行讀取并執(zhí)行修改。
原始APDL命令如下:
finish
/clear
/prep7
et,1,beam3
E=2.1e11
nuxy=0.3
mp,ex,1,E
mp,prxy,1,nuxy
h=2
w=0.2
f=1000
a=10
b=15
r,1,w*h,w*h**3/12,h
k,1
k,2,a
k,3,a,-b
*do,i,1,2
l,i,i+1
*enddo
lesize,all,0.5
lmesh,all
dk,3,,,,,all
fk,1,fy,f
fk,2,fx,-f
allsel,all
/solu
Solve
相應的Fortran程序如下:
1.program main
2.implicit none
3.integer counter
4.integer::lines=30
5.character(len=30)::APDLfile(50)
6.character(len=10)::f_force="2000"
7.open(unit=10,file="E:\CAE_Tech\78&
8.&利用Fortran讀取并修改ANSYS的APDL\APDL.txt")
9.open(unit=11,file="E:\CAE_Tech\&
10.&78利用Fortran讀取并修改ANSYS的APDL\Fortran_APDL.txt")
11.do counter=1,lines,1
12. read(10,"(A30)")APDLfile(counter)
13. if(counter==11) then
14. APDLfile(counter)="f="http://f_force
15. write(11,"(1X,A30)")APDLfile(counter)
16. else
17. write(11,"(A30)")APDLfile(counter)
18. end if
19.end do
20.close(10,status="KEEP")
21.close(11,status="KEEP")
22.end program main
第二行:聲明變量需要通過聲明才能調(diào)用
第三行:定義整型變量counter
第四行:定義整型變量并且賦值,表示APDL文件的總行數(shù)
第五行:聲明字符串變量數(shù)組,字符串寬度為30,一共50組
第六行:聲明字符變量
第七、八行:打開文件夾下的APDL文件,文件代號為10
第九、十行:打開輸出后的文件,如果文件已經(jīng)存在則覆蓋,如果不存在則新建
第十一、十二行:循環(huán)讀取APDL文件的內(nèi)容,分別賦值給字符變量數(shù)組APDLfile
第十三到十八行:將讀取的APDL字符串輸出到Fortran_APDL.txt文件中,并且修改其中的第十一行
第二十、二十一行:關閉文件,釋放資源
最終生成的Fortran_APDL.txt文件內(nèi)容如下:
finish
/clear
/prep7
et,1,beam3
E=2.1e11
nuxy=0.3
mp,ex,1,E
mp,prxy,1,nuxy
h=2
w=0.2
f=2000
a=10
b=15
r,1,w*h,w*h**3/12,h
k,1
k,2,a
k,3,a,-b
*do,i,1,2
l,i,i+1
*enddo
lesize,all,0.5
lmesh,all
dk,3,,,,,all
fk,1,fy,f
fk,2,fx,-f
allsel,all
/solu
solve
紅色標記部分為可以在Fortran語句中改變變量f_force來改變。
相關標簽搜索:利用Fortran讀取并修改ANSYS的APDL Ansys有限元培訓 Ansys workbench培訓 ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術咨詢 ansys基礎知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓