FEA P1: Static/frequency analysis of a PIPE-network

! Created by Tamas Poloskei, 2018. February 20.

Units: [N],[N],[N/m^2=Pa],[N/m]

D=200E-3
T=10E-3
R=300E-3
ALPHA=1.2E-5
THEATA=150
S=50E3
P=6e6
G=9.81
RHO=7800
ELAST=2E11
NU=0.3
M=100

! Second moment of inertias
THETAX=0
THETAY=0
THETAZ=0
! Modified inertias, uncomment it for the modified problem
! THETAX=2.25
! THETAY=2.25
! THETAZ=4.56

SETUP: element types, real constants, materials, sections

KEYWKEYW, Keyword, KEY Sets a keyword used by the GUI for context filtering (GUI).,PR_STRUC,1! Preferences / Structural
/PREP7/PREP7 Enters the model creation preprocessor.
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,1,PIPE16
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,2,PIPE17
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,3,PIPE18
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,4,COMBIN14
ETET, ITYPE, Ename, KOP1, KOP2, KOP3, KOP4, KOP5, KOP6, INOPR Defines a local element type from the element library.,5,MASS21
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,1,DD = 200E-3 ,TT = 10E-3
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,2,DD = 200E-3 ,TT = 10E-3 ,1,,,,DD = 200E-3 ,TT = 10E-3 ,1,,,DD = 200E-3 ,TT = 10E-3 ,1
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,3,DD = 200E-3 ,TT = 10E-3 ,RR = 300E-3
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,4,SS = 50E3
RR, NSET, R1, R2, R3, R4, R5, R6 Defines the element real constants.,5,MM = 100,MM = 100,MM = 100,THETAXTHETAX = 0,THETAYTHETAY = 0,THETAZTHETAZ = 0
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.,1,0
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,EX,1,,ELASTELAST = 2E11
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,PRXY,1,,NUNU = 0.3
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.,1,0
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,DENS,1,,RHORHO = 7800
MPTEMPMPTEMP, SLOC, T1, T2, T3, T4, T5, T6 Defines a temperature table for material properties.,1,0
UIMPUIMP, MAT, Lab1, Lab2, Lab3, VAL1, VAL2, VAL3 Defines constant material properties (GUI).,1,REFT
MPDATAMPDATA, Lab, MAT, SLOC, C1, C2, C3, C4, C5, C6 Defines property data to be associated with the temperature table.,ALPX,1,,ALPHAALPHA = 1.2E-5

