I think this is correct, except I think the final END is not needed. I am currently using 2 subroutines (SDVINI and UMAT) and the syntax is:
SUBROUTINE SDVINI(STATEV,COORDS,NSTATV,NCRDS,NOEL,NPT,
LAYER,KSPT)
.........
RETURN
END
SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD...
Abaqus uses the material property based on which value of the vield variable is closest to the actual result. In your case, the actual value of strain is compared to the value of your field variables. It checks which value in your table is closest to the strain value obtained. Once this value is...
I havent worked with the material orientation in the remodelling algorithm. I assumed that the material (bone) was isotropic. The Youngs modulus changed simply on the total strain for each element. If you think that it might be the issue, first try it without your definition of material...
Thank you all for the suggestions. I did move the cup towards the mold and now I get a solution.
The only problem now is that the cup just won't move downwards into the shape of the mold. I use relatively large forces (4000N) and the materials that I use have a very low E-modulus and low yield...
I think that if you have a compilation error you have to check for syntax-mistakes in your code. Anaqus creates a log file in your directory (jobname.log). If you open that it tells you the reason why it can not compile (e.g. EXAMPLE.for(27): error #5082: Syntax error, found ')' when expecting...
I dont have a lot of experience with penetration errors, but I think it is best to indeed first do the simulation without your remodelling algorithm, then see if you obtain some errors. Try to solve the errors first before you apply your algorithm. For the remodelling I wrote a Fortran code...
Hi all,
I am trying to do a simulation where a hemispherical cup is placed into some sort of mold (see image). The cup is a little bit larger than the mold, and the idea is that the cup is forced in and deforms (in the plasticity range). Then, by friction the cup should stay in place. My goal...
I recently did a similar simulation where the bone-properties were adapted according to the strain in the bone. For this I used an UMAT subroutine with fortran. It worked for me to use c3D4 tet elements. The model also included an implant and I had no problems with the contact between the...
Thanks all for the suggestions.
Although the shape is very straightforward, I could not recreate it in abaqus due to the exact orientation and location of the part, which is critical.
I followed the suggestion of Mustaine and checked in the manual how scripting works. It was indeed exactly...
Greetings all,
Maybe a simple question, but I can't find an easy way to do this. I have a hemisphere with a diameter of 25mm, see attachment. My goal is to define two elemen sets to apply different material properties. One element sets that contains all inner elements (let's say all the...
Hi all,
I did a simulation with Abaqus with 1 static load. After the simulation has been completed, my disk space decreased from 91gb of free space to 79 gb of free space. I checked the sizes of all the files (odb, cae etc) and they do not even add up to 0.5 gb. When I deleted all the files, I...
I'm quite positive you cannot do this, since the poissons ratio is by definition the ratio between the tranvserse train to the axial strain. You need both values in order to be able to calculate the ratio.
Hi all,
The problem is solved, so I thought I write a solution here if anyone else in running into these problems.
Apparently some keywords are not recognized by Abaqus CAE from the input file, even if the expression might be correct. I downloaded an example input file for an SDVINI subroutine...
Hi all,
I wrote an UMAT-subroutine which changes the Youngs modulus based on the strains. It seems to work, except that the initial conditions of the state variables are not read by Abaqus. I defined the intial state variables in an SDVINI subroutine, and put this in the same file as the UMAT...
