-
Notifications
You must be signed in to change notification settings - Fork 0
/
StrainsEpsilon.m
98 lines (96 loc) · 5.28 KB
/
StrainsEpsilon.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
function out1 = StrainsEpsilon(x1_1,x2_1,x3_1,x4_1,x1_2,x2_2,x3_2,x4_2,u1_1,u2_1,u3_1,u4_1,u1_2,u2_2,u3_2,u4_2)
%STRAINSEPSILON
% OUT1 = STRAINSEPSILON(X1_1,X2_1,X3_1,X4_1,X1_2,X2_2,X3_2,X4_2,U1_1,U2_1,U3_1,U4_1,U1_2,U2_2,U3_2,U4_2)
% This function was generated by the Symbolic Math Toolbox version 8.1.
% 25-Jun-2020 19:56:07
t2 = sqrt(3.0);
t3 = t2.*(1.0./1.2e1);
t4 = t3+1.0./4.0;
t5 = t2.*x1_2.*(1.0./3.0);
t6 = t2.*x3_2.*(1.0./3.0);
t7 = x1_1.*x2_2;
t8 = x1_2.*x4_1;
t9 = x2_1.*x3_2;
t10 = x3_1.*x4_2;
t11 = t2.*x1_1.*x2_2.*(1.0./3.0);
t12 = t2.*x1_2.*x4_1.*(1.0./3.0);
t13 = t2.*x2_2.*x3_1.*(1.0./3.0);
t14 = t2.*x2_1.*x4_2.*(2.0./3.0);
t15 = t2.*x3_2.*x4_1.*(1.0./3.0);
t19 = x1_2.*x2_1;
t20 = x1_1.*x4_2;
t21 = x2_2.*x3_1;
t22 = x3_2.*x4_1;
t23 = t2.*x1_2.*x2_1.*(1.0./3.0);
t24 = t2.*x1_1.*x4_2.*(1.0./3.0);
t25 = t2.*x2_1.*x3_2.*(1.0./3.0);
t26 = t2.*x2_2.*x4_1.*(2.0./3.0);
t27 = t2.*x3_1.*x4_2.*(1.0./3.0);
t16 = t7+t8+t9+t10+t11+t12+t13+t14+t15-t19-t20-t21-t22-t23-t24-t25-t26-t27;
t17 = 1.0./t16;
t28 = t2.*x2_2.*(1.0./3.0);
t29 = t2.*x4_2.*(1.0./3.0);
t18 = t5+t6-t28-t29+x1_2-x2_2-x3_2+x4_2;
t30 = t3-1.0./4.0;
t31 = t5+t6-t28-t29+x1_2+x2_2-x3_2-x4_2;
t32 = t4.*t17.*t18.*2.0;
t33 = t2.*x1_2.*x3_1.*(2.0./3.0);
t36 = t2.*x1_1.*x3_2.*(2.0./3.0);
t34 = t7+t8+t9+t10+t11-t12-t13+t15-t19-t20-t21-t22-t23+t24+t25-t27+t33-t36;
t35 = 1.0./t34;
t37 = t5+t6-t28-t29-x1_2-x2_2+x3_2+x4_2;
t38 = t4.*t35.*t37.*2.0;
t39 = t7+t8+t9+t10-t11-t12-t13-t14-t15-t19-t20-t21-t22+t23+t24+t25+t26+t27;
t40 = 1.0./t39;
t41 = t30.*t37.*t40.*2.0;
t42 = t5+t6-t28-t29-x1_2+x2_2+x3_2-x4_2;
t43 = t4.*t37.*t40.*2.0;
t44 = t7+t8+t9+t10-t11+t12+t13-t15-t19-t20-t21-t22+t23-t24-t25+t27-t33+t36;
t45 = 1.0./t44;
t46 = t4.*t31.*t45.*2.0;
t47 = t4.*t42.*t45.*2.0;
t48 = t2.*x1_1.*(1.0./3.0);
t49 = t2.*x3_1.*(1.0./3.0);
t51 = t2.*x2_1.*(1.0./3.0);
t52 = t2.*x4_1.*(1.0./3.0);
t50 = t48+t49-t51-t52+x1_1-x2_1-x3_1+x4_1;
t53 = t48+t49-t51-t52+x1_1+x2_1-x3_1-x4_1;
t54 = t4.*t17.*t50.*2.0;
t55 = t48+t49-t51-t52-x1_1-x2_1+x3_1+x4_1;
t56 = t4.*t35.*t55.*2.0;
t57 = t30.*t40.*t55.*2.0;
t58 = t48+t49-t51-t52-x1_1+x2_1+x3_1-x4_1;
t59 = t4.*t45.*t53.*2.0;
t60 = t4.*t45.*t58.*2.0;
t62 = t2.*(1.0./3.0);
t61 = t62+1.0;
t63 = t62-1.0;
t64 = t61.^2;
t65 = t63.^2;
t66 = t4.*t17.*t53.*2.0;
t67 = t17.*t30.*t50.*2.0;
t68 = t4.*t17.*t31.*2.0;
t69 = t17.*t18.*t30.*2.0;
t70 = t54-t17.*t30.*t53.*2.0;
t72 = t17.*t30.*t31.*2.0;
t71 = t69-t72;
t73 = t30.*t35.*t50.*2.0;
t74 = t4.*t35.*t50.*2.0;
t75 = t18.*t30.*t35.*2.0;
t78 = t30.*t35.*t55.*2.0;
t76 = t74-t78;
t77 = t4.*t18.*t35.*2.0;
t80 = t30.*t35.*t37.*2.0;
t79 = t77-t80;
t81 = t30.*t40.*t58.*2.0;
t82 = t4.*t40.*t55.*2.0;
t83 = t41-t4.*t40.*t42.*2.0;
t84 = t30.*t45.*t53.*2.0;
t87 = t30.*t45.*t58.*2.0;
t85 = t84-t87;
t86 = t30.*t31.*t45.*2.0;
t90 = t30.*t42.*t45.*2.0;
t88 = t86-t90;
t89 = t59-t60;
t91 = t46-t47;
out1 = reshape([t71.*u3_1-u4_1.*(t32-t17.*t30.*t31.*2.0)+u2_1.*(t68-t17.*t18.*t30.*2.0)+u1_1.*(t32-t4.*t17.*(t5+t6+x1_2+x2_2-x3_2-x4_2-t2.*x2_2.*(1.0./3.0)-t2.*x4_2.*(1.0./3.0)).*2.0),t70.*u4_2-u2_2.*(t66-t17.*t30.*t50.*2.0)-u3_2.*(t67-t17.*t30.*t53.*2.0)-u1_2.*(t54-t4.*t17.*(t48+t49+x1_1+x2_1-x3_1-x4_1-t2.*x2_1.*(1.0./3.0)-t2.*x4_1.*(1.0./3.0)).*2.0),-(t64.*u1_1.*(1.0./4.0)+t65.*u3_1.*(1.0./4.0)-t61.*t63.*u2_1.*(1.0./4.0)-t61.*t63.*u4_1.*(1.0./4.0))./(t64.*x1_1.*(1.0./4.0)+t65.*x3_1.*(1.0./4.0)-t61.*t63.*x2_1.*(1.0./4.0)-t61.*t63.*x4_1.*(1.0./4.0)),-t71.*u3_2-t70.*u4_1+u3_1.*(t67-t17.*t30.*t53.*2.0)-u1_2.*(t32-t68)+u1_1.*(t54-t66)+u4_2.*(t32-t72)+u2_1.*(t66-t67)-u2_2.*(t68-t69),0.0,0.0,-t79.*u3_1-u2_1.*(t38-t4.*t18.*t35.*2.0)+u1_1.*(t38-t18.*t30.*t35.*2.0)+u4_1.*(t75-t30.*t35.*t37.*2.0),t76.*u3_2+u2_2.*(t56-t4.*t35.*t50.*2.0)-u1_2.*(t56-t30.*t35.*t50.*2.0)-u4_2.*(t73-t30.*t35.*t55.*2.0),-(t64.*u2_1.*(1.0./4.0)+t65.*u4_1.*(1.0./4.0)-t61.*t63.*u1_1.*(1.0./4.0)-t61.*t63.*u3_1.*(1.0./4.0))./(t64.*x2_1.*(1.0./4.0)+t65.*x4_1.*(1.0./4.0)-t61.*t63.*x1_1.*(1.0./4.0)-t61.*t63.*x3_1.*(1.0./4.0)),-t76.*u3_1+t79.*u3_2-u1_2.*(t38-t75)+u2_2.*(t38-t77)+u1_1.*(t56-t73)-u2_1.*(t56-t74)+u4_1.*(t73-t78)-u4_2.*(t75-t80),0.0,0.0,t83.*u2_1-u3_1.*(t43-t4.*t40.*t42.*2.0)-u1_1.*(t41-t30.*t40.*t42.*2.0)+u4_1.*(t43-t30.*t40.*t42.*2.0),-u2_2.*(t57-t4.*t40.*t58.*2.0)+u1_2.*(t57-t30.*t40.*t58.*2.0)+u3_2.*(t82-t4.*t40.*(t48+t49-t51-t52-x1_1+x2_1+x3_1-x4_1).*2.0)+u4_2.*(t81-t4.*t40.*(t48+t49-t51-t52-x1_1-x2_1+x3_1+x4_1).*2.0),-(t65.*u1_1.*(1.0./4.0)+t64.*u3_1.*(1.0./4.0)-t61.*t63.*u2_1.*(1.0./4.0)-t61.*t63.*u4_1.*(1.0./4.0))./(t65.*x1_1.*(1.0./4.0)+t64.*x3_1.*(1.0./4.0)-t61.*t63.*x2_1.*(1.0./4.0)-t61.*t63.*x4_1.*(1.0./4.0)),-t83.*u2_2-u3_1.*(t82-t4.*t40.*t58.*2.0)-u1_1.*(t57-t81)-u4_1.*(t81-t82)+u3_2.*(t43-t4.*t40.*(t5+t6-t28-t29-x1_2+x2_2+x3_2-x4_2).*2.0)+u1_2.*(t41-t30.*t40.*(t5+t6-t28-t29-x1_2+x2_2+x3_2-x4_2).*2.0)+u2_1.*(t57-t4.*t40.*(t48+t49-t51-t52-x1_1+x2_1+x3_1-x4_1).*2.0)+u4_2.*(t30.*t40.*t42.*2.0-t4.*t40.*(t5+t6-t28-t29-x1_2-x2_2+x3_2+x4_2).*2.0),0.0,0.0,-t88.*u2_1-t91.*u4_1-u1_1.*(t47-t30.*t31.*t45.*2.0)+u3_1.*(t46-t30.*t42.*t45.*2.0),t85.*u2_2+t89.*u4_2+u1_2.*(t60-t30.*t45.*t53.*2.0)-u3_2.*(t59-t30.*t45.*t58.*2.0),-(t65.*u2_1.*(1.0./4.0)+t64.*u4_1.*(1.0./4.0)-t61.*t63.*u1_1.*(1.0./4.0)-t61.*t63.*u3_1.*(1.0./4.0))./(t65.*x2_1.*(1.0./4.0)+t64.*x4_1.*(1.0./4.0)-t61.*t63.*x1_1.*(1.0./4.0)-t61.*t63.*x3_1.*(1.0./4.0)),-t85.*u2_1+t88.*u2_2-t89.*u4_1+t91.*u4_2+u1_2.*(t47-t86)-u1_1.*(t60-t84)-u3_2.*(t46-t90)+u3_1.*(t59-t87),0.0,0.0],[6,4]);