!6.3  用殼體單元分析受均布荷載作用的固支圓板大撓度彎曲
! 本程序來源于邢靜忠等編著的《ANSYS分析實例與工程應用》，機械工業出版社（2004年）
FINISH
/CLEAR, NOSTART
! (1) 設置工程參數,分析類型,單元類型,實常數和材料參數
/FileName, EX6.3
/PREP7
/TITLE, EX6.3, LARGE DEFLECTION BENDING OF A CIRCULAR PLATE
C***    THEORY OF PLATES AND SHELLS, TIMOSHENKO, P. 401, EQ. 232
ANTYPE,STATIC
NLGEOM,ON
SSTIF,ON
ET,1,SHELL51,,,,,,,1
R,1,.0025
MP,EX,1,2E11
MP,NUXY,1,0.3
! (2)定義節點和單元
N,1
N,11,0.25
FILL
E,1,2
EGEN,10,1,1
! (3)設置求解選項,位移約束和荷載
OUTPR,,LAST
NSEL,S,LOC,X,0
DSYM,SYMM,X
NSEL,ALL
D,11,ALL
SFE,1,1,PRES,,6585.2
*REPEAT,10,1
FINISH
! (4)進入求解模塊后,開始求解
/SOLU
SOLVE
FINISH
! (5)在后處理模塊中,顯示計算結果
/POST1
*GET,UY,NODE,1,U,Y
*status,parm
! (6)定義數組,將計算結果保存到數組變量,然后輸出到文件
*DIM,LABEL,CHAR,1,2
*DIM,VALUE,,1,3
LABEL(1,1) = 'Y DISP '
LABEL(1,2) = ' m'
*VFILL,VALUE(1,1),DATA,-0.00125
*VFILL,VALUE(1,2),DATA,UY
*VFILL,VALUE(1,3),DATA,ABS(UY/0.00125)
/COM
/OUT, EX6_3, out
/COM,-------------------  EX6.3 RESULTS COMPARISON -------------
/COM,
/COM,   LOAD STEP 3   |   TARGET   |   ANSYS   |   RATIO
/COM,
*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
(1X,A8,A8,'   ',F10.5,'  ',F10.5,'   ',1F5.3)
/COM,-----------------------------------------------------------
/OUT
FINISH
*LIST, EX6_3, out