*Heading
 _Scenario_Mechanical_
** ---------------------------------------------------------------
** Use this input deck to compute the in-situ fail stress of
** the constituents in the multiscale model.
** The fail stresses of the composite are applied in 6 load cases,
** the in-situ fail stresses are corresponding constituent stresses (S_fiber and S_matrix) in the field output.
** ---------------------------------------------------------------
*Preprint, echo=NO, model=no, history=NO, contact=NO
** ----------------------------------------------------------------
** 
** PART INSTANCE: Part-1-1
** 
*Node
      1,           1.,           1.,           1.
      2,           1.,           0.,           1.
      3,           1.,           1.,           0.
      4,           1.,           0.,           0.
      5,           0.,           1.,           1.
      6,           0.,           0.,           1.
      7,           0.,           1.,           0.
      8,           0.,           0.,           0.
*Element, type=C3D8R
1, 5, 6, 8, 7, 1, 2, 4, 3
*Nset, nset=Part-1-1_Set-1, generate
 1,  8,  1
*Elset, elset=Part-1-1_Set-1
 1,
** Section: Section-1
*Solid Section, elset=Part-1-1_Set-1, material=composite,orientation=ori
,
*orientation,name=ori
1,0,0,0,1,0
3,0
*System
*Node
      9,           0.,           0.,           0.
     10,           0.,           0.,           0.
*Nset, nset=RP-Normal
 9,
*Nset, nset=RP-Shear
10,
*Elset, elset=RVE_AG_ALL_ORIG_ELEMENTS
 1,
*Nset, nset=RVE_AG_NegX, generate
 5,  8,  1
*Nset, nset=RVE_AG_NegX_RVE_AGSORT, unsorted
 5, 6, 7, 8
*Nset, nset=RVE_AG_NegY, generate
 2,  8,  2
*Nset, nset=RVE_AG_NegY_RVE_AGSORT, unsorted
 6, 8
*Nset, nset=RVE_AG_NegZ
 3, 4, 7, 8
*Nset, nset=RVE_AG_NegZ_RVE_AGSORT, unsorted
 8,
*Nset, nset=RVE_AG_PinNode
 8,
*Nset, nset=RVE_AG_PosX, generate
 1,  4,  1
*Nset, nset=RVE_AG_PosX_RVE_AGSORT, unsorted
 1, 2, 3, 4
*Nset, nset=RVE_AG_PosY, generate
 1,  7,  2
*Nset, nset=RVE_AG_PosY_RVE_AGSORT, unsorted
 5, 7
*Nset, nset=RVE_AG_PosY_RVE_AG_RED
 5, 7
*Nset, nset=RVE_AG_PosZ
 1, 2, 5, 6
*Nset, nset=RVE_AG_PosZ_RVE_AGSORT, unsorted
 6,
*Nset, nset=RVE_AG_PosZ_RVE_AG_RED
 6,
