{"id":370956,"date":"2024-07-08T14:36:22","date_gmt":"2024-07-08T14:36:22","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/problems-with-orientation-angle-when-using-smeared-reinforcements\/"},"modified":"2024-07-08T14:36:22","modified_gmt":"2024-07-08T14:36:22","slug":"problems-with-orientation-angle-when-using-smeared-reinforcements","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/problems-with-orientation-angle-when-using-smeared-reinforcements\/","title":{"rendered":"Problems with orientation angle, when using smeared Reinforcements"},"content":{"rendered":"

Hi all,<\/p>\n

I use the following APDL-Code to model a quarter torus with fiber reinforcements within t3D continuum elements (SOLID186). The orientation of the reinforcements (Smeared modelling approach with REINF 265 Elements) should vary over the cross-section from 5 to 90 degrees. Somehow ANSYS gives me problems using the EREINF command, so that all elements adopt the orientation of the last defined reinforcing elements. In this case 90 degrees. Can anyone help me? Does this approach work at all?<\/p>\n

The background of the question is, that in future a toroidal helix should be modelled, to simulate fibre-reinforced pressurized hoses in which the winding angle changes over the arc length.<\/p>\n

!====================================================================== <\/div>\n
!==================        Start of Input Deck        ================= <\/div>\n
!====================================================================== <\/div>\n
fini<\/div>\n
\/clear,all  <\/div>\n
\/prep7<\/div>\n
!====================================================================== <\/div>\n
!===                          Parameter                             ===<\/div>\n
!====================================================================== <\/div>\n
*SET,ro,24\/2  <\/div>\n
*SET,Do,ro*2<\/div>\n
*SET,ri,20\/2 <\/div>\n
*SET,Di,ri*2<\/div>\n
*SET,t,ro-ri<\/div>\n
*SET,rm,(ro+ri)\/2<\/div>\n
*SET,rk,50<\/div>\n
 <\/div>\n
*SET,Emodul,750<\/div>\n
*SET,nu,0.4 <\/div>\n
*SET,pressure,1<\/div>\n
*SET,Emodul_faser,Emodul*100 <\/div>\n
*SET,nu_faser,nu <\/div>\n
*SET,wickelwinkel,60 <\/div>\n
*SET,wickelwinkel_rad,wickelwinkel*3.14159265359\/180 <\/div>\n
 <\/div>\n
 <\/div>\n
!====================================================================== <\/div>\n
!===                         Materials            === <\/div>\n
!====================================================================== <\/div>\n
mp,  ex,1,Emodul  ! Base Material<\/div>\n
mp,  ex,2,Emodul_faser ! Reinforcing Material<\/div>\n
mp,prxy,1,nu<\/div>\n
mp,prxy,2,nu_faser<\/div>\n
!====================================================================== <\/div>\n
!===                                    Geometry                    === <\/div>\n
!====================================================================== <\/div>\n
 <\/div>\n
k,1,rk,0,0<\/div>\n
k,2,0,rk,0 <\/div>\n
k,3,0,0,0<\/div>\n
LARC,1,2,3,rk<\/div>\n
 <\/div>\n
wprota,0,90,0<\/div>\n
cyl4,rk,0,ri,0,ro,180<\/div>\n
wprota,0,-90,0<\/div>\n
 <\/div>\n
vdrag,1,,,,,,1<\/div>\n
!====================================================================== <\/div>\n
!===                      Elements               === <\/div>\n
!====================================================================== <\/div>\n
et,1,SOLID186<\/div>\n
mat,1<\/div>\n
!====================================================================== <\/div>\n
!===                 Meshing                                        === <\/div>\n
!====================================================================== <\/div>\n
 <\/div>\n
esize,2.2 <\/div>\n
vmesh, 1<\/div>\n
 <\/div>\n
!====================================================================== <\/div>\n
!===                  SYMMETRY-BC                === <\/div>\n
!======================================================================<\/div>\n
local,11,3,0,0,0,0,0,0,rk ! Erzeugt ein lokales Koordinatensystem – Toruskoordinaten mit einem Radius rk<\/div>\n
 <\/div>\n
 <\/div>\n
! SYMMETRIE_RB POLOIDAL<\/div>\n
asel,s,,,3<\/div>\n
asel,a,,,5<\/div>\n
nsla,s,1<\/div>\n
nrotat,all<\/div>\n
d,all,uz,0<\/div>\n
ALLSEL,ALL<\/div>\n
 <\/div>\n
!====================================================================== <\/div>\n
!===                  BC                                    === <\/div>\n
!====================================================================== <\/div>\n
Csys,11 ! Cylindrical with the global Cartesian z as the axis of rotation<\/div>\n
asel,s,,,1<\/div>\n
nsla,s,1<\/div>\n
nrotat,all<\/div>\n
 <\/div>\n
d,all,ux,0<\/div>\n
d,all,uy,0<\/div>\n
d,all,uz,0<\/div>\n
 <\/div>\n
ALLSEL,ALL<\/div>\n
 <\/div>\n
!====================================================================== <\/div>\n
!===                        REINFORCEMENTS                          ===<\/div>\n
!====================================================================== <\/div>\n
 <\/div>\n
*DO,i,1,18,1<\/div>\n
*SET,j,5*i<\/div>\n
*SET,wickelwinkel,j<\/div>\n
sect,1,REINF,SMEAR<\/div>\n
secd,2,0.5,1.0,,wickelwinkel,ELEF,1,0.5<\/div>\n
esel,all<\/div>\n
*SET,startElem,i<\/div>\n
*SET,endElem,793+i<\/div>\n
esel,s,,,startElem,endElem,18<\/div>\n
EREINF<\/div>\n
*ENDDO<\/div>\n

 <\/p>\n

Thank you very much for your help!<\/p>\n","protected":false},"template":"","class_list":["post-370956","topic","type-topic","status-closed","hentry","topic-tag-composites","topic-tag-orientation-of-fibers-2","topic-tag-reinforcement"],"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