# Solution #analisis= predefined_solutions.simple_static_linear(feProblem) solver = predefined_solutions.SolutionProcedure() analisis = solver.simpleLagrangeStaticLinear(feProblem) #analisis= solver.simpleTransformationStaticLinear(feProblem) result = analisis.analyze(1) nodes.calculateNodalReactions(False, 1e-7) reactionNode10 = n10.getReaction ratio1 = reactionNode10.Norm() svdReactionNodes = nodal_reactions.getReactionFromNodes( nodes, "UVWRxRyRz", elem.getNodes.getExternalNodes) actionNode10 = xc.Vector(loadOnDOFs) actionNode10Norm = actionNode10.Norm() svdAction = nodal_reactions.getSlidingVectorsSystemfromSlidingVector( "UVWRxRyRz", n10.get3dCoo, actionNode10) svdResid = svdReactionNodes + svdAction ratio2 = svdResid.getResultant().getModulo() / actionNode10Norm ratio3 = svdResid.getMoment().getModulo() / actionNode10Norm reactionNode20 = n20.getReaction ratio4 = reactionNode20.Norm() actionNode20 = xc.Vector([0, 0, 0, 0, 0, 0]) actionNode20Norm = actionNode20.Norm() # print "svdAction= ", svdAction # print "svdReactionNodes= ", svdReactionNodes # print "svdResid= ", svdResid # print "ratio1= ", ratio1 # print "ratio2= ", ratio2 # print "ratio3= ", ratio3
#We add the load case to domain. lPatterns.addToDomain("0") # Solution solver= predefined_solutions.SolutionProcedure() analisis= solver.simpleTransformationStaticLinear(feProblem) result= analisis.analyze(1) nodes.calculateNodalReactions(False,1e-7) reactionNode10= n10.getReaction ratio1= reactionNode10.Norm() svdReactionNodes= nodal_reactions.getReactionFromNodes(nodes,"UVWRxRyRz",elem.getNodes.getExternalNodes) actionNode10= xc.Vector(loadOnDOFs) actionNode10Norm= actionNode10.Norm() svdAction= nodal_reactions.getSlidingVectorsSystemfromSlidingVector("UVWRxRyRz",n10.get3dCoo,actionNode10) svdResid= svdReactionNodes+svdAction ratio2= svdResid.getResultant().getModulo()/actionNode10Norm ratio3= svdResid.getMoment().getModulo()/actionNode10Norm # print "svdAction= ", svdAction # print "svdReactionNodes= ", svdReactionNodes # print "svdResid= ", svdResid # print "ratio1= ", ratio1 # print "ratio2= ", ratio2 # print "ratio3= ", ratio3 # print "RN2= ", RN2 # print "RN3= ", RN3 # print "RN4= ", RN4 import os