FINISH $/CLEAR $/PREP7

R=0.04 $ L=6 $ B=1.5 $ H=0.14 $ DIS=0.06 $ D1=0.08 $ D2=0.03 $ L2=1.2
H1=0.03 $ H2=0.022      
g=4.45e3 $ q=2e3 $ Q0=(1.2*g+1.4*q)
Q0=8.5e3
FC= 19.1e6 $ FT= 1.71e6
ESIZE1= 0.2 $ ESIZE2= 0.05 $  ESIZE3= 0.025     
ET,1, MESH200,6

ET,2, SOLID65 $ KEYOPT,2,7,0   
ET,3, LINK180
MPTEMP,1,0           
MPDATA,EX,1,,3.25e10  
MPDATA,PRXY,1,,0.2   
MPDATA,DENS,1,,2500
TB,MISO,1,1,11,0   
TBTEMP,0
TBPT,,0.0002,3.25e10*2e-4 $ TBPT,,0.0004,FC* 0.42 $ TBPT,,0.0006,FC* 0.51  ! 0.36->0.42
TBPT,,0.0008,FC* 0.64 $ TBPT,,0.0010,FC* 0.75 $ TBPT,,0.0012,FC* 0.84
TBPT,,0.0014,FC* 0.91 $ TBPT,,0.0016,FC* 0.96 $ TBPT,,0.0018,FC* 0.99
TBPT,,0.0020,FC $ TBPT,,0.0033,FC            
TB,CONC,1,1,9,                                
TBTEMP,0
TBDATA,,0.5,0.95,FT,-1,,                        
MPTEMP,1,0                        
MPDATA,EX,2,,2.05e11                   
MPDATA,PRXY,2,,0.3                
MPDATA,DENS,2,,7800
TB,BISO,2,1,2,                    
TBTEMP,0
TBDATA,,1110e6,2e9,,,,                 
MPTEMP,1,0 $ UIMP,2,REFT,,,0 $ MPDATA,ALPX,2,,1.2e-5
       
PI= ACOS( - 1)                        ! 定义PI
R,1,0.25*PI* 0.005*0.005              
R,2,0.25*PI* 0.010*0.010
R,3,0.25*PI* 0.005*0.005 *7

BLC4,,,B*4,H,L
WPCSYS, -1 $ CYL4,R+D1,R+H1,R,,,,L

VGEN,10,2, , ,DIS+R*2, , , ,0
*DO,I,1,10,1                 
VGEN,4,I+1, , ,B, , , ,0   
*ENDDO

*DO,I,1,10*4,1                 
VSBV,ALL,I+1              
*ENDDO

WPCSYS, -1 $ WPOFF,,H2 $ WPROTA,,90 $ VSBW, ALL   
WPCSYS, -1 $ *DO,I,1,4,1 $ WPOFF,,,L2 $ VSBW, ALL $ *ENDDO  

WPCSYS, -1 $ WPROTA,,,90                           
*DO,I,1,4,1
WPOFF,,,D2 $ VSBW, ALL
WPOFF,,,D1+R*2+DIS*0.5-D2 $ VSBW, ALL
WPOFF,,,R*2+DIS $ VSBW, ALL
WPOFF,,,R*2+DIS $ VSBW, ALL
WPOFF,,,R*2+DIS $ VSBW, ALL
WPOFF,,,R*2+DIS $ VSBW, ALL
WPOFF,,,R*2+DIS $ VSBW, ALL
WPOFF,,,R*2+DIS $ VSBW, ALL
WPOFF,,,R*2+DIS $ VSBW, ALL
WPOFF,,,D1+R*2+DIS*0.5-D2 $ VSBW, ALL
WPOFF,,,D2
*ENDDO

WPCSYS, -1
SELTOL,1.0E-12     


ESIZE,ESIZE3
LSEL,S,TAN1,Z,-1 $ LSEL,A,TAN1,Z, 1 $ CM, LINE1,LINE $ LESIZE,ALL,ESIZE1
! ASEL,S,LOC,Z,0 $ MESHAPE,0,2D $ MSHKEY,0 $ SMRT,5 $ AMESH,ALL,,1
VSEL,ALL $ VSWEEP , ALL
! ALLSEL $ VCLEAR,ALL


LSEL,S,TAN1,Z $ LSEL,R,LOC,Y,H2 $ LSEL,R,LOC,Z,L2,L2*4 $ CM,HYGJ,LINE $ LATT,2,1,3 $ LESIZE,ALL,ESIZE1

LSEL,S,TAN1,Z,1 $ LSEL,A,TAN1,Z,-1 $ LSEL,R,LOC,Y,H2 $ LSEL,R,LOC,X,D2,B*4-D2 $ CM,ZYGJ,LINE $ LATT,2,1,3 $ LESIZE,ALL,ESIZE1

CMSEL,S,HYGJ $ CMSEL,A,ZYGJ $ CM,GJ,LINE
LMESH ,ALL $ LSEL,ALL


TREF,0,            
NCNV,2                              

/SOLU

LSEL,S,LOC,Y,0 $ LSEL,R,LOC,Z,    0 $ DL, ALL,,UX $ DL, ALL,,UY $ DL, ALL,,UZ
LSEL,S,LOC,Y,0 $ LSEL,R,LOC,Z,  B*4 $ DL, ALL,,UX $ DL, ALL,,UY $ DL, ALL,,UZ
LSEL,S,LOC,Y,0 $ LSEL,R,LOC,X,    0 $ DL, ALL,,UX $ DL, ALL,,UY $ DL, ALL,,UZ
LSEL,S,LOC,Y,0 $ LSEL,R,LOC,X,  B*4 $ DL, ALL,,UX $ DL, ALL,,UY $ DL, ALL,,UZ
BFL,HYGJ,TEMP,-458 $ BFL,ZYGJ,TEMP,-458     

ASEL,S,LOC,Y,H $ SFA,ALL,1, PRES ,Q0 $ ALLSEL, ALL
ACEL,,9.8
ANTYPE,0 $ NSUBST, 100 $ OUTRES ,ALL, ALL $ AUTOS,ON
NEQIT,50                              
CNVTOL,U,,0. 05                       
SOLVE
/POST1
PLDISP                 
PLNSOL,S,1 @wuming524，

我这里看不到有没有正确上传附件，所以楼上把代码附上了。把代码复制到txt里面，另存为.lgw文件，可以用ansys直接读取。
迭代运行一段时间后，可能会提示位移过大 @wuming524

上面是报错提示：
the value of ux at node 27584 is greater than the current limit. This generally indicates rigid body motion as a result of an unconstrained model. Verify that your model is properly constrained. @wuming524

上图是迭代过程，一直挺合理，到一定程度，突然挠度无穷大，程序直接崩溃 @wuming524

约束条件有问题，检查模型