/COM,ANSYS MEDIA REL. 2024R2 (05/10/2024) REF. VERIF. MANUAL: REL. 2024R2
/VERIFY,VMP09-T52-187
/TITLE,VMP09-T52-187, SIMPLY-SUPPORTED SOLID SQUARE PLATE
/COM,   SELECTED BENCHMARKS FOR NATURAL FREQUENCY ANALYSIS
/COM,   JUNE 1987, TEST 52

/PREP7
N1=4
N2=1
ET,1,187
MP,EX,1,200E9
MP,NUXY,1,.3
MP,DENS,1,8000
K,1
K,2,10
K,3,10,10
K,4,,10
KGEN,2,1,4,1,,,1
L,1,5,N2
*REP,4,1,1
V,1,2,3,4,5,6,7,8
ESIZE,,N1
VMESH,1
NSEL,S,LOC,Y,0
NSEL,R,LOC,Z,0
D,ALL,UZ,0
NSEL,S,LOC,X,0
NSEL,R,LOC,Z,0
D,ALL,UZ,0
NSEL,S,LOC,Y,10
NSEL,R,LOC,Z,0
D,ALL,UZ,0
NSEL,S,LOC,X,10
NSEL,R,LOC,Z,0
D,ALL,UZ,0
NSEL,ALL
FINISH
/out,scratch
/SOLVE
ANTYPE,MODAL
MODOPT,LANB,20
MXPAND,10
SOLVE
FINISH
/out
/POST1
*DIM,FREQ,ARRAY,10
*GET,FREQ(1,1),MODE,4,FREQ
*GET,FREQ(2,1),MODE,5,FREQ
*GET,FREQ(3,1),MODE,6,FREQ
*GET,FREQ(4,1),MODE,7,FREQ
*GET,FREQ(5,1),MODE,8,FREQ
*GET,FREQ(6,1),MODE,9,FREQ
*GET,FREQ(7,1),MODE,10,FREQ
R1=FREQ(1,1)/45.897
R2=FREQ(2,1)/109.44
R3=FREQ(3,1)/109.44
R4=FREQ(4,1)/167.89
R5=FREQ(5,1)/193.59
R6=FREQ(6,1)/206.19
R7=FREQ(7,1)/206.19
*DIM,VALUE,,7,2
*DIM,LABEL,CHAR,10
*DIM,RATIO,,7,1
LABEL(1) = 'FREQ 4','FREQ 5','FREQ 6','FREQ 7','FREQ 8','FREQ 9','FREQ 10'
*VFILL,VALUE(1,1),DATA,45.897,109.44,109.44,167.89,193.59,206.19,206.19
*VFILL,VALUE(1,2),DATA,FREQ(1,1),FREQ(2,1),FREQ(3,1),FREQ(4,1),FREQ(5,1),FREQ(6,1),FREQ(7,1)
*VFILL,RATIO(1,1),DATA,R1,R2,R3,R4,R5,R6,R7
FINISH

/PREP7
ET,1,187
LUMPM,ON
FINISH
/out,scratch
/SOLU
SOLVE
FINISH
/out
/POST1
*DIM,FREQ1,ARRAY,10
*GET,FREQ1(1,1),MODE,4,FREQ
*GET,FREQ1(2,1),MODE,5,FREQ
*GET,FREQ1(3,1),MODE,6,FREQ
*GET,FREQ1(4,1),MODE,7,FREQ
*GET,FREQ1(5,1),MODE,8,FREQ
*GET,FREQ1(6,1),MODE,9,FREQ
*GET,FREQ1(7,1),MODE,10,FREQ
R8=FREQ1(1,1)/45.897
R9=FREQ1(2,1)/109.44
R10=FREQ1(3,1)/109.44
R11=FREQ1(4,1)/167.89
R12=FREQ1(5,1)/193.59
R13=FREQ1(6,1)/206.19
R14=FREQ1(7,1)/206.19
*DIM,VALUE1,,7,2
*DIM,LABEL1,CHAR,10
*DIM,RATIO1,,7,1
LABEL1(1) = 'FREQ 4','FREQ 5','FREQ 6','FREQ 7','FREQ 8','FREQ 9','FREQ 10'
*VFILL,VALUE1(1,1),DATA,45.897,109.44,109.44,167.89,193.59,206.19,206.19
*VFILL,VALUE1(1,2),DATA,FREQ1(1,1),FREQ1(2,1),FREQ1(3,1),FREQ1(4,1),FREQ1(5,1),FREQ1(6,1), FREQ1(7,1)
*VFILL,RATIO1(1,1),DATA,R8,R9,R10,R11,R12,R13,R14
/COM,
/COM,------------ VMP09-T52-187 RESULTS COMPARISON ----------------
/COM,
/COM,         |   TARGET   |   Mechanical APDL   |   RATIO
/COM,
*VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),RATIO(1,1)
(1X,A8,'   ',F10.3,'  ',F15.3,'   ',1F8.3)
/COM
/COM,---------------- LUMPED MASS RESULTS --------------------
/COM
/COM,         |   TARGET   |   Mechanical APDL   |   RATIO
/COM,
*VWRITE,LABEL1(1),VALUE1(1,1),VALUE1(1,2),RATIO1(1,1)
(1X,A8,'   ',F10.3,'  ',F15.3,'   ',1F8.3)
FINISH

/POST26
*DIM,LABEL3,CHAR,5
*DIM,VALUE2,,5,3
*DIM,LABEL4,CHAR,5
LABEL3(1) = ' MODE',' MODE',' MODE',' MODE',' MODE'
LABEL4(1) = '4','5,6','7','8','9,10'
*VFILL,VALUE2(1,1),DATA,FREQ(1,1),FREQ(3,1),FREQ(4,1),FREQ(5,1),FREQ(7,1),
*VFILL,VALUE2(1,2),DATA,R1,R3,R4,R5,R7
*DIM,LABEL2,CHAR,5
*DIM,LABEL5,CHAR,5
LABEL2(1) = 'vmp09-t5','vmp09-t5','vmp09-t5','vmp09-t5','vmp09-t5'
LABEL5(1) = '2-187','2-187','2-187','2-187','2-187'

/OUT,vmp09-t52-187,vrt
/COM
/COM,------------------- VMP09-T52 RESULTS COMPARISON ---------------------
/COM,
/COM,            |   Mechanical APDL   |   RATIO  |         INPUT         |
/COM,
/COM, SOLID187
*VWRITE,LABEL3(1),LABEL4(1),VALUE2(1,1),VALUE2(1,2),LABEL2(1),LABEL5(1)
(1X,A5,A5,'      ',F13.3,'    ',F9.4,'        ',A8,A5)
/COM,
/COM,----------------------------------------------------------------
/OUT

FINISH
*LIST,vmp09-t52-187,vrt