# HEAD Example of 3-point testing (short section). Units: MN and m #--------------------------------------------------------------------------- # Control data # Maxit Ndim Isolvr Npoint Iprint Conr Gac Istres CONTROL 500 3 2 16 01 1.e-9 0.0 stressfree #CONTROL 500 3 2 16 01 1.e-9 0.0 restart 561 # Matat Alpha1 Alpha2 Alpha DYNCONT 1 0.0 0.09 -0.05 # T Dt Dtvi Dt0 Type StepType Iterco Itcrit Maxit Maxdiv Conr TIMECO 0.5 0.5 0.5 10.0 STATIC auto NONE all 30 5 1e-5 TIMECO 1.0 0.025 0.5 10.0 dynamic auto none all 15 5 1E-5 #TIMECO 2.0 0.025 0.5 10.0 dynamic auto none all 15 5 1E-5 #TIMECO 30.0 0.01 0.25 10.0 dynamic auto none all 15 5 1E-5 # Mode Factor Result Visres Integration 1 Sigma-xx-ax Sigma-xx #--------------------------------------------------------------------------- ## Nocoor input # #The centre # Type Node Xcor Ycor Zcor Nocoor Coordinates 1 -0.6 0 0 41 0.6 0 0 #The core radial nodes # Type Node Xcor Ycor Zcor Nocoor Coordinates 4001 -0.6 0 0 4040 0.6 0 0 #The tape radial nodes # Type Node Xcor Ycor Zcor Nocoor Coordinates 5001 -0.6 0 0 5040 0.6 0 0 # #The outer layer radial nodes # Type Node Xcor Ycor Zcor Nocoor Coordinates 6001 -0.6 0 0 6040 0.6 0 0 # # # Tensile layer 1 # Type X0 Y0 Z0 Beta1 Beta2 Beta3 R Node Xcor Theta Nocoor Polar 0.0 0.0 0.0 0.0 0.0 0.0 70.3e-3 10001 -0.6 0.6*11.11 10041 0.6 -0.6*11.11 # N Nodeinc Xinc Thetainc Repeat 16 41 0.0 2*3.14159/16 # # Tensile layer 2 # Type X0 Y0 Z0 Beta1 Beta2 Beta3 R Node Xcor Theta Nocoor Polar 0.0 0.0 0.0 0.0 0.0 0.0 74.3e-3 20001 -0.6 -0.6*10.52 20041 0.6 0.6*10.52 # N Nodeinc Xinc Thetainc Repeat 16 41 0.0 2*3.14159/16 # #------------------------------------------------------------------------ # Creating additionl nodes at identical coordinates as start and end of the helices, # keep fixed in all dofs # Will later be used for attachment of springs toward helix nodes to set up "end-fittings" # Nocoor Polar 0.0 0.0 0.0 0.0 0.0 0.0 70.3e-3 100101 -0.6 0.6*11.11 Repeat 16 1 0.0 2*3.14159/16 Nocoor Polar 0.0 0.0 0.0 0.0 0.0 0.0 70.3e-3 100201 0.6 -0.6*11.11 Repeat 16 1 0.0 2*3.14159/16 # Nocoor Polar 0.0 0.0 0.0 0.0 0.0 0.0 74.3e-3 100301 -0.6 -0.6*10.52 Repeat 16 1 0.0 2*3.14159/16 Nocoor Polar 0.0 0.0 0.0 0.0 0.0 0.0 74.3e-3 100401 0.6 0.6*10.52 Repeat 16 1 0.0 2*3.14159/16 # BONCON GLOBAL 100101 1 REPEAT 16 1 BONCON GLOBAL 100101 2 REPEAT 16 1 BONCON GLOBAL 100101 3 REPEAT 16 1 BONCON GLOBAL 100101 4 REPEAT 16 1 BONCON GLOBAL 100101 5 REPEAT 16 1 BONCON GLOBAL 100101 6 REPEAT 16 1 # BONCON GLOBAL 100201 1 REPEAT 16 1 BONCON GLOBAL 100201 2 REPEAT 16 1 BONCON GLOBAL 100201 3 REPEAT 16 1 BONCON GLOBAL 100201 4 REPEAT 16 1 BONCON GLOBAL 100201 5 REPEAT 16 1 BONCON GLOBAL 100201 6 REPEAT 16 1 # BONCON GLOBAL 100301 1 REPEAT 16 1 BONCON GLOBAL 100301 2 REPEAT 16 1 BONCON GLOBAL 100301 3 REPEAT 16 1 BONCON GLOBAL 100301 4 REPEAT 16 1 BONCON GLOBAL 100301 5 REPEAT 16 1 BONCON GLOBAL 100301 6 REPEAT 16 1 # BONCON GLOBAL 100401 1 REPEAT 16 1 BONCON GLOBAL 100401 2 REPEAT 16 1 BONCON GLOBAL 100401 3 REPEAT 16 1 BONCON GLOBAL 100401 4 REPEAT 16 1 BONCON GLOBAL 100401 5 REPEAT 16 1 BONCON GLOBAL 100401 6 REPEAT 16 1 # #------------------------------------------------------------------------ ### The core layer # Elgroup Elty Matname Elid Node1 Node2 Node3 Node4 Elcon core hshear363 plast_PVDF 3001 1 2 4001 # N Nelinc Nodinc Repeat 40 1 1 # # Elgroup Elty Matname Elid Node1 Node2 Node3 Node4 Elcon pressarm hshear363 steel 4001 1 2 4001 # N Nelinc Nodinc Repeat 40 1 1 ## Tensile Layers # # Tensile Layer 1 # Elgroup Elty Flexcrossname Elid Node1 Node2 Node3 Node4 Elcon tenslayer1 hshear353 steel 10001 1 2 10001 10002 10016 1 2 10616 10617 # N Nelinc Nodinc Repeat 40 16 1 ### The tape layer # Elgroup Elty Matname Elid Node1 Node2 Node3 Node4 Elcon tape hshear363 tapeMat 5001 1 2 5001 # N Nelinc Nodinc Repeat 40 1 1 # # Tensile Layer 2 # Elgroup Elty Flexcrossname Elid Node1 Node2 Node3 Node4 Elcon tenslayer2 hshear353 steel 20001 1 2 20001 20002 20016 1 2 20616 20617 # N Nelinc Nodinc Repeat 40 16 1 # ### outer layer # Elgroup Elty Matname Elid Node1 Node2 Node3 Node4 Elcon outer hshear363 plast_PVDF 6001 1 2 6001 # N Nelinc Nodinc Repeat 40 1 1 # ### outer layer - dummy only for contact with rollers # Elgroup Elty Matname Elid Node1 Node2 Node3 Node4 Elcon couter pipe31 dummymat 7001 1 2 # N Nelinc Nodinc Repeat 40 1 1 # ### Contact Layers # # Contact Layer1-between core and tensile layer 1 # Elgroup Elty Flexcrossname Elid Node1 Node2 Node3 Node4 Elcon contlayer1 hcont463 contmat 30001 4001 10001 10002 30016 4001 10616 10617 # N Nelinc Nodinc repeat 40 16 1 # Contact Layer1-between tensile layer 1 and tape # Elgroup Elty Flexcrossname Elid Node1 Node2 Node3 Node4 Elcon contlayer2 hcont463 contmat 40001 10001 10002 5001 40016 10616 10617 5001 # N Nelinc Nodinc repeat 40 16 1 # # Contact Layer 2-between tape and tensile layer 2 # Elgroup Elty Flexcrossname Elid Node1 Node2 Node3 Node4 Elcon contlayer3 hcont463 contmat 50001 5001 20001 20002 50016 5001 20616 20617 # N Nelinc Nodinc repeat 40 16 1 # Contact Layer 2-between outer layer and tensile layer 2 # Elgroup Elty Flexcrossname Elid Node1 Node2 Node3 Node4 Elcon contlayer4 hcont463 contmat 60001 20001 20002 6001 60016 20616 20617 6001 # N Nelinc Nodinc repeat 40 16 1 # #--------------------------------------------------------------------------- # Orient input # # The core layer # Type Elid X Y Z Elorient Coordinates 4001 0 1e3 0 4040 0 1e3 0 # The pressarm layer # Type Elid X Y Z Elorient Coordinates 3001 0 1e3 0 3040 0 1e3 0 # The tape layer # Type Elid X Y Z Elorient Coordinates 5001 0 1e3 0 5040 0 1e3 0 # The outer layer # Type Elid X Y Z Elorient Coordinates 6001 0 1e3 0 6040 0 1e3 0 Elorient Coordinates 7001 0 1e3 0 7040 0 1e3 0 # # Tensile Layer 1 # Type Elid X Y Z Elorient Coordinates 10001 0 1e3 0 10016 0 1e3 0 # N Inc Xinc Yinc Zinc repeat 40 16 0 0 0 # # Tensile Layer 2 # Type Elid X Y Z Elorient Coordinates 20001 0 1e3 0 20016 0 1e3 0 # N Inc Xinc Yinc Zinc repeat 40 16 0 0 0 # # Contact Layer1 # Type Elid X Y Z Elorient Eulerangle 30001 0 0 0 30640 0 0 0 # # Type Elid X Y Z Elorient Eulerangle 40001 0 0 0 40640 0 0 0 # # Contact Layer2 # Type Elid X Y Z Elorient Eulerangle 50001 0 0 0 50640 0 0 0 # Type Elid X Y Z Elorient Eulerangle 60001 0 0 0 60640 0 0 0 # ## Contact interfaces # GroupName MasterName SlaveNamei IS1 ISn IStx ISty IStz MaxIT IGap # 1 means that friction is independent (tape between layers) else (isotropic model) CONTINT contlayer1 core tenslayer1 1 3 5.0 0 0 60 0 CONTINT contlayer2 tenslayer1 tape 1 3 5.0 0 0 60 0 CONTINT contlayer3 tape tenslayer2 1 3 5.0 0 0 60 0 CONTINT contlayer4 tenslayer2 outer 1 3 5.0 0 0 60 0 # #------------------------------------------------------------------ # # Element property input # # =3 turn off both # ELPROP ELGRP SHEARHELIX GEOTYPE W TH WD WS SCALEFACT [PHIST GHIST AXISYM] ELPROP core shearhelix tube 58.35e-3 5.1e-3 0.2e-6 0.0 1.0 # # ELPROP ELGRP SHEAR2HELIX GEOTYPE W TH WD WS RAD ALFA SCALEFACT [PHIST GHIST ] ELPROP pressarm shear2helix rectangle 15.38e-3 5.1512e-3 0.2e-6 0.0 64.10e-3 1.532626 1.0 # # ELPROP ELGRP SHEARHELIX GEOTYPE W TH WD WS SCALEFACT [PHIST GHIST AXISYM] ELPROP tenslayer1 shearhelix rectangle 5e-3 2e-3 0.0785e-6 0.0 3.81 300 300 -10 ELPROP tape shearhelix tube 72.30e-3 2.0e-3 0.45e-6 0.0 1.0 ELPROP tenslayer2 shearhelix rectangle 5e-3 2e-3 0.0785e-6 0.0 4.06 300 300 -10 # ELPROP ELGRP SHEAR2HELIX GEOTYPE W TH WD WS RAD ALFA SCALEFACT [PHIST GHIST ] ELPROP outer shear2helix rectangle 15.38e-3 2e-3 0.2e-6 0.0 76.3e-3 1.532626 1.0 # name eptyp rad th CDr Cdt CMr CMt wd ws ODp ODw rks ELPROP couter pipe 76.3e-3 2e-3 0.0 0.0 0.0 0.0 0.02e-6 0.0 2*76.3e-3 2*76.3e-3 0.5 # # EleGroupName EleType Gap0 TuneTime AUTOMNPC AutoSearch Scalefact ELPROP contlayer1 layercontact D D D D 3.81 ELPROP contlayer2 layercontact D D D D 3.81 ELPROP contlayer3 layercontact D D D D 4.06 ELPROP contlayer4 layercontact D D D D 4.06 #----------------------------------------------------------------------- # Boundary condition data relating to the pipe internal structure # # Type NodId DOF BONCON GLOBAL 4001 2 repeat 40 1 BONCON GLOBAL 4001 3 repeat 40 1 BONCON GLOBAL 5001 2 repeat 40 1 BONCON GLOBAL 5001 3 repeat 40 1 BONCON GLOBAL 6001 2 repeat 40 1 BONCON GLOBAL 6001 3 repeat 40 1 # BONCON GLOBAL 10001 4 repeat 656 1 BONCON GLOBAL 20001 4 repeat 656 1 # #----------------------------------------------- # Constraint input - no ovalization of tensile wires CONSTR CONEQ GLOBAL 10002 3 0.0 10001 3 1.0 repeat 655 1 0 CONSTR CONEQ GLOBAL 20002 3 0.0 20001 3 1.0 repeat 655 1 0 # #----------------------------------------------------------------------- # Boundary condition data relating to the global behaviiour # # Type NodId DOF # end 1 # # mid node # # end 2 # #----------------------------------------------------------- #Load input data # cload applied at the same node as strong spring springmat1. # vertical displacement will be governed by the spring stiffness # # Hist Dir NODE Load CLOAD 100 3 21 1.5e7 # Preshist Gravhist PELOAD 300 300 # Hist Elnr1 P1 Elnr2 P2 PILOAD 200 3001 20 3040 20 # #----------------------------------------------------------- # History data # cload # NO T1 T2 FAC1 THIST_R 100 0.0 5.0 rampcos 0.0 5.0 10.0 rampcos 1.0 10.0 20.0 rampcos -1.0 20.0 30.0 rampcos 1.0 # # internal pressure # NO T1 FAC1 THIST 200 0 0.1 4 1.0 # gravity THIST 300 0 1.0 200 1.0 # #------------------------------------------------------------------------------ # Material data # Matname Mtype Poiss Density Talfa Tecond Heatc EM GM Etrans MATERIAL steel elastic 0.3 7850e-6 11.7e-6 50 800 2.1e5 8.076e4 2.1e5 # # Matname Mtype Poiss Talfa Tecond Heatc Beta EA EIY EIZ GIT EM GM Density Etrans MATERIAL dummymat linear 0.0 11.7e-6 50 800 0 1e-10 1e-10 1e-10 1e-10 1 1 MATERIAL tapeMat elastic 0.35 1E-6 11.7e-6 2.0 50 179 MATERIAL plast_PVDF elastic 0.33 1E-6 11.7e-6 2.0 50 869 # # Contact with viscous effects # tag::fricvisc[] #MATERIAL MNAME FRICVISC KSTICK KTRANS KSURF MU C1 C2 A1 #MATERIAL contmat fricvisc 1e3 1e3 1.4e3 0.0 6.70E1 0.0 1.0 # end::fricvisc[] # # Contact with coloumb friction # tag::coulomb[] # Mname Type MuXY XYName ZName #MATERIAL contmat isocontact 0.2 bellx bellz # # name type alfa eps sig #MATERIAL bellx epcurve 1 0 0 # 0.0001 1.0 # 1000 1.02 #MATERIAL bellz hycurve -1000 -1.4e12 # 1000 1.4e12 # end::coulomb[] # # Contact with new friction model, constant slope # 12 EA sin(alfa)^2/R^2 = 2000 # # tag::friction[] #MATERIAL MNAME FRICONTACT TYPE MUS MUD KEL KEL2 RATIO KDYN C1 C2 DIM kz MATERIAL contmat fricontact 6 0.2 0.2 2e3 2e3 2e-15 10000.0E-6 1.0 1.0 -2 5e5 # end::friction[] #--------------------------------------------------------------------------- # test rig 2 pairs of rollers + spring for displacement in the middle #--------------------------------------------------------------------------- # rollers # tag::testrignodes[] NOCOOR COORDINATES 100000 0.00000 0.00000 0.00000 NOCOOR COORDINATES 100001 -0.40000 0.00000 0.12800 NOCOOR COORDINATES 100002 -0.40000 0.00000 -0.12800 NOCOOR COORDINATES 100003 0.40000 0.00000 0.12800 NOCOOR COORDINATES 100004 0.40000 0.00000 -0.12800 # BONCON GLOBAL 100000 1 REPEAT 5 1 BONCON GLOBAL 100000 2 REPEAT 5 1 BONCON GLOBAL 100000 3 REPEAT 5 1 BONCON GLOBAL 100000 4 REPEAT 5 1 BONCON GLOBAL 100000 5 REPEAT 5 1 BONCON GLOBAL 100000 6 REPEAT 5 1 # end::testrignodes[] # # tag::rollers[] ELCON spool01_01 cont164 rollermat1 100001 100001 ELCON spool01_02 cont164 rollermat1 100002 100002 ELCON spool02_01 cont164 rollermat1 100003 100003 ELCON spool02_02 cont164 rollermat1 100004 100004 # ELORIENT EULERANGLES 100001 0 0 0 REPEAT 4 1 0 0 0 # ELECC radius 100001 1 2 0.0 -0.3 0.0 0.0 0.3 0.0 ELECC radius 100002 1 2 0.0 -0.3 0.0 0.0 0.3 0.0 ELECC radius 100003 1 2 0.0 -0.3 0.0 0.0 0.3 0.0 ELECC radius 100004 1 2 0.0 -0.3 0.0 0.0 0.3 0.0 # # diam ELPROP spool01_01 roller 0.1 ELPROP spool01_02 roller 0.1 ELPROP spool02_01 roller 0.1 ELPROP spool02_02 roller 0.1 # CONTINT spool01_01 spool01_01 couter 7001 7040 1000.0 1000.0 0.0 40 -1 CONTINT spool02_01 spool02_01 couter 7001 7040 1000.0 1000.0 0.0 40 -1 CONTINT spool01_02 spool01_02 couter 7001 7040 1000.0 1000.0 0.0 40 -1 CONTINT spool02_02 spool02_02 couter 7001 7040 1000.0 1000.0 0.0 40 -1 # MATERIAl rollx hycurve -10 -20e3 10 20e3 # MATERIAl rolly hycurve -10 -20e3 10 20e3 # MATERIAl rollz1 hycurve -1 -1e3 1 1e3 # MATERIAL rollermat1 contact 0.3 0.3 rollx rolly rollz1 # end::rollers[] # # spring for application of displacement in pipe center # tag::loadspring[] ELCON lspring1 spring137 lspringmat 100000 100000 21 # ELORIENT eulerangle 100000 0.0 0.0 0.0 ELPROP lspring1 genspring 0.0 0.0 0.0 0.0 0.0 0.0 # MATERIAL lspringmat genspring springmat1 springmat1 springmat1 springmat1 springmat1 springmat1 MATERIAl springmat1 hycurve -10 -1e10 10 1e10 MATERIAL dummy hycurve -10 -0.1 10 0.1 # end::loadspring[] # # tag::endfittings[] ELCON endspring11 spring137 lspringmat2 100101 10001 100101 100116 10616 100116 ELCON endspring12 spring137 lspringmat2 100201 10041 100201 100216 10656 100216 ELCON endspring21 spring137 lspringmat2 100301 20001 100301 100316 20616 100316 ELCON endspring22 spring137 lspringmat2 100401 20041 100401 100416 20656 100416 ELORIENT eulerangle 100101 0.0 0.0 0.0 100116 0.0 0.0 0.0 ELORIENT eulerangle 100201 0.0 0.0 0.0 100216 0.0 0.0 0.0 ELORIENT eulerangle 100301 0.0 0.0 0.0 100316 0.0 0.0 0.0 ELORIENT eulerangle 100401 0.0 0.0 0.0 100416 0.0 0.0 0.0 ELPROP endspring11 genspring 0.0 0.0 1000.0 1000.0 0.0 0.0 1 ELPROP endspring12 genspring 0.0 0.0 1000.0 1000.0 0.0 0.0 1 ELPROP endspring21 genspring 0.0 0.0 1000.0 1000.0 0.0 0.0 1 ELPROP endspring22 genspring 0.0 0.0 1000.0 1000.0 0.0 0.0 1 # MATERIAL lspringmat2 genspring springmat2 springmat2 springmat2 springmat2 springmat2 springmat2 # MATERIAl springmat2 hycurve -10 -1e8 10 1e8 # # end::endfittings[]