FEAP User Forum

FEAP => Input File Issues => Topic started by: Yaakov on September 24, 2016, 03:48:04 PM

Title: Tie commnad
Post by: Yaakov on September 24, 2016, 03:48:04 PM

Dear FEAP Team & Users,

I want to define a hybrid material structure in FEAP (e.g. ITEST), and I use  tie command, tie mate 1 2 (e.g.). It works, but If I try command tie mate 2 1. the result would be wrong.

which influence has the "direction" of materials in TIE command, e.g. 1->2, 2->1, 1->3, 3->1 ... ?

and how can I find global number of nodes in FEAP, in my example, I use 4 elements, 4x4 nodes, global number of node should be 9.

Best Regards,
Yaakov

Title: Re: Tie commnad
Post by: Prof. R.L. Taylor on September 25, 2016, 10:11:56 AM

You are right!  It is wrong, we will work to find what is happening.  It would seem one could check the sequence and just use min/max to always get 1 -> 2 ; but first we need to make sure that will fix all cases.

Title: Re: Tie commnad
Post by: Yaakov on September 25, 2016, 11:42:48 AM

Dear Prof. Taylor,

Thanks very much for your answer! such problem I could find in other TIE command, e.g. TIE REGI 1 2 (2 1)...

Is it possible, that  I could find the global number of node?

Title: Re: Tie commnad
Post by: Prof. R.L. Taylor on September 25, 2016, 02:36:57 PM

I made the following changes in pcontr for now
        elseif(pcomp(titl(16:19),'regi',4)) then
          l1 = max(    0,int(td(1)))
          l2 = min(mxreg,int(td(2)))
          l3 = min(l1,l2)
          l4 = max(l1,l2)
and for 'mate'
        elseif(pcomp(titl(16:19),'mate',4)) then
          l1 = max(     1,int(td(1)))
          l2 = min(nummat,max(1,int(td(2))))
          l3 = min(l1,l2)
          l4 = max(l1,l2)
          l1 = 1

I will look at it more later to see if there is another way to ensure things work correctly.

Title: Re: Tie commnad
Post by: Yaakov on September 26, 2016, 02:27:53 AM

Dear Prof. Taylor,
Thanks very much for your detail answers!