HEAD 4"WS-FPS2000 p=20MPa RUC Thinw. HS353-363 HC453-463 Oldfricmodel #--------------------------------------------------------------------------- # Control data # # maxit ndim isolvr npoint ipri conr gacc iproc Control 500 3 2 16 01 1.e-9 0.0 stressfree #restart 130 # # # Matat Alpha1 Alpha2 Alpha Dyncont 1 0.0 0.09 -0.05 # # 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 (radial dof fixed for stability) 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 # core (outside flexbody until tensile1= 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 # flexwear2 (radial dof fixed for stability because of direct contact vis hcont453) 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 # 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 #------------------------------------------------------------------------------------- # # 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 hshear363 mypipe 1001 1 2 1001 repeat 2 1 1 Elcon flexbarrier hshear363 mypipe 2001 1 2 2001 repeat 2 1 1 Elcon flexloc hshear363 mypipe 3001 1 2 2001 repeat 2 1 1 Elcon flexwear1 hshear363 mypipe 4001 1 2 2001 repeat 2 1 1 Elcon flextensile1 hshear353 mypipe 30001 1 2 30001 30002 30061 1 2 30181 30182 repeat 2 61 1 Elcon flexwear2 hshear363 mypipe 6001 1 2 6001 repeat 2 1 1 Elcon flextensile2 hshear353 mypipe 40001 1 2 40001 40002 40065 1 2 40193 40194 repeat 2 65 1 Elcon flextape1 hshear363 mypipe 8001 1 2 8001 repeat 2 1 1 Elcon flexshield hshear363 mypipe 9001 1 2 8001 repeat 2 1 1 # end::elcons[] # Contact layers # # tag::elconc[] Elcon flexwear1contact hcont463 mypipe 53001 2001 30001 30002 53061 2001 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 # end::elconc[] # #--------------------------------------------------------------------------- # 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 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 #-------------------------------------------------------------------------------- # # 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 2-2-flexbarrier 3-3-flexloc 4-4-flexwear1-flexwear1contact 5-5-flexwear2inwardcontact-flextensile1 6-6-flexwear2 7-7-flextape1contact-flextensile2 8-8-flextape1 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 #------------------------------------------------------------------------------------------------------------ # By activating this option via elprop you will obtain a constant stick stiffness that do not vary with contact pressure # If not activation the shearmodel option the stick stiffness is based (as before ) on a given stick displacement that is set at the time = timeini # Then the stiffness will vary with contact pressure # tag:shearmodel[] #elprop flexwear1contact shearmodel d d d #elprop flexwear2inwardcontact shearmodel d d d #elprop flextape1contact shearmodel d d d # end::shearmodel[] # #----------------------------------------------------------------------- # 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 4 repeat 183 1 BONCON gLObAL 40001 4 repeat 195 1 #--------------------------------------------------------------- BONCON gLObAL 1001 1 repeat 2 1 BONCON gLObAL 1001 2 repeat 2 1 BONCON gLObAL 1001 3 repeat 2 1 BONCON gLObAL 2001 2 repeat 2 1 BONCON gLObAL 2001 3 repeat 2 1 BONCON gLObAL 6001 1 repeat 2 1 BONCON gLObAL 6001 2 repeat 2 1 BONCON gLObAL 6001 3 repeat 2 1 BONCON gLObAL 8001 2 repeat 2 1 BONCON gLObAL 8001 3 repeat 2 1 #------------------------------------------------------------------------------- # tag::constr[] # Constraint equations # 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 alle 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. tensile 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[] # #-------------------------------------------------------------------------------------------------------------------- # Loads and sequence #---------------------------------------------------- TIMECO 1.0 0.1 0.10 100.0 static auto none all 500 5 1e-5 TIMECO 1.1 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 1.0 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