Differenze

Queste sono le differenze tra la revisione selezionata e la versione attuale della pagina.

--- wikipaom2017:051.020.000 [2017/03/31 08:34] – creata ebertocchi
+++ wikipaom2017:051.020.000 [2017/03/31 11:43] (versione attuale) – ebertocchi
@@ Linea 1: / Linea 1: @@
+<code>
+c=============================== calcolo matrice Jacobiana trasposta ===
+      subroutine defJt(xieta,xord,Jt)
+      implicit none
+      include 'dim.for'
+      double precision xieta(ncrd),xord(ncrd,nnode),Jt(ncrd,ncrd),
+     +dN(ncrd,nnode)
+      integer inode,icrd,jcrd
+c     azzero la matrice Jt
+      call dzero(Jt,ncrd*ncrd)
+c     accumulo i contributi relativi alle varie funzioni di forma nodali
+      do inode=1,nnode
+       call defdN(xieta,dN)
+c      scorro sulle righe della matrice jacobiana, ossia sull'indice
+c      di coordinata fisica 1:'x', 2:'y', 3:'z'
+       do icrd=1,ncrd
+c       scorro sulle colonne della matrice jacobiana, ossia sull'indice
+c       di coordinata naturale 1:'xi', 2:'eta', 3:'zeta'
+        do jcrd=1,ncrd
+         Jt(icrd,jcrd) = Jt(icrd,jcrd) + xord(icrd,inode)*dN(jcrd,inode)
+        enddo
+       enddo
+      enddo
+c     fatto.
+      return
+      end subroutine
+</code>