FEAP User Forum

FEAP => General questions => Topic started by: Dr. Ahmad Chihadeh on August 30, 2024, 02:59:39 AM

Title: Viscoelastic material
Post by: Dr. Ahmad Chihadeh on August 30, 2024, 02:59:39 AM

Dear Colleagues,

Following the derivations of the Linear Viscoelasticity in the Version 8.6 Theory Manual and looking at subroutine viscoe

the formula
qi(i,n) = exp_n*qi(i,n) + dq_n*(ee(i) - en(i))

shouldnt be
qi(i,n) = mu_n*exp_n*qi(i,n) + dq_n*(ee(i) - en(i))
?

as eqn. 9.37 in the manual
s(t) = 2G[mu_0 e(t) + mu_1 h(t)]

Best regards
Ahmad

Title: Re: Viscoelastic material
Post by: Prof. R.L. Taylor on August 30, 2024, 03:31:08 PM

I believe it is correct as coded. If you multiplied as you state the left and right would be inconsistent by the mu value.

Title: Re: Viscoelastic material
Post by: Dr. Ahmad Chihadeh on September 01, 2024, 05:10:32 AM

but dq_n is also multiplied by mu_n

do n = 1,nv
        mu_n  = d(2*n+49)
        dtau  = dt/d(2*n+50)

        dq_n = mu_n * hvisc(dtau,exp_n)
        gfac = gfac + dq_n
        mu_0 = mu_0 + mu_n

!       Update history and compute viscoelastic deviatoric stress
        do i = 1,ntm
          qi(i,n) = exp_n*qi(i,n) + dq_n*(ee(i) - en(i))
          sig(i)  = sig(i) + qi(i,n)
        end do ! i
      end do ! n

Title: Re: Viscoelastic material
Post by: Prof. S. Govindjee on September 01, 2024, 03:18:53 PM

I tried a simple creep test computation and did not get the result I was expecting.  I will have a further look at the coding to see if there is an error or if my expectation was flawed.  Attached is my input file, where I was expecting the displacement to double upon relaxation.

[PS: There is a printing error in elements/material/inmate.f, whereby the volumetric viscoelastic parameters are not being printed to the output file.]