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...
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
Thank you so much for such a brilliant tutorial and videos indeed!
@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
Жыл бұрын
@@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
Жыл бұрын
@@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
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
Жыл бұрын
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
Жыл бұрын
@@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
Жыл бұрын
@@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
Жыл бұрын
@@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
Жыл бұрын
@@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.
Hello, how can i send email for you?
@luisantos090
Жыл бұрын
Hi Hasan, You can contact me on LinkedIn if you'd like: www.linkedin.com/in/luisantos090
@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
Жыл бұрын
@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