search

TarletonGroup / CrystalPlasticity

CP UMAT and CZM UEL for Abaqus
138 stars 64 forks source link

Question about slipmodel=2 #57

Closed louis96108 closed 3 months ago

louis96108 commented 3 months ago

I've tried to use slipmodel2 in two ways, but I'm having problems with both. The model I used is from‘Example-Polycrystal with PROPS’ (1) I modified the slip model in case2 and changed it to slipmodel=2 with the corresponding parameter slipparam(1) = 2.54d6 slipparam(2) = 0.85，slipparam(3) = 1.36，slipparam(4) = 1.6d-19 (2) I changed the Material ID into 7, and try to use this material that inherently contains slipmodel2 However, in both cases there were problems and I didn't know how to fix them. wrong

EralpDemir commented 3 months ago

Hi, I have tested these example cases with the latest version. They work, I added the updated files.

When using PROPS, these material properties shall be updated for all of the grains in the mesh, which is somewhat tedious.

I am sending a single crystal case with an element that uses materialID-7 for the given properties. I am also sending the usermaterial.f file that contains these values. I ran it and it worked for those material properties.

I have added an if statement to the slip law for slip that is negligibly small. Please use the latest version (v.20) in the Github or in the attached.

Let me know if you have any other issues. I close this issue this with this.

Best, Eralp

v2.20.zip

louis96108 commented 3 months ago

Dear Eralp Demir, Thank you very much for replying to me and updating the code. With the new code, the single crystal model with MaterialID=7 can be fully calculated.However, there are still problems with the use of double exponent law, I added a custom material (case(13)) to the usermaterials.f file, using the double exponent law (slipmodel=2) slip model and the Voce type hardening (hardeningmodel=1) hardening model, the relevant parameters are derived from the reference and the example on Github. The rest of the files have not changed, and the calculations have not converged. I provided the file in the attachment and I don't know if the different slip models, hardening models in the code can be arbitrarily combined to get a custom material. I'm look forward to hearing from you. Wish you a nice day!

Best regards, Siyu Lu

At 2024-05-20 17:15:05, "Eralp Demir" @.***> wrote:

Hi, I have tested these example cases with the latest version. They work, I added the updated files.

When using PROPS, these material properties shall be updated for all of the grains in the mesh, which is somewhat tedious.

I am sending a single crystal case with an element that uses materialID-7 for the given properties. I am also sending the usermaterial.f file that contains these values. I ran it and it worked for those material properties.

I have added an if statement to the slip law for slip that is negligibly small. Please use the latest version (v.20) in the Github or in the attached.

Let me know if you have any other issues. I close this issue this with this.

Best, Eralp

v2.20.zip

— Reply to this email directly, view it on GitHub, or unsubscribe. You are receiving this because you authored the thread.Message ID: @.***>