Coding Ogden hyperelastic model - Finite Element Modeling
https://polymerfem.com/community/finite-element-modeling-aa/coding-ogden-hyperelastic-model/
PolymerFEM.com Discussion Boarden-USFri, 19 Apr 2024 06:37:04 +0000wpForo60
https://polymerfem.com/community/finite-element-modeling-aa/coding-ogden-hyperelastic-model/#post-25291
Mon, 21 Nov 2016 19:35:29 +0000Finite Element ModelingJorgenhttps://polymerfem.com/community/finite-element-modeling-aa/coding-ogden-hyperelastic-model/#post-25291Coding Ogden hyperelastic model
https://polymerfem.com/community/finite-element-modeling-aa/coding-ogden-hyperelastic-model/#post-29561
Wed, 05 Oct 2016 15:58:02 +0000 I am trying to implement Ogden hyperelastic model (N=1)(I want to modify it later) in Abaqus. The way in which I am coding isFrom degradNew and stretchNew get F and U. Hence get R = F * U^-1Compute b = J**(-2/3)*(F * F^T)Compute eigen values and eigen vector of bCompute lambda as sqrt of eigen values of bCompute principal Cauchy stress as s(k,k) = 2/J*mu/alpha*( lamda_k**alpha - 1/3(lamda_1**alpha + lamda_2**alpha + lamda_3**alpha) + 2/D(J-1)Rotate back the principal Cauchy stress, S = as eigenVec * S * eigenVec^TGet back to corotational frame as S = R^T * S * TMy results are not the same as what Abaqus gets. I am not sure what I am doing wrong. Any pointersS]]>Finite Element Modelinggamblepath41https://polymerfem.com/community/finite-element-modeling-aa/coding-ogden-hyperelastic-model/#post-29561