Not sure

Author: ZibraMax

Examples/safarrppp.py Examples/safarrancho7-2.pyExamples/safarrppp.py renamed to Examples/safarrancho7-2.py

@@ -1,6 +1,6 @@
 if __name__ == '__main__':
     from FEM import NewtonTotalLagrangian, MGDCM, ContinumTotalLagrangian, QuadMembraneLinear, Geometry3D
-    # from FEM.Elements.ContinumElements import QuadShellLinearReddy as QuadShellLinear
+    from FEM.Elements.ContinumElements import QuadShellLinearReddy as QuadShellLinear
     from FEM.Elements.ContinumElements import QuadShellLinear
     import matplotlib.pyplot as plt
     import numpy as np

Examples/safarrancho7-3.py

@@ -98,3 +98,14 @@ def cm(_E):
             uij_num = JINV @ uij_num.T
             uij_num = uij_num.T
             a = 0
+    for w in W:
+        JS, _JS = e.get_jacobians(w)
+        FS = e.calculate_deformation_gradients(w)
+        thetas = e.derivatives_transformation()
+        for i, z in enumerate(e.Z):
+            JINV = _JS[i]
+            dpx = e.calculate_dpxs(w, i, JINV)
+            BNL = e.calculate_BNL(dpx)
+            UIJ = FS[i] - np.eye(3)
+            theta = thetas[i]
+            a = 0

src/FEM/Elements/ContinumElements/Shells_reddy.py

@@ -105,15 +105,15 @@ def calculate_A(self, F):
         dudx = F+np.eye(3)
         A = np.zeros((6, 9))
         A[0, 0] = dudx[0, 0]
-        A[0, 3] = dudx[0, 1]
-        A[0, 6] = dudx[0, 2]
+        A[0, 3] = dudx[1, 0]
+        A[0, 6] = dudx[2, 0]
 
-        A[1, 1] = dudx[1, 0]
+        A[1, 1] = dudx[0, 1]
         A[1, 4] = dudx[1, 1]
-        A[1, 7] = dudx[1, 2]
+        A[1, 7] = dudx[2, 1]
 
-        A[2, 2] = dudx[2, 0]
-        A[2, 5] = dudx[2, 1]
+        A[2, 2] = dudx[0, 2]
+        A[2, 5] = dudx[1, 2]
         A[2, 8] = dudx[2, 2]
 
         fila = np.array([dudx[0, 1], dudx[0, 0], 0, dudx[1, 1],

src/FEM/Elements/ContinumElements/Shells_rem.py

@@ -311,6 +311,22 @@ def calculate_deformation_gradients(self, w):
             FS.append(F)
         return FS
 
+    def derivatives_transformation(self):
+        thetas = []
+        e1s, e2s, e3s = self.calculate_e1_e2(deformed=True)
+        for k in range(len(self._p)):
+            phik = self._p[k]
+            theta = np.zeros((3, 3))
+            for i in range(len(self.coords)):
+                e1 = e1s[i]
+                e2 = e2s[i]
+                e3 = e3s[i]
+                theta[:, 0] += e1*phik[i]
+                theta[:, 1] += e2*phik[i]
+                theta[:, 2] += e3*phik[i]
+            thetas.append(theta)
+        return thetas
+
     def elementMatrices(self) -> None:
         """Calculate element matrices and vectors.
 
@@ -325,6 +341,7 @@ def elementMatrices(self) -> None:
         t_z, t_w = np.polynomial.legendre.leggauss(self.n_gauss_thickness)
         Ke = 0.0
         Fe = 0.0
+        THETAS = self.derivatives_transformation()
         for z_t, w_t in zip(t_z, t_w):
             JS, _JS = self.get_jacobians(z_t)
             FS = self.calculate_deformation_gradients(z_t)