/prep7
ET,1,LINK10
!找形時(shí)采用很小的彈性模量
R,1,0.001468,0.9,
MP,EX,1,0.60551e9
MP,NUXY,1,0.3
!定義節(jié)點(diǎn)
*do,i,1,5
*do,j,1,2*i-1
k=(i-1)*(i-1)+j
n,k,(j-i)*9.15,-36.6+(i-1)*9.15,0
*enddo
*enddo
*do,i,1,4
*do,j,1,2*i-1
k=(i-1)*(i-1)+j+25
n,k,(j-i)*9.15,36.6-(i-1)*9.15,0
*enddo
*enddo
!定義單元
*do,i,1,4
*do,j,1,2*i
k=i*i+j
e,k,k+1
*enddo
*enddo
*do,i,1,3
*do,j,1,2*i
k=i*i+j+25
e,k,k+1
*enddo
*enddo
*do,i,1,4
*do,j,1,2*i-1
k=(i-1)*(i-1)+j
e,k,k+2*i
*enddo
*enddo
*do,i,1,3
*do,j,1,2*i-1
k=(i-1)*(i-1)+j+25
e,k,k+2*i
*enddo
*enddo
*do,i,18,24
e,i,i+17
*enddo
!施加位移約束
d,1,ux,0
d,1,uy,0
d,1,uz,3.66
d,2,ux,0
d,2,uy,0
d,2,uz,1.83
d,4,ux,0
d,4,uy,0
d,4,uz,1.83
d,5,ux,0
d,5,uy,0
d,5,uz,0
d,9,ux,0
d,9,uy,0
d,9,uz,0
d,10,ux,0
d,10,uy,0
d,10,uz,-1.83
d,16,ux,0
d,16,uy,0
d,16,uz,-1.83
d,17,ux,0
d,17,uy,0
d,17,uz,-3.66
d,25,ux,0
d,25,uy,0
d,25,uz,-3.66
d,35,ux,0
d,35,uy,0
d,35,uz,-1.83
d,41,ux,0
d,41,uy,0
d,41,uz,-1.83
d,30,ux,0
d,30,uy,0
d,30,uz,0
d,34,ux,0
d,34,uy,0
d,34,uz,0
d,27,ux,0
d,27,uy,0
d,27,uz,1.83
d,29,ux,0
d,29,uy,0
d,29,uz,1.83
d,26,ux,0
d,26,uy,0
d,26,uz,3.66
!設(shè)定時(shí)間步
TIME,1
AUTOTS,0
NSUBST,20, , ,1
KBC,0
!輸出每個(gè)子步的結(jié)果
OUTPR,BASIC,ALL,
OUTRES,ALL,ALL,
!設(shè)定大變形，應(yīng)力剛化
NLGEOM,1
SSTIF,ON
FINISH
!求截
/SOLU
/STATUS,SOLU
SOLVE
FINISH
!后處理
/POST1
PRNSOL,DOF,
PRESOL,SMISC,1
PLESOL,SMISC,1
PLNSOL,U,X,0,1
finish
/PREP7
!更新幾何形狀
!a.rst為計(jì)算結(jié)果文件名，最后一個(gè)為目錄
!這兩個(gè)參數(shù)應(yīng)根據(jù)你的計(jì)算情況定
UPGEOM,1,LAST,LAST,a,rst,E:\JZD\1\
!彈性模量恢復(fù)為真值
R,1,0.001468,0.0027248,
MP,EX,1,2.0e11
!重新施加位移約束
d,1,uz,0
d,2,uz,0
d,4,uz,0
d,5,uz,0
d,9,uz,0
d,10,uz,0
d,16,uz,0
d,17,uz,0
d,25,uz,0
d,35,uz,0
d,41,uz,0
d,30,uz,0
d,34,uz,0
d,27,uz,0
d,29,uz,0
d,26,uz,0
!求截
/SOLU
/STATUS,SOLU
SOLVE
FINISH
/POST1
!觀察施加荷載之前結(jié)構(gòu)的位移
!可見(jiàn)位移基本為零,預(yù)張力基本不變
PLNSOL,U,Z,0,1
PLESOL,SMISC,1
FINISH
!寫(xiě)荷載工況文件01
/PREP7
LSWRITE,01,
!施加節(jié)點(diǎn)荷載
*do,i,1,41
f,i,fz,-167445
*enddo
!設(shè)置第2荷載步
TIME,2
AUTOTS,0
NSUBST,20, , ,1
KBC,0
LSWRITE,02,
!求解
FINISH
/SOLU
LSSOLVE,1,2,1,
后處理
/POST26
!可以看到，該結(jié)果和書(shū)中結(jié)果一致
NSOL,2,21,U,Z,
PLVAR,2, , , , , , , , , ,
ESOL,3,33,1,F,Y,
PLVAR,3, , , , , , , , , ,