[转帖]线性屈曲分析实例

AaronSpark

/PREP7<BR>ET,1,SHELL63<BR>MP,EX,1,2.1E11<BR>MP,PRXY,1,0.3<BR>R,1,0.001<BR><BR>CSYS,2<BR>N,10000,10,,90<BR>*DO,I,1,72<BR>N,I,10,5*(I-1),85<BR>*ENDDO<BR><BR>NGEN,33,72,1,72,1,,,-5<BR>*DO,I,1,71<BR>E,10000,I,I+1<BR>*ENDDO<BR>E,10000,72,1<BR><BR>*DO,J,1,32<BR>*DO,I,1+(J-1)*72,(J-1)*72+71<BR>E,I,I+72,I+73,I+1<BR>*ENDDO<BR>E,72+(J-1)*72,144+(J-1)*72,(J-1)*72+73,1+(J-1)*72<BR>*ENDDO<BR><BR>NUMMRG,ALL<BR>NUMCMP,ALL<BR><BR>*GET,E_TOTAL,ELEM,,COUNT<BR>*DIM,E_CENT,ARRAY,E_TOTAL,3<BR>*DO,I,1,E_TOTAL<BR>*GET,E_CENT(I,1),ELEM,I,CENT,X<BR>*GET,E_CENT(I,2),ELEM,I,CENT,Y<BR>*GET,E_CENT(I,3),ELEM,I,CENT,Z<BR>*ENDDO<BR><BR>*DIM,PRS,ARRAY,E_TOTAL,1<BR><BR>*DO,I,1,E_TOTAL<BR>PRS(I)=SIN(E_CENT(I,2)/45*ATAN(1))*COS((E_CENT(I,3))/45*ATAN(1))<BR>SFE,I,,PRES,,PRS(I)<BR>*ENDDO<BR>NSEL,S,LOC,Z,-75<BR>D,ALL,ALL<BR>ALLSEL,ALL<BR><BR>CSYS,0<BR><BR>FINI<BR>/SOLU<BR><BR>PSTRES,ON<BR><BR>SOLVE<BR>FINI<BR>/SOLU<BR>ANTYPE,BUCKLE<BR>BUCOPT,LANB,10,<BR>SOLVE<BR>FINI<BR>/SOLU<BR>MXPAND,10,,,1<BR>SOLVE