yade-users team mailing list archive

[Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>]

> I checked the new functor with simulation inspection and detected the
> first problem.
>
> I set CapillaryPressure=10000 (see line 56 in the attachment), but in
> simulation inspection CapillaryPressure=0. The same with betan and
> betas (lines 60-62).
>
> So, if you have the time, please have a look at the new functor code
> and the script. Furthermore I will need more help, when we want to
> create a new law functor.
>

In Ip2_FrictMat_FrictMat_MindlinCapillaryPhys::go, you only duplicated
the FrictPhys code apparently, which means capillary pressure and
Mindlin-specific data is never set.
Have a look at what Ip2_FrictMat_FrictMat_MindlinPhys is doing, you have
to do the same + also set capillary pressure.

You will not have to write any law functor. The existing Mindlin laws
and capillary law will do their jobs correctly, as soon as the
interaction physics contains the correct data.

Bruno

p.s. private mail -> understood. There is maybe a "attach" button in
launchpad, not sure...





> Do I have to be in the yade-dev group to submit code?
>
> Christian
>
> P.S. I write private mail to you, because I can not attach files, when
> I post a question.
>
>
>
> Zitat von Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>:
>
>> (Please avoid private mails and keep posting to yade-users Christian)
>>
>> If everything looks correct, then you have to validate the results and
>> submit your code as candidate for commit. :)
>> It is worth adding combined Mindlin and capillary law to Yade.
>>
>> Bruno
>>
>> On 02/11/11 12:42, Christian Jakob wrote:
>>> Thank you very very much. I created following new files (see
>>> attachment):
>>>
>>> MindlinCapillaryPhys.hpp
>>> Ip2_FrictMat_FrictMat_MindlinCapillaryPhys.hpp
>>> and the corresponding .cpp's
>>>
>>> scons was able to compile yade with the new files and my test script
>>> is running.
>>>
>>> Did I forget something? Is there anything else to do?
>>>
>>> Christian.
>>>
>>> Zitat von Bruno Chareyre <bruno.chareyre@xxxxxxxxxxx>:
>>>
>>>> Hi,
>>>> You need to create the type of interaction that will contain all the
>>>> data for Mindlin + capillary.
>>>> In the case of FrictPhys, we added the capillary data by using
>>>> inheritance:
>>>> class CapillaryPhys : public FrictPhys (CapillaryPhys.hpp)
>>>>
>>>> OTOH, MindlinPhys also inherits from FrictPhys:
>>>> class MindlinPhys: public FrictPhys (HertzMindlin.hpp)
>>>>
>>>> The best thing to do is to derive the new type from MindlinPhys:
>>>> class MindlinCapillaryPhys: public MindlinPhys
>>>>
>>>> Then the functor you just created will have to instantiate this new
>>>> type
>>>> of interaction.
>>>> Make sure that Ip2_FrictMat_FrictMat_MindlinCapillaryPhys
>>>> reproduces the
>>>> behaviour of Ip2_FrictMat_FrictMat_MindlinPhys for the solid part, as
>>>> found in HertzMindlin.cpp (lines 40 to 96). For that matter, you will
>>>> probably be inclined to duplicate the code, but we can avoid that
>>>> in the
>>>> end when everything else is solved.
>>>>
>>>> Bruno
>>>>
>>>> On 01/11/11 10:08, Christian Jakob wrote:
>>>>> Dear Mr. Chareyre,
>>>>>
>>>>> Since two hours I am looking at the source code of
>>>>> Ip2_FrictMat_FrictMat_CapillaryPhys and I really do not know what to
>>>>> do ...
>>>>>
>>>>> First I created a Ip2_FrictMat_FrictMat_MindlinCapillaryPhys.hpp
>>>>> and a
>>>>> Ip2_FrictMat_FrictMat_MindlinCapillaryPhys.cpp by simply copy and
>>>>> rename the Ip2_FrictMat_FrictMat_CapillaryPhys code. Then I replaced
>>>>> the cross references, so that we have a
>>>>> Ip2_FrictMat_FrictMat_MindlinCapillaryPhys, that is doing the same as
>>>>> Ip2_FrictMat_FrictMat_CapillaryPhys.
>>>>>
>>>>> I do not understand the structure of YADE. Can you give me a hint,
>>>>> what there is to do in the next step?
>>>>>
>>>>> Christian.
>>>>>
>>>>
>>>>
>>>> -- 
>>>> _______________
>>>> Bruno Chareyre
>>>> Associate Professor
>>>> ENSE³ - Grenoble INP
>>>> 11, rue des Mathématiques
>>>> BP 46
>>>> 38402 St Martin d'Hères, France
>>>> Tél : +33 4 56 52 86 21
>>>> Fax : +33 4 76 82 70 43
>>>> ________________
>>>>
>>>>
>>>> _______________________________________________
>>>> 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
>>>>
>>>
>>
>>
>> -- 
>> _______________
>> Bruno Chareyre
>> Associate Professor
>> ENSE³ - Grenoble INP
>> 11, rue des Mathématiques
>> BP 46
>> 38402 St Martin d'Hères, France
>> Tél : +33 4 56 52 86 21
>> Fax : +33 4 76 82 70 43
>> ________________
>>
>>
>> _______________________________________________
>> 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
>>
>


-- 
_______________
Bruno Chareyre
Associate Professor
ENSE³ - Grenoble INP
11, rue des Mathématiques
BP 46
38402 St Martin d'Hères, France
Tél : +33 4 56 52 86 21
Fax : +33 4 76 82 70 43
________________