! Set angle of view
/RGB/RGB, Kywrd, PRED, PGRN, PBLU, N1, N2, NINC, NCNTR Specifies the RGB color values for indices and contours.,INDEX,100,100,100,0
/RGB/RGB, Kywrd, PRED, PGRN, PBLU, N1, N2, NINC, NCNTR Specifies the RGB color values for indices and contours.,INDEX,0,0,0,15
/VIEW/VIEW, WN, XV, YV, ZV Defines the viewing direction for the display.,1,1,1,1
/ANG/ANGLE, WN, THETA, Axis, KINCR Rotates the display about an axis.,1
/REP/REPLOT, Label Automatically reissues the last display command for convenience.,FAST
! Create the keypoints
KK, NPT, X, Y, Z Defines a keypoint.,1,8500E-3,-5000E-3,-1000E-3
KK, NPT, X, Y, Z Defines a keypoint.,2,8500E-3,-5000E-3,-300E-3
KK, NPT, X, Y, Z Defines a keypoint.,3,8500E-3,-4700E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,4,8500E-3,-3500E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,5,8500E-3,-2300E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,6,8500E-3,-2000E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,7,8500E-3,-1700E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,8,8500E-3,-500E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,9,8500E-3,1700E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,10,8200E-3,2000E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,11,7500E-3,2000E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,12,3300E-3,2000E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,13,3000E-3,1700E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,14,3000E-3,1000E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,15,3000E-3,300E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,16,2700E-3,0,0
KK, NPT, X, Y, Z Defines a keypoint.,17,0,0,0
KK, NPT, X, Y, Z Defines a keypoint.,18,8500E-3,-2000E-3,300E-3
KK, NPT, X, Y, Z Defines a keypoint.,19,8500E-3,-2000E-3,3000E-3
KK, NPT, X, Y, Z Defines a keypoint.,20,7500E-3,2000E-3,-1000E-3
KK, NPT, X, Y, Z Defines a keypoint.,21,5000E-3,2000E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,22,8500E-3,-4700E-3,-300E-3
KK, NPT, X, Y, Z Defines a keypoint.,23,8200E-3,1700E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,24,3300E-3,1700E-3,0
KK, NPT, X, Y, Z Defines a keypoint.,25,2700E-3,300E-3,0
! Create the lines
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,1,2
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,3,4
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,4,5
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,7,8
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,8,9
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,10,11
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,11,21
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,21,12
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,13,14
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,14,15
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,16,17
LL, P1, P2, NDIV, SPACE, XV1, YV1, ZV1, XV2, YV2, ZV2 Defines a line between two keypoints.,18,19
! Meshing setup
LESIZELESIZE, NL1, SIZE, ANGSIZ, NDIV, SPACE, KFORC, LAYER1, LAYER2, KYNDIV Specifies the divisions and spacing ratio on unmeshed lines.,ALL,100E-3
! Mesh straight lines with PIPE16 element
TYPETYPE, ITYPE Sets the element type attribute pointer.,1
REALREAL, NSET Sets the element real constant set attribute pointer.,1
LMESHLMESH, NL1, NL2, NINC Generates nodes and line elements along lines.,ALL
! Create further nodes
NKPTNKPT, NODE, NPT Defines a node at an existing keypoint location.,,22
NKPTNKPT, NODE, NPT Defines a node at an existing keypoint location.,,6
NKPTNKPT, NODE, NPT Defines a node at an existing keypoint location.,,23
NKPTNKPT, NODE, NPT Defines a node at an existing keypoint location.,,20
NKPTNKPT, NODE, NPT Defines a node at an existing keypoint location.,,24
NKPTNKPT, NODE, NPT Defines a node at an existing keypoint location.,,25
! Create PIPE17 T element manually
TYPETYPE, ITYPE Sets the element type attribute pointer.,2
REALREAL, NSET Sets the element real constant set attribute pointer.,2
KI = 5
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KJ = 6
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KJ)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEJ,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KK = 7
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KK)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KK)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KK)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEK,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KL = 18
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KL)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KL)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KL)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEL,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,NODEI,NODEJ,NODEK,NODEL
! Create PIPE18 elbow elements
TYPETYPE, ITYPE Sets the element type attribute pointer.,3
REALREAL, NSET Sets the element real constant set attribute pointer.,3
KI = 2
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KJ = 3
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KJ)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEJ,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KK = 22
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KK)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KK)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KK)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEK,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,NODEI,NODEJ,NODEK
KI = 9
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KJ = 10
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KJ)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEJ,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KK = 23
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KK)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KK)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KK)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEK,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,NODEI,NODEJ,NODEK
KI = 12
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KJ = 13
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KJ)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEJ,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KK = 24
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KK)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KK)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KK)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEK,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,NODEI,NODEJ,NODEK
KI = 15
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KJ = 16
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KJ)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEJ,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KK = 25
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KK)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KK)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KK)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEK,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,NODEI,NODEJ,NODEK
! Create COMBIN14 spring element
TYPETYPE, ITYPE Sets the element type attribute pointer.,4
REALREAL, NSET Sets the element real constant set attribute pointer.,4
KI = 11
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
KJ = 20
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KJ)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KJ)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEJ,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,NODEI,NODEJ
! Create MASS21 structural mass element
TYPETYPE, ITYPE Sets the element type attribute pointer.,5
REALREAL, NSET Sets the element real constant set attribute pointer.,5
KI = 8
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
EE, I, J, K, L, M, N, O, P Defines an element by node connectivity.,NODEI
! Define kinematic constraints
KI = 1
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,NODEI,ALL
KI = 17
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,NODEI,ALL
KI = 19
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,NODEI,ALL
KI = 20
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,SS = 50E3 ,LOC,X,KX(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Y,KY(KI)
NSELNSEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of nodes.,RR = 300E-3,LOC,Z,KZ(KI)
*GET*GET, Par, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM Retrieves a value and stores it as a scalar parameter or part of an array parameter.,NODEI,NODE,0,NUM,MAX
ALLSELALLSEL, LabT, Entity Selects all entities with a single command.
DD, Node, Lab, VALUE, VALUE2, NEND, NINC, Lab2, Lab3, Lab4, Lab5, Lab6 Defines degree-of-freedom constraints at nodes.,NODEI,ALL

Static analysis

Load case 1 - internal pressure

! Select the pipe elements, REAL constants 1, 2 and 3
ESELESEL, Type, Item, Comp, VMIN, VMAX, VINC, KABS Selects a subset of elements.,SS = 50E3 ,REAL,,1,3,1
! Apply internal pressure
SFESFE, Elem, LKEY, Lab, KVAL, VAL1, VAL2, VAL3, VAL4 Specifies surface loads on elements.,ALL,1,PRES,,PP = 6e6
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,1

Load case 2 - own weight

! Delete internal pressure
SFEDELESFEDELE, ELEM, LKEY, Lab Deletes surface loads from elements.,ALL,1,PRES
ACELACEL, ACEL_X, ACEL_Y, ACEL_Z Specifies the linear acceleration of the global Cartesian reference frame for the analysis.,0,0,GG = 9.81
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,2

Load case 3 - thermal expansion

! Delete gravity field
ACELACEL, ACEL_X, ACEL_Y, ACEL_Z Specifies the linear acceleration of the global Cartesian reference frame for the analysis.,0,0,0
TUNIFTUNIF, TEMP Assigns a uniform temperature to all nodes.,20
BFLBFL, Line, Lab, VAL1, VAL2, VAL3, VAL4 Defines a body force load on a line.,ALL,TEMP,THETA
LSWRITELSWRITE, LSNUM Writes load and load step option data to a file.,3
! Add thickness and shape of element
/RGB/RGB, Kywrd, PRED, PGRN, PBLU, N1, N2, NINC, NCNTR Specifies the RGB color values for indices and contours.,INDEX,100,100,100,0
/RGB/RGB, Kywrd, PRED, PGRN, PBLU, N1, N2, NINC, NCNTR Specifies the RGB color values for indices and contours.,INDEX,0,0,0,15
/VIEW/VIEW, WN, XV, YV, ZV Defines the viewing direction for the display.,1,1,1,1
/ANG/ANGLE, WN, THETA, Axis, KINCR Rotates the display about an axis.,1
/ESHAPE/ESHAPE, SCALE, KEY Displays elements with shapes determined from the real constants or section definition.,1.5
/EFACET/EFACET, NUM Specifies the number of facets per element edge for PowerGraphics displays.,1
FINISHFINISH Exits normally from a processor.
! Solution
/SOL/SOLU Enters the solution processor.
LSSOLVELSSOLVE, LSMIN, LSMAX, LSINC Reads and solves multiple load steps.,1,3,1
! Post processing
/POST1/POST1 Enters the database results postprocessor.
LCDEFLCDEF, LCNO, LSTEP, SBSTEP, KIMG Creates a load case from a set of results on a results file.,1,1,Last
LCDEFLCDEF, LCNO, LSTEP, SBSTEP, KIMG Creates a load case from a set of results on a results file.,2,2,Last
LCDEFLCDEF, LCNO, LSTEP, SBSTEP, KIMG Creates a load case from a set of results on a results file.,3,3,Last
! element tables
FINISHFINISH Exits normally from a processor.

Modal analysis

/SOL/SOLU Enters the solution processor.
ANTYPEANTYPE, Antype, Status, LDSTEP, SUBSTEP, Action Specifies the analysis type and restart status.,2
MODOPTMODOPT, Method, NMODE, FREQB, FREQE, Cpxmod, Nrmkey, ModType, BlockSize, --, --, Scalekey Specifies modal analysis options.,LANB,10
EQSLVEQSLV, Lab, TOLER, MULT, --, KeepFile Specifies the type of equation solver.,SPAR
MXPANDMXPAND, NMODE, FREQB, FREQE, Elcalc, SIGNIF, MSUPkey, ModeSelMethod Specifies the number of modes to expand and write for a modal or buckling analysis.,10,,,1
MODOPTMODOPT, Method, NMODE, FREQB, FREQE, Cpxmod, Nrmkey, ModType, BlockSize, --, --, Scalekey Specifies modal analysis options.,LANB,10,,,,OFF
/STATUS/STATUS, Lab Lists the status of items for the run.,SOLU
SOLVESOLVE, Action Starts a solution.
FINISHFINISH Exits normally from a processor.
/POST1/POST1 Enters the database results postprocessor.
SETSET, Lstep, Sbstep, Fact, KIMG, TIME, ANGLE, NSET, ORDER Defines the data set to be read from the results file.,LIST