** connector friction tests for JOIN+ROTATION, JOIN+CARDAN, JOIN+Flex ** -CUSTOMIZED + PRESTRESS ** friction results should be identical between the three connector elements **------------------------------------------------ *parameter rota = 0.8 Rs = 0.20 *preprint,model=yes *Node,nset = all 101,0,0,0 102,0,0,0 201,1,0,0 202,1,0,0 301,2,0,0 302,2,0,0 401,3,0,0 402,3,0,0 *orientation,name=global 1,0,0, 0,1,0 3,0.0 *orientation,name=global2 1,0,0, 0,1,0 3,0.0 **------------------------------------------------- *element,type=mass,elset=mass 10001,101 10002,102 20001,201 20002,202 30001,301 30002,302 40001,401 40002,402 *mass,elset=mass 1.0, *element,type=rotaryI,elset=rotaryI1 110001,101 110002,102 *rotary inertia,elset=rotaryI1 1.0,1.1,1.2 *element,type=rotaryI,elset=rotaryI2 120001,201 120002,202 *rotary inertia,elset=rotaryI2 2.0,2.1,2.2 *element,type=rotaryI,elset=rotaryI3 130001,301 130002,302 *rotary inertia,elset=rotaryI3 3.0,3.1,3.2 *element,type=rotaryI,elset=rotaryI4 140001,401 140002,402 *rotary inertia,elset=rotaryI4 4.0,4.1,4.2 **------------------------------------------------- *element,type=conn3d2,elset=ballCustomRota 1001,101,102 *connector section,elset=ballCustomRota,behavior=Rota1 join,rotation global,global2 **-------------- *element,type=conn3d2,elset=ballCustomCard 2001,201,202 *connector section,elset=ballCustomCard,behavior=Rota2 join,cardan global,global2 **-------------- *element,type=conn3d2,elset=ballCustomFlex 3001,301,302 *connector section,elset=ballCustomFlex,behavior=Rota3 join,flexion global,global2 **-------------- *element,type=conn3d2,elset=ballCustomEuler 4001,401,402 *connector section,elset=ballCustomEuler,behavior=Rota4 join,euler global,global2 **-------------- *connector behavior,name=Rota1 *connector friction 1.25 *connector potential 4, 5, 6, *friction,elastic slip=1.d-8 0.1, *connector behavior,name=Rota2 *connector friction 1.25 *connector potential 4, 5, 6, *friction,elastic slip=1.d-8 0.1, *connector behavior,name=Rota3 *connector friction 1.25 *connector potential 4, 5, 6, *friction,elastic slip=1.d-8 0.1, *connector behavior,name=Rota4 *connector friction 1.25 *connector potential 4, 5, 6, *friction,elastic slip=1.d-8 0.1, **-------------------------------------------------- *boundary 101,1,6 201,1,6 301,1,6 401,1,6 *elset,elset=conn ballCustomRota,ballCustomCard, ballCustomFlex,ballCustomEuler *PSD-DEFINITION, NAME=WHITE, TYPE=FORCE 1.0,0.0,0.0 1.0,0.0,60.0 *SPECTRUM, TYPE=DISPLACEMENT, NAME=RESPD, INPUT=conn_quake_dis.inp **-------------------------------------------------- *step,nlgeom *static 0.01,1.0,0.01,0.01 *boundary 102,4,6, 202,4,6, 302,4,6, 402,4,6, *OUTPUT, FIELD, FREQ=0 *OUTPUT,HISTORY,freq=2 *NODE OUTPUT,nset=all U,RF *ELEMENT OUTPUT,elset=conn CUR1,CUR2,CUR3, CASUR1,CASUR2,CASUR3,CASUC, CSM1,CSM2,CSFC, CNFC, *node print,freq=1000 u, rf, *el file, elset=conn,freq=1000 CU,CASU,CEF,CSF,CTF,CNF,CSF *el print,freq=100,sum=no, tot=no,elset=conn CU, CASU, CEF, CTF, CNF, CSF, CIVC, *end step **-------------------------------------------------- *step *frequency 12,,,-1 *boundary, op=new 101,1,6 201,1,6 301,1,6 401,1,6 *OUTPUT, FIELD, freq=0 *NODE OUTPUT U, *el print,sum=no,tot=no,elset=conn CUR1,CUR2,CUR3, CSM1,CSM2,CSM3,CSFC *el file,elset=conn CSF, *end step *** **-------------------------------------------------- *step *RANDOM RESPONSE 12.0,60.0,20,3,0 *SELECT EIGENMODES,generate 1,12,1 *MODAL DAMPING, MODAL=DIRECT 1,12,0.02 *CONNECTOR LOAD, LOAD CASE=1 conn,6,-1.0 *CONNECTOR LOAD, LOAD CASE=2 conn,6,-1.0 ** *CORRELATION, PSD=WHITE, TYPE=UNCORRELATED, COMPLEX=NO 1,1 2,1 *OUTPUT, FIELD, FREQ=0 *OUTPUT,HISTORY,freq=2 *NODE OUTPUT,nset=all U1, *ELEMENT OUTPUT,elset=conn CUR1,CUR2,CUR3, CSM1,CSM2,CSM3, CNFC,RCNFC *el print,elset=conn CUR1,CTM1,CTM2,CTM3,CNFC,CSFC,RCNFC,RCSFC *end step **-------------------------------------------------- *STEP *RESPONSE SPECTRUM, COMP=SRSS, SUM=SRSS RESPD, 1.0, 0.0, 0.0, 1.0 ** *MODAL DAMPING, MODAL=DIRECT 1,12,0.02 *SELECT EIGENMODES,generate 1,12,1 *OUTPUT, FIELD, FREQ=0 *OUTPUT,HISTORY,Freq=0 *NODE OUTPUT,nset=all U1, *ELEMENT OUTPUT,elset=conn CUR1,CUR2,CUR3, CSM1,CSM2,CSM3,CSFC, CNFC, ** *el print,elset=conn CUR1,CTM1,CTM2,CTM3,CNFC,CSFC, *END STEP