Abaqus python script for writing to text file

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Abaqus python script for writing to text file

Abaqus Users mailing list
Hi,
I have some abaqus models for modeling soil structure interaction of buried pipelines. I have created the model using Thick Pipe (PIPE31) elements and also using Elbow (Elbow31B) elements. the reason for the two approaches is to see which one is better suited for the analysis. i created some python script files that access the odb and basically extract stresses (S11, S22, etc.) and nodal displacements to a text file. i then import the text files into excel and plot the behavior. When we use either the Pipe or Elbow elements, there are numerous points on the inner fiber or outer fiber that we can output the stresses. for pipe section i use: odb.steps['apply_temperature'].frames[5].fieldOutputs['S'].getSubset(region=PipeElements,position=CENTROID,elementType='PIPE31') for elbow section i use: odb.steps['apply_temperature'].frames[5].fieldOutputs['S'].getSubset(region=PipeElements,position=CENTROID,elementType='ELBOW31B') in either of the two methods above i can extract the stresses, S11, S22, S33, S12 and in excel i have written macros that sieves thru the data, sort, and extract the data. 1) Is there a way to extract the max tension or compression for each element? or max envelope? say for longitudinal stresses (S11) for the inner and outer fiber? 2) can we extract the mises stress? envelope of the inner or outer fiber? 3) what are the lines or changes that need to be made? i changed fieldOutputs['S'] to fieldOutputs['Mises'] as per manual and it gives an error. currently i painfully calculate the mises stresses and then seieve the data, sort and extract for plotting.
I can share the python script file if necessary.. any suggestion / advise is most welcome.


-Milind
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

RE: Abaqus python script for writing to text file

Abaqus Users mailing list
You just need to take advantage of the tools available in Python (and/or NumPy).  For example:

step = odb.steps[‘apply_temperature’]
frame = step.frames[5]
field = frame.fieldOutputs[‘S’].getSubset(region=…)
s11Max = max([x.data[0] for x in field.values])

Regards

Dave Lindeman
Lead Research Specialist
SEMS Research Laboratory
3M Center 235-3G-08
St. Paul, MN 55144
651-733-6383

From: [hidden email] [mailto:[hidden email]]
Sent: Tuesday, May 02, 2017 8:45 AM
To: [hidden email]
Subject: [EXTERNAL] [Abaqus] Abaqus python script for writing to text file



Hi,

I have some abaqus models for modeling soil structure interaction of buried pipelines. I have created the model using Thick Pipe (PIPE31) elements and also using Elbow (Elbow31B) elements. the reason for the two approaches is to see which one is better suited for the analysis. i created some python script files that access the odb and basically extract stresses (S11, S22, etc.) and nodal displacements to a text file. i then import the text files into excel and plot the behavior. When we use either the Pipe or Elbow elements, there are numerous points on the inner fiber or outer fiber that we can output the stresses. for pipe section i use: odb.steps['apply_temperature'].frames[5].fieldOutputs['S'].getSubset(region=PipeElements,position=CENTROID,elementType='PIPE31') for elbow section i ! use: odb.steps['apply_temperature'].frames[5].fieldOutputs['S'].getSubset(region=PipeElements,position=CENTROID,elementType='ELBOW31B') in either of the two methods above i can extract the stresses, S11, S22, S33, S12 and in excel i have written macros that sieves thru the data, sort, and extract the data. 1) Is there a way to extract the max tension or compression for each element? or max envelope? say for longitudinal stresses (S11) for the inner and outer fiber? 2) can we extract the mises stress? envelope of the inner or outer fiber? 3) what are the lines or changes that need to be made? i changed fieldOutputs['S'] to fieldOutputs['Mises'] as per manual and it gives an error. currently i painfully calculate the mises stresses and then seieve the data, sort and extract for plotting.
I can share the python script file if necessary.. any suggestion / advise is most welcome.

-Milin! d





3M security scanners have not detected any malicious content in this message.
Click here<https://spam.mmm.com:443/pem/pages/digestProcess/digestProcess.jsf?content=aedaaa864ecbae9483c6fa90572e3f0356ac76184387f507251011ddfb5d8327c6c4ff5fb22097b592aab29dac82584f30f60c161a797f8d8d8a621478a2f25250f75f3a4b570e12125b08ef0dd8e5d1dbc7a045714a146255659761471276b985c737a3334ed90b1ff78053eaa8c27c6dc37ee21f7f3ae2a65d7f79d67ffbd5ff53b47a8a233403> to report this email as spam







[Non-text portions of this message have been removed]


Loading...