FEAP User Forum

FEAP => General questions => Topic started by: Yaakov on May 05, 2017, 02:31:01 AM

Title: About use of command "ORTH VECT"
Post by: Yaakov on May 05, 2017, 02:31:01 AM

Dear All,

I would ask, how to use the command "ORTH VECT"?

until now I can't find the way to get the results with this command "ORTH VECT"...

Thanks very much for your answer!

Kind Regards,
Yaakov

Title: Re: About use of command "ORTH VECT"
Post by: FEAP_Admin on May 05, 2017, 08:20:49 AM

The ordering of the words is the other way around.
From elements/material/inmate.f

Code: [Select]

c       VECTor ORTHotropic

        elseif(ietype.ne.6 .and. pcomp(text(1),'vect',4)) then

c         Orthotropic orientation VECTors: V_1 & V_2

          if(pcomp(text(2),'orth',4)) then
            d(242) = 1.0d0
            do i = 1,6
              d(i+242) = ev(i)
            end do ! i
         endif

The command looks for 6 numbers.  The first three are V1 and the second 3 are V2, stored in d(243:248).  d(242) is the flag indicating that they have been set.  Note that FEAP will orthonormalize what ever you enter.  So what gets stored is V1/|V1| and [V2 - V2.V1/|V1|]/|V2 - V2.V1/|V1||.  The usage is

Code: [Select]

vect orth v1x v1y v1z v2x v2y v2z


Title: Re: About use of command "ORTH VECT"
Post by: Yaakov on May 05, 2017, 03:51:47 PM

Dear FEAP Administrator,

Thanks very much for your answer!

I give the following values in input file, it seems not working... more details in attachments IEX3A, which I change FEAP own examples for orthotropic cases.

Code: [Select]

  ELAStic ORTH 210000 170000 290000 0.3 0.2 0.27 16500 26500 36500
    VECT ORTH 0.88698698639195 0.46179441959743  0 -0.46179441959743  0.88698698639195 0
And I check feap resource code from elements/material/small/ estrsd.f, I couldn't good understand about the theory part of orthonormalization process in FEAP. e.g.     call pushr2,    call pushr4 ...

Code: [Select]

      if(nint(d(242)).eq.1) then
        do i = 1,3
          ro(i,1) = d(242+i)
          ro(i,2) = d(245+i)
        end do ! i
        call triad(ro)
        psi = 0.0d0
        call dmat2d(d,psi,dm,betm)      ! Transform thermal stress
        call pushr2(ro,betm,beta,1.0d0)
        call tranr4(ro,ro,tt,.false.)
        call pushr4(tt,tt,dd,dm,1.d0)   ! Transform moduli
As far as I know, VECT ORTH is a new command in FEAP 8.4 Version.

I would ask, is it possible that you can add a benchmark Input file for VECT ORTH problems in our FEAP forum?

Many thanks for your support!

Kind regards,
Yaakov