I would like to model the growth phenomenon in the cardiac muscle according to the following equations below. the code I have written is the one where I calculate the growth tensor to form the elastic tensor.
But I have a quadrature error.
Has anyone ever modelled this phenomenon before? Please help.
where theta is the growth factor, lamda_e the elastic stretch and the lamda_h is the hydrostatic stretch.
here is the code function
def growthfunction (f_0, s_0, n_0, d,v):
t = np.linspace(0,1,4)
I = Identity(d)
for j in range(N):
F_gt=F_g.T FF_g=inner(f_0,(F_g.T*F_g)*f_0) lamda_g=sqrt(FF_g) #lamda_e=(1/lamda_g)*(f_0*(F.T/O_0)*(F/O_0)*f_0)^0.5 O_1=O_0+dt*((1/gam)*(((O_max-O_0)/(O_max-O_min))**lam)*(lamda_g-lamda_h)) F_g=O_1*I O_0=O_1 print ('\n growth steps ' ,j)