** 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
alpha=0.0
**matrix fail stresses
S22t_m=45
S22c_m=134
S33t_m=3.5
S33c_m=10.6
S12_m=33
S23_m=27
**fiber fail stresses
S11t_f=4013
S11c_f=3210
S12_f=61.4
**degradation
sdegm=1e-6
sdegf=1e-6
Emd=Em*sdegm
E1d=E1*sdegf
E2d=E2*sdegf
G12d=G12*sdegf
G23d=G23*sdegf
Gf_mt= 1.0
Gf_mc= 1.0
Gf_ft= 12.5
Gf_fc= 12.5
etam=etaf
*parameter table type, name = "matrixdam_params",parameters=10
FLOAT,,"Tensile stress limit in fiber direction"
FLOAT,,"Compressive stress limit in fiber direction"
FLOAT,,"Tensile stress limit in transverse direction"
FLOAT,,"Compressive stress limit in transverse direction"
FLOAT,,"Longitudinal shear stress limit"
FLOAT,,"Transverse shear stress limit"
FLOAT,,"Tensile fracture energy in transverse direction"
FLOAT,,"Compression fracture energy in transverse direction"
FLOAT,,"Viscosity parameter"
FLOAT,,"Maximum degradation"
*parameter table type, name = "fiberdam_params",parameters=8
FLOAT,,"Tensile stress limit in fiber direction"
FLOAT,,"Compressive stress limit in fiber direction"
FLOAT,,"Longitudinal shear stress limit"
FLOAT,,"Alpha"
FLOAT,,"Tensile fracture energy in 1-direction"
FLOAT,,"Compression fracture energy in 1-direction"
FLOAT,,"Viscosity parameter"
FLOAT,,"Maximum degradation"
*Material, name=epoxy
*Elastic,dependencies=1
<Em>, 0.41,0,0
<Emd>,0.41,0,1
*userdefinedfield
*depvar
4,
1,"DAMAGET","Tensile damage variable"
2,"DAMAGEC","Compression damage variable"
3,"DAMAGET_DV","Regularized damage variable in tension"
4,"DAMAGEC_DV","Regularized damage variable in compression"
*parameter table,type="matrixdam_params"
<S22t_m>,<S22c_m>,<S33t_m>,<S33c_m>,<S12_m>,<S23_m>,<Gf_mt>,<Gf_mc>,
<etam>,<sdegm>
*Material, name=fiber_glass
*Elastic,type=engineering constants,dependencies=1
<E1>,<E2>,<E2>,<v12>,<v12>,<v23>,<G12>,<G12>,
<G23>,0,0
<E1d>,<E2d>,<E2d>,<v12>,<v12>,<v23>,<G12d>,<G12d>,
<G23d>,0,1
*userdefinedfield
*depvar
4,
1,"DAMAGET","Tensile damage variable"
2,"DAMAGEC","Compression damage variable"
3,"DAMAGET_DV","Regularized damage variable in tension"
4,"DAMAGEC_DV","Regularized damage variable in compression"
*parameter table,type="fiberdam_params"
<S11t_f>,<S11c_f>,<S12_f>,<alpha>,<Gf_ft>,<Gf_fc>,<etaf>,<sdegf>
*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