Python for Abaqus - 1.4: Post-Processing - SS I-beams

Final script (watch until video 1.5 to understand the full thing):
github.com/luisantos090/Abaqu...
Follow this link to avoid issues with Notepad ++
stackoverflow.com/questions/4...

Пікірлер: 14

  • @shervinsafaeifaegh1540
    @shervinsafaeifaegh1540 Жыл бұрын

    Thank you so much for such a brilliant tutorial and videos indeed!

  • @luisantos090

    @luisantos090

    Жыл бұрын

    Thank you Shervin. If you have any recommendations, let me know. I feel the post processing was the hardest part to demonstrate but I'm happy you enjoyed it

  • @shervinsafaeifaegh1540

    @shervinsafaeifaegh1540

    Жыл бұрын

    @@luisantos090 Definitely yes. But the procedure you mentioned in the video is really helpful. As some suggestions, it would be more than great if you extend the application of this procedure to some more expertized fields like the iteration in Eigen buckling analysis of a member under different shape factors or imperfections, XFEM, and propagation of cracks induced by low cycle fatigue recognized with python script. Anyway, thank you so much for distributing your knowledge for free!

  • @luisantos090

    @luisantos090

    Жыл бұрын

    @@shervinsafaeifaegh1540 thank you for that. I will show the buckling analysis and the imperfections in future videos. Regarding cracks and XFEM, it's not my field but I hope you can solve it nonetheless

  • @azizseddiki6310
    @azizseddiki6310 Жыл бұрын

    Thank you very much for this video tutorial, I have been able to solve many doubts I had about this topic. I would like to find out if there is any automatic way to detect the maximum value of the vector "values" and not look for it manually and put the value. That is, instead of putting: print(odb.steps['Step-1'].frames[1].fieldOutputs['S'].values[1]). put something like this: print(odb.steps['Step-1'].frames[1].fieldOutputs['S'].values[X].mises) X being the one that has the maximum value of mises Is there a way to save this data for each run to be able to compare it with the other runs? Thank you again

  • @luisantos090

    @luisantos090

    Жыл бұрын

    Hi Aziz. Thank you for your comment, that means a lot. There is indeed such command, I'll send it to you in 30 mins or so ;)

  • @azizseddiki6310

    @azizseddiki6310

    Жыл бұрын

    @@luisantos090 Thank you very much for your quick response, and for the effort. What I want is to save the maximum value of the Von Mises voltage for each run, either in an excel table or any other file, or just return those values to me via Abaqus command box.

  • @luisantos090

    @luisantos090

    Жыл бұрын

    @@azizseddiki6310 To save the data on an excel file, I use the following script (you will need to adapt it a bit but it should work: opFile = NewFolder+'/'+'DatainExcel.csv' try: opFileU = open(opFile,'w') opFileU.write("%10s,%10s "%('Force', 'U2') ) except IOError: PrintToScreen('cannot open opFILE') exit(0) for i in range(NrOfSteps): displacement = odb.steps['Step-1'].frames[i].fieldOutputs['U'].values[23].data[1] force = RefPointOutput[i][1] Displacements_midspan.append(-displacement) Applied_forces.append(-force/1000.0) opFileU.write("%10f,%10f " % (force, displacement)) opFileU.close()

  • @luisantos090

    @luisantos090

    Жыл бұрын

    @@azizseddiki6310 To get maximum values for displacements (you will need to dig to find the equivalent for von Mises voltage) I use the following code: Disps = odb.steps['Step-1'].frames[i].fieldOutputs['U'].getSubset(position = NODAL).bulkDataBlocks[0].data followed by: np.max(Disps,axis=0)[0] This will give me the maximum in the x-direction. You can do "[1]" or "[2]" for the y- and z-directions, respectively. Also, for this to work you need, at the top of your file, to "import numpy as np". I hope this helps :)

  • @azizseddiki6310

    @azizseddiki6310

    Жыл бұрын

    @@luisantos090 Thank you very much for your efforts to find a solution. I will try your instructions and give you feedback since I have tested it. I appreciate the time spent to solve my doubt. Thank you.

  • @hasanh3968
    @hasanh3968 Жыл бұрын

    Hello, how can i send email for you?

  • @luisantos090

    @luisantos090

    Жыл бұрын

    Hi Hasan, You can contact me on LinkedIn if you'd like: www.linkedin.com/in/luisantos090

  • @hasanh3968

    @hasanh3968

    Жыл бұрын

    ​@@luisantos090 Hi Luis, Thank you. there is a problem that i could not send message for you. I want to get temperature, y-coordinate and pressure for some nodes. As usual when i get node number in the visualisation we see the node label. But the data for each node on their node list index number. what is your solution for this problem? actually node list index number and node label are different.

  • @luisantos090

    @luisantos090

    Жыл бұрын

    @Hasan H try to make the index equal to the node number minus 1. So, for example, if you want to know the output for node 2367, the index should be 2366. Give it a try and let me know if it doesn't work