用SHELL单元做工字梁求侧压力的极限强度中遇到的问题

chenzhimei

用SHELL单元做工字梁求侧压力时的极限强度:简支的边界条件在ANSYS中具体应怎么实现?(向上为Y,向右为X,由纸向外为Z,0为受约束) 我用的是PLATE:T[1 0 1] , R[1 0 0] WEB:T[1 1 0] R[0 1 0] FLANGE:T[1 0 1] , R[1 0 0] 但我计算出的值比理论值大太多.我还试了只在质点上加边界条件,但不稳定.请大家好心帮我看看我的方法上有没有错呀.我这东西做了太久了,一直找不到原因, 没的办法,苦恼呀...谢谢
/PREP7  
ET,1,SHELL181   
R,1,19.5, , , , , ,
RMORE, , , , , , ,  
R,2,10,0,0,0,0,0,   
RMORE,0,0,0,0,0,0,  
R,3,22,0,0,0,0,0,   
RMORE,0,0,0,0,0,0,  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,205800
MPDATA,PRXY,1,,0.3  
TB,BISO,1,1,2,  
TBTEMP,0
TBDATA,,315,2058,,,,
WPSTYLE,,,,,,,,1
WPSTYLE,,,,,,,,0
WPSTYLE,,,,,,,,1
wprot,0,90,0
RECTNG,0,5100/2,-455,455,   
wprot,0,90,0
RECTNG,0,5100/2,0,19.5/2+22/2+610,  
wpoff,0,(19.5/2+22/2+610),0
wprot,0,-90,0   
RECTNG,0,5100/2,175/2,-175/2,   
FLST,2,3,5,ORDE,2   
FITEM,2,1   
FITEM,2,-3  
AOVLAP,P51X
FLST,5,6,4,ORDE,6   
FITEM,5,1   
FITEM,5,3   
FITEM,5,5   
FITEM,5,7   
FITEM,5,9   
FITEM,5,11  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,100/2, , , , ,1  
FLST,5,4,4,ORDE,2   
FITEM,5,13  
FITEM,5,-16
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,9, , , , ,1  
FLST,5,4,4,ORDE,2   
FITEM,5,17  
FITEM,5,-20
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,2, , , , ,1  
FLST,5,2,4,ORDE,2   
FITEM,5,6   
FITEM,5,8   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
LESIZE,_Y1, , ,12, , , , ,1
MSHKEY,0
FLST,5,5,5,ORDE,3   
FITEM,5,2   
FITEM,5,4   
FITEM,5,-7  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
FINISH  
/PREP7  
FINISH  
/SOL
FLST,2,7,4,ORDE,7   
FITEM,2,1   
FITEM,2,3   
FITEM,2,6   
FITEM,2,14  
FITEM,2,-15
FITEM,2,17  
FITEM,2,19  
DL,P51X, ,SYMM  
FLST,2,2,4,ORDE,2   
FITEM,2,13  
FITEM,2,16  
/GO
DL,P51X, ,UY,   
FLST,2,2,4,ORDE,2   
FITEM,2,13  
FITEM,2,16  
/GO
DL,P51X, ,ROTY,
FLST,2,2,4,ORDE,2   
FITEM,2,13  
FITEM,2,16  
/GO
DL,P51X, ,ROTX,
FLST,2,11,1,ORDE,2  
FITEM,2,114
FITEM,2,-124
/GO
D,P51X, , , , , ,UZ,ROTX,ROTZ, , ,  
FINISH  
/PREP7  
FLST,4,35,1,ORDE,14
FITEM,4,1   
FITEM,4,64  
FITEM,4,114
FITEM,4,-124
FITEM,4,664
FITEM,4,723
FITEM,4,-730
FITEM,4,1132
FITEM,4,1182
FITEM,4,-1189   
FITEM,4,1582
FITEM,4,1634
FITEM,4,1686
FITEM,4,1736
CP,1,UX,P51X
FINISH  
/SOL
FLST,2,2,5,ORDE,2   
FITEM,2,4   
FITEM,2,-5  
/GO
SFA,P51X,1,PRES,2   
ANTYPE,0
NLGEOM,1
NSUBST,100,200,50   
OUTRES,ERASE
OUTRES,ALL,1
NEQIT,5
RESCONTRL,DEFINE,ALL,1,1