!********************************************************************************
!			這是一個真正的梁單元計算簡支梁的例子			*
!										*
!				2000.3.8完成					*
!********************************************************************************
!先定義接點（其中第四點是為定義梁截面方向而設的）
/PREP7  
K, ,,,, 
K, ,400,,,  
K, ,800,,,  
K, ,1400,100,,   	!在梁‘體’外，顯示接點誤差較小，如1400為400則誤差很大。
!連線
LSTR,       1,       2  
LSTR,       2,       3  
!*  
!定義單元類型
ET,1,BEAM188
!*  
!定義材料屬性
UIMP,1,EX, , ,1e5,  
UIMP,1,NUXY, , ,0.3,  
!
!*定義截面形狀（MM>SECTIONS>-BEAM-COMMON SECTIONS) 
!
SECTYPE,   1, BEAM, RECT,   
SECOFFSET, CENT 
SECDATA,50,20,0,0,0,0,0,0,0,0   
!
!
!定義梁截面方向
!
FLST,5,2,4,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
CM,_Y,LINE  
LSEL, , , ,P51X 
CM,_Y1,LINE 
CMSEL,S,_Y  
!*  
!*  
CMSEL,S,_Y1 
LATT,1, ,1, ,       4, ,1   
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
!分網
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,1   
FITEM,5,-2  
CM,_Y,LINE  
LSEL, , , ,P51X 
CM,_Y1,LINE 
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,50, , , , ,1 
CMDELE,_Y   
CMDELE,_Y1  
!*  
FLST,2,2,4,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
LMESH,P51X  
/UI,MESH,OFF
!*  
!*  
!顯示梁截面 
!*  
!*  
!*  
/SHRINK,0   
/ESHAPE,1.0 
/EFACET,1   
/RATIO,1,1,1
/REPLOT 
!*  
!*  
!*  
!加約束邊界條件  
!
!* 約束  
!
FINISH  
/SOLU   
FLST,2,1,3,ORDE,1   
FITEM,2,1   
!*  
/GO 
DK,P51X, , , ,0,UX,UY,UZ, , , , 
FLST,2,1,3,ORDE,1   
FITEM,2,3   
!*  
/GO 
DK,P51X, , , ,0,UY,UZ, , , , ,  
!
!加載
!
FLST,2,1,3,ORDE,1   
FITEM,2,2   
!*  
/GO 
FK,P51X,FY,-1000
!
!求解
!
/STATUS,SOLU
SOLVE   
!
!顯示應力分布云圖
!
FINISH  
/POST1  
/EFACE,1
AVPRIN,0,0, 
!*  
PLNSOL,S,1,0,1