*HEADING CONN_XPL_FORCE_3D: forced, steady-state oscillation [CONN3D2] Steady-state damped motion: u(t) = p_0/(2 xi k) sin(Omega*t) = 0.48 sin(2t) u(0) = 0.0, v(0) = p_0 Omega/(2 xi k) f(t) = p_0 cos(Omega*t) Omega = w_n = 2.0 xi = c/(2 m w_n) = 5/48 p_0 = k/10 ************************************************* ** create a rigid base for the connector elements ************************************************* *NODE,NSET=PLATE 101, 0., 0., 1. 112, 0., 11., 1. 312, 2., 11., 1. 314, 2., 13., 1. 514, 4., 13., 1. 518, 4., 17., 1. 718, 6., 17., 1. 701, 6., 0., 1. *NGEN,NSET=E1 101,112,1 *NGEN,NSET=E2 112,312,100 *NGEN,NSET=E3 312,314,1 *NGEN,NSET=E4 314,514,100 *NGEN,NSET=E5 514,518,1 *NGEN,NSET=E6 518,718,100 *NGEN,NSET=E7 701,718,1 *NGEN,NSET=E8 101,701,100 *NGEN,NSET=INTERIOR 201,212,1 301,314,1 401,414,1 501,518,1 601,618,1 *ELEMENT,TYPE=S4R,ELSET=PLATE 101, 101,201,202,102 301, 301,401,402,302 501, 501,601,602,502 *ELGEN,ELSET=PLATE 101,11,1,1,2,100,100 301,13,1,1,2,100,100 501,17,1,1,2,100,100 *SHELL SECTION,ELSET=PLATE,MATERIAL=MAT 1.0, *MATERIAL,NAME=MAT *ELASTIC 1.0,0.0 *DENSITY 1.0, *NODE,NSET=REF 1000, 0.,0.,0. *RIGID BODY, REFNODE=1000, ELSET=PLATE ************************************************* *NODE,NSET=TOP 1303, 2., 2., 2. 1311, 2.,10., 2. 1503, 4., 2., 2. 1513, 4.,12., 2. 1703, 6., 2., 2. 1717, 6.,16., 2. *NGEN,NSET=TOP 1303,1311,2 1503,1513,2 1703,1717,2 ******************************************* ** put mass and rotary inertia at top nodes ******************************************* *ELEMENT, TYPE=MASS, ELSET=MASS 1011,1303 1021,1503 1031,1703 *ELGEN, ELSET=MASS 1011,5,2,1 1021,6,2,1 1031,8,2,1 *MASS, ELSET=MASS 12.0 *ELEMENT, TYPE=ROTARYI, ELSET=ROTI 2011,1303 2021,1503 2031,1703 *ELGEN, ELSET=ROTI 2011,5,2,1 2021,6,2,1 2031,8,2,1 *ROTARY INERTIA, ELSET=ROTI 120.0, 120.0, 120.0 ******************************** ** define the connector elements ******************************** *ELEMENT, TYPE=CONN3D2, ELSET=AX 11, 202, 1303 *ELEMENT, TYPE=CONN3D2, ELSET=CAR 12, 204, 1305 *NODE ** Special definition of node 1307 for radial-thrust ** so that the radial direction is properly defined ** 1307 = (2,6,2) + (0,-1,1)/sqrt(2) 1307, 2.0, 5.29289322, 2.70710678 *ELEMENT, TYPE=CONN3D2, ELSET=RTH 13, 206, 1307 *ELEMENT, TYPE=CONN3D2, ELSET=SLP 14, 208, 1309 *ELEMENT, TYPE=CONN3D2, ELSET=SLT 15, 210, 1311 *ELEMENT, TYPE=CONN3D2, ELSET=CDN 21, 402, 1503 *ELEMENT, TYPE=CONN3D2, ELSET=EUL 22, 404, 1505 *ELEMENT, TYPE=CONN3D2, ELSET=FLT 23, 406, 1507 *ELEMENT, TYPE=CONN3D2, ELSET=REV 24, 408, 1509 *ELEMENT, TYPE=CONN3D2, ELSET=ROT 25, 410, 1511 *ELEMENT, TYPE=CONN3D2, ELSET=UNI 26, 412, 1513 *ELEMENT, TYPE=CONN3D2, ELSET=AX-EUL 31, 602, 1703 *ELEMENT, TYPE=CONN3D2, ELSET=CAR-CDN 32, 604, 1705 *ELEMENT, TYPE=CONN3D2, ELSET=CAR-ROT 33, 606, 1707 *ELEMENT, TYPE=CONN3D2, ELSET=JN-FLT 34, 608, 1709 *NODE ** define a special node for link so that bc's can ** be defined to keep the node from moving 1711, 5., 9., 2. *BOUNDARY 1711, 1,2 *ELEMENT, TYPE=CONN3D2, ELSET=LNK-UNI 35, 610, 1711 *NODE ** Special definition of node 1713 for radial-thrust ** so that the radial direction is properly defined ** 1713 = (6,12,2) + (0,-1,1)/sqrt(2) 1713, 6.0, 11.29289322, 2.70710678 *ELEMENT, TYPE=CONN3D2, ELSET=RTH-CV 36, 612, 1713 *ELEMENT, TYPE=CONN3D2, ELSET=SLP-REV 37, 614, 1715 *ELEMENT, TYPE=CONN3D2, ELSET=SLT-ALN 38, 617, 1717 *ORIENTATION, NAME=GLOBAL ** Global directions 1.,0.,0., 0.,1.,0. *ORIENTATION, NAME=LDIRECTION ** orientation used by most connectors 1.,1.,1., 0.,-1.,1. *ORIENTATION, NAME=ORI_RAD-THRUST ** 1-2-3 for radial-thrust is 2-3-1 for others. 0.,-1.,1., 2.,-1.,-1. *ORIENTATION, NAME=ORI_UNI ** 1-2-3 for universal is 2-1-(-3) for others. 0.,-1.,1., 1.,1.,1. *ORIENTATION, NAME=L45 1.,0.,0., 0.,1.,1. ***** *CONNECTOR SECTION, ELSET=AX, BEHAVIOR=CB_AX AXIAL, *CONNECTOR BEHAVIOR, NAME=CB_AX *CONNECTOR ELASTICITY, COMPONENT=1 48.0, *CONNECTOR DAMPING, COMPONENT=1 5.0, *CONNECTOR FRICTION,COMPONENT=1 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=CAR, BEHAVIOR=CB_CAR CARTESIAN, LDIRECTION, *CONNECTOR BEHAVIOR, NAME=CB_CAR *CONNECTOR ELASTICITY, COMPONENT=1 48.0, *CONNECTOR DAMPING, COMPONENT=1 5.0, *CONNECTOR ELASTICITY, COMPONENT=2 48.0, *CONNECTOR DAMPING, COMPONENT=2 5.0, *CONNECTOR ELASTICITY, COMPONENT=3 48.0, *CONNECTOR DAMPING, COMPONENT=3 5.0, *CONNECTOR FRICTION,COMP=1 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=2 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=3 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=RTH, BEHAVIOR=CB_RTH RADIAL-THRUST, ORI_RAD-THRUST, *CONNECTOR BEHAVIOR, NAME=CB_RTH *CONNECTOR ELASTICITY, COMPONENT=1 48.0, *CONNECTOR DAMPING, COMPONENT=1 5.0, *CONNECTOR ELASTICITY, COMPONENT=3 48.0, *CONNECTOR DAMPING, COMPONENT=3 5.0, *CONNECTOR FRICTION,COMP=1 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=3 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=SLP, BEHAVIOR=CB_SLP SLIDE-PLANE, LDIRECTION, *CONNECTOR BEHAVIOR, NAME=CB_SLP *CONNECTOR ELASTICITY, COMPONENT=2 48.0, *CONNECTOR DAMPING, COMPONENT=2 5.0, *CONNECTOR ELASTICITY, COMPONENT=3 48.0, *CONNECTOR DAMPING, COMPONENT=3 5.0, *CONNECTOR FRICTION,COMP=2 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=3 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=SLT, BEHAVIOR=CB_SLT SLOT, LDIRECTION, *CONNECTOR BEHAVIOR, NAME=CB_SLT *CONNECTOR ELASTICITY, COMPONENT=1 48.0, *CONNECTOR DAMPING, COMPONENT=1 5.0, *CONNECTOR FRICTION,COMP=1 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=CDN, BEHAVIOR=CB_CDN ** Make CARDAN behave like UNIVERSAL CARDAN, ORI_UNI, *CONNECTOR BEHAVIOR, NAME=CB_CDN *CONNECTOR ELASTICITY, COMPONENT=4 480.0, *CONNECTOR DAMPING, COMPONENT=4 50.0, *CONNECTOR ELASTICITY, COMPONENT=5 1.0E6, *CONNECTOR ELASTICITY, COMPONENT=6 480.0, *CONNECTOR DAMPING, COMPONENT=6 50.0, *CONNECTOR FRICTION,COMP=4 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=5 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=6 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=EUL, BEHAVIOR=CB_EUL EULER, GLOBAL,L45 *CONNECTOR BEHAVIOR, NAME=CB_EUL *CONNECTOR ELASTICITY, COMPONENT=4 480.0, *CONNECTOR DAMPING, COMPONENT=4 50.0, *CONNECTOR ELASTICITY, COMPONENT=5 480.E2, *CONNECTOR DAMPING, COMPONENT=5 50.E2, *CONNECTOR ELASTICITY, COMPONENT=6 480.E2, *CONNECTOR DAMPING, COMPONENT=6 50.E2, *CONNECTOR FRICTION,COMP=4 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=5 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=6 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=FLT, BEHAVIOR=CB_FLT FLEXION-TORSION, GLOBAL,L45 *CONNECTOR BEHAVIOR, NAME=CB_FLT *CONNECTOR ELASTICITY, COMPONENT=4 480.E2, *CONNECTOR DAMPING, COMPONENT=4 50.E2, *CONNECTOR ELASTICITY, COMPONENT=5 480.0, *CONNECTOR DAMPING, COMPONENT=5 50.0, *CONNECTOR ELASTICITY, COMPONENT=6 480.d5, *CONNECTOR FRICTION,COMP=4 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=5 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=6 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=REV, BEHAVIOR=CB_REV REVOLUTE, LDIRECTION, *CONNECTOR BEHAVIOR, NAME=CB_REV *CONNECTOR ELASTICITY, COMPONENT=4 480.0, *CONNECTOR DAMPING, COMPONENT=4 50.0, *CONNECTOR FRICTION,COMP=4 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=ROT, BEHAVIOR=CB_ROT ROTATION, LDIRECTION, *CONNECTOR BEHAVIOR, NAME=CB_ROT *CONNECTOR ELASTICITY, COMPONENT=4 480.0, *CONNECTOR DAMPING, COMPONENT=4 50.0, *CONNECTOR ELASTICITY, COMPONENT=5 480.0, *CONNECTOR DAMPING, COMPONENT=5 50.0, *CONNECTOR ELASTICITY, COMPONENT=6 480.0, *CONNECTOR DAMPING, COMPONENT=6 50.0, *CONNECTOR FRICTION,COMP=4 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=5 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=6 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=UNI, BEHAVIOR=CB_UNI UNIVERSAL, ORI_UNI, *CONNECTOR BEHAVIOR, NAME=CB_UNI *CONNECTOR ELASTICITY, COMPONENT=4 480.0, *CONNECTOR DAMPING, COMPONENT=4 50.0, *CONNECTOR ELASTICITY, COMPONENT=6 480.0, *CONNECTOR DAMPING, COMPONENT=6 50.0, *CONNECTOR FRICTION,COMP=4 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=6 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=AX-EUL, BEHAVIOR=CB_AX-EUL EULER, AXIAL GLOBAL, L45 *CONNECTOR BEHAVIOR, NAME=CB_AX-EUL *CONNECTOR ELASTICITY, COMPONENT=1 48.0, *CONNECTOR DAMPING, COMPONENT=1 5.0, *CONNECTOR ELASTICITY, COMPONENT=4 480.0, *CONNECTOR DAMPING, COMPONENT=4 50.0, *CONNECTOR ELASTICITY, COMPONENT=5 480.E2, *CONNECTOR DAMPING, COMPONENT=5 50.E2, *CONNECTOR ELASTICITY, COMPONENT=6 480.E2, *CONNECTOR DAMPING, COMPONENT=6 50.E2, *CONNECTOR FRICTION,COMP=1 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=4 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=5 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=6 10.0, *friction 0.2, ***** *CONNECTOR SECTION, ELSET=CAR-CDN, BEHAVIOR=CB_CAR-CDN CARTESIAN,CARDAN ORI_UNI, *CONNECTOR BEHAVIOR, NAME=CB_CAR-CDN *CONNECTOR ELASTICITY, COMPONENT=1 48.0, *CONNECTOR DAMPING, COMPONENT=1 5.0, *CONNECTOR ELASTICITY, COMPONENT=2 48.0, *CONNECTOR DAMPING, COMPONENT=2 5.0, *CONNECTOR ELASTICITY, COMPONENT=3 48.0, *CONNECTOR DAMPING, COMPONENT=3 5.0, *CONNECTOR ELASTICITY, COMPONENT=4 480.0, *CONNECTOR DAMPING, COMPONENT=4 50.0, *CONNECTOR ELASTICITY, COMPONENT=5 1.0E6, *CONNECTOR ELASTICITY, COMPONENT=6 480.0, *CONNECTOR DAMPING, COMPONENT=6 50.0, *CONNECTOR FRICTION,COMP=1 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=2 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=3 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=4 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=5 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=6 10.0, *friction 0.2, ***** *CONNECTOR SECTION, ELSET=CAR-ROT, BEHAVIOR=CB_CAR-ROT CARTESIAN,ROTATION LDIRECTION, *CONNECTOR BEHAVIOR, NAME=CB_CAR-ROT *CONNECTOR ELASTICITY, COMPONENT=1 48.0, *CONNECTOR DAMPING, COMPONENT=1 5.0, *CONNECTOR ELASTICITY, COMPONENT=2 48.0, *CONNECTOR DAMPING, COMPONENT=2 5.0, *CONNECTOR ELASTICITY, COMPONENT=3 48.0, *CONNECTOR DAMPING, COMPONENT=3 5.0, *CONNECTOR ELASTICITY, COMPONENT=4 480.0, *CONNECTOR DAMPING, COMPONENT=4 50.0, *CONNECTOR ELASTICITY, COMPONENT=5 480.0, *CONNECTOR DAMPING, COMPONENT=5 50.0, *CONNECTOR ELASTICITY, COMPONENT=6 480.0, *CONNECTOR DAMPING, COMPONENT=6 50.0, *CONNECTOR FRICTION,COMP=1 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=2 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=3 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=4 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=5 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=6 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=JN-FLT, BEHAVIOR=CB_JN_FLT JOIN, FLEXION-TORSION GLOBAL,L45 *CONNECTOR BEHAVIOR, NAME=CB_JN_FLT *CONNECTOR ELASTICITY, COMPONENT=4 480.E5, *CONNECTOR DAMPING, COMPONENT=4 50.E2, *CONNECTOR ELASTICITY, COMPONENT=5 480.0, *CONNECTOR DAMPING, COMPONENT=5 50.0, *CONNECTOR ELASTICITY, COMPONENT=6 480.d5, *CONNECTOR FRICTION,COMP=4 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=5 10.0, *friction 0.2, *CONNECTOR FRICTION,COMP=6 10.0, *friction 0.2, ***** *CONNECTOR SECTION, ELSET=LNK-UNI, BEHAVIOR=CB_LNK-UNI LINK,UNIVERSAL ORI_UNI, *CONNECTOR BEHAVIOR, NAME=CB_LNK-UNI *CONNECTOR ELASTICITY, COMPONENT=4 480.0, *CONNECTOR DAMPING, COMPONENT=4 50.0, *CONNECTOR ELASTICITY, COMPONENT=6 480.0, *CONNECTOR DAMPING, COMPONENT=6 50.0, *CONNECTOR FRICTION,COMP=4 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=6 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=RTH-CV, BEHAVIOR=CB_RTH_CV RADIAL-THRUST,CONSTANT VELOCITY ORI_RAD-THRUST, *CONNECTOR BEHAVIOR, NAME=CB_RTH_CV *CONNECTOR ELASTICITY, COMPONENT=1 48.0, *CONNECTOR DAMPING, COMPONENT=1 5.0, *CONNECTOR ELASTICITY, COMPONENT=3 48.0, *CONNECTOR DAMPING, COMPONENT=3 5.0, *CONNECTOR FRICTION,COMP=1 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=3 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=SLP-REV, BEHAVIOR=CB_SLP-REV SLIDE-PLANE,REVOLUTE LDIRECTION, *CONNECTOR BEHAVIOR, NAME=CB_SLP-REV *CONNECTOR ELASTICITY, COMPONENT=2 48.0, *CONNECTOR DAMPING, COMPONENT=2 5.0, *CONNECTOR ELASTICITY, COMPONENT=3 48.0, *CONNECTOR DAMPING, COMPONENT=3 5.0, *CONNECTOR ELASTICITY, COMPONENT=4 480.0, *CONNECTOR DAMPING, COMPONENT=4 50.0, *CONNECTOR FRICTION,COMP=2 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=3 10.0 *friction 0.2, *CONNECTOR FRICTION,COMP=4 10.0 *friction 0.2, ***** *CONNECTOR SECTION, ELSET=SLT-ALN, BEHAVIOR=CB_SLT-ALN SLOT,ALIGN LDIRECTION, *CONNECTOR BEHAVIOR, NAME=CB_SLT-ALN *CONNECTOR ELASTICITY, COMPONENT=1 48.0, *CONNECTOR DAMPING, COMPONENT=1 5.0, *CONNECTOR FRICTION,COMP=1 10.0 *friction 0.2, ******************************************** ** Define elsets for elements with available ** components in the directions indicated: ******************************************** *ELSET,ELSET=CONN1 AX, CAR, RTH, SLT AX-EUL, CAR-CDN, CAR-ROT, RTH-CV, SLT-ALN *ELSET,ELSET=CONN2 CAR, SLP CAR-CDN, CAR-ROT, SLP-REV *ELSET,ELSET=CONN3 CAR, RTH, SLP CAR-CDN, CAR-ROT, RTH-CV, SLP-REV *ELSET,ELSET=CONN4 CDN, EUL, REV, ROT, UNI AX-EUL, CAR-CDN, CAR-ROT, LNK-UNI, SLP-REV *ELSET,ELSET=CONN5 FLT, ROT CAR-ROT, JN-FLT *ELSET,ELSET=CONN6 CDN, ROT, UNI CAR-CDN, CAR-ROT, LNK-UNI *ELSET,ELSET=CONNALL CONN1,CONN2,CONN3,CONN4,CONN5,CONN6 *BOUNDARY REF, 1,6 *AMPLITUDE, NAME=COSINE, DEF=PERIODIC 1, 2.0, 0.0, 0.0 1.0, 0.0 *NSET,NSET=D100 ** initial velocity in local (1,0,0) direction. ** AXIAL, SLOT, AX-EUL, SLT-ALN 1303, 1311, 1703, 1717 *NSET,NSET=D111 ** initial velocity in local (1,1,1) direction. ** CARTESIAN, CAR-ROT 1305, 1707 *NSET,NSET=D101 ** initial velocity in RADIAL-THRUST local (1,0,1) direction. ** RADIAL-THRUST, RTH-CV 1307, 1713 *NSET,NSET=D011 ** initial velocity in local (0,1,1) direction. ** SLIDE-PLANE, SLP-REV 1309, 1715 *NSET,NSET=DUNI ** initial velocity for ORI_UNI orientation ** CAR-CDN 1705, *NSET,NSET=R100 ** initial angular velocity in local (1,0,0) direction. ** REVOLUTE, SLP-REV 1509, 1715 *NSET,NSET=R101 ** UNIVERSAL, CARDAN, CAR-CDN, LNK-UNI 1503, 1513, 1705, 1711 *NSET,NSET=R111 ** initial angular velocity in local (1,1,1) direction. ** ROTATION, CAR-ROT 1511, 1707 *NSET,NSET=REUL ** initial angular velocity in global (0,0,1) direction. ** EULER, AX-EUL 1505, 1703 *NSET,NSET=RFLT ** initial angular velocity in global (0,-1,1) direction. ** FLT, JN-FLT 1507, 1709 *INITIAL CONDITIONS, TYPE=VELOCITY ** each displacement component has initial velocity 0.96 ** in all the available component of relative motion. D100,1, 0.5542563 D100,2, 0.5542563 D100,3, 0.5542563 ** D111,1, 1.33809298 D111,2,-0.51648461 D111,3, 0.84116041 ** D101,1, 0.5542563 D101,2,-0.12456625 D101,3, 1.23307877 ** D011,1, 0.78383672 D011,2,-1.07074087 D011,3, 0.28690415 ** DUNI,1,-0.22958046 DUNI,2, 0.26735211 DUNI,3, 1.62499713 ** R100,4, 0.5542563 R100,5, 0.5542563 R100,6, 0.5542563 ** R111,4, 1.33809298 R111,5,-0.51648461 R111,6, 0.84116041 ** R101,4, -0.7838367176 R101,5, -0.2869041514 R101,6, 1.070740869 ** REUL,4, 0.0 REUL,5, 0.0 REUL,6, 0.96 ** RFLT,4, 0.0 RFLT,5,-0.67882251 RFLT,6, 0.67882251 ************************************ *ELSET,ELSET=QA_TEST_CONNALL CONNALL, ** *STEP,NLGEOM=YES Force steady-state vibration. Every dof should give the same answer. *DYNAMIC, EXPLICIT, DIRECT 0.0075, 8.0 *CONNECTOR LOAD, AMP=COSINE CONN1,1, 4.8 CONN2,2, 4.8 CONN3,3, 4.8 CONN4,4, 48.0 CONN5,5, 48.0 CONN6,6, 48.0 *FILE OUTPUT, NUM=4 *EL FILE, ELSET=CONNALL CP,CU *OUTPUT,HISTORY, FREQ=1 *ELEMENT OUTPUT, ELSET=CONNALL CP,CU,CV,CA,CSF,CNF *NODE OUTPUT, NSET=TOP U,V,A *OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=2 ***OUTPUT,HISTORY,VARIABLE=PRESELECT *END STEP