*Heading *Node 1, 0.5, 0.5, 1. 2, 0.5, -0.5, 1. 3, 0.5, -0.5, 0. 4, 0.5, 0.5, 0. 5, -0.5, -0.5, 1. 6, -0.5, -0.5, 0. 7, -0.5, 0.5, 0. 8, -0.5, 0.5, 1. 9, 0., 0., 0.5 *Elset, elset=cube1 1,2 *Elset, elset=cube2 3,4 *Elset, elset=cube3 5,6 *Elset, elset=cube4 7,8 *Elset, elset=cube5 9,10 *Elset, elset=cube6 11,12 *Element, type=C3D4, elset=cube 1, 9, 6, 8, 7 2, 9, 3, 6, 4 3, 9, 4, 6, 7 4, 9, 5, 8, 6 5, 9, 3, 5, 6 6, 9, 5, 2, 8 7, 9, 2, 5, 3 8, 9, 4, 2, 3 9, 9, 1, 2, 4 10, 9, 8, 2, 1 11, 9, 7, 1, 4 12, 9, 8, 1, 7 *Distribution Table, name=OrientationTable coord3d, coord3d *Distribution, name=distOri, table=orientationTable, location=element , 1., 0., 0., 0., 1., 0. cube1, 1., 1., 1., -1., 1.,-1. cube2, 1., 1.,-1., -1., 1., 1. cube3, -1., 1.,-1., 1., 1., 1. cube4, -1.,-1., 1., 1.,-1.,-1. cube5, -1., 1., 1., 1., 1.,-1. cube6, -1.,-1.,-1., 1.,-1., 1. *Distribution Table, name=AngleTable angle *Distribution, name=distAngle, table=AngleTable, location=element ,0 cube1, 0.0 cube2, 1.0 cube3, 2.0 cube4, 3.0 cube5, 4.0 cube6, -1.0 *Orientation, name=ori distOri 3, distAngle *Solid Section, elset=cube, material=Material-1, orientation=ori 1., *Nset, nset=Fix 1, 2, 3, 5, 6 *Nset, nset=Load 4, 7, 8 *Material, name=Material-1 *Elastic, type=ENGINEERING CONSTANTS 3e+06, 4e+06, 5e+06, 0.1, 0.2, 0.3,300000.,400000. 500000., *Density 0.750E-3 *Elset, elset=eall, generate 1, 12, 1 *psd-definition, name=white, type=force 1.0,0.1,0.0 1.0,0.1,60.0 *spectrum, name=spec,type=acceleration 0., 0., 0. 1.e5, 10., 0. 1.e5, 100., 0. *Boundary Fix, 1,3 ** ** -------------------------------------------------- ** *step, name=step-1 frequency analysis *frequency 5, *El print s, *Node print u, *Output, field *Node output u, *Element output s, *end step ** ** -------------------------------------------------- ** *step, name=step-2 SSD mode-based *steady state dynamics, interval=eigenfrequency 60, *El print, freq=999 s, *Node print, freq=0 *Output, field, op=new *Element output s, *cload load, 2, 1.e5 load, 3, 1.e5 *select eigenmodes, generate 1, 5 *modal damping, rayleigh 1, 5, 0.5 *end step ** ** -------------------------------------------------- ** *step, name=step-3 modal dynamic *modal dynamic 0.1,1. *cload load, 2, 1.e5 load, 3, 1.e5 *select eigenmodes, generate 1, 5 *modal damping, rayleigh 1, 5, 0.5 *end step ** ** -------------------------------------------------- ** *step, name=step-4 SSD direct *steady state dynamics, direct, interval=eigenfrequency 60., *cload load, 2, 100. *end step ** ** -------------------------------------------------- ** *step, name=step-5 SSD subspace *steady state dynamics, subspace projection, interval=eigenfrequency 60., *cload load, 2, 100. *end step ** ** -------------------------------------------------- ** *step, name=step-6 random response *random response 60., *select eigenmodes, generate 1, 5 *cload, load case=1 load, 2, 100. *cload, load case=2 load, 2, 100. *correlation, psd=white, type=uncorrelated, complex=no 1,1 2,1 *end step ** ** -------------------------------------------------- ** *step, name=step-7 response spectrum *response spectrum, sum=abs, comp=alg spec,0.,1.,0.,1. *select eigenmodes, generate 1, 5 *end step ** ** -------------------------------------------------- ** *step, nlgeom=yes, name=step-8 static *static 1.0, 1.0 *Boundary load, 1, 1, 0.025 load, 2, 2, 0.050 *El print, freq=999 s, *Node print, freq=999 u, *Output, field, op=new *Node output u, *Element output s, *end step ** ** -------------------------------------------------- ** *step, perturbation, name=step-9 static *static *El print, freq=999 s, *Node print, freq=0 *Output, field, op=new *Element output s, *load case, name=load1 *dload eall, bx, 1.e3 *end load case *load case, name=load2 *dload eall, grav, 1.e6, 0.,1.,0. *end load case *load case, name=load3 *dload eall, cent, 100.,0.,0.,0.,1.,0. *end load case *load case, name=load4 *dload eall, centrif, 1.e6,0.,0.,0.,1.,0. *end load case *end step