g_f if(solid2.sl11>=0,solid2.sl11/Sigmats1_f,-solid2.sl11/Sigmacs1_f)-1 "Failure criterion of fiber"
sf_f if(solid2.sl11>=0,Sigmats1_f/solid2.sl11,Sigmacs1_f/-solid2.sl11) "Safety factor of fiber failure criterion"
g_m solid2.mises^2/(Sigmats_m*Sigmacs_m)+(1/Sigmats_m-1/Sigmacs_m)*solid2.I1s-1 "Failure criterion of matrix"
a solid2.mises^2/(Sigmats_m*Sigmacs_m) "Coefficient of quadratic equation"
b (1/Sigmats_m-1/Sigmats_m)*solid2.I1s "Coefficient of quadratic equation"
c -1 "Coefficient of quadratic equation"
r1 (-b+sqrt(b^2-4*a*c))/(2*a) "First root"
r2 (-b-sqrt(b^2-4*a*c))/(2*a) "Second root"
sf_m r1*(r1>0)*(r2<=0)+r2*(r1<=0)*(r2>0)+min(r1,r2)*(r1>0)*(r2>0) "Safety factor of matrix failure criterion"
Tax solid2.sxx*solid2.nx+solid2.sxy*solid2.ny+solid2.sxz*solid2.nz "Traction, x-component"
Tay solid2.sxy*solid2.nx+solid2.syy*solid2.ny+solid2.syz*solid2.nz "Traction, y-component"
Taz solid2.sxz*solid2.nx+solid2.syz*solid2.ny+solid2.szz*solid2.nz "Traction, z-component"
Ta1 solid2.sysT11*Tax+solid2.sysT21*Tay+solid2.sysT31*Taz "Traction, tangential component"
Ta2 solid2.sysT12*Tax+solid2.sysT22*Tay+solid2.sysT32*Taz "Traction, circumferential component"
Ta3 solid2.sysT13*Tax+solid2.sysT23*Tay+solid2.sysT33*Taz "Traction, normal component"
g_i (if(Ta3>=0,Ta3,0)/Sigman_i)^2+(Ta1/Sigmas_i)^2-1 "Failure criterion of interface"
sf_i 1/sqrt(g_i+1) "Safety factor of interface failure criterion"
eXX comp2.at2(X0,Y0,Z0,comp2.lshell.atxd1(xd,comp2.lshell.eXX)) "Elastic strain, XX-component from layered shell interface"
eYY comp2.at2(X0,Y0,Z0,comp2.lshell.atxd1(xd,comp2.lshell.eYY)) "Elastic strain, YY-component from layered shell interface"
eZZ comp2.at2(X0,Y0,Z0,comp2.lshell.atxd1(xd,comp2.lshell.eZZ)) "Elastic strain, ZZ-component from layered shell interface"
eXY comp2.at2(X0,Y0,Z0,comp2.lshell.atxd1(xd,comp2.lshell.eXY)) "Elastic strain, XY-component from layered shell interface"
eYZ comp2.at2(X0,Y0,Z0,comp2.lshell.atxd1(xd,comp2.lshell.eYZ)) "Elastic strain, YZ-component from layered shell interface"
eXZ comp2.at2(X0,Y0,Z0,comp2.lshell.atxd1(xd,comp2.lshell.eXZ)) "Elastic strain, XZ-component from layered shell interface"