Dear all,
Actually, I am wondering how can I write specific data to a text file via Umat as. I need this data for further purposes?
Many thanks for your help
Dear Ali Nassr,
You can write the data to text file using the following commands:
*************************************************************
open(105,file='D:\Ranjith\My_work\tt.txt',position='append')
WRITE(105,*)Z, E1
close(105)
*************************************************************
The above lines are self explanatory :-) I'm writing two variable Z and E1 to the text file tt.txt.
Thanks & Regards,
Ranjith.
Dear Ali Nassr,
You can write the data to text file using the following commands:
*************************************************************
open(105,file='D:\Ranjith\My_work\tt.txt',position='append')
WRITE(105,*)Z, E1
close(105)
*************************************************************
The above lines are self explanatory :-) I'm writing two variable Z and E1 to the text file tt.txt.
Thanks & Regards,
Ranjith.
Hi Ali,
I usually do that in Abaqus viewer during post processing, via Python script. Before running my model, I specify I want it to save the Sate Variables (SDV's) and the I read them and export them to a txt. file.
Here's the code I use to export Mises Stress, Triaxiality and some state variables from a particular Integration point of a particular element in the region o interest in my model (frame-by-frame in the load-step):
for frame in frames:
# getting frame ID and time at this frame
frameID = frame.frameId
time = frame.frameValue
# getting quantities at the integration point closest to the origin:
# Stresses = frame.fieldOutputs['Snodal'].getSubset(region=neckCenter)
# ... Stresses:
Stresses = frame.fieldOutputs['S'].getSubset(region=e0)
q = Stresses.values[indSecPoint-1].mises
p = Stresses.values[indSecPoint-1].press
if q==0:
triax = 0.0
else:
triax = -p/q
# ... State variables:
# ...... Matrix material eqiv. plastic strain:
State1 = frame.fieldOutputs['SDV1'].getSubset(region=e0)
PEmatrix = State1.values[indSecPoint-1].data
# ...... Void Ratio:
State2 = frame.fieldOutputs['SDV2'].getSubset(region=e0)
voidRatio = State2.values[indSecPoint-1].data
# getting diameter reduction at neck
NeckDisp = frame.fieldOutputs['U'].getSubset(region=neckDiam)
DiamReduc = -2*NeckDisp.values[0].data[0]
# getting Reaction Force of the test
Reaction = frame.fieldOutputs['RF'].getSubset(region=RP)
RF = Reaction.values[0].data[1]
f.write( '%i \t %f \t %f \t %f \t %f \t %f \t %f \n' %
(frameID, time, DiamReduc, RF, triax, PEmatrix, voidRatio) )
f.close()
Regards,
Read the Utility from the following link:
http://www.maths.cam.ac.uk/computing/software/abaqus_docs/docs/v6.12/books/sub/default.htm
I know this is a 3 year old forum but I hope some one can help me regarding a similar problem in writing UMAT.
I want to write stress values in the .dat or .msg file of Abaqus I am writing it as:
WRITE(6,*) (STRESS(i),i=1,NDI)
but it is not showing results
Hi, I recently needed to write my subroutine results into a file, I used the following script :
open(105,file='D:\tt.txt',position='append')
WRITE(105,*)x,dist
close(105)
So for you it should be :
open(105,file='D:\tt.txt',position='append')
WRITE(105,*) (STRESS(i),i=1,NDI)
close(105)
Try and tell me if it worked.
Have a good day,
Sam
Guys, Thanks for your information. I also want to know how I can write the output just for centroid of each element. That command writes it for all the integration points.
There exists a default channel (6) for the "*.dat " file, where * stands for the job-name of your project on ABAQUS CAE. So, a statement on UMAT composed using FORTRAN language, viz.,
write(6,*) "statement to print"
will append "statement to print" in the "*.dat" file upon running.
- Badges
- Science topic