← Back to team overview

yade-dev team mailing list archive

Re: yade-dev Digest, Vol 16, Issue 5

 

Get it, thank you so much , Mr. Bruno .

Regards.
Kan.


On 3/6/08, yade-dev-request@xxxxxxxxxxxxxxxx <
yade-dev-request@xxxxxxxxxxxxxxxx> wrote:
>
> Send yade-dev mailing list submissions to
>        yade-dev@xxxxxxxxxxxxxxxx
>
> To subscribe or unsubscribe via the World Wide Web, visit
>        https://lists.berlios.de/mailman/listinfo/yade-dev
> or, via email, send a message with subject or body 'help' to
>        yade-dev-request@xxxxxxxxxxxxxxxx
>
> You can reach the person managing the list at
>        yade-dev-owner@xxxxxxxxxxxxxxxx
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of yade-dev digest..."
>
>
> Today's Topics:
>
>   1. Could I ask the method to calculate the shear forces? (kan)
>   2. (no subject)
>   3. Re: Could I ask the method to calculate the shear forces?
>      (Bruno Chareyre)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 5 Mar 2008 10:30:28 -0600
> From: kan <nicessgg@xxxxxxxxx>
> Subject: [yade-dev] Could I ask the method to calculate the shear
>        forces?
> To: yade-dev@xxxxxxxxxxxxxxxx
> Message-ID:
>        <45213f870803050830g47e237hed5a66af7a104597@xxxxxxxxxxxxxx>
> Content-Type: text/plain; charset="utf-8"
>
> Hi, dear all,
>
> When I develop my code to calculate the shear force I meet some problem.
>
> ------------------------------
>
> Message: 2
> Message-ID: <mailman.60.1204801280.17254.yade-dev@xxxxxxxxxxxxxxxx>
>
> Deta_Shear Force = ks*deta_shear_displacement
> while the deta_shear_displacement is the deta displacement at the shear
> direction at each contact.
> my problem now is that, how to calculate this displacement ?
>
> I know for normal displacement at each contact, at each time step, when we
> get the resultant displacement during this time step,
> then we get the new position for each particle, so we can get the distance
> between the center of the particles, which is the normal distance. so we
> can
> get the normal displacement and use the
> NormalForce=kn*normal_displacement.
>
> but how can I calculate the deta shear displacement? it seems that I need
> to
> store the old position of each particles, which means a doubled storage of
> the memory and lower speed.
>
> could you please give me a hint, and which file is it in to calculate the
> shear force or the resultant forces?
>
> Thanks for your nice information in advance.
>
> Kan
>
> ------=_Part_27553_5946959.1204734628668
> Content-Type: text/html; charset=UTF-8
> Content-Transfer-Encoding: 7bit
> Content-Disposition: inline
>
> <div>Hi, dear all,</div>
> <div>&nbsp;</div>
> <div>When I develop my code to calculate the shear force I meet some
> problem.</div>
> <div>From the literature, we know that</div>
> <div>Deta_Shear Force = ks*deta_shear_displacement</div>
> <div>while the deta_shear_displacement is the deta displacement at the
> shear direction&nbsp;at each contact. </div>
> <div>my&nbsp;problem now is that,&nbsp;how to calculate this displacement
> ?</div>
> <div>&nbsp;</div>
> <div>I know for normal displacement at each contact, at each time step,
> when we get the resultant displacement during this time step, </div>
> <div>then we get the new position for each particle, so we can get the
> distance between the center of the particles, which is the normal distance.
> so we can get the normal displacement and use the
> NormalForce=kn*normal_displacement.</div>
>
> <div>&nbsp;</div>
> <div>but how can I calculate the deta shear displacement? it seems that I
> need to store the old position of each particles, which means a doubled
> storage of the memory and lower speed.</div>
> <div>&nbsp;</div>
> <div>could you please give me a hint, and which file is it in to calculate
> the shear force or the resultant forces?</div>
> <div>&nbsp;</div>
> <div>Thanks for your nice information in advance.</div>
> <div>&nbsp;</div>
> <div>Kan</div>
> <div>&nbsp;</div>
>
> ------=_Part_27553_5946959.1204734628668--
>
>
> ------------------------------
>
> Message: 3
> Date: Wed, 05 Mar 2008 19:19:13 +0100
> From: Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>
> Subject: Re: [yade-dev] Could I ask the method to calculate the shear
>        forces?
> To: Yade Development Group <yade-dev@xxxxxxxxxxxxxxxx>
> Message-ID: <47CEE421.8090701@xxxxxxxxxxx>
> Content-Type: text/plain; charset=UTF-8; format=flowed
>
>
>
>
> kan a ?crit :
> > Hi, dear all,
> >
> > When I develop my code to calculate the shear force I meet some problem.
> > From the literature, we know that
> > Deta_Shear Force = ks*deta_shear_displacement
> > while the deta_shear_displacement is the deta displacement at the
> > shear direction at each contact.
> > my problem now is that, how to calculate this displacement ?
> >
> > I know for normal displacement at each contact, at each time step,
> > when we get the resultant displacement during this time step, then we
> > get the new position for each particle, so we can get the distance
> > between the center of the particles, which is the normal distance. so
> > we can get the normal displacement and use the
> > NormalForce=kn*normal_displacement.
> Its a matter of vocabulary.
> I think you make a confusion between "displacement" and "position". In
> the normal direction, you get at each time-step the "POSITION" (i.e.
> interpenetration), not the displacement.
> So that NormalForce=kn*interpenetration.
>
> >
> > but how can I calculate the deta shear displacement? it seems that I
> > need to store the old position of each particles, which means a
> > doubled storage of the memory and lower speed.
> shear_displacement = (velocity1 - velocity2 + r1*rotational_velocity1 -
> r2*rotational_velocity2)*dt
>
> This IS a displacement.
>
> So that shear_force = shear_force + shear_displacement*ks
>
> No need to store initial orientations.
>
> Check ElasticContactLaw.
>
> Bruno
>
> >
> > ------------------------------------------------------------------------
> >
> > _______________________________________________
> > yade-dev mailing list
> > yade-dev@xxxxxxxxxxxxxxxx
> > https://lists.berlios.de/mailman/listinfo/yade-dev
> >
>
>
> --
>
> _______________
> Chareyre Bruno
> Maitre de conference
>
> Institut National Polytechnique de Grenoble
> Laboratoire 3S (Soils Solids Structures) - bureau E145
> BP 53 - 38041, Grenoble cedex 9 - France
> T?l : 33 4 56 52 86 21
> Fax : 33 4 76 82 70 43
> ________________
>
>
>
> ------------------------------
>
> _______________________________________________
> yade-dev mailing list
> yade-dev@xxxxxxxxxxxxxxxx
> https://lists.berlios.de/mailman/listinfo/yade-dev
>
>
> End of yade-dev Digest, Vol 16, Issue 5
> ***************************************
>
_______________________________________________
yade-dev mailing list
yade-dev@xxxxxxxxxxxxxxxx
https://lists.berlios.de/mailman/listinfo/yade-dev