너희들 이게 몬지 알아..??

/AUX15
IOPTN,IGES,NODEFEAT
IGESIN,'600lb_100','IGS',''
FINI

C*** DEFINE PARAMETER

THETA    = 45  ! FULL OPEN ANGLE
R      = 25 ! OPENING RATIO
XC     = -44  ! HINGE POINT
YC     = 90

D = 92 ! INLET DIAMETER
L = 305 ! VALVE LENGTH

I_VELOCITY = -114
O_PRESSURE = 8.9
F_DENSITY = 667.1
F_VISCOSITY = 0.000078
fix_1 = 30
fix_2 = 31

/FILENAME,600#_100_%R%_Rev
/TITLE, SWING CHECK VALVE 600lb 100

C***
C*** 1. Build the model of the entire domain:
C***

/PREP7

NUMCMP,ALL

ET,1,141
ET,2,42

LOCAL,11,2,XC,YC,0
AGEN,,2,,,,-THETA,,,,1
AGEN,,2,,,,THETA*R/100,,,,1

RECT,-(D*10+L),-(L/2),-(D/2),(D/2)
RECT,(L/2),(D*10+L),-(D/2),(D/2)

AOVLAP,ALL

NUMCMP,ALL

LESIZE,10,,,50,-10,,,,0
LESIZE,20,,,50,-10,,,,0
LESIZE,32,,,50,-10,,,,0
LESIZE,33,,,50,-10,,,,0
LESIZE,34,,,100,10,,,,0
LESIZE,35,,,100,10,,,,0
LESIZE,36,,,100,10,,,,0
LESIZE,37,,,100,10,,,,0
AMESH,2,3

CSYS,0
LSEL,S,LOC,X,-L+0.0001,L-0.0001
LESIZE,ALL,1,,0,0,4,0,-2
ALLS

AMESH,4

TYPE,2 $MAT,2
ESIZE,2
AMESH,1
ALLS

ARSCALE,ALL,,,0.001,0.001,1,,0,1

C***
C*** 2. Create Physics Environment for the Fluid
C***

et,1,141
et,2,0 ! Disk becomes the Null Element

/com, CFD Solution Control
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TURB,1
FLDATA2,ITER,EXEC,400,
FLDATA5,OUTP,SUMF,10,
/com, CFD Property Information
FLDATA7,PROT,DENS,CONSTANT
FLDATA8,NOMI,DENS,F_DENSITY,
FLDATA7,PROT,VISC,CONSTANT
FLDATA8,NOMI,VISC,F_VISCOSITY,

LSEL,s,LOC,X,-(L/2)*0.001
LSEL,A,LOC,X,(L/2)*0.001
LSEL,A,LOC,X,(D*10+L)*0.001
LSEL,A,LOC,X,-(D*10+L)*0.001
LSEL,INVE

DL,all, ,VX,0.0,1
DL,all, ,VY,0.0,1

LSEL,S,LOC,X,(D*10+L)*0.001
DL,all, ,VX,I_VELOCITY,1
DL,all, ,VY,0.0,1

lsel,s,loc,x,-(D*10+L)*0.001
DL,all, ,PRES,0.0,1
alls

/TITLE,FLUID ANALYSIS
PHYSICS,WRITE,FLUID,600#100_%R%_FLUID
PHYSICS,CLEAR

C***
C*** 3.CREATE PHYSICS ENVIRONMENT FOR THE STRUCTURE
C***

/PREP7
ET,1,0 ! THE NULL ELEMENT FOR THE FLUIDE
ET,2,42 ! DISC ELEMENT
MP,EX,2,195E9
MP,PRXY,2,0.3

DK,FIX_1,ALL
DK,FIX_2,ALL

/TITLE,STRUCTURE ANALYSIS
PHYSICS,WRITE,STRUC,600#_100_%R%_STURC
PHYSICS,CLEAR
FINI

SAVE

/SOLU
PHYSICS,READ,FLUID ! READ IN FLUID
SOLVE
FINI

/SOLU
!ANTYPE,STAT,REST
PHYSICS,READ,STRUC
LDREAD,PRES,LAST,,,,,rfl ! APPLY PRESSURE SURFACE LOAD FROM FLOTRAN
SOLVE
FINI

SAVE

/POST1
PHYSICS,READ,STRUC
SET,LAST

KSEL,S,,,FIX_1
KSEL,A,,,FIX_2

NSLK,R

*get,in_1_point,node,0,num,max
*get,in_2_point,node,0,num,min

*get,in_1_x,node,in_1_point,rf,fx,
*get,in_1_y,node,in_1_point,rf,fy,
*get,in_2_x,node,in_2_point,rf,fx,
*get,in_2_y,node,in_2_point,rf,fy,

/out,600#_100_%R%_REACTION,lis
/COM,- Reaction force for 600#_100_%r%_REV
/COM,
/COM, %r%% open
/COM,
/COM, At node No. %in_1_point%, Reaction force of 1st point FX = %in_1_x% N
/COM, At node No. %in_1_point%, Reaction force of 1st point FY = %in_1_y% N
/COM
/COM, At node No. %in_2_point%, Reaction force of 2st point FX = %in_2_x% N
/COM, At node No. %in_2_point%, Reaction force of 2st point FY = %in_2_y% N
/OUT
fini

save

/POST1

PHYSICS,READ,FLUID
set,last

in_x = (d*10+l)
in_n = node(in_x,0,0)

*get,in_p,node,in_n,pres

/out,600#_100_%R%_REV,lis
/COM,- Differential Pressure for 600#_100_%R%_REV
/COM,
/COM, %r%% open
/com, At node No. %in_n%, PSID = %in_p% Pa
/OUT
fini

SAVE

나도 몰라...ㅡㅡㅋ