Eliasmp
Marine/Ocean
- Jan 8, 2015
- 1
Hello,
I'm modelling a composite cylinder along the x axis with solid186 elements.
My problem is that I can't get the element coordinate system that I want in order to correctly model the layup of the cylinder.
I want the x axis of the elements to be parallel to the axis of the cylinder (which is in advance parallel to the global cartesian x axis).
I've already tried the EORIENT command with all the available directions but I did not get the result I wanted. Actually both Y and Z axes can be parallel to the axial direction but not X axis.
I will be glad to hear that there is a way to achieve this because I'm in the final stage of my Thesis and need to complete some runs with a solid model before finish.
Thanks in advance.
I'm also posting my apdl code:
/prep7
!* --------
!* GEOMETRY
!* --------
!*
/units,si
*SET,L , 0.692 !* Length in m
*SET,Di , 0.250 !* Internal diameter in m
*SET,t , 0.005 !* Thickness in m
*SET,R , Di/2 !* Internal radius
*SET,Do, Di+(2*t) !External diameter
*SET,pi, acos(-1)!*
!*
!*
!*---------------
!*DEFINE GEOMETRY
!*---------------
csys,0
wprota,,,90 !Workplane Rotation so that the Axis of the Cylinder coincides with the global X-Axis
CYL4,0,0,Di/2, ,Do/2, ,L
wprota,,,-90
!* --------
!* ELEMENTS
!* --------
!csys,0
ET,1,SOLID186 !element type
!*
KEYOPT,1,2,0 !Uniform reduced integration (default)
KEYOPT,1,3,1 !Layered Structural Solid
KEYOPT,1,6,0 !Use pure displacement formulation (default)
KEYOPT,1,8,1 !Store top and bottom data for all layers.
!*
CSYS,6
Clocal,11,1
ESYS,11 !Sets the element coordinate system attribute pointer.
MP,EX,1,37.04*1e9 !Defines property data to be associated with the temperature table.
MP,EY,1,15.04*1e9
MP,EZ,1,15.04*1e9
MP,PRXY,1,0.28
MP,PRYZ,1,0.28
MP,PRXZ,1,0.28
MP,GXY,1,5.5*1e9
MP,GYZ,1,2.75*1e9
MP,GXZ,1,5.5*1e9
!*
!* ------
!* LAYERS
!* ------
!*
*SET,NL , 12 !* Number of layers
*SET,tl , t/NL !* Thickness of the layer
*SET,AG , 45 !* Angle
*SET,AM , -45 !* Angle (-)
!* ---------
!* SECTIONS
!* ---------
sectype,1,shell
secoffset,bot
Wangle=AG !Composite Layers
*do,Layer,1,NL,1
secdata,tl,1, Wangle,9
Wangle=-Wangle
*enddo
!*
!* ---------------
!* MESH DEFINITION
!* ---------------
!*
*SET,A , 60 !* Axial line mesh
*SET,C , 16 !* Circumferential Mesh
!*
!Longitudinal Lines
lesize,17,,,A,,,,,1 !Specifies the divisions and spacing ratio on unmeshed lines.
lesize,18,,,A,,,,,1
lesize,19,,,A,,,,,1
lesize,20,,,A,,,,,1
!Circumferential Lines
lesize,1,,,C,,,,,1 !Specifies the divisions and spacing ratio on unmeshed lines.
lesize,2,,,C,,,,,1
lesize,3,,,C,,,,,1
lesize,4,,,C,,,,,1
lesize,5,,,C,,,,,1
lesize,6,,,C,,,,,1
lesize,7,,,C,,,,,1
lesize,8,,,C,,,,,1
lesize,9,,,C,,,,,1
lesize,10,,,C,,,,,1
lesize,11,,,C,,,,,1
lesize,12,,,C,,,,,1
lesize,13,,,C,,,,,1
lesize,14,,,C,,,,,1
lesize,15,,,C,,,,,1
lesize,16,,,C,,,,,1
vsweep,all !volume sweep
/PSYMB,ESYS,1
/REPLOT
EORIENT,1,NEGX,
/REPLOT
_________________
Elias Bilalis
Student at National Technical University of Athens
I'm modelling a composite cylinder along the x axis with solid186 elements.
My problem is that I can't get the element coordinate system that I want in order to correctly model the layup of the cylinder.
I want the x axis of the elements to be parallel to the axis of the cylinder (which is in advance parallel to the global cartesian x axis).
I've already tried the EORIENT command with all the available directions but I did not get the result I wanted. Actually both Y and Z axes can be parallel to the axial direction but not X axis.
I will be glad to hear that there is a way to achieve this because I'm in the final stage of my Thesis and need to complete some runs with a solid model before finish.
Thanks in advance.
I'm also posting my apdl code:
/prep7
!* --------
!* GEOMETRY
!* --------
!*
/units,si
*SET,L , 0.692 !* Length in m
*SET,Di , 0.250 !* Internal diameter in m
*SET,t , 0.005 !* Thickness in m
*SET,R , Di/2 !* Internal radius
*SET,Do, Di+(2*t) !External diameter
*SET,pi, acos(-1)!*
!*
!*
!*---------------
!*DEFINE GEOMETRY
!*---------------
csys,0
wprota,,,90 !Workplane Rotation so that the Axis of the Cylinder coincides with the global X-Axis
CYL4,0,0,Di/2, ,Do/2, ,L
wprota,,,-90
!* --------
!* ELEMENTS
!* --------
!csys,0
ET,1,SOLID186 !element type
!*
KEYOPT,1,2,0 !Uniform reduced integration (default)
KEYOPT,1,3,1 !Layered Structural Solid
KEYOPT,1,6,0 !Use pure displacement formulation (default)
KEYOPT,1,8,1 !Store top and bottom data for all layers.
!*
CSYS,6
Clocal,11,1
ESYS,11 !Sets the element coordinate system attribute pointer.
MP,EX,1,37.04*1e9 !Defines property data to be associated with the temperature table.
MP,EY,1,15.04*1e9
MP,EZ,1,15.04*1e9
MP,PRXY,1,0.28
MP,PRYZ,1,0.28
MP,PRXZ,1,0.28
MP,GXY,1,5.5*1e9
MP,GYZ,1,2.75*1e9
MP,GXZ,1,5.5*1e9
!*
!* ------
!* LAYERS
!* ------
!*
*SET,NL , 12 !* Number of layers
*SET,tl , t/NL !* Thickness of the layer
*SET,AG , 45 !* Angle
*SET,AM , -45 !* Angle (-)
!* ---------
!* SECTIONS
!* ---------
sectype,1,shell
secoffset,bot
Wangle=AG !Composite Layers
*do,Layer,1,NL,1
secdata,tl,1, Wangle,9
Wangle=-Wangle
*enddo
!*
!* ---------------
!* MESH DEFINITION
!* ---------------
!*
*SET,A , 60 !* Axial line mesh
*SET,C , 16 !* Circumferential Mesh
!*
!Longitudinal Lines
lesize,17,,,A,,,,,1 !Specifies the divisions and spacing ratio on unmeshed lines.
lesize,18,,,A,,,,,1
lesize,19,,,A,,,,,1
lesize,20,,,A,,,,,1
!Circumferential Lines
lesize,1,,,C,,,,,1 !Specifies the divisions and spacing ratio on unmeshed lines.
lesize,2,,,C,,,,,1
lesize,3,,,C,,,,,1
lesize,4,,,C,,,,,1
lesize,5,,,C,,,,,1
lesize,6,,,C,,,,,1
lesize,7,,,C,,,,,1
lesize,8,,,C,,,,,1
lesize,9,,,C,,,,,1
lesize,10,,,C,,,,,1
lesize,11,,,C,,,,,1
lesize,12,,,C,,,,,1
lesize,13,,,C,,,,,1
lesize,14,,,C,,,,,1
lesize,15,,,C,,,,,1
lesize,16,,,C,,,,,1
vsweep,all !volume sweep
/PSYMB,ESYS,1
/REPLOT
EORIENT,1,NEGX,
/REPLOT
_________________
Elias Bilalis
Student at National Technical University of Athens