*HEADING
BARREL VAULT  UGENS  4 X 4 MESH  S8R
*NODE,SYSTEM=C,NSET=CORNER
 1,300.,0.,0.,
 9,300.,,300.
 81,300.,40.,
 89,300.,40.,300.
*NGEN,LINE=C,NSET=WALL
 1,81,10
*NGEN,LINE=C,NSET=SYMM
 9,89,10,,,,300.
*NFILL,NSET=NALL
WALL,SYMM,8,1
*NCOPY, OLD SET=NALL, NEW SET=NALL1,CHANGE NUMBER=100,SHIFT
0.,0.,350.

*NSET,NSET=WALL,GEN
101,181,10
*NSET,NSET=SYMM,GEN
109,189,10,,,,650.
*NSET,NSET=CORNER
101,109,181,189
*NSET,NSET=APEX,GENERATE
1,9,1
101,109,1
*ELEMENT,TYPE=S8R,ELSET=ROOF1
 1,1,21,23,3,11,22,13,2
*ELGEN,ELSET=ALL
 1,4,20,4, 4,2,1
*ELCOPY,ELEMENT SHIFT=100,OLD SET=ALL, NEW SET=ALL1, SHIFT NODES=100
*ELSET,ELSET=FEW
 1,2,3,4
101,102,103,104
*BOUNDARY
 APEX,YSYMM
 SYMM,ZSYMM
 WALL,1,2
*SHELL GENERAL SECTION,ELSET=ALL,DENSITY=1.0,USER,VARIABLES=60,PROPERTIES=24
 3.0, 
 9.0E6, 0.0, 9.0E6, 0.0, 0.0, 4.5E6, 0.0, 0.0
 0.0, 6.75E6, 0.0, 0.0, 0.0, 0.0, 6.75E6, 0.0
 0.0, 0.0, 0.0, 0.0, 3.375E6,    70.,  1E-6, 0.
** last 3 items on line above:  ZERO, ALPHA, NU
*TRANSVERSE SHEAR STIFFNESS
 3.75E+06, 3.75E+06
*HOURGLASS STIFFNESS
 7500.0, 24998.0, 5625.0
*SHELL GENERAL SECTION,ELSET=ALL1,DENSITY=1.0,USER,VARIABLES=60,PROPERTIES=21,
IPROPERTIES=3, UNSYMM
 3.0, 
 9.0E6, 0.0, 9.0E6, 0.0, 0.0, 4.5E6, 0.0, 0.0
 0.0, 6.75E6, 0.0, 0.0, 0.0, 0.0, 6.75E6, 0.0
 0.0, 0.0, 0.0, 0.0, 3.375E6,    70,  1E+6, 0
** last 3 items on line above:  ZERO, ALPHA, NU
*TRANSVERSE SHEAR STIFFNESS
 3.75E+06, 3.75E+06
*HOURGLASS STIFFNESS
 7500.0, 24998.0, 5625.0
************
*INITIAL CONDITIONS,TYPE=TEMPERATURE
NALL,70.
NALL1,70.
********************
*STEP
*FREQUENCY,EIGENSOLVER=SUBSPACE
 4, 
*ELPRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE
*END STEP
************
*STEP,NLGEOM
 SELF-WEIGHT
*STATIC
 0.1,1.0
*TEMPERATURE,OP=NEW
*DLOAD
 ALL,BX,-1.0
 ALL1,BX,-1.0
*ENERGY PRINT,FREQ=1
*EL PRINT,ELSET=FEW
 SF, 
 SE, 
 TSHR,
 TEMP, 
*EL PRINT,ELSET=FEW,POSITION=AVERAGED AT NODES
 SF, 
 SE, 
 TSHR,
 TEMP, 
*NODE PRINT,NSET=CORNER
 U, 
 RF, 
*EL FILE,ELSET=FEW
 SF, 
 SE, 
 TSHR, 
 TEMP, 
 SDV,
*NODE FILE,NSET=CORNER
 U, 
 RF, 
*END STEP
************
*STEP
*BUCKLE
 3, 
*DLOAD
 ALL,BX,-1.2
 ALL1,BX,-1.2
*ELPRINT,FREQ=0
*NODE PRINT,FREQ=0
*MODAL FILE
*output,field,freq=0
*END STEP
************
*STEP,PERTURBATION
THERMAL LOADING
*STATIC
*TEMPERATURE
NALL,120.
NALL1,120.
*ENERGY PRINT
*EL PRINT,ELSET=FEW
 SF, 
 SE, 
 TSHR,
 TEMP,
*EL PRINT,ELSET=FEW,POSITION=AVERAGED AT NODES
 SF, 
 SE, 
 TSHR,
 TEMP,
*NODE PRINT,NSET=CORNER
 U, 
 RF, 
*EL FILE,ELSET=FEW
 SF, 
 SE, 
 TSHR, 
 TEMP, 
*NODE FILE,NSET=CORNER
 U, 
 RF, 
*END STEP
************
*STEP
THERMAL LOADING 2
*STATIC
*TEMPERATURE
NALL,120.
NALL1,120.
*ENERGY PRINT
*EL PRINT,ELSET=FEW
 SF, 
 SE, 
 TSHR,
 TEMP,
*EL PRINT,ELSET=FEW,POSITION=AVERAGED AT NODES
 SF, 
 SE, 
 TSHR,
 TEMP,
*NODE PRINT,NSET=CORNER
 U, 
 RF, 
*EL FILE,ELSET=FEW
 SF, 
 SE, 
 TSHR, 
 TEMP, 
 SDV,
*NODE FILE,NSET=CORNER
 U, 
 RF, 
*END STEP