search

TarletonGroup / CrystalPlasticity

CP UMAT and CZM UEL for Abaqus
152 stars 72 forks source link

Issue caused by units #64

Closed louis96108 closed 4 months ago

louis96108 commented 4 months ago

Hi, In the globalvariables.f, the Boltzmann constant is set to 1.380649d-23 [m2 kg s-2 K-1]. And in hardening.f, line 298, the codes for calculating k2 are as follows: k2 = k1Xb/g*(1.-KBT/D/b*3dlog(pdot/pdot0)) In the answers to other questions, the use of μm and MPa was mentioned. So, after putting the parameter from the Document into the code: 图片 I feel that the part of the calculation that the Boltzmann constant is involved in is problematic due to the unit mismatch 图片 I look forward to hearing from you, thank you very much

Best regards！ Louis

EralpDemir commented 4 months ago

Hi, The units of the entries should be consistent. Thanks