# HEAD Example of wire break (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 # Matat Alpha1 Alpha2 Alpha DYNCONT 1 0.0 0.09 -0.05 # # T Dt Dtvi Dt0 Type StepType Iterco Itcrit Maxit Maxdiv Conr #TIMECO 10.0 0.005 0.5 10.0 dynamic auto none all 20 5 1E-5 TIMECO 0.005 0.005 0.005 10.0 dynamic auto none all 20 5 1E-5 # Mode Factor Result Visres Integration 1 Sigma-xx-ax #--------------------------------------------------------------------------- ## 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 # #------------------------------------------------------------------------ ### 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 # ### 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 # # 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 CONTINT contlayer1 core tenslayer1 1 3 4.0 4.0 0 60 0 CONTINT contlayer2 tenslayer1 tape 1 3 4.0 4.0 0 60 0 CONTINT contlayer3 tape tenslayer2 1 3 4.0 4.0 0 60 0 CONTINT contlayer4 tenslayer2 outer 1 3 4.0 4.0 0 60 0 # #------------------------------------------------------------------ # # Element property input # # 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] #tag::scaleproperty[] ELPROP tenslayer1 shearhelix rectangle 5e-3 2e-3 0.0785e-6 0.0 3.81 300 300 -10 ELPROP tenslayer2 shearhelix rectangle 5e-3 2e-3 0.0785e-6 0.0 4.06 300 300 -10 #end::scaleproperty[] ELPROP tape shearhelix tube 72.30e-3 2.0e-3 0.45e-6 0.0 1.0 # 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 # #tag::scalecontact[] # 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 #end::scalecontact[] #----------------------------------------------------------------------- # 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 # #tag::armbonc[] BONCON GLOBAL 10001 1 repeat 16 41 BONCON GLOBAL 10001 2 repeat 16 41 BONCON GLOBAL 10001 3 repeat 16 41 BONCON GLOBAL 10001 4 repeat 16 41 # BONCON GLOBAL 10041 1 repeat 16 41 BONCON GLOBAL 10041 2 repeat 16 41 BONCON GLOBAL 10041 3 repeat 16 41 BONCON GLOBAL 10041 4 repeat 16 41 # BONCON GLOBAL 20001 1 repeat 16 41 BONCON GLOBAL 20001 2 repeat 16 41 BONCON GLOBAL 20001 3 repeat 16 41 BONCON GLOBAL 20001 4 repeat 16 41 # BONCON GLOBAL 20041 1 repeat 16 41 BONCON GLOBAL 20041 2 repeat 16 41 BONCON GLOBAL 20041 3 repeat 16 41 BONCON GLOBAL 20041 4 repeat 16 41 #end::armbonc[] #----------------------------------------------------------------------- # Boundary condition data relating to the global behaviiour # #tag::globon[] # Type NodId DOF # end 1 BONCON GLOBAL 1 1 BONCON GLOBAL 1 2 BONCON GLOBAL 1 3 BONCON GLOBAL 1 4 BONCON GLOBAL 1 5 BONCON GLOBAL 1 6 #end::globon[] # # end 2 # #----------------------------------------------------------- #Load input data # #tag::cload[] # Hist Dir NODE Load CLOAD 100 1 41 0.25 # # Hist Elnr1 P1 Elnr2 P2 PILOAD 200 3001 2 3040 2 # #end::cload[] # # Preshist Gravhist PELOAD 300 300 #tag::break[] # elstart elend hist ELHIST 10321 10321 500 # THIST 500 0.0 1.0 4.0 1.0 6.0 -1.0 10.0 -1.0 #end::break[] #----------------------------------------------------------- # History data # cload # NO T FAC1 THIST 100 0.0 0.0 4.0 0.5 6.0 0.5 10.0 1.0 # # internal pressure # NO T1 FAC1 THIST 200 0 0.01 4 1.0 # gravity THIST 300 0 0.0 200 0.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 new friction model, constant slope # 12 EA sin(alfa)^2/R^2 = 2000 # #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