normal contact force

[Liam]

Hello all,

I am doing 3-D contact Problem, and the friction cases are not fulfilled yet. That is why I need to code it by myself. I am wondering that which subroutine stores the contact normal force (Fn) in 3-D cases.

Therefore, I am very looking forward to all you guys' replying. Many many many thanks!!!

Best regards

Liam

[fheinrich]

Dear Liam,

were you able to implement the friction-case for 3-D contact problems?

In thread http://feap.berkeley.edu/forum/index.php?topic=675.msg1743#msg1743 it seems
as if friction was partially implemented for certain cases. However, I am not sure if this is true
or whether it was/is not implemented in FEAP Ver. 8.2/8.4 at all?

Best regards

[Prof. R.L. Taylor]

In contact elements the frictional case depends on the parameter 'iffric' which allows you to follow where friction is being implemented.  In the case of 3-d it is apparently in the array ch3 with the parameter p3(15).  Similarly, the normal force is in ch3 with parameter p3(14).

Other variables are similarly defined in cnts3d.f file using the arrays w1(*) and w3(*) for storage of names.