[cig-commits] r15956 - short/3D/PyLith/branches/pylith-friction/libsrc/feassemble

[brad at geodynamics.org]

Author: brad
Date: 2009-11-11 17:26:26 -0800 (Wed, 11 Nov 2009)
New Revision: 15956

Added:
   short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/jacobian3d_lgdeform.wxm
Log:
Updated Jacobian calculation for large deformations.

Added: short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/jacobian3d_lgdeform.wxm
===================================================================
--- short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/jacobian3d_lgdeform.wxm	                        (rev 0)
+++ short/3D/PyLith/branches/pylith-friction/libsrc/feassemble/jacobian3d_lgdeform.wxm	2009-11-12 01:26:26 UTC (rev 15956)
@@ -0,0 +1,76 @@
+/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
+/* [ Created with wxMaxima version 0.8.3a ] */
+
+/* [wxMaxima: input   start ] */
+Bnli: matrix([Nip,0,0],
+[Niq,0,0],
+[Nir,0,0],
+[0,Nip,0],
+[0,Niq,0],
+[0,Nir,0],
+[0,0,Nip],
+[0,0,Niq],
+[0,0,Nir]);
+/* [wxMaxima: input   end   ] */
+
+/* [wxMaxima: input   start ] */
+Bnlj: matrix([Njp,0,0],
+[Njq,0,0],
+[Njr,0,0],
+[0,Njp,0],
+[0,Njq,0],
+[0,Njr,0],
+[0,0,Njp],
+[0,0,Njq],
+[0,0,Njr]);
+/* [wxMaxima: input   end   ] */
+
+/* [wxMaxima: input   start ] */
+S: matrix([s11,s12,s13,0,0,0,0,0,0],
+[s12,s22,s23,0,0,0,0,0,0],
+[s13,s23,s33,0,0,0,0,0,0],
+[0,0,0,s11,s12,s13,0,0,0],
+[0,0,0,s12,s22,s23,0,0,0],
+[0,0,0,s13,s23,s33,0,0,0],
+[0,0,0,0,0,0,s11,s12,s13],
+[0,0,0,0,0,0,s12,s22,s23],
+[0,0,0,0,0,0,s13,s23,s33]);
+/* [wxMaxima: input   end   ] */
+
+/* [wxMaxima: input   start ] */
+x: transpose(Bnli) . S . Bnlj;
+/* [wxMaxima: input   end   ] */
+
+/* [wxMaxima: input   start ] */
+Bi: matrix([Nip*(1+l11), l21*Nip, l31*Nip],
+[l12*Niq,(1+l22)*Niq, l32*Niq],
+[l13*Nir, l23*Nir, (1+l33)*Nir],
+[(1+l11)*Niq+l12*Nip,(1+l22)*Nip+l21*Niq,l31*Niq+l32*Nip],
+[l12*Nir+l13*Niq, (1+l22)*Nir+l23*Niq, (1+l33)*Niq+l32*Nir],
+[(1+l11)*Nir+l13*Nip, l21*Nir+l23*Nip, (1+l33)*Nip+l31*Nir]);
+/* [wxMaxima: input   end   ] */
+
+/* [wxMaxima: input   start ] */
+Bj: matrix([Njp*(1+l11), l21*Njp, l31*Njp],
+[l12*Njq,(1+l22)*Njq, l32*Njq],
+[l13*Njr, l23*Njr, (1+l33)*Njr],
+[(1+l11)*Njq+l12*Njp,(1+l22)*Njp+l21*Njq,l31*Njq+l32*Njp],
+[l12*Njr+l13*Njq, (1+l22)*Njr+l23*Njq, (1+l33)*Njq+l32*Njr],
+[(1+l11)*Njr+l13*Njp, l21*Njr+l23*Njp, (1+l33)*Njp+l31*Njr]);
+/* [wxMaxima: input   end   ] */
+
+/* [wxMaxima: input   start ] */
+C: matrix([C1111, C1122, C1133, C1112, C1123, C1113],
+[C1122, C2222, C2233, C2212, C2223, C2213],
+[C1133, C2233, C3333, C3312, C3323, C3313],
+[C1112, C2212, C3312, C1212, C1223, C1213],
+[C1123, C2223, C3323, C1223, C2323, C2313],
+[C1113, C2213, C3313, C1213, C2313, C1313]);
+/* [wxMaxima: input   end   ] */
+
+/* [wxMaxima: input   start ] */
+y: transpose(Bi) . C . Bj;
+/* [wxMaxima: input   end   ] */
+
+/* Maxima can't load/batch files which end with a comment! */
+"Created with wxMaxima"$