Differenze

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

--- wikipaom2017:051:010:000 [2017/03/31 13:48] – ebertocchi
+++ wikipaom2017:051:010:000 [2017/04/19 20:03] (versione attuale) – [valutazione delle funzioni di forma] ebertocchi
@@ Linea 1: / Linea 1: @@
+{{ :wikipaom2017:050:010:derivate_funzionidiforma_iso4.wxmx |}}
+====== Isoparametrico 4 nodi ======
+===== Funzioni di forma =====
+==== valutazione delle funzioni di forma ====
+<code>
+c234567---------------------------------------------------------------72
+c================== calcolo funzioni di forma per l'isoparam. 4 nodi ===
+      subroutine defN(xieta, N)
+      implicit none
+      include 'dim.for'
+      double precision xieta(ncrd), N(nnode)
+      N(1) = (1d0-xieta(1))*(1d0-xieta(2))/4d0
+      N(2) = (1d0+xieta(1))*(1d0-xieta(2))/4d0
+      N(3) = (1d0+xieta(1))*(1d0+xieta(2))/4d0
+      N(4) = (1d0-xieta(1))*(1d0+xieta(2))/4d0
+      return
+      end subroutine
+</code>
+==== Derivate direzionali delle stesse nelle coordinate naturali ====
+<code>
+c234567---------------------------------------------------------------72
+c========= calcolo derivate funzioni di forma per l'isoparam. 4 nodi ===
+      subroutine defdN(xieta,dN)
+      implicit none
+      include 'dim.for'
+      double precision xieta(ncrd), dN(ncrd,nnode)
+      dN(1,1) = -(      1d0 -xieta(2))/4.d+0
+      dN(1,2) =  (      1d0 -xieta(2))/4.d+0
+      dN(1,3) =  ( xieta(2) +     1d0)/4.d+0
+      dN(1,4) = -( xieta(2) +     1d0)/4.d+0
+      dN(2,1) = -(      1d0 -xieta(1))/4.d+0
+      dN(2,2) = -( xieta(1) +1d0     )/4.d+0
+      dN(2,3) =  ( xieta(1) +1d0     )/4.d+0
+      dN(2,4) =  (      1d0 -xieta(1))/4.d+0
+      return
+      end
+</code>