簡介

一開始想做一個繩索單元懸掛重物的模擬，其實想法很簡單，但是在實現過程中卻出現了很大的問題，糾結了很久，初步歸結為：繩索單元在垂直其單元方向上受力，會導致其產生很大的變形，從而導致其不收斂。因此專門設計了用繩索單元拉伸樑的中點的一個案例，經過模擬可以求解，程式碼和結果如下。

APDL程式碼

FINISH
/CLEAR
/PREP7
!建立單元beam188
ET,1,BEAM188
!設定beam的section
SECTYPE,   1, BEAM, CTUBE, BEAM1, 0
SECOFFSET, CENT
SECDATA,0.01,0.015,0,0,0,0,0,0,0,0,0,0


ET,2,LINK180            !選擇link180單元

SECTYPE,2,LINK, ,link1        !建立Sections，link，編號為1，名字為link1
SECDATA,7.8E-5,              !定義介面尺寸，10mm直徑為
SECCONTROL,0,1              !定義屬性為拉伸，壓縮為2，拉壓為0

!材料1 鋼
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.1E11
MPDATA,PRXY,1,,0.269
MPDATA,DENS,1,,7850

!畫五個點
FLST,3,5,8  
FITEM,3,1,0,0   
FITEM,3,0,0,0   
FITEM,3,0,1,0   
FITEM,3,1,1,0   
FITEM,3,2,1,0   
K, ,P51X
!將五個點連成直線
LSTR,       1,       2  
LSTR,       2,       3  
LSTR,       3,       4  
LSTR,       4,       5  
LSTR,       1,       4  


FLST,5,4,4,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
CM,_Y,LINE  
LSEL, , , ,P51X 
CM,_Y1,LINE 
CMSEL,S,_Y  
CMSEL,S,_Y1 
LATT,1, ,1, , , ,1  !設定直線的單元屬性為beam
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1    
CM,_Y,LINE 

LSEL, , , ,       5 
CM,_Y1,LINE 
CMSEL,S,_Y  
CMSEL,S,_Y1 
LATT,1, ,2, , , ,2  !設定豎直線的單元屬性為link
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
FLST,5,5,4,ORDE,2   
FITEM,5,1   
FITEM,5,-5  
CM,_Y,LINE  
LSEL, , , ,P51X 
CM,_Y1,LINE 
CMSEL,,_Y    
LESIZE,_Y1, , ,10, , , , ,1 !十等分
!*  
FLST,2,5,4,ORDE,2   
FITEM,2,1   
FITEM,2,-5  
LMESH,P51X  
FINISH  


!設定載荷
FLST,2,1,3,ORDE,1   
FITEM,2,2   
/GO 
DK,P51X, , , ,0,ALL, , , , , ,  
FLST,2,1,3,ORDE,1   
FITEM,2,5   

/GO 
FK,P51X,FY,100  
FINISH  
/SOL
ANTYPE,0
NLGEOM,1!開啟大變形



!求解
/STATUS,SOLU
SOLVE   
FINISH  

!後處理，顯示結果
/POST1  
SET,FIRST   
!*  
/EFACET,1   
PLNSOL, U,SUM, 0,1.0
 

結果圖