Hello:
I am quite new to both abaqus and fortran programming. I am working on a project where I am tying phase field to fracture propagation where the void phase is taken into the kinematics of the crack propagation. I will not go too much into detail!
I have attached 3 files with my question here:
a) umat_sdvini3.f (my user subroutine file that has both sdvini and umat in it) b) Job-User-mod-9.inp (the .inp input file that I created with the help of Abaqus CAE) c) Job-User-15.msg (the file that contains the error message!)
When I am telling Abaqus to run the job with the user subroutine, the job gets aborted for some reason! In the previous job runs (11 to 14), when in the subroutine file, the if/else statement for the state variable was commented out, everything ran fine!
With the second if/else statement uncommented in the attached fortran file, there is an error and the job is getting aborted!
The fortran 77 code for my user subroutine: (it is the second if/else statement command block not running properly! I thought I could access the values of the solution dependent state variables via coords(1) which is the x-coordinates!
Could any of you please help me with what it could be? Please!
Thank you,
Mousumi
Hello Mousumi Ghosh
Why are you still using Fortran 77 program?
The programming style of fortran 77 is dated and cumbersome. Please upgrade to either Fortran 90 or Fortran 95. You can obtain free Fortran 90 or 95 from
http://force.lepsch.com/2009/05/downloads.html
Hello:
I wish Abaqus allowed their user subroutines to be written in fortran free format, but, unfortunately, that is not the case. That is why I am forced to write the subroutines in fortran77 language!
1- nstatev is not determined in your code. Do not mix up, ntens are prescribed within the abaqus and are input and do not need to be determined, but nstatev are your input and you must determine its dimension both within the umat and your above subroutine. By defining the value of nstatev as a parameter or by directly entering it : i.e. statev(9) or whatever.
2- in fortran 77 you cannot define arquments after executive expressions. So move all definitions like :
real f_i=1.0e-10
real w_bar=0.25d0
real walpha=2.d0
to above of the first executive espression: i.e. above of the
ddsdde = 0.d0
3- You have a mistype:
statev(1)=state(1) + kinc
it should be maybe
statev(1)=statev(1) + kinc
4- not recommended to use END DO in F77
recommended is:
DO "0000"
"0000" CONTINUE
format
5- The second if is wrong
comparing a Decimal value with an integer brings error in F77.
change it to this format:
if(coords(1).le.0.0 .and. coords(2).eq.0.0) then ....
------------------------
I hope these notes due to your code attached could help
Sincerely
Pejman
- Badges
- Science topic