本帖最后由 骨灰级浪子 于 2011-5-10 10:49 编辑
FINISH
/CLEAR
/FILNAME,expansive analysis
/TITLE,expansive analysis
/PLOPTS,DATE,0
/PREP7
/COM,THERMAL
ET,1,PLANE55
MP,KXX,1,1E-5 ! PERMEABILITY
MP,KXX,2,1E-7
MP,DENS,1,2277
MP,DENS,2,2500
MP,C,1,0.0015
MP,C,2,0.0025
K,1,0,0
K,2,8.2,0
K,3,0,2
K,4,4,2
K,5,4.2,2
K,6,8.2,2
K,7,0,3.2
K,8,4,3.2
K,9,4.2,3.2
K,10,8.2,3.2
L,1,2
L,2,6
L,1,3
L,3,4
L,4,5
L,5,6
L,4,8
L,3,7
L,7,8
L,5,9
L,9,10
L,10,6
L,8,9
LESIZE,ALL,0.1
AL,1,2,3,4,5,6
AL,4,7,9,8
AL,6,10,11,12
AL,5,7,13,10
AGLUE,1,2
AGLUE,1,3
AGLUE,1,4
AGLUE,4,2
AGLUE,4,3
LSEL,S,LINE,,4
LSEL,A,LINE,,5
LSEL,A,LINE,,6
LCCAT,ALL
ALLSEL,ALL
/REPLOT
MSHAPE,0,2D
MSHKEY,1
ASEL,,,,4
AATT,1,,1
AMESH,ALL ! MESH AREAS
ALLSEL,ALL
ASEL,,,,2
ASEL,A,,,3
AATT,2,,1
AMESH,ALL ! MESH AREAS
ALLSEL,ALL
ASEL,,,,1
AATT,2,,1
AMESH,ALL ! MESH AREAS
ALLSEL,ALL
NUMMRG,NODE
NUMMRG,ELEM
NUMMRG,KP
LSEL,S,LINE,,13
NSLL,S,1
D,ALL,TEMP,0.1718
ALLSEL,ALL
LSEL,S,LINE,,11
LSEL,A,LINE,,9
NSLL,S,1
D,ALL,TEMP,0
ALLSEL,ALL
LSEL,S,,,4
LSEL,A,LINE,,5
LSEL,A,LINE,,6
NSLL,S,1
D,ALL,TEMP,0.2118
ALLSEL,ALL
FINISH
/SOLU
ANTYPE,STATIC
SOLVE
SAVE
/POST1
/DEVICE,VECTOR,0
PLNSOL,TEMP
SAVE
ALLSEL,ALL
/REPLOT
/COM,Structural
/PREP7
ETCHG,TTS !热-结构单元转换
MP,ALPX,1,0.4784
MP,ALPY,1,0.0797
MP,ALPX,2,0
MP,ALPY,2,0
MP,EX,1,1.0E9
MP,NUXY,1,0.3
MP,EX,2,20.2E9
MP,NUXY,2,0.26
TB,DP,1
TBDATA,1,0.16e6,30,30
TB,DP,2
TBDATA,1,0.34e6,36.4,36.4
ALLSEL,ALL
/DEVICE,VECTOR,0
/REPLOT
MSHAPE,0,2D
MSHKEY,1
LSEL,,,,2
LSEL,A,,,3
LSEL,A,,,8
LSEL,A,,,12
NSLL,S,1
D,ALL,UX,0
ALLSEL,ALL
LSEL,S,,,1
NSLL,S,1
D,ALL,ALL,0
ALLSEL,ALL
/PSF,PRES,NORM,2,0,1
ACEL,0,9.8,0
FINISH
/SOLU
TREF,0
LDREAD,TEMP,,,,,,rth
/PBF,TEMP,,1
ANTYPE,STATIC
NLGEOM,ON
TIME,360
AUTOTS,ON
NSUBST,60,500,10
OUTRES,ALL,ALL
SOLVE
SAVE
/POST1
!/VIEW,1,-1,1,1
SET,LAST ! 读入最后荷载步的计算结果
PLDISP,2 ! 显示变形图(保留原结构轮廓)
PLNSOL,S,X ! 在图形窗口显示节点上的X方向正应力
PLNSOL,S,Y ! 在图形窗口显示节点上的Y方向正应力
PLNSOL,S,XY ! 在图形窗口显示节点上的XY方向的剪应力
PLNSOL,S,XZ
PLNSOL,S,YZ
PLNSOL,S,1
PLNSOL,S,2
PLNSOL,S,3
PLNSOL,S,EQV
上述命令流为热应力耦合 间接法求解
应力求解时提示 Solution not converged at time 0.72 run terminated 计算收敛不了 请指教! |