e_damp.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	e_damp (co, nk, konl, lakonl, p1, p2, omx, bodyfx, nbody, s, sm, ff, nelem, elcon, nelcon, rhcon, nrhcon, alcon, nalcon, alzero, ielmat, ielorien, norien, orab, ntmat_, t0, t1, ithermal, vold, iperturb, nelemload, sideload, xload, nload, idist, sti, stx, iexpl, plicon, nplicon, plkcon, nplkcon, xstiff, npmat_, dtime, matname, mi, ncmat_, ttime, time, istep, iinc, nmethod)

Function/Subroutine Documentation

e_damp()

subroutine e_damp	(	real*8, dimension(3,*)	co,
		integer	nk,
		integer, dimension(20)	konl,
		character*8	lakonl,
		real*8, dimension(3)	p1,
		real*8, dimension(3)	p2,
		real*8	omx,
		real*8, dimension(3)	bodyfx,
		integer	nbody,
		real*8, dimension(60,60)	s,
		real*8, dimension(60,60)	sm,
		real*8, dimension(60)	ff,
		integer	nelem,
		real*8, dimension(0:ncmat_,ntmat_,*)	elcon,
		integer, dimension(2,*)	nelcon,
		real*8, dimension(0:1,ntmat_,*)	rhcon,
		integer, dimension(*)	nrhcon,
		real*8, dimension(0:6,ntmat_,*)	alcon,
		integer, dimension(2,*)	nalcon,
		real*8, dimension(*)	alzero,
		integer, dimension(mi(3),*)	ielmat,
		integer, dimension(mi(3),*)	ielorien,
		integer	norien,
		real*8, dimension(7,*)	orab,
		integer	ntmat_,
		real*8, dimension(*)	t0,
		real*8, dimension(*)	t1,
		integer	ithermal,
		real*8, dimension(0:mi(2),*)	vold,
		integer	iperturb,
		integer, dimension(2,*)	nelemload,
		character*20, dimension(*)	sideload,
		real*8, dimension(2,*)	xload,
		integer	nload,
		integer	idist,
		real*8, dimension(6,mi(1),*)	sti,
		real*8, dimension(6,mi(1),*)	stx,
		integer	iexpl,
		real*8, dimension(0:2*npmat_,ntmat_,*)	plicon,
		integer, dimension(0:ntmat_,*)	nplicon,
		real*8, dimension(0:2*npmat_,ntmat_,*)	plkcon,
		integer, dimension(0:ntmat_,*)	nplkcon,
		real*8, dimension(27,mi(1),*)	xstiff,
		integer	npmat_,
		real*8	dtime,
		character*80, dimension(*)	matname,
		integer, dimension(*)	mi,
		integer	ncmat_,
		real*8	ttime,
		real*8	time,
		integer	istep,
		integer	iinc,
		integer	nmethod
	)

!
!     computation of the element matrix and rhs for the element with
!     the topology in konl
!
!     ff: rhs without temperature and eigenstress contribution
!
      implicit none
!
      character*8 lakonl
      character*20 sideload(*)
      character*80 matname(*)
!
      integer konl(20),nelemload(2,*),nk,nbody,nelem,
     &  ithermal,iperturb,nload,idist,i,j,i1,i2,mi(*),
     &  nelcon(2,*),nrhcon(*),nalcon(2,*),ielmat(mi(3),*),
     &  ielorien(mi(3),*),ntmat_,nope,norien,iexpl,kode,imat,ncmat_,
     &  istep,iinc,nmethod,nplicon(0:ntmat_,*),nplkcon(0:ntmat_,*),
     &  npmat_
!
      real*8 co(3,*),xl(3,20),
     &  s(60,60),p1(3),p2(3),bodyfx(3),ff(60),
     &  elcon(0:ncmat_,ntmat_,*),
     &  rhcon(0:1,ntmat_,*),alcon(0:6,ntmat_,*),alzero(*),orab(7,*),
     &  t0(*),t1(*),voldl(3,20),vold(0:mi(2),*),xload(2,*),omx,
     &  sm(60,60),
     &  sti(6,mi(1),*),stx(6,mi(1),*),t0l,t1l,elas(21),elconloc(21)
!
      real*8 plicon(0:2*npmat_,ntmat_,*),plkcon(0:2*npmat_,ntmat_,*),
     &  xstiff(27,mi(1),*),dtime,ttime,time
!
      imat=ielmat(1,nelem)
!
      read(lakonl(8:8),'(i1)') nope
      nope=nope+1
!
!     computation of the coordinates of the local nodes
!
      do i=1,nope
        do j=1,3
          xl(j,i)=co(j,konl(i))
        enddo
      enddo
!
!     displacements for 2nd order static and modal theory
!
      if(iperturb.ne.0) then
         do i1=1,nope
            do i2=1,3
               voldl(i2,i1)=vold(i2,konl(i1))
            enddo
         enddo
      endif
!
!     initialisation of s
!
      do i=1,3*nope
        do j=1,3*nope
          s(i,j)=0.d0
        enddo
      enddo
!
!     calculating the stiffness matrix for the contact spring elements
!
         if(lakonl(7:7).eq.'A') then
            kode=nelcon(1,imat)
            t0l=0.d0
            t1l=0.d0
            if(ithermal.eq.1) then
               t0l=(t0(konl(1))+t0(konl(2)))/2.d0
               t1l=(t1(konl(1))+t1(konl(2)))/2.d0
            elseif(ithermal.ge.2) then
               t0l=(t0(konl(1))+t0(konl(2)))/2.d0
               t1l=(vold(0,konl(1))+vold(0,konl(2)))/2.d0
            endif
         endif
         call dashdamp(xl,elas,konl,voldl,s,imat,elcon,nelcon,
     &      ncmat_,ntmat_,nope,lakonl,t0l,t1l,kode,elconloc,plicon,
     &      nplicon,npmat_,iperturb,time,nmethod)
         return
      return