!********************************************************************************
!										*
!		BEAMUNIT是驗證長度，力和應力及E，MU等單位的小程序		*
!				2000.3.8修改完成				*
!			梁長=400MM，梁高=20MM，梁寬=50MM			*
!				  E=1E5，MU=0.3					*
!********************************************************************************
!
/prep7  
ET,1,SOLID95
!*  

!*  
UIMP,1,EX, , ,1e5,   
UIMP,1,PRXY, , ,.3, 
R,1,H, , , , , ,        !Define Real Constant(DELTA)
!*  
!   
*SET,B1,50              !wide1
*SET,B2,50              !wide2
*SET,L,400              !length
*SET,H,20               !height
*SET,P,1000
V=H*L*(B1+B2)/2
!*  
!BLOCK,0,L,B1,B2,H,, 
K, ,0,0,0,              !KP1
K, ,0,0,B1,             !KP2
K, ,L,0,0,              !KP3
K, ,L,0,B2,             !KP4
K, ,0,H,0,              !KP5
K, ,0,H,B1,             !KP6
K, ,L,H,0,              !KP7
K, ,L,H,B2,             !KP8
!
LSTR,       Kp(0,0,0),       Kp(0,0,B1)      !line1 ( 1,       2 )
LSTR,       Kp(0,0,B1),      Kp(L,0,B2)      !line2 ( 2,       4 )
LSTR,       Kp(L,0,B2),      Kp(L,0,0)       !line3 ( 4,       3 )
LSTR,       Kp(L,0,0),       Kp(0,0,0)       !line4 ( 3,       1 )
LSTR,       Kp(0,0,0),       Kp(0,H,0)       !line5 ( 1,       5 )
LSTR,       Kp(0,H,0),       Kp(0,H,B1)      !line6 ( 5,       6 )
LSTR,       Kp(0,H,B1),      Kp(L,H,B2)      !line7 ( 6,       8 )
LSTR,       Kp(L,H,B2),      Kp(L,H,0)       !line8 ( 8,       7 )
LSTR,       Kp(L,H,0),       Kp(0,H,0)       !line9 ( 7,       5 )
LSTR,       Kp(0,0,B1),      Kp(0,H,B1)      !line10( 2,       6 )
LSTR,       Kp(L,0,0),       Kp(L,H,0)       !line11( 3,       7 )
LSTR,       Kp(L,0,B2),      Kp(L,H,B2)      !line12( 4,       8 )
!
!Create areas       
!1
FLST,5,4,4,ORDE,4   
FITEM,5,1   
FITEM,5,5   
FITEM,5,6  
FITEM,5,10  
LSEL,R, , ,P51X 
AL,all 
!2
FLST,5,4,4,ORDE,4   
FITEM,5,3   
FITEM,5,8   
FITEM,5,11  
FITEM,5,12  
LSEL,R, , ,P51X 
LPLOT   
AL,all 
!3
FLST,5,4,4,ORDE,4   
FITEM,5,6   
FITEM,5,7   
FITEM,5,8  
FITEM,5,9  
LSEL,R, , ,P51X 
LPLOT   
AL,all 
!4
FLST,5,4,4,ORDE,4   
FITEM,5,1   
FITEM,5,2   
FITEM,5,3  
FITEM,5,4  
LSEL,R, , ,P51X 
LPLOT   
AL,all 
!5
FLST,5,4,4,ORDE,4   
FITEM,5,4   
FITEM,5,5   
FITEM,5,9  
FITEM,5,11  
LSEL,R, , ,P51X 
LPLOT   
AL,all 
!6
FLST,5,4,4,ORDE,4   
FITEM,5,2   
FITEM,5,7   
FITEM,5,10  
FITEM,5,12  
LSEL,R, , ,P51X 
LPLOT   
AL,all 
!
!Create Volume
!
allsel,all
!FLST,2,6,5,ORDE,2   
VA,all                                      !VOLUME HAS BEEN CREATED
!   
!Select line KP5-KP6
!ALLSEL,ALL
!LSEL,S,LOC,X,0                              !Select lines located x=0  
!LSEL,R,LOC,Y,H                              !Select lines located Y=H
!LESIZE,all,,,2,1.0,1.0,,                    !devide into 4 sections
!
!*  
!Select line KP2-KP6
!ALLSEL,ALL
!LSEL,S,LOC,X,0                              !Select lines located x=0  
!LSEL,R,LOC,Z,B1                             !Select lines located Z=H  
!LESIZE,all,,,2,1.0,1.0,,                    !devide into 4 sections
!   
!* 
!ALLSEL,ALL
!LSEL,S,LOC,z,0                               !Select lines located Y=0  
!LSEL,R,LOC,y,H                               !Select lines located Z=H   
!LESIZE,all,,,5,1.0,1.0,,                    !devide into 10 sections   
!
!
!*  
MSHAPE,0,3D 
MSHKEY,1
!*  
CM,_Y,VOLU  
VSEL, , , ,       1 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDEL,_Y
CMDEL,_Y1   
CMDEL,_Y2   
!
!
!Mesh completed!
!
!以上拷貝等強度梁優(yōu)化小程序（E=1E5，MU=0.3）
!======================================================================
!以下加邊界條件和載荷
!
!三接點UY=0
!
FINISH  
/SOLU   
FLST,2,3,1,ORDE,3   
FITEM,2,1   
FITEM,2,-2  
FITEM,2,4   
/GO 
D,P51X, , , , , ,UY, , , , ,
!
!
!一面對稱
DA,       2,SYMM
!
!加載P=1000N
!
FLST,2,1,1,ORDE,1   
FITEM,2,33  
!*  
/GO 
F,P51X,FY,-P/2 
/STATUS,SOLU
SOLVE   
FINISH  
/POST1  

繪制應力分布云圖
AVPRIN,0,0, 
!*  
PLNSOL,S,1,0,1