HEAD 4"WS-FPS2000 p=20MPa RUC Thickw. HS353-364 HC453-463-464 Newfricmodel #--------------------------------------------------------------------------- # Control data # # maxit ndim isolvr npoint ipri conr gacc iproc Control 500 3 2 16 01 1.e-9 0.0 stressfree # # # Matat Alpha1 Alpha2 Alpha Dyncont 1 0.0 0.09 -0.05 # tag::dynres[] dynres_i 1 30015 3 1 dynres_i 1 30030 3 1 # end::dynres[] Visres Integration 1 Sigma-xx-ax Sigma-xx Sigma-yy Sigma-zz sigma-xx-my sigma-xx-mz Vconfor-z Vcondis-z #----------------------------------------------------------------------------------------------------------------- # Nocoor input # # core no x y z # Nocoor Coordinates 1 0.0 0.0 0.0 3 5.6536037E-01*2/61 0.0 0.0 # flexbody no x y z Nocoor Coordinates 1001 5.6536037E-01*2/61/4 0.0 0.0 1002 5.6536037E-01*2/61*3/4 0.0 0.0 1003 5.6536037E-01*2/61/4 0.0 0.0 1004 5.6536037E-01*2/61*3/4 0.0 0.0 # flexbarrier no x y z Nocoor Coordinates 2001 5.6536037E-01*2/61/4 0.0 0.0 2002 5.6536037E-01*2/61*3/4 0.0 0.0 2003 5.6536037E-01*2/61/4 0.0 0.0 2004 5.6536037E-01*2/61*3/4 0.0 0.0 # flexloc no x y z Nocoor Coordinates 3001 5.6536037E-01*2/61/4 0.0 0.0 3002 5.6536037E-01*2/61*3/4 0.0 0.0 3003 5.6536037E-01*2/61/4 0.0 0.0 3004 5.6536037E-01*2/61*3/4 0.0 0.0 # flexwear1 no x y z Nocoor Coordinates 4001 5.6536037E-01*2/61/4 0.0 0.0 4002 5.6536037E-01*2/61*3/4 0.0 0.0 4003 5.6536037E-01*2/61/4 0.0 0.0 4004 5.6536037E-01*2/61*3/4 0.0 0.0 # flextensile1 no x y z # helix node system # flexwear2 no x y z Nocoor Coordinates 6001 5.6536037E-01*2/61/4 0.0 0.0 6002 5.6536037E-01*2/61*3/4 0.0 0.0 6003 5.6536037E-01*2/61/4 0.0 0.0 6004 5.6536037E-01*2/61*3/4 0.0 0.0 # Loc node dir BONCON GLOBAL 6001 1 BONCON GLOBAL 6002 1 BONCON GLOBAL 6003 1 BONCON GLOBAL 6004 1 # flexflextensile2 no x y z # helix node system # flextape1 no x y z Nocoor Coordinates 8001 5.6536037E-01*2/61/4 0.0 0.0 8002 5.6536037E-01*2/61*3/4 0.0 0.0 8003 5.6536037E-01*2/61/4 0.0 0.0 8004 5.6536037E-01*2/61*3/4 0.0 0.0 # flexshield no x y z Nocoor Coordinates 9001 5.6536037E-01*2/61/4 0.0 0.0 9002 5.6536037E-01*2/61*3/4 0.0 0.0 9003 5.6536037E-01*2/61/4 0.0 0.0 9004 5.6536037E-01*2/61*3/4 0.0 0.0 #------------------------------------------------------------------------------------- # # flextensile1 # no x0 y0 z0 b1 b2 b3 R node xcor theta Nocoor Polar 0.0 0.0 0.0 0.0 0.0 0.0 0.0703 30001 0.00 6.28318530/61 30003 5.653603652E-01*2/61 -6.28318530/61 # N Nodeinc Xinc Thetainc Repeat 61 3 0.0 6.28318530/61 # flextensile2 # no x0 y0 z0 b1 b2 b3 R node xcor theta Nocoor Polar 0.0 0.0 0.0 0.0 0.0 0.0 0.0743 40001 0.00 -6.28318530/65 40003 5.653603652E-01*2/65 6.28318530/65 # N Nodeinc Xinc Thetainc Repeat 65 3 0.0 6.28318530/65 # # #------------------------------------------------------------------------ # Elcon input # # # Structural layers # # tag::elcons[] # group elty material no n1 n2 n3 n4 Elcon flexbody hshear364 mypipe 1001 1 2 1001 1003 repeat 2 1 1 Elcon flexbarrier hshear364 mypipe 2001 1 2 2001 2003 repeat 2 1 1 Elcon flexloc hshear364 mypipe 3001 1 2 3001 3003 repeat 2 1 1 Elcon flexwear1 hshear364 mypipe 4001 1 2 4001 4003 repeat 2 1 1 Elcon flextensile1 hshear353t mypipe 30001 1 2 30001 30002 30061 1 2 30181 30182 repeat 2 61 1 Elcon flexwear2 hshear364 mypipe 6001 1 2 6001 6003 repeat 2 1 1 Elcon flextensile2 hshear353t mypipe 40001 1 2 40001 40002 40065 1 2 40193 40194 repeat 2 65 1 Elcon flextape1 hshear364 mypipe 8001 1 2 8001 8003 repeat 2 1 1 Elcon flexshield hshear364 mypipe 9001 1 2 9001 9003 repeat 2 1 1 # end::elcons[] # # Contact layers # # tag::elconc1[] Elcon flexbodycontact hcont464 mypipe 50001 1003 2001 repeat 2 1 1 Elcon flexbarriercontact hcont464 mypipe 51001 2003 3001 repeat 2 1 1 Elcon flexloccontact hcont464 mypipe 52001 3003 4001 repeat 2 1 1 Elcon flexshieldcontact hcont464 mypipe 57001 8003 9001 repeat 2 1 1 # end::elconc1[] # tag::elconc2[] Elcon flexwear1contact hcont463 mypipe 53001 4003 30001 30002 53061 4003 30181 30182 repeat 2 61 1 Elcon flexwear2inwardcontact hcont453 mypipe 54001 30001 30002 40001 40002 54061 30181 30182 40181 40182 repeat 2 61 1 Elcon flextape1contact hcont463 mypipe 56001 40001 40002 8001 56065 40193 40194 8001 repeat 2 65 1 # elprop flexwear1contact shearmodel d d d elprop flexwear2inwardcontact shearmodel d d d elprop flextape1contact shearmodel d d d # end::elconc2[] #--------------------------------------------------------------------------- # Orient input # # structural els # no x y z Elorient Coordinates 1001 0 1e3 0 1002 0 1e3 0 Elorient Coordinates 2001 0 1e3 0 2002 0 1e3 0 Elorient Coordinates 3001 0 1e3 0 3002 0 1e3 0 Elorient Coordinates 4001 0 1e3 0 4002 0 1e3 0 Elorient Coordinates 6001 0 1e3 0 6002 0 1e3 0 Elorient Coordinates 8001 0 1e3 0 8002 0 1e3 0 Elorient Coordinates 9001 0 1e3 0 9002 0 1e3 0 Elorient Coordinates 30001 0 1e3 0 30122 0 1e3 0 Elorient Coordinates 40001 0 1e3 0 40130 0 1e3 0 # # contact elorient eulerangle 50001 0.0 0.0 0.0 50002 0.0 0.0 0.0 elorient eulerangle 51001 0.0 0.0 0.0 51002 0.0 0.0 0.0 elorient eulerangle 52001 0.0 0.0 0.0 52002 0.0 0.0 0.0 elorient eulerangle 53001 0.0 0.0 0.0 53122 0.0 0.0 0.0 elorient eulerangle 54001 0.0 0.0 0.0 54122 0.0 0.0 0.0 #elorient eulerangle 55001 0.0 0.0 0.0 # 55130 0.0 0.0 0.0 elorient eulerangle 56001 0.0 0.0 0.0 56130 0.0 0.0 0.0 elorient eulerangle 57001 0.0 0.0 0.0 57002 0.0 0.0 0.0 #-------------------------------------------------------------------------------- # # tag::flexcross[] # name type ID Timeini ilaint ilaext ielbfl fimod content nelgr CROSSECTION mypipe 353FLEXCROSS 101.6e-3 1.0 2 9 1 0 1000 9 1-1-flexbody-flexbodycontact 2-2-flexbarrier-flexbarriercontact 3-3-flexloc-flexloccontact 4-4-flexwear1-flexwear1contact 5-5-flexwear2inwardcontact-flextensile1 6-6-flexwear2 7-7-flextape1contact-flextensile2 8-8-flextape1-flexshieldcontact 9-9-flexshield #FLEXBODY A=5*12.7*.563 = 36.0, FLEXLOK A=6.4*14.1*.878 = 79.23 # CTYPE TH matname FRIC LAYANG RNUM TEMP nlmat CCODE CFATFL AREA IT INY IKS WIDTH CARC 5.0e-3 steel_316 0.15 87.828 1 0.0 none MANUAL NONE 36.0e-6 0.000e+00 0.000e+00 0.000e+00 7.20e-3 THER 5.1e-3 plast_PVDF 0.15 0.000 0 0.0 none NONE NONE 0.00 0.000e+00 0.000e+00 0.000e+00 0.00 ZETA 6.4e-3 steel_110 0.15 87.813 1 0.0 none MANUAL NONE 79.23e-6 0.000e+00 0.000e+00 0.000e+00 12.38e-3 THER 2.0e-3 plast_PA11 0.15 0.000 0 0.0 none NONE NONE 0.00 0.000e+00 0.000e+00 0.000e+00 0.00 # 1st armour TENS 2.00e-3 steel_190 0.15 -38 61 0.0 none FLEXTENSILE NONE 0.00 0.000e+00 0.000e+00 0.000e+00 0.00 THER 2.00e-3 plast_PA11 0.15 0.000 0 0.0 none NONE NONE 0.00 0.000e+00 0.000e+00 0.000e+00 0.00 # 2nd armour TENS 2.00e-3 steel_190 0.15 37.8 65 0.0 none FLEXTENSILE NONE 0.00 0.000e+00 0.000e+00 0.000e+00 0.00 THER 0.3e-3 rubber 0.15 0.000 0 0.0 none NONE NONE 0.00 0.000e+00 0.000e+00 0.000e+00 0.00 THER 6.0e-3 rubber 0.15 0.000 0 0.0 none NONE NONE 0.00 0.000e+00 0.000e+00 0.000e+00 0.00 # end::flexcross[] # #CROSS-SECTION BOUNDARY DATA # NAME type X0 Y0 CCURV P1 P2 P3 P4 NINTER ICODE # CROSSGEOM TENS-FLEXTENSILE BFLEX 0 0 S 5.0e-3 0.0000 0.0000 0.0000 30 1 S 2.0e-3 90.0000 0.0000 0.0000 20 1 S 5.0e-3 180.0000 0.0000 0.0000 30 1 S 2.0e-3 270.0000 0.0000 0.0000 20 1 # # #-------------------------------------------------------------------- # # name type diameter inside #ELPROP bscontact bellmouth 300.0 1 # #----------------------------------------------------------------------- # Boundary condition data # Loc node dir BONCON GLOBAL 2 1 BONCON GLOBAL 2 2 BONCON GLOBAL 2 3 BONCON GLOBAL 2 4 BONCON GLOBAL 2 5 BONCON GLOBAL 2 6 #------------------------------------------not for bending CONSTR CONEQ GLOBAL 1 1 0.0 3 1 -1.0 #------------------------------------------------------------------------------------------- # fix the relative disp at ends - now we are using the thickwalled option and dof 4 must be free #BONCON gLObAL 30001 2 repeat 183 1 #BONCON gLObAL 40001 2 repeat 195 1 #--------------------------------------------------------------- BONCON gLObAL 1001 2 repeat 4 1 BONCON gLObAL 1001 3 repeat 4 1 BONCON gLObAL 2001 2 repeat 4 1 BONCON gLObAL 2001 3 repeat 4 1 BONCON gLObAL 3001 2 repeat 4 1 BONCON gLObAL 3001 3 repeat 4 1 BONCON gLObAL 4001 2 repeat 4 1 BONCON gLObAL 4001 3 repeat 4 1 BONCON gLObAL 6001 2 repeat 4 1 BONCON gLObAL 6001 3 repeat 4 1 BONCON gLObAL 8001 2 repeat 4 1 BONCON gLObAL 8001 3 repeat 4 1 BONCON gLObAL 9001 2 repeat 4 1 BONCON gLObAL 9001 3 repeat 4 1 #------------------------------------------------------------------------------- # Constraint equations # # tag::constr[] # type slave dof master dof CONSTR CONEQ GLOBAL 30002 3 0.0 30001 3 1.0 repeat 182 1 0 # CONSTR CONEQ GLOBAL 30001 1 0.0 30005 1 1.0 repeat 60 3 3 CONSTR CONEQ GLOBAL 30001 2 0.0 30005 2 1.0 repeat 60 3 3 CONSTR CONEQ GLOBAL 30001 5 0.0 30005 5 1.0 repeat 60 3 3 CONSTR CONEQ GLOBAL 30001 6 0.0 30005 6 1.0 repeat 60 3 3 # CONSTR CONEQ GLOBAL 30181 1 0.0 30002 1 1.0 CONSTR CONEQ GLOBAL 30181 2 0.0 30002 2 1.0 CONSTR CONEQ GLOBAL 30181 5 0.0 30002 5 1.0 CONSTR CONEQ GLOBAL 30181 6 0.0 30002 6 1.0 # CONSTR CONEQ GLOBAL 30006 1 0.0 30002 1 1.0 repeat 60 3 3 CONSTR CONEQ GLOBAL 30006 2 0.0 30002 2 1.0 repeat 60 3 3 CONSTR CONEQ GLOBAL 30006 5 0.0 30002 5 1.0 repeat 60 3 3 CONSTR CONEQ GLOBAL 30006 6 0.0 30002 6 1.0 repeat 60 3 3 # CONSTR CONEQ GLOBAL 30003 1 0.0 30182 1 1.0 CONSTR CONEQ GLOBAL 30003 2 0.0 30182 2 1.0 CONSTR CONEQ GLOBAL 30003 5 0.0 30182 5 1.0 CONSTR CONEQ GLOBAL 30003 6 0.0 30182 6 1.0 # # end::constr[] #------------------------------------------------------------------------ CONSTR CONEQ GLOBAL 40003 1 0.0 40005 1 1.0 repeat 64 3 3 CONSTR CONEQ GLOBAL 40003 2 0.0 40005 2 1.0 repeat 64 3 3 CONSTR CONEQ GLOBAL 40003 5 0.0 40005 5 1.0 repeat 64 3 3 CONSTR CONEQ GLOBAL 40003 6 0.0 40005 6 1.0 repeat 64 3 3 # CONSTR CONEQ GLOBAL 40004 1 0.0 40002 1 1.0 repeat 64 3 3 CONSTR CONEQ GLOBAL 40004 2 0.0 40002 2 1.0 repeat 64 3 3 CONSTR CONEQ GLOBAL 40004 5 0.0 40002 5 1.0 repeat 64 3 3 CONSTR CONEQ GLOBAL 40004 6 0.0 40002 6 1.0 repeat 64 3 3 # CONSTR CONEQ GLOBAL 40001 1 0.0 40194 1 1.0 CONSTR CONEQ GLOBAL 40001 2 0.0 40194 2 1.0 CONSTR CONEQ GLOBAL 40001 5 0.0 40194 5 1.0 CONSTR CONEQ GLOBAL 40001 6 0.0 40194 6 1.0 # CONSTR CONEQ GLOBAL 40195 1 0.0 40002 1 1.0 CONSTR CONEQ GLOBAL 40195 2 0.0 40002 2 1.0 CONSTR CONEQ GLOBAL 40195 5 0.0 40002 5 1.0 CONSTR CONEQ GLOBAL 40195 6 0.0 40002 6 1.0 #------------------------------------------------------------------------- # CONSTR CONEQ GLOBAL 40002 3 0.0 40001 3 1.0 repeat 194 1 0 # #------------------------------------------------------------------------------ # Material data # # name type poiss talfa tecond heatc beta ea eiy eiz git em gm den Etrans MATERIAL steel_316 linear 0.3 11.7e-6 2.0 50 800 0 0 0 0 2.00E11 0.8e11 8E3 2.00E11 MATERIAL steel_110 linear 0.3 11.7e-6 2.0 50 800 0 0 0 0 1.96E11 0.8e11 8E3 1.96E11 MATERIAL steel_190 linear 0.3 11.7e-6 2.0 50 800 0 0 0 0 2.07E11 0.8e11 8E3 2.07E11 # MATERIAL plast_PVDF linear 0.33 11.7e-6 2.0 50 800 0 0 0 0 869e6 326.0e6 1E3 869e6 MATERIAL plast_PA11 linear 0.35 11.7e-6 2.0 50 800 0 0 0 0 179e6 66.0e6 1E3 179e6 MATERIAL glass_fil linear 0.35 11.7e-6 2.0 50 800 0 0 0 0 179e6 66.0e6 1E3 179e6 MATERIAL rubber linear 0.40 11.7e-6 2.0 50 800 0 0 0 0 155e6 63.0e6 1E3 155e6 #------------------------------------------------------------------------------------- # tag::symspring[] # extra node and element to ensure symmetry of the tensile armour behaviour # additional node Nocoor Coordinates 20001 0.0 0.0 0.0 # node fixed in all dofs BONCON gLObAL 20001 1 BONCON gLObAL 20001 2 BONCON gLObAL 20001 3 BONCON gLObAL 20001 4 BONCON gLObAL 20001 5 BONCON gLObAL 20001 6 # two spring elements, one pr. layer # group elty flexcrossname no n1 n2 n3 n4 Elcon symspring spring137 springmat 20001 20001 30002 20002 20001 40002 Elorient eulerangle 20001 0.0 0.0 0.0 20002 0.0 0.0 0.0 # turnofftransformation ELPROP symspring genspring 0 0 0 0 0 0 1 # MATERIAL springmat genspring belly belly zero zero zero zero # name type alfa eps sig # MATERIAL belly hycurve -1000 -765e9 1000 765e9 # MATERIAL zero hycurve -1000 0 1000 0 # end::symspring[] #----------------------------------------------------------------i---------------------------------------------------- # Loads and sequence #---------------------------------------------------- TIMECO 1.0 0.1 0.10 100.0 static auto none all 100 5 1e-7 TIMECO 2.0 0.01 0.10 100.0 dynamic auto none all 20 5 1e-5 #TIMECO 26.0 0.01 0.10 100.0 dynamic auto none all 20 5 1e-5 #--------------------------------------------- # load input CLOAD 450 1 1 -1000.0 CLOAD 450 1 3 1000.0 # PELOAD 100 100 # to be treated as external pressure for carcsss PILOAD 200 -1001 1e6 -1002 1e6 # then internal pressure PILOAD 200 2001 1e6 2002 1e6 # hist dir no1 r1 no2 r2 n m inistr 350 1 9001 0.106 9002 0.106 # # constr pdisp global 1 5 0.0009268 300 # constr pdisp global 3 5 -0.0009268 300 # tload 400 30001 1.0 30122 1.0 #--------------------------------------------------------------------------------------- # mass THIST 100 0.0 1.0 0.5 1.0 # piload THIST_R 200 0.0 1.0 rampcos 20.0 # pdisp THIST_R 300 0.0 1.0 rampcos 0.0 1.0 13.0 rampcos 0.2 13.0 26.0 rampcos 0.0 # inistr outer sheath THIST_R 350 0.0 1.0 rampcos 0.0 # tload THIST_R 400 0.0 1.0 rampcos 0.0 1.0 3.0 rampcos 0.0 THIST_R 450 0.0 1.0 rampcos 0.0 1.0 2.0 rampcos 0.0 2.0 3.0 rampcos 0.0