How do I perform Load Case combination of principal stresses using DPF (Python Result) in Mechanical
Ayush Kumar
Member, Moderator, Employee Posts: 457
✭✭✭✭
How do I perform Load Case combination of principal stresses using DPF (Python Result) in Mechanical?
Example: LC1 LC2
User defined results Solution Combination in Mechanical does a combination on the calculated results eg. Principal or SEQV, which is not the right way. The combination should be done on the stress tensor level and then Principal or SEQV should be calculated based on the new stress tensor. How do I do this using DPF and plot the results as contour on the body in Mechanical?
Tagged:
2
Answers

You can use the following code to plot load case combination on S1 (Maximum Principal Stress), the time scoping needs to be adjusted as per your analysis.
For S2 and S3 use
p_inv.outputs.field_eig_2
andp_inv.outputs.field_eig_3
respectively.def post_started_301(sender, analysis):# Do not edit this line define_dpf_workflow(analysis) def define_dpf_workflow(analysis): import mech_dpf import Ans.DataProcessing as dpf dataSource = dpf.DataSources(analysis.ResultFileName) # Instantiate Stress Tensor stress_tensor = dpf.operators.result.stress() # For LC1 timeScop = dpf.Scoping() timeScop.Ids = [1, 2] stress_tensor.inputs.time_scoping.Connect(timeScop) stress_tensor.inputs.data_sources.Connect(dataSource) # Scale LC2 to 1 stress_tensor_lc2_sc = dpf.operators.math.scale(field=stress_tensor.outputs.fields_container.GetData()[1], ponderation=1.0) # Add load cases stress_tensor_combi = dpf.operators.math.add(fieldA=stress_tensor.outputs.fields_container.GetData()[0], fieldB=stress_tensor_lc2_sc) # Principal Invariants  Principal Stresses p_inv = dpf.operators.invariant.principal_invariants() p_inv.inputs.field.Connect(stress_tensor_combi) # Plot Combined S1 or Maximum Principal Stress dpf_workflow = dpf.Workflow() dpf_workflow.Add(p_inv) # Plot contour dpf_workflow.SetOutputContour(p_inv.outputs.field_eig_1) wf_id = dpf_workflow.Record('wf_id', True) this.WorkflowId = dpf_workflow.GetRecordedId()
0