** connector friction tests for JOIN+ROTATION, JOIN+CARDAN, JOIN+Flex ** -CUSTOMIZED + PRESTRESS ** friction results should be identical between the three connector elements **------------------------------------------------ *preprint,model=yes *parameter Rs = 0.20 *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 *orientation,name=global 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=rotaryI 110001,101 110002,102 120001,201 120002,202 130001,301 130002,302 140001,401 140002,402 *rotary inertia,elset=rotaryI 1.0,1.0,1.0 **------------------------------------------------- *element,type=conn3d2,elset=ballCustomRota 1001,101,102 *connector section,elset=ballCustomRota,behavior=fricCustomRota join,rotation global, **-------------- *connector behavior,name=fricCustomRota *connector derived component,name=norm123 1,2,3 1.0,1.0,1.0 *connector friction,contact force=norm123 1.25,0.0 1.20,1.0 *Connector potential 4, 5, 6, *friction,elastic slip=1.d-8 0.15, **------ *connector elastic,comp=4 1.0, *connector elastic,comp=5 1.0, *connector elastic,comp=6 1.0, **------------------------------------------------- *element,type=conn3d2,elset=ballCustomCard 2001,201,202 *connector section,elset=ballCustomCard,behavior=fricCustomCard join,cardan global, **-------------- *connector behavior,name=fricCustomCard *connector derived component,name=norm123 1,2,3 1.0,1.0,1.0 *connector friction, contact force=norm123 1.25, 0.0 1.20, 1.0 *Connector potential 4, 5, 6, *friction,elastic slip=1.d-8 0.15, **------ *connector elastic,comp=4 1.0, *connector elastic,comp=5 1.0, *connector elastic,comp=6 1.0, **------------------------------------------------- *element,type=conn3d2,elset=ballCustomFlex 3001,301,302 *connector section,elset=ballCustomFlex,behavior=fricCustomFlex join,flexion global, **-------------- *connector behavior,name=fricCustomFlex *connector derived component,name=norm123 1,2,3 1.0,1.0,1.0 *connector friction, contact force=norm123 1.25, 0.0 1.20, 1.0 *Connector potential 4, 5, 6, *friction,elastic slip=1.d-8 0.15, **------ *connector elastic,comp=4 1.0, *connector elastic,comp=5 1.0, *connector elastic,comp=6 1.0, **-------------------------------------------------- *boundary 101,1,6 201,1,6 301,1,6 *elset,elset=conn ballCustomRota, ballCustomCard, ballCustomFlex *amplitude,name=backforth 0,0, 0.6, 0.6, 1.0,0.2 *amplitude,name=step 0,1,1,1 **-------------------------------------------------- *step,nlgeom *static 0.01,1.0,0.01,0.01 *boundary,amplitude=backforth 102,4,6,1.0 202,4,6,1.0 302,4,6,1.0 *cload,amplitude=step **load that creates constraint forces 102,1,10.0 102,2,20.0 102,3,30.0 **----------- 202,1,10.0 202,2,20.0 202,3,30.0 **----------- 302,1,10.0 302,2,20.0 302,3,30.0 **----------- *OUTPUT, FIELD, FREQ=0 *OUTPUT,HISTORY,freq=10 *NODE OUTPUT,nset=all U,RF *ELEMENT OUTPUT,elset=conn CUR1,CUR2,CUR3, CASUR1,CASUR2,CASUR3,CASUC, CSM1,CSM2,CSM3,CSFC, CTF1,CTF2,CTF3, 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,nlgeom *dynamic 0.01,1.0,0.01,0.01 *boundary,op=new 101,1,6 201,1,6 301,1,6 *cload,amplitude=step **load that creates constraint forces 102,1,10.0 102,2,20.0 102,3,30.0 **----------- 202,1,10.0 202,2,20.0 202,3,30.0 **----------- 302,1,10.0 302,2,20.0 302,3,30.0 **load that creates motion *cload 102,4,10.0 102,5,20.0 102,6,30.0 **----------- 202,4,10.0 202,5,20.0 202,6,30.0 **----------- 302,4,10.0 302,5,20.0 302,6,30.0 **-------------------------------------------------- *OUTPUT, FIELD, FREQ=0 *OUTPUT,HISTORY,freq=10 *NODE OUTPUT,nset=all U,RF *ELEMENT OUTPUT,elset=conn CUR1,CUR2,CUR3, CASUR1,CASUR2,CASUR3,CASUC, CSM1,CSM2,CSM3,CSFC, CTF1,CTF2,CTF3, 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