*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 10, -0.25, 0.25, 0.75 11, -0.5, 0., 0.5 12, -0.25, -0.25, 0.25 13, -0.5, 0., 0. 14, -0.25, 0.25, 0.25 15, -0.5, 0.5, 0.5 16, 0., -0.5, 0. 17, 0.25, -0.25, 0.25 18, 0.5, 0., 0. 19, 0.25, 0.25, 0.25 20, 0., 0., 0. 21, 0., 0.5, 0. 22, -0.5, 0., 1. 23, -0.25, -0.25, 0.75 24, -0.5, -0.5, 0.5 25, 0., -0.5, 0.5 26, 0.25, -0.25, 0.75 27, 0., -0.5, 1. 28, 0., 0., 1. 29, 0.5, -0.5, 0.5 30, 0.5, 0., 0.5 31, 0.5, 0., 1. 32, 0.25, 0.25, 0.75 33, 0.5, 0.5, 0.5 34, 0., 0.5, 1. 35, 0., 0.5, 0.5 *Element, type=C3D10, elset=cube1 1, 9, 6, 8, 7, 12, 11, 10, 14, 13, 15 2, 9, 3, 6, 4, 17, 16, 12, 19, 18, 20 *Solid Section, elset=cube1, material=Material-1, orientation=ori1 1., *Orientation, name=ori1 1.,1.,1., -1.,1.,-1. *Element, type=C3D10, elset=cube2 3, 9, 4, 6, 7, 19, 20, 12, 14, 21, 13 4, 9, 5, 8, 6, 23, 22, 10, 12, 24, 11 *Solid Section, elset=cube2, material=Material-2, orientation=ori2 1., *Orientation, name=ori2 1.,1.,-1., -1.,1.,1. 2, 1.0 *Element, type=C3D10, elset=cube3 5, 9, 3, 5, 6, 17, 25, 23, 12, 16, 24 6, 9, 5, 2, 8, 23, 27, 26, 10, 22, 28 *Solid Section, elset=cube3, material=Material-3, orientation=ori3 1., *Orientation, name=ori3 -1.,1.,-1., 1.,1.,1. 2, 2.0 *Element, type=C3D10, elset=cube4 7, 9, 2, 5, 3, 26, 27, 23, 17, 29, 25 8, 9, 4, 2, 3, 19, 30, 26, 17, 18, 29 *Solid Section, elset=cube4, material=Material-4, orientation=ori4 1., *Orientation, name=ori4 -1.,-1.,1., 1.,-1.,-1. 2, 3.0 *Element, type=C3D10, elset=cube5 9, 9, 1, 2, 4, 32, 31, 26, 19, 33, 30 10, 9, 8, 2, 1, 10, 28, 26, 32, 34, 31 *Solid Section, elset=cube5, material=Material-5, orientation=ori5 1., *Orientation, name=ori5 -1.,1.,1., 1.,1.,-1. 2, -2.0 *Element, type=C3D10, elset=cube6 11, 9, 7, 1, 4, 14, 35, 32, 19, 21, 33 12, 9, 8, 1, 7, 10, 34, 32, 14, 15, 35 *Solid Section, elset=cube6, material=Material-6, orientation=ori6 1., *Orientation, name=ori6 -1.,-1.,-1., 1.,-1.,1. 2, -1.0 *Nset, nset=Fix 1, 2, 3, 5, 6, 16, 24, 25, 27, 29 *Nset, nset=Load 4, 7, 8, 15, 21, 33, 34, 35 ** ** MATERIALS ** *Material, name=Material-1 *Density 0.751E-3 *Elastic, type=ORTHOTROPIC 3.38e6,0.888e6,4.74e6,1.45e6,2.06e6,6.25e6,2.99e5,3.99e5, 4.99e5 *Plastic 2.0e4, 0. 2.2e4, .0009 2.2e4, .0029 *Material, name=Material-2 *Density 0.752E-3 *Elastic, type=ORTHOTROPIC 3.39e6,0.889e6,4.75e6,1.46e6,2.07e6,6.26e6,3.00e5,4.00e5, 5.00e5 *Plastic 2.0e4, 0. 2.2e4, .0009 2.2e4, .0029 *Material, name=Material-3 *Density 0.753E-3 *Elastic, type=ORTHOTROPIC 3.40e6,0.890e6,4.76e6,1.47e6,2.08e6,6.27e6,3.01e5,4.01e5, 5.01e5 *Plastic 2.0e4, 0. 2.2e4, .0009 2.2e4, .0029 *Material, name=Material-4 *Density 0.754E-3 *Elastic, type=ORTHOTROPIC 3.41e6,0.891e6,4.77e6,1.48e6,2.09e6,6.28e6,3.02e5,4.02e5, 5.02e5 *Plastic 2.0e4, 0. 2.2e4, .0009 2.2e4, .0029 *Material, name=Material-5 *Density 0.755E-3 *Elastic, type=ORTHOTROPIC 3.42e6,0.892e6,4.78e6,1.49e6,2.10e6,6.29e6,3.03e5,4.03e5, 5.03e5 *Plastic 2.0e4, 0. 2.2e4, .0009 2.2e4, .0029 *Material, name=Material-6 *Density 0.756E-3 *Elastic, type=ORTHOTROPIC 3.43e6,0.893e6,4.79e6,1.50e6,2.11e6,6.30e6,3.04e5,4.04e5 5.04e5 *Plastic 2.0e4, 0. 2.2e4, .0009 2.2e4, .0029 *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, peeq, *Node print, freq=999 u, *Output, field, op=new *Node output u, *Element output s,peeq *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