** Constraint: RVE_AG_1_RVE_AG_NegX_RVE_AGSORT_to_RVE_AG_PosX_RVE_AGSORT
*Equation
3
RVE_AG_PosX_RVE_AGSORT, 1, -1.
RVE_AG_NegX_RVE_AGSORT, 1, 1.
RP-Normal, 1, 1.
** Constraint: RVE_AG_1_RVE_AG_NegY_RVE_AGSORT_to_RVE_AG_PosY_RVE_AGSORT
*Equation
3
RVE_AG_PosY_RVE_AGSORT, 1, -1.
RVE_AG_NegY_RVE_AGSORT, 1, 1.
RP-Shear, 1, 1.
** Constraint: RVE_AG_1_RVE_AG_NegZ_RVE_AGSORT_to_RVE_AG_PosZ_RVE_AGSORT
*Equation
3
RVE_AG_PosZ_RVE_AGSORT, 1, -1.
RVE_AG_NegZ_RVE_AGSORT, 1, 1.
RP-Shear, 2, 1.
** Constraint: RVE_AG_2_RVE_AG_NegX_RVE_AGSORT_to_RVE_AG_PosX_RVE_AGSORT
*Equation
3
RVE_AG_PosX_RVE_AGSORT, 2, -1.
RVE_AG_NegX_RVE_AGSORT, 2, 1.
RP-Shear, 1, 1.
** Constraint: RVE_AG_2_RVE_AG_NegY_RVE_AGSORT_to_RVE_AG_PosY_RVE_AGSORT
*Equation
3
RVE_AG_PosY_RVE_AGSORT, 2, -1.
RVE_AG_NegY_RVE_AGSORT, 2, 1.
RP-Normal, 2, 1.
** Constraint: RVE_AG_2_RVE_AG_NegZ_RVE_AGSORT_to_RVE_AG_PosZ_RVE_AGSORT
*Equation
3
RVE_AG_PosZ_RVE_AGSORT, 2, -1.
RVE_AG_NegZ_RVE_AGSORT, 2, 1.
RP-Shear, 3, 1.
** Constraint: RVE_AG_3_RVE_AG_NegX_RVE_AGSORT_to_RVE_AG_PosX_RVE_AGSORT
*Equation
3
RVE_AG_PosX_RVE_AGSORT, 3, -1.
RVE_AG_NegX_RVE_AGSORT, 3, 1.
RP-Shear, 2, 1.
** Constraint: RVE_AG_3_RVE_AG_NegY_RVE_AGSORT_to_RVE_AG_PosY_RVE_AGSORT
*Equation
3
RVE_AG_PosY_RVE_AGSORT, 3, -1.
RVE_AG_NegY_RVE_AGSORT, 3, 1.
RP-Shear, 3, 1.
** Constraint: RVE_AG_3_RVE_AG_NegZ_RVE_AGSORT_to_RVE_AG_PosZ_RVE_AGSORT
*Equation
3
RVE_AG_PosZ_RVE_AGSORT, 3, -1.
RVE_AG_NegZ_RVE_AGSORT, 3, 1.
RP-Normal, 3, 1.
** 
** MATERIALS
** 
** units mpa
*parameter
vf=0.6
Em=4888.965569
E1=88366.522
E2=14155.981
G12=24489.72344
v12=0.28
v23=0.516
G23=E2*0.5/(1.0+v23)
**damage parameters
Xt=2500
Xc=2000
Yt=50
Yc=150
S_L=50
S_T=33
*Material, name=epoxy
*Elastic
<Em>, 0.41
*Material, name=fiber_glass
*Elastic,type=engineering constants
<E1>,<E2>,<E2>,<v12>,<v12>,<v23>,<G12>,<G12>,
<G23>
*material, name=composite
*mean field homogenization,formulation=mt
*constituent,material=epoxy,type=matrix,name=matrix
*constituent,material=fiber_glass,type=inclusion,shape=CYLINDER,name=fiber,direction=fixed
<vf>,, 1,0,0
** 
** BOUNDARY CONDITIONS
** 
** Name: PinnedNode Type: Symmetry/Antisymmetry/Encastre
*Boundary
RVE_AG_PinNode, PINNED
*Parameter
S12=2*S_L
S23=2*S_T
** ----------------------------------------------------------------
** 
** STEP: Load-Initial
** 
*Step, name=Load-Initial, nlgeom=NO, perturbation
RP_stress
*Static
** 
** OUTPUT REQUESTS
** 
** FIELD OUTPUT: F-Output-1
** 
*Output, field, variable=PRESELECT
*Element Output, directions=YES,micro
S,E
** 
** LOAD CASES
** 
*Load Case, name=Tensile-S11
*Cload
RP-Normal, 1, <Xt>
*End Load Case
*Load Case, name=Compression-S11
*Cload
RP-Normal, 1, -<Xc>
*End Load Case
*Load Case, name=Tensile-S22
*Cload
RP-Normal, 2, <Yt>
*End Load Case
*Load Case, name=Compression-S22
*Cload
RP-Normal, 2, -<Yc>
*End Load Case
*Load Case, name=Shear-S12
*Cload
RP-Shear, 1, <S12>
*End Load Case
*Load Case, name=Shear-S23
*Cload
RP-Shear, 3, <S23>
*End Load Case
*End Step