← Back to team overview

yade-dev team mailing list archive

Re: ViscoElastic test

 

Hi,

One last thing : in the 3rd case (young modulus
 case
), it would be useful to determine
 cn from E1,E2,massR and en. I know the formula for normal cn, but not for
cs. Do you have an idea ?

Thanks,

François.


2014-04-03 10:15 GMT+02:00 François <francois.kneib@xxxxxxxxx>:

> 2014-04-03 10:00 GMT+02:00 Anton Gladky <gladky.anton@xxxxxxxxx>:
>
> 2014-04-03 9:51 GMT+02:00 François <francois.kneib@xxxxxxxxx>:
>> > It will for sure break my scripts, as I used massMultiply=0. But your
>> > implementation is part of my "todo list", which I can remove now :-)
>> > You have made huge changes at once, so I need a little bit of time to
>> check
>> > everything. What I can see now is that you use "isnormal(r)" function
>> into
>> > the Ip2, which returns false if r==0. It will be a problem, for example
>> if
>> > we want to set cn, but cs=0 (same matter for en and et).
>>
>> I just wanted to check, that it is not(isnan). Thank for notice,  will
>> change to
>> not(isnan).
>>
> Yes that's what I supposed.
>
>
>>
>> I think you can  just remove massmultiply from your scripts and it should
>> work
>> fine.
>>
> I will re-check everything and see whether I am consistent with the
> current implementation. You can delete the "massmultiply" argument in the
> "contactParameterCalculation" function, as it's not used any more.
>
>> Anton
>>
>
> François
>
>

Follow ups

References