Elmer Discussion Forum

I get an output now, but unfortunately not the right values: mustatic: 0.29999999999999999 mudynamic: 0.20000000000000001 NodeLoad: 0.0000000000000000 DetermineContactSet: Creating fields out of normal and stick contact sets Fstatic: 4.5454039417394682E-322 Fdynamic: 6.9529800193766885E-310 the fric...

So little update: I deleted the line "Stick Contact Global = Logical True" and now the function is called from the beginning on. But my output is still not there. I even added the line CALL Info('TangentContactSet', 'Hallo_outofif') outside of the IF-loop to see if the If-condition is mayb...

That's what I thought. So right now I am trying to output Fdynamic and Fstatic directly to the terminal. It is not working yet though. Are you well versed in Fortran90? The code can be found in SolverUtils.F90, Subroutine TangentContactSet. I added lines to IF( FrictionContact .AND. & ListGetLog...

I changed "Normal-Tangential Velocity = Logical True" to "Normal-Tangential Displacement = Logical True" in both the Master and the Slave BC and deleted "Normal-Tangential Velocity Rotate = Logical False" and now the cube does not fall through the floor anymore and cont...

Thanks again for your effort and sorry for the late answer, I had to call in sick. So I tried the same thing and got the same output. I am really confused by this, because once contact is recognized, usually the non-penetration-condition of contact (contact gap >= 0) should hold, shouldn't it? Howev...

Thanks!

Thank you so much for your effort!

I changed everything according to your explanations, but now somehow no more contact is recognized and the contact surface does not hold the cube up, the cube penetrates the floor and just falls through it. I attached two pictures.

Thank you! That makes sense. But wouldn't it be Stress Bodyforce 2 = $ -9.81 * 7870.0 * "Volume" (Normal Force = 9,81 * mass instead of density) I don't know the measurements of the beam, but in my case (the cube is 1cm³) I would put Stress Bodyforce 3 = $ -9.81 * 7810.0*10e-6. Do you agre...

Do you mean the normal force of the cube? I didn't define it, because I assumed ELMER would calculate it automatically. In which section would you define it?

Hello Kevin,

thanks for your reply! I've seen this test case before. Unfortunately It didn't really help me since they do not use "Apply Contact BCs" and "Friction Contact".

Sincerely,
Maike