!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
