Beam results, dpf, Mechanical

Options

Is there any way to access beam results using Python?

Tagged:

Best Answer

  • Rohith Patchigolla
    Rohith Patchigolla Member, Moderator, Employee Posts: 115
    First Answer First Anniversary Name Dropper Solution Developer Community of Practice Member
    Answer ✓
    Options

    Hello @Joshua Fillion

    Here is a simple script using DPF that you can use in Mechanical Scripting Console to print beam reactions (X axis. This can be extended to Y, Z and Total as well) for all the beam connections in the Tree. You can also integrate this script into "Python Code" object with "Target Callback" --> "After Post" in Mechanical to export this data into excel (for example).

    import os
    import wbjn
    
    analysis = DataModel.AnalysisList[0]
    first_solution = analysis.Solution
    solver_data = first_solution.SolverData    
    allBeams = DataModel.GetObjectsByType(DataModelObjectCategory.Beam)
    for beam in allBeams:
        beamData = solver_data.GetObjectData(beam)
        beamMatid = beamData.MaterialId
        beamElem = solver_data.ElementIdsByMaterialId(str(beamMatid))
    
        import mech_dpf
        import Ans.DataProcessing as dpf
        #Result Data    
        dataSource = dpf.DataSources(analysis.Solution.ResultFilePath)
    
        # Read mesh in results file
        mesh_op = dpf.operators.mesh.mesh_provider() 
        mesh_op.inputs.data_sources.Connect(dataSource)
        my_mesh = mesh_op.outputs.mesh.GetData()
    
        #time scoping
        timelist = dpf.operators.metadata.time_freq_provider(data_sources=dataSource).outputs.gettime_freq_support().TimeFreqs.Data
        my_time_scoping = range(1, len(timelist)+1)
        #scoping
        my_scoping = dpf.Scoping()
        my_scoping.Location = "Elemental"
        my_scoping.Ids = beamElem
    
        time_data = []
        Contact_TotalReaction = []
        for setnum in range(len(timelist)):
            CX = dpf.operators.result.mapdl.smisc(time_scoping=my_time_scoping,mesh=my_mesh,data_sources=dataSource,item_index = 1,mesh_scoping=my_scoping).outputs.fields_container.GetData()[setnum]
            #print timelist[setnum],',',(CX.Data)*sum(CX.Data) + sum(CY.Data)*sum(CY.Data) + sum(CZ.Data)*sum(CZ.Data))
            #reactionTotal = sqrt(sum(CX.Data)*sum(CX.Data) + sum(CY.Data)*sum(CY.Data) + sum(CZ.Data)*sum(CZ.Data))
            time_data.append(timelist[setnum])
            Contact_TotalReaction.append(CX.Data) 
    print(Contact_TotalReaction)
    

Answers

  • Chris Harrold
    Chris Harrold Member, Administrator, Employee Posts: 129
    First Answer First Comment First Anniversary Ansys Employee
    admin
    Options

    @Ramdane or someone in @AKD-Scripting-Team might have some thoughts here?

  • Joshua Fillion
    Joshua Fillion Member Posts: 2
    First Comment
    Options

    Hey, thanks guys. At the time I ended up utilizing data written out to .csv files for the elements of interest via Ansys scripting and then calling a dpf python script separately for plotting via fields factory. Somewhat bubblegum and duct tape, but it got the job done too!