FEAP User Forum

FEAP => Programming => Topic started by: Liam on April 07, 2015, 07:40:25 AM

Title: normal contact force
Post by: Liam on April 07, 2015, 07:40:25 AM

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

Title: Re: normal contact force
Post by: fheinrich on September 15, 2015, 02:17:22 AM

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

Title: Re: normal contact force
Post by: Prof. R.L. Taylor on September 15, 2015, 06:09:29 AM

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.