Hi,
I am using FEAP8.3.19 and not sure if you correct this in FEAP8.4. Around the line 2114 in inmate.f, you coded:
c Output Regular Mooney-Rivlin
elseif(imat.eq.9) then
if(incomp) then
bulk = 0.0d0
d(21) = 0.0d0
else
bulk = e1/(1.d0 - nu12*2.d0)/3.d0
d(21) = bulk - 2.d0/3.d0*g12
endif
g12 = e1/(1.d0 + nu12)/2.d0
......
in which you calculated g12 after calculating d(21). This leads to d(21) = bulk because g12 is still 0 when implementing d(21) = bulk -2/3*g12. Apparently this is wrong because the first Lame parameter d(21) and the bulk modulus should certainly be different.
Best,
Ziyu
