Analysis of Load Optimization in Solid Rocket Motor Propellant Grain with Pressure Cure

Algorithm 1

“for step in o.steps.values():

  for frame in step.frames:

    fopstress = frame.fieldOutputs[‘S’]

    stress = fopstress.getSubset(region = \

    propellant_e)

    stressOutputExists = TRUE

    for stressValue in stress.values:

     if (not maxValue or stressValue.mises\

       > maxValue.mises):

       maxValue = stressValue

       maxStep,maxFrame = step, frame

elementL = maxValue.elementLabel

  integration_point = maxValue.integrationPoint”