{"id":376048,"date":"2024-08-02T19:11:05","date_gmt":"2024-08-02T19:11:05","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/reply\/376048\/"},"modified":"2024-08-02T19:11:05","modified_gmt":"2024-08-02T19:11:05","slug":"376048","status":"publish","type":"reply","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/reply\/376048\/","title":{"rendered":"Reply To: Reinforcing section (steel wires and rubber) details for multilayer spiral pipe"},"content":{"rendered":"

<p>Volume 1 appears to be subdivided (sliced) in such a way that there is a fairly sharp angle between the plane of the cut and the inner cylindrical surface. For what reason is the the geometry subdivided in this way? You mentioned contact. Is there contact between one of the surfaces on the right end of the reinforced pipe, with the pipe subjected to an axial load? Did you end up meshing with tetrahedra? If the subdivision was made to clearly define contact surface, I think you can do away with it… the pressure distribution from the contacting body (if you model it explicity) should be fairly well represented.<\/p><p>There are multiple options for defining smeared reinforcing. I believe that in the option you are using, the location and orientation of the fibers is defined with respect to the order of nodes in each element:<\/p><p><\/p><p>If you have a tetrahedral mesh, with random node number order within each element and multiple layers of elements through the thickness of the pipe wall, you might end up with a corresponspondingly random reinforcing locations and orientations.<\/p><p>I can offer you the following (copy the lines below into a text file and read the text file into MAPDL with the \/INPUT command). Try two different values of parameter r_reinf (0.001 & 0.0005). The calculated results should differ substantially. I have the luxury of using very simple geometry (not sliced as yours is). Even with this “unfair advantage”, I had to resort to a “trick” to align my SOLID185 elements (initially define them to be layered so I could use EORIENT to renumber the elements in such a way as to get consistent REINF265 orientation in the cylinder wall). I added comments to explain the numerous other settings I used in the setup.<\/p><p>Another option would be to create MESH200 elements and assign the reinf section properties to them.<\/p><p> <\/p><p>I hope this helps!<\/p><p>Bill<\/p><p> <\/p><p> <\/p><div>fini<\/div><div>\/cle<\/div><div> <\/div><div>\/sys,del file*.png<\/div><div> <\/div><div>\/vie,1,1,1,1<\/div><div>\/vup,1,z<\/div><div> <\/div><div> <\/div><div>C********************************************************<\/div><div>C*** PARAMETERS<\/div><div>C********************************************************<\/div><div>pi=acos(-1)<\/div><div> <\/div><div>r1=0.010 ! CYLINDER INNER RADIUS<\/div><div>r2=0.020 ! CYLINDER OUTER RADIUS<\/div><div>l=0.01 ! CLINDER LENGTH<\/div><div> <\/div><div>r_reinf=0.001 ! REINFORCING FIBER RADIUS (COMPARE 0.001 & 0.0005)<\/div><div>a=pi*r_reinf**2 ! CORRESPONDING REINF CROSS SECTION AREA<\/div><div>facenum=3 ! FACE NUMBER FOR “ELEF” REINF SECTION PARAMETER<\/div><div>s=0.0025 ! REINF SPACING<\/div><div>thta=30 ! MAGNITUDE OF REINF FIBER ORIENTATION ANGLE<\/div><div> <\/div><div>E_base=0.05e9 ! BASE ELEMENT ELASTIC MODULUS (~RUBBER)<\/div><div>nu_base=0.20 ! BASE ELEMENET POISSON’S<\/div><div> <\/div><div>E_reinf=2e11 ! REINF ELASTIC MODULUS<\/div><div>nu_reinf=0.20 ! REINF POISSON’S<\/div><div> <\/div><div>u_top=0.001 ! MAGNITUDE OF AXIAL DISPLACEMENT IMPOSED ON CYLINDER IN TRIAL SOLVE<\/div><div>f_top=2e3 ! MAGNITUDE OF AXIAL FORCE IMPOSED ON CYLINDER IN TRIAL SOLVE<\/div><div> <\/div><div> <\/div><div>\/title,R_fiber=%r_reinf%, S_fiber=%s%, E_fiber=%E_reinf%<\/div><div> <\/div><div>C********************************************************<\/div><div>C*** GEOMETRY & BASE ELEMENT MESH<\/div><div>C********************************************************<\/div><div>wpcs,-1,0 ! CYLINDRICAL ELEMENT COORDINATE SYSTEM<\/div><div>cswp,11,1<\/div><div> <\/div><div>\/prep7 ! GEOMETRY (FOUR 90 DEG SECTORS)<\/div><div>cyli,r1,r2,0,l,0,90<\/div><div>cyli,r1,r2,0,l,90,180<\/div><div>cyli,r1,r2,0,l,180,270<\/div><div>cyli,r1,r2,0,l,270,360<\/div><div>numm,kp ! MERGE COINCIDENT KEYPOINTS<\/div><div>vatt,1,1,1,11 ! ASSIGN ATTRIBUTE NUMBERS (MAT, REAL, TYPE, ESYS)<\/div><div> <\/div><div>et,1,185,,,1 ! LAYERED SOLID185 ELEMENT TYPE<\/div><div>mp,ex,1,E_base ! ELASTIC MODULUS<\/div><div>mp,nuxy,1,nu_base ! POISSON’S<\/div><div>vmes,all ! MESH VOLUMES<\/div><div> <\/div><div> <\/div><div>C********************************************************<\/div><div>C*** REORIENT BASE ELEMENTS<\/div><div>C********************************************************<\/div><div>eorient,lysl,posz ! RENUMBER SO THAT FACE 1 IS PARALLEL +Z ESYS AXIS<\/div><div> <\/div><div> <\/div><div>C********************************************************<\/div><div>C*** REINFORCING ELEMENT PROPERTIES<\/div><div>C********************************************************<\/div><div>mp,ex,2,E_reinf ! ELASTIC MODULUS<\/div><div>mp,nuxy,2,nu_reinf ! POISSON’S<\/div><div> <\/div><div>sect,2,reinf,smear ! REINF SECTION FOR INNER RADIUS<\/div><div>secd,2,A,s,,thta,elef,facenum,0.5<\/div><div>seccontrol,0,1,3 ! REMOVE BASE MATERIAL IN DOMAIN OCCUPIED BY REINFORCING, INCLUDE TRANSVERSE SHEAR & BENDING<\/div><div> <\/div><div>sect,3,reinf,smear ! REINF SECTION FOR OUTER RADIUS<\/div><div>secd,2,A,s,,-thta,elef,facenum,0.5<\/div><div>seccontrol,0,1,3 ! REMOVE BASE MATERIAL IN DOMAIN OCCUPIED BY REINFORCING, INCLUDE TRANSVERSE SHEAR & BENDING<\/div><div> <\/div><div> <\/div><div>C********************************************************<\/div><div>C*** CREATE REINFORCING ELEMENTS<\/div><div>C********************************************************<\/div><div>csys,11 ! DETERMINE INNER\/OUTER RADII<\/div><div>*get,rmin,node,,mnloc,x<\/div><div>*get,rmax,node,,mxloc,x<\/div><div> <\/div><div>nsel,s,loc,x,rmin ! REINFORCE BASE ELEMENTS ON INNER RADIUS<\/div><div>esln<\/div><div>esel,r,type,,1<\/div><div>nsle<\/div><div>secn,2<\/div><div>ereinf<\/div><div> <\/div><div>nsel,s,loc,x,rmax ! REINFORCE BASE ELEMENTS ON OUTER RADIUS<\/div><div>esln<\/div><div>esel,r,type,,1<\/div><div>nsle<\/div><div>secn,3<\/div><div>ereinf<\/div><div> <\/div><div> <\/div><div>C********************************************************<\/div><div>C*** MAKE BASE ELEMENT HOMOGENEOUS, PLOT ELEMENTS<\/div><div>C********************************************************<\/div><div>et,1,185 ! CHANGE BASE ELEMENT TYPE FROM LAYERED TO HOMOGENEOUS OPTION<\/div><div> <\/div><div>esel,s,type,,1 ! MAKE BASE ELEMENTS TRANSLUCENT<\/div><div>\/trlcy,elem,0.9<\/div><div>alls<\/div><div>\/psy,layr,-1 ! DISPLAY LAYERS (REINFORCING)<\/div><div>\/esh,1<\/div><div>\/dev,vect,0 ! RASTER FILL DISPLAY<\/div><div>eplo ! PLOT ELEMENTS<\/div><div>\/sho,png $eplo $\/sho,close $\/wait,2<\/div><div> <\/div><div>esel,s,type,,2 ! SELECT REINFORCING ELEMENTS<\/div><div>\/dev,vect,1 ! WIRE FRAME DISPLAY<\/div><div>eplo ! PLOT ELEMENTS<\/div><div>\/sho,png,,1 $eplo $\/sho,close $\/wait,2<\/div><div> <\/div><div>\/dev,vect,0 ! REVERT TO RASTER FILL<\/div><div>esel,s,type,,1 ! MAKE BASE ELEMENTS OPAQUE<\/div><div>\/trlcy,elem,0<\/div><div>alls<\/div><div> <\/div><div>fini<\/div><div> <\/div><div> <\/div><div>C********************************************************<\/div><div>C*** TRIAL SOLVE<\/div><div>C********************************************************<\/div><div>\/solu<\/div><div> <\/div><div>nlgeom,on ! MAY BE NECESSARY TO INVOKE TENSION-ONLY REINF BEHAVIOR<\/div><div> <\/div><div>autots,off<\/div><div>outr,all,all<\/div><div> <\/div><div>nsel,s,loc,z ! FIX BASE<\/div><div>d,all,all<\/div><div> <\/div><div>nsel,s,loc,z,l<\/div><div>cp,1,uz,all<\/div><div>nd_top=ndnext(0)<\/div><div> <\/div><div>alls<\/div><div> <\/div><div>!d,nd_top,uz,u_top ! TENSION ON TOP SURFACE<\/div><div>f,nd_top,fz,f_top ! TENSION ON TOP SURFACE<\/div><div>nsub,5<\/div><div>solv<\/div><div> <\/div><div>!d,nd_top,uz,-u_top ! COMPRESSION ON TOP SURFACE<\/div><div>f,nd_top,fz,-f_top ! COMPRESSION ON TOP SURFACE<\/div><div>nsub,10<\/div><div>solv<\/div><div> <\/div><div>fini<\/div><div> <\/div><div> <\/div><div>C********************************************************<\/div><div>C*** POST PROCESSING<\/div><div>C********************************************************<\/div><div> <\/div><div>\/post26<\/div><div>!rfor,2,nd_top,f,z<\/div><div>nsol,2,nd_top,u,z<\/div><div>\/axl,y,Axial Displacement<\/div><div>plva,2<\/div><div>\/sho,png $plva,2 $\/sho,close $\/wait,2<\/div><div> <\/div><div> <\/div><div>\/post1<\/div><div>set,1<\/div><div> <\/div><div>esel,s,ename,,265<\/div><div>ples,s,x<\/div><div>\/sho,png $ples,s,x $\/sho,close $\/wait,2<\/div><div> <\/div><div>esel,s,ename,,185<\/div><div>rsys,1<\/div><div>plns,s,z<\/div><div>\/sho,png $plns,s,z $\/sho,close $\/wait,2<\/div><div> <\/div><\/p>\n","protected":false},"template":"","class_list":["post-376048","reply","type-reply","status-publish","hentry"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t