finish ! Finishes the last action/setting /CLEAR,nostart ! Close old analysis and start new one, clear everything !-------------------------------------------------------------------------------- ! Initiate pre-processing !-------------------------------------------------------------------------------- /prep7 PI = ACOS(-1) t = 2 ! thickness b = 5 ! width l = 10 ! length ET,1,SOLSH190 ! Element type, Reference number, Element name KEYOPT,1,2,1 ! 1 for Include enhanced transverse-shear strains KEYOPT,1,6,0 ! 0 for Use pure displacement formulation (default) KEYOPT,1,8,0 ! For multilayer elements, store data for top and bottom for all layers (1) ET,2,SOLID185 ! Element type, Reference number, Element name KEYOPT,2,2,2 ! 2 for Enhanced strain formulation against shear and volume locking KEYOPT,2,3,1 ! 0 for homogenious, 1 for layered solid KEYOPT,2,6,0 ! 0 for Use pure displacement formulation (default) KEYOPT,2,8,0 ! For multilayer elements, store data for top and bottom for all layers (1) ! Material Properties (MAT = 1) MPTEMP,1,0 MPDATA,EX,1,,135000 ! Tensile E-Modulus in x (1)-direction [MPa] MPDATA,EY,1,,8500 ! Tensile E-Modulus in y (2)-direction [MPa] MPDATA,EZ,1,,8500 ! Tensile E-Modulus in z (3)-direction [MPa] MPDATA,PRXY,1,,0.35 ! Poisson's ratio xy outer cylinder [] MPDATA,PRYZ,1,,0.35 ! Poisson's ratio yz outer cylinder [] MPDATA,PRXZ,1,,0.35 ! Poisson's ratio xz outer cylinder [] MPDATA,GXY,1,,4200 ! Shear Modulus in xy-direction [MPa] MPDATA,GYZ,1,,4200 ! Shear Modulus in yz-direction [MPa] MPDATA,GXZ,1,,4200 ! Shear Modulus in xz-direction [MPa] MPDATA,DENS,1,,1.59e-9 ! Density UD-ply [t/mm^3] UIMP,1,REFT,,,0 ! Reference temperature = 0 degree !MPDATA,ALPX,1,,300e-12 !MPDATA,ALPY,1,,-100e-12 !MPDATA,ALPZ,1,,-100e-12 MPDATA,ALPX,1,,0 MPDATA,ALPY,1,,0 MPDATA,ALPZ,1,,0 ! Material properties (MAT = 2) !----------------------------------------------------------------------- MPTEMP,1,0 MPDATA,EX,2,,30000 ! E-Modulus in x-direction [MPa] MPDATA,EY,2,,15500 ! E-Modulus in y-direction [MPa] MPDATA,EZ,2,,15500 ! E-Modulus in z-direction [MPa] MPDATA,PRXY,2,,0.4 ! Poisson's ratio xy outer cylinder [] MPDATA,PRYZ,2,,0.35 ! Poisson's ratio yz outer cylinder [] MPDATA,PRXZ,2,,0.4 ! Poisson's ratio xz outer cylinder [] MPDATA,GXY,2,,10.7 ! Shear Modulus in xy-direction [MPa] MPDATA,GYZ,2,,5.7 ! Shear Modulus in yz-direction [MPa] MPDATA,GXZ,2,,10.7 ! Shear Modulus in xz-direction [MPa] MPDATA,DENS,2,,4.7e-9 ! Density [t/mm^3] UIMP,2,REFT,,,0 ! Reference temperature = 0 degree MPDATA,ALPX,2,,418e-09 ! CTE analogy to d_33 in 1/K longitudinal MPDATA,ALPY,2,,-198e-12 ! CTE analogy to d_31 in 1/K transversal MPDATA,ALPZ,2,,-198e-12 ! CTE analogy to d_31 in 1/K transversal ! Local cylindric coordinate system !---------------------------------------------------------------------------- LOCAL,11,1,0,0,0 ! Defining a new local cylindrical coordinate system, origin 0,0,0 !-------------------------------------------------------------------------------- ! Definition of geometry / cylindrical structure !-------------------------------------------------------------------------------- MSHKEY,1 ! Specifies whether free (0) or mapped meshing (1) is used CYLIND, b,b+t, 0,l,-45,45 ! Cylindric geometry defined by parameters csys,11 ! Activating coordinate system 11 esys,11 sect,1,shell,, ! Definition of layup with "section" command secdata,t/2,1,0.0,3 secdata,t/2,1,90,3 secoffset,MID seccontrol,0,0,0,0,1,1,1 VEORIENT,1,KP,8,5 ! Specifies orientation of z-axis, x and y are determined by right-hand-rule LESIZE,8,,,3 ! Specifies divisions and spacing ratio on unmeshed lines LESIZE,7,1 ! Specifies divisions and spacing ratio on unmeshed lines LESIZE,12,1 ! Specifies divisions and spacing ratio on unmeshed lines VATT,1,,1 ! Associated element attributes with selected, unmeshed volumes VMESH,1 ! Generates nodes and volume elements nsel,all ! Select all nodes cm,inner_cylinder_nodes,node ! Groups outer cylinder nodes to a component nsel,none ! Unselect all nodes MSHKEY,1 ! Specifies whether free (0) or mapped meshing (1) is used CYLIND, b+t,b+t+t/2, 0,l,-45,45 ! Cylindric geometry defined by parameters csys,11 ! Activating coordinate system 11 esys,11 sect,2,shell,, ! Definition of layup with "section" command secdata,t/2,2,0.0,3 secoffset,MID seccontrol,0,0,0,0,1,1,1 VEORIENT,2,KP,16,13 ! Specifies orientation of z-axis, x and y are determined by right-hand-rule LESIZE,20,,,1 ! Specifies divisions and spacing ratio on unmeshed lines LESIZE,19,0.5 ! Specifies divisions and spacing ratio on unmeshed lines LESIZE,24,0.5 ! Specifies divisions and spacing ratio on unmeshed lines VATT,2,,2 ! Associated element attributes with selected, unmeshed volumes VMESH,2 ! Generates nodes and volume elements nsel,all ! Select all nodes cmsel,u,inner_cylinder_nodes cm,outer_cylinder_nodes,node ! Groups outer cylinder nodes to a component nsel,none ! Unselect all nodes ! Boundaries and node components csys,11 ! Cylindric coordinate system cmsel,s,inner_cylinder_nodes ! Select only all nodes of inner cylinder nsel,r,loc,x,b,b+t ! Select all nodes along x-axis within given boundaries nsel,r,loc,z,0,0.01 cm,inner_cylinder_fixation_nodes,node ! Groups inner cylinder fixation nodes to a component nsel,none ! Unselect all nodes csys,11 ! Cylindric coordinate system cmsel,s,inner_cylinder_nodes ! Select only all nodes of inner cylinder nsel,r,loc,x,b,b+t ! Select nodes through thickness nsel,r,loc,y,-45,-45+0.01 ! Reselect nodes depending on angle cm,inner_cylinder_sym_nodes_1,node ! Groups symmetry nodes to a component nsel,none ! Unselect all nodes csys,11 ! Cylindric coordinate system cmsel,s,inner_cylinder_nodes ! Select only all nodes of inner cylinder nsel,r,loc,x,b,b+t ! Select nodes through thickness nsel,r,loc,y,45-0.01,45 ! Reselect nodes depending on angle cm,inner_cylinder_sym_nodes_2,node ! Groups symmetry nodes to a component nsel,none ! Unselect all nodes csys,11 ! Cylindric coordinate system cmsel,s,outer_cylinder_nodes ! Select only all nodes of inner cylinder nsel,r,loc,x,b+t,b+t+t/2 ! Select all nodes along x-axis within given boundaries nsel,r,loc,z,0,0.01 cm,outer_cylinder_fixation_nodes,node ! Groups inner cylinder fixation nodes to a component nsel,none ! Unselect all nodes csys,11 ! Cylindric coordinate system cmsel,s,outer_cylinder_nodes ! Select only all nodes of inner cylinder nsel,r,loc,x,b+t,b+t+t/2 ! Select nodes through thickness nsel,r,loc,y,-45,-45+0.01 ! Reselect nodes depending on angle cm,outer_cylinder_sym_nodes_1,node ! Groups symmetry nodes to a component nsel,none ! Unselect all nodes csys,11 ! Cylindric coordinate system cmsel,s,outer_cylinder_nodes ! Select only all nodes of inner cylinder nsel,r,loc,x,b+t,b+t+t/2 ! Select nodes through thickness nsel,r,loc,y,45-0.01,45 ! Reselect nodes depending on angle cm,outer_cylinder_sym_nodes_2,node ! Groups symmetry nodes to a component nsel,none ! Unselect all nodes ! Rotating nodes into cylindrical coordinate system nsel,all ! Select all nodes csys,11 ! Cylindric coordinate system nrotat,all ! Rotate all nodes ! Load Step FINISH ! Exits normally from a processor ! Static analysis !---------------------------------------------------------------------------- /SOLU ! Enters the solution processor ANTYPE,0 ! A new analysis type is specified ! Boundary conditions !---------------------------------------------------------------------------- csys,11 ! Select cylindric coordinate system allsel ! BCs !---------------------------------------------------------------------------- D,inner_cylinder_sym_nodes_1,uy,0 D,inner_cylinder_sym_nodes_2,uy,0 D,inner_cylinder_fixation_nodes,all,0 D,outer_cylinder_sym_nodes_1,uy,0 D,outer_cylinder_sym_nodes_2,uy,0 D,outer_cylinder_sym_nodes_1,uz,0 D,outer_cylinder_sym_nodes_2,uz,0 ! Temperature Load !---------------------------------------------------------------------------- tref,0 bfunif,temp,1500/0.0005 ! Solution, static analysis !---------------------------------------------------------------------------- NLGEOM,OFF ! Large-deflections allowed or not OUTRES,LOCI,all ! Integration point locations written to database OUTRES,all,all ! All solution item except LOCI and SVAR written to database nsub,3,50,3 ! Number of increments within substep (first value), max (second value) and min (third value) *GET,wall_start,ACTIVE,0,TIME,WALL ! Start time of solution /STATUS,SOLU SOLVE ! Solve the system FINISH *GET,wall_end,ACTIVE,0,TIME,WALL solving_time = wall_end - wall_start ! Get solving time ! Postprocessing !---------------------------------------------------------------------------- /Post1 PLNSOL,u,sum,1,1.0