if(abs(My)>1): epsilonZPos= scc.getSection().getStrain(0.0,1.0) fibraEpsZMax= prop_statistics.getItemWithMaxProp(scc.getSection().getFibers(),"getMaterial.getStrain") epsZMax= fibraEpsZMax.getMaterial().getStrain() yEpsZMax= fibraEpsZMax.getPos().x zEpsZMax= fibraEpsZMax.getPos().y epsMy= scc.getSectionDeformationByName("defMy")#defMy Mz= scc.getStressResultantComponent("Mz") MzTot+= Mz if(abs(Mz)>1): epsilonYPos= scc.getSection().getStrain(1.0,0.0) fibraEpsYMax= prop_statistics.getItemWithMaxProp(scc.getSection().getFibers(),"getMaterial.getStrain") epsYMax= fibraEpsYMax.getMaterial().getStrain() yEpsYMax= fibraEpsYMax.getPos().x zEpsYMax= fibraEpsYMax.getPos().y fibraEpsYMin= prop_statistics.getItemWithMinProp(scc.getSection().getFibers(),"getMaterial.getStrain")#IMinProp("getMaterial.getStrain") epsYMin= fibraEpsYMin.getMaterial().getStrain() yEpsYMin= fibraEpsYMin.getPos().x zEpsYMin= fibraEpsYMin.getPos().y epsilonYPos= scc.getStrain(yEpsYMin,zEpsYMin) epsilonYPosTeor= -Mz/(E*Iz)*yEpsYMin epsMz= scc.getSectionDeformationByName("defMz")#defMz lPatterns.removeFromDomain(key) dom= preprocessor.getDomain dom.revertToStart() ratio1= abs((epsilonYPos-epsilonYPosTeor)/epsilonYPosTeor) ratio2= abs((epsYMin+epsilonYPosTeor)/epsilonYPosTeor) ratio3= abs((epsilonZPos-5*F/(E*Iy))) ratio4= abs((epsilonYPos-6*F/(E*Iz)))