下面以一個(gè)彈塑性立方塊求解為例.
基本建模
/prep7
Et,1,solid45
Mp,ex,1,3e4
Mp,nuxy,1,0.3
TB,BKIN,1,1,2,1
TBDATA,,20,10,,,,
Block,,10,,10,,10
Esize,1
Vmesh,1
Asel,s,loc,y,0
Da,all,all
!求解部分共三種方法比較
方法1:采用do循環(huán)批處理
*DO,I,1,2
/SOLU
TIME,I
nsel,s,loc,y,10
D,all,uy,-1*I
allsel,all
Solve
*ENDDO
I=1時(shí): 時(shí)間TIME從0 算到1
I=2時(shí): 時(shí)間TIME從0算到2 ,都是基于初始構(gòu)型計(jì)算的,不具有對(duì)前一步計(jì)算結(jié)果的繼承性,相當(dāng)于時(shí)間函數(shù)只有一個(gè)點(diǎn)(1,2)
最終循環(huán)計(jì)算的結(jié)果是number step= 1 ; time=2
方法2:step by step求解步
ALLSEL,ALL
/SOLU
TIME,1 !一個(gè)時(shí)間對(duì)應(yīng)一個(gè)荷載步,后一步是在前一步繼承計(jì)算的
NSEL,S,LOC,Y,10
D,all,uy,-1
allsel,all
Solve
TIME,2
NSEL,S,LOC,Y,10
D,all,uy,-2
allsel,all
Solve
該方法是在前一步計(jì)算結(jié)果的基礎(chǔ)上繼承計(jì)算的,時(shí)間函數(shù)上有兩個(gè)點(diǎn)(1,1),(2,2)
最終計(jì)算的結(jié)果是NUMBER STEP=2, TIME=2
!如果TIME=1計(jì)算順利完畢,采用重啟動(dòng)是多余的,如果中途結(jié)束,需要重啟動(dòng)把TIME=1算完
方法3:
用lswire文件的方法,設(shè)置荷載步(在線性計(jì)算時(shí)相當(dāng)于荷載工況)和時(shí)間,此方法對(duì)前一步結(jié)果有
繼承性,計(jì)算結(jié)果同方法2
方法4:
對(duì)第一種方法,要具有繼承性,采取重啟動(dòng)的方法:
/solu
RESCONTRL,DEFINE,ALL,1,1 !保存每一子步計(jì)算結(jié)果,用于重啟動(dòng)
*DO,I,1,2
*IF,I,GT,1,THEN
ANTYPE,,REST !重啟動(dòng)以后,變量I恢復(fù)到上一步了,新的荷載步需要I=I+1,重啟動(dòng)文件將變量也保存了
!I=I+1
*ENDIF
TIME,I
nsel,s,loc,y,10
D,all,uy,-1*I
allsel,all
Solve
SAVE
*ENDDO
計(jì)算結(jié)果同方法2
方法5:是對(duì)方法4的一點(diǎn)補(bǔ)充,如果不明示TIME
/solu
RESCONTRL,DEFINE,ALL,1,1 !保存每一子步計(jì)算結(jié)果,用于重啟動(dòng)
*DO,I,1,2
*IF,I,GT,1,THEN
ANTYPE,,REST,I-1,LAST,0 !重啟動(dòng)以后,變量I恢復(fù)到上一步了,新的荷載步需要I=I+1,重啟動(dòng)文件將變量也保存了
I=I+1
*ENDIF
nsel,s,loc,y,10
D,all,uy,-1*I
allsel,all
Solve
*ENDDO
計(jì)算結(jié)果,time的值自動(dòng)+1,計(jì)算結(jié)果同上
注意:1. 每個(gè)時(shí)間輸入的荷載,是荷載-時(shí)間曲線上的點(diǎn),也就是荷載水平,不是荷載增量,這一點(diǎn)與ADINA相同
2. 后繼時(shí)間點(diǎn)一定1,3要大于前面的時(shí)間點(diǎn),如果小于,則系統(tǒng)在前面時(shí)間自動(dòng)+1
3. NUMBER STEP(在線性分析中相當(dāng)于荷載工況)的數(shù)目,是通過solve來確定的,每執(zhí)行solve一次,就會(huì)有一個(gè)記錄,與時(shí)間沒什么關(guān)系
4. 關(guān)于*do循環(huán)與Fortran語法的區(qū)別
*do,i,1,3
i=i-1 !這里的變量不參與循環(huán),開始就把循環(huán)次數(shù)算完了,沒有遞歸
*enddo
!即當(dāng):i=1, i=i-1=0,此時(shí)i=0沒有返回循環(huán)結(jié)構(gòu)*DO中
! i=2, i=i-1=1
! i=3, i=i-1=2
但在FORTRAN中,這種寫法會(huì)出現(xiàn)錯(cuò)誤Error: An assignment to a DO variable within a DO body is invalid. [I],循環(huán)體中不能出現(xiàn)
相關(guān)文章
- 2021-09-08BIM技術(shù)叢書Revit軟件應(yīng)用系列Autodesk Revit族詳解 [
- 2021-09-08全國(guó)專業(yè)技術(shù)人員計(jì)算機(jī)應(yīng)用能力考試用書 AutoCAD2004
- 2021-09-08EXCEL在工作中的應(yīng)用 制表、數(shù)據(jù)處理及宏應(yīng)用PDF下載
- 2021-08-30從零開始AutoCAD 2014中文版機(jī)械制圖基礎(chǔ)培訓(xùn)教程 [李
- 2021-08-30從零開始AutoCAD 2014中文版建筑制圖基礎(chǔ)培訓(xùn)教程 [朱
- 2021-08-30電氣CAD實(shí)例教程AutoCAD 2010中文版 [左昉 等編著] 20
- 2021-08-30電影風(fēng)暴2:Maya影像實(shí)拍與三維合成攻略PDF下載
- 2021-08-30高等院校藝術(shù)設(shè)計(jì)案例教程中文版AutoCAD 建筑設(shè)計(jì)案例
- 2021-08-29環(huán)境藝術(shù)制圖AutoCAD [徐幼光 編著] 2013年P(guān)DF下載
- 2021-08-29機(jī)械A(chǔ)utoCAD 項(xiàng)目教程 第3版 [繆希偉 主編] 2012年P(guān)DF