← Back to team overview

yade-users team mailing list archive

Re: [Question #175394]: numerical instability and small time steps

 

Question #175394 on Yade changed:
https://answers.launchpad.net/yade/+question/175394

Chiara Modenese proposed the following answer:
On 29 November 2011 09:25, Christian Jakob <
question175394@xxxxxxxxxxxxxxxxxxxxx> wrote:

> Question #175394 on Yade changed:
> https://answers.launchpad.net/yade/+question/175394
>
> Christian Jakob posted a new comment:
> It is not double-checked, but quad-checked (maybe even more).
>
Did you check the computation of the time step in PFC in case of viscous
damping and Hertzian law? This has to be quad-checked. In Yade there is no
such a computation of time step that considers viscous damping at the
contact (AFAIK).

C.

>
> YADE input file:
>
> #GENERAL SETTINGS:
> initial_porosity                = 0.4
> shear_modulus                   = 3e7
> poisson_ratio                   = 0.3
> young_modulus                   = 2*shear_modulus*(1+poisson_ratio)
> friction_coeff                  = 2
> angle                                   = math.atan(friction_coeff)
> rho_p                                   = 2650  #density of particles
>
> #SET DAMPING CONSTANTS:
> local_damping                   = 0.01
> viscous_normal                  = 0.021
> viscous_shear                   = 0.8*viscous_normal
>
> PFC input file:
>
>        ;GENERAL SETTINGS:
>        initial_porosity                = 0.4
>        shear_modulus           = 3e7
>        poisson_ratio           = 0.3
>        friction_coeff          = 2
>        rho_p                           = 2650  ;density of particles
>
>        ;SET DAMPING CONSTANTS:
>        local_damping           = 0.01
>        viscous_normal          = 0.021
>        viscous_shear           = 0.8*viscous_normal
>
> The only thing I could imagine, that is wrong, is
>
> young_modulus                   = 2*shear_modulus*(1+poisson_ratio)
>
> but wikipedia tells me this is ok...
>
> --
> You received this question notification because you are a member of
> yade-users, which is an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.