Welcome to the forums. Please post in English or French.

You are not logged in.

#1 2019-06-04 11:32:53

bhattarai
Member
Registered: 2014-01-31
Posts: 136

Calculation of right Cauchy-Green strain [Solved]

Dear Developers and experts,

I needed the right Cauchy-Green strain tensor for the introduction of a new material model.
I found nmgeom.F90 calculates the deformation gradient. But when I make an output (AS CALL  FLUSH) of
the deformation gradient [F] matrix and cauchy green strain tensor as C = transpose (F)*F, they are always an identity matrix at each time step with proper iteration. I was wondering why it is not updated with increased load and obvious deformation at the end of the simulation. Can anyone hint on this?

Also, I created a separate file to compute the matrix with following procedure:

! updated strain {E} in vector form
    eps(1)=deps(1)+epsm(1)
    eps(2)=deps(2)+epsm(2)
    eps(3)=deps(3)+epsm(3)
    eps(4)=deps(4)+epsm(4)
    eps(5)=deps(5)+epsm(5)
    eps(6)=deps(6)+epsm(6)

!modifying shear terms of the strain {E}
    do i = 1, 3
        epstot(i)=eps(i)
        epstot(3+i)=eps(3+i)/sqrt(2.d0)
    end do

! symmetric right cauchy green tensor {C} = 2{E} + {I} in vector form
        c11 = 2.d0*epstot(1)+1.d0
        c22 = 2.d0*epstot(2)+1.d0
        c33 = 2.d0*epstot(3)+1.d0
        c12 = 2.d0*epstot(4)
        c13 = 2.d0*epstot(5)
        c23 = 2.d0*epstot(6)

Here the components of the matrix is different and not unit matrix. can anyone provide
some light how to verify if the calculation of the {C} is correct or not?

Many thanks.
Bhattarai

Last edited by bhattarai (2019-06-25 11:04:04)

Offline

#2 2019-06-04 14:01:37

bhattarai
Member
Registered: 2014-01-31
Posts: 136

Re: Calculation of right Cauchy-Green strain [Solved]

Part of the time increment is attached here.


Attachments:
Screenshot from 2019-06-04 14-58-37.png.tar.gz, Size: 96.04 KiB, Downloads: 39

Offline

#3 2019-06-11 10:02:44

bhattarai
Member
Registered: 2014-01-31
Posts: 136

Re: Calculation of right Cauchy-Green strain [Solved]

Hello again,

I am wondering if my question does not fit to the Code_Aster development part of the forum. Can anyone put light on my question please?

Thank you again in advance.

Bests,
Bhattarai

Offline

#4 2019-06-12 19:39:48

AsterO'dactyle
Administrator
Registered: 2007-11-29
Posts: 252

Re: Calculation of right Cauchy-Green strain [Solved]

Hello,

If you want to compute large strains (gradient F or Cauchy-Green tensor) you have to set "large" flag to .true. (parameter of nmgeom)

In code_aster, when you set DEFORMATION='GROT_GDEP' or 'SIMO_MIEHE', it's the case.

Where are you in the code ?


Code_Asterの開発者

Offline

#5 2019-06-24 14:34:30

bhattarai
Member
Registered: 2014-01-31
Posts: 136

Re: Calculation of right Cauchy-Green strain [Solved]

Hi AsterO'dactyle,

My apologies for late reply. I set the parameter grand to .true. in geom.F90 and in my implementation file as,

       call nmgeom(ndim, nno, axi, .true., geom,&
                  kpg, ipoids, ivf, idfde, depl,&
                  ldfdi, poids, dfdi, f, eps,&
                  r)

but got an error:

Error: Type mismatch in argument ‘grand’ at (1); passed LOGICAL(8) to LOGICAL(1)

Bhattarai

Last edited by bhattarai (2019-06-25 11:03:49)

Offline

#6 2019-06-24 20:37:51

AsterO'dactyle
Administrator
Registered: 2007-11-29
Posts: 252

Re: Calculation of right Cauchy-Green strain [Solved]

Hello,

aster_logical is boolean(kind=1)

So, set .true._1 instead of .true.


Code_Asterの開発者

Offline

#7 2019-06-25 11:03:31

bhattarai
Member
Registered: 2014-01-31
Posts: 136

Re: Calculation of right Cauchy-Green strain [Solved]

Thank you AsterO'dactyle. The error is solved now, but got another problem which I shall post creating a new thread if I am not successful to solve it by myself.

Bhattarai

Offline

Board footer