-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathStrainsDifference.m
198 lines (196 loc) · 6.75 KB
/
StrainsDifference.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
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
function out1 = StrainsDifference(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)
%STRAINSDIFFERENCE
% OUT1 = STRAINSDIFFERENCE(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:57:59
t2 = sqrt(3.0);
t3 = t2.*(1.0./1.2e1);
t4 = t3+1.0./4.0;
t5 = t2.*x1_1.*(1.0./3.0);
t6 = t2.*x3_1.*(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_1.*(1.0./3.0);
t29 = t2.*x4_1.*(1.0./3.0);
t18 = t5+t6-t28-t29+x1_1-x2_1-x3_1+x4_1;
t30 = t3-1.0./4.0;
t31 = t5+t6-t28-t29+x1_1+x2_1-x3_1-x4_1;
t32 = t4.*t17.*t18.*2.0;
t71 = t4.*t17.*t31.*2.0;
t72 = t17.*t18.*t30.*2.0;
t73 = t17.*t30.*t31.*2.0;
t74 = t72-t73;
t133 = t74.*u3_1;
t134 = u4_1.*(t32-t73);
t33 = t133-t134+u2_1.*(t71-t17.*t18.*t30.*2.0)+u1_1.*(t32-t4.*t17.*(t5+t6+x1_1+x2_1-x3_1-x4_1-t2.*x2_1.*(1.0./3.0)-t2.*x4_1.*(1.0./3.0)).*2.0);
t34 = t2.*x1_2.*(1.0./3.0);
t35 = t2.*x3_2.*(1.0./3.0);
t37 = t2.*x2_2.*(1.0./3.0);
t38 = t2.*x4_2.*(1.0./3.0);
t36 = t34+t35-t37-t38+x1_2-x2_2-x3_2+x4_2;
t39 = t34+t35-t37-t38+x1_2+x2_2-x3_2-x4_2;
t40 = t4.*t17.*t36.*2.0;
t76 = t4.*t17.*t39.*2.0;
t77 = t17.*t30.*t36.*2.0;
t78 = t17.*t30.*t39.*2.0;
t79 = t77-t78;
t141 = t79.*u3_1;
t142 = u4_1.*(t40-t78);
t41 = t141-t142+u2_1.*(t76-t17.*t30.*t36.*2.0)+u1_1.*(t40-t4.*t17.*(t34+t35+x1_2+x2_2-x3_2-x4_2-t2.*x2_2.*(1.0./3.0)-t2.*x4_2.*(1.0./3.0)).*2.0);
t42 = t2.*x1_2.*x3_1.*(2.0./3.0);
t45 = t2.*x1_1.*x3_2.*(2.0./3.0);
t43 = t7+t8+t9+t10+t11-t12-t13+t15-t19-t20-t21-t22-t23+t24+t25-t27+t42-t45;
t44 = 1.0./t43;
t46 = t5+t6-t28-t29-x1_1-x2_1+x3_1+x4_1;
t47 = t4.*t18.*t44.*2.0;
t48 = t18.*t30.*t44.*2.0;
t81 = t4.*t44.*t46.*2.0;
t82 = t48-t81;
t83 = t47-t81;
t84 = t30.*t44.*t46.*2.0;
t85 = t47-t84;
t86 = t48-t84;
t147 = t82.*u1_1;
t148 = t83.*u2_1;
t149 = t85.*u3_1;
t150 = t86.*u4_1;
t49 = t147-t148+t149-t150;
t50 = t34+t35-t37-t38-x1_2-x2_2+x3_2+x4_2;
t51 = t4.*t44.*t50.*2.0;
t88 = t30.*t36.*t44.*2.0;
t89 = t4.*t36.*t44.*2.0;
t90 = t30.*t44.*t50.*2.0;
t91 = t89-t90;
t92 = t88-t90;
t156 = t91.*u3_1;
t157 = t92.*u4_1;
t52 = t156-t157+u2_1.*(t51-t4.*t36.*t44.*2.0)-u1_1.*(t51-t30.*t36.*t44.*2.0);
t53 = t7+t8+t9+t10-t11-t12-t13-t14-t15-t19-t20-t21-t22+t23+t24+t25+t26+t27;
t54 = 1.0./t53;
t55 = t30.*t46.*t54.*2.0;
t56 = t5+t6-t28-t29-x1_1+x2_1+x3_1-x4_1;
t57 = t4.*t46.*t54.*2.0;
t94 = t30.*t54.*t56.*2.0;
t95 = t55-t94;
t96 = t4.*t54.*t56.*2.0;
t97 = t57-t96;
t98 = t57-t94;
t162 = t95.*u1_1;
t164 = t98.*u4_1;
t58 = -t162+t164-t97.*u3_1+u2_1.*(t55-t96);
t59 = t30.*t50.*t54.*2.0;
t60 = t34+t35-t37-t38-x1_2+x2_2+x3_2-x4_2;
t61 = t4.*t50.*t54.*2.0;
t100 = t30.*t54.*t60.*2.0;
t101 = t59-t100;
t102 = t4.*t54.*t60.*2.0;
t103 = t61-t102;
t104 = t61-t100;
t169 = t101.*u1_1;
t171 = t104.*u4_1;
t62 = -t169+t171-t103.*u3_1+u2_1.*(t59-t102);
t63 = t7+t8+t9+t10-t11+t12+t13-t15-t19-t20-t21-t22+t23-t24-t25+t27-t42+t45;
t64 = 1.0./t63;
t65 = t4.*t31.*t64.*2.0;
t66 = t4.*t56.*t64.*2.0;
t106 = t30.*t31.*t64.*2.0;
t107 = t30.*t56.*t64.*2.0;
t108 = t106-t107;
t109 = t65-t107;
t110 = t65-t66;
t177 = t108.*u2_1;
t178 = t109.*u3_1;
t179 = t110.*u4_1;
t67 = t177-t178+t179+u1_1.*(t66-t30.*t31.*t64.*2.0);
t68 = t4.*t39.*t64.*2.0;
t69 = t4.*t60.*t64.*2.0;
t112 = t30.*t39.*t64.*2.0;
t113 = t30.*t60.*t64.*2.0;
t114 = t112-t113;
t115 = t68-t113;
t116 = t68-t69;
t185 = t114.*u2_1;
t186 = t115.*u3_1;
t187 = t116.*u4_1;
t70 = t185-t186+t187+u1_1.*(t69-t30.*t39.*t64.*2.0);
t131 = t32-t71;
t132 = t71-t72;
t135 = t131.*u1_2;
t136 = t132.*u2_2;
t137 = t74.*u3_2;
t138 = u4_2.*(t32-t73);
t75 = t135+t136+t137-t138;
t139 = t40-t76;
t140 = t76-t77;
t143 = t139.*u1_2;
t144 = t140.*u2_2;
t145 = t79.*u3_2;
t146 = u4_2.*(t40-t78);
t80 = t143+t144+t145-t146;
t151 = t82.*u1_2;
t152 = t83.*u2_2;
t153 = t85.*u3_2;
t154 = t86.*u4_2;
t87 = t151-t152+t153-t154;
t155 = t51-t88;
t158 = t155.*u1_2;
t159 = u2_2.*(t51-t89);
t160 = t91.*u3_2;
t161 = t92.*u4_2;
t93 = -t158+t159+t160-t161;
t163 = t4.*t54.*(t5+t6-t28-t29-x1_1+x2_1+x3_1-x4_1).*2.0;
t165 = t95.*u1_2;
t166 = u2_2.*(t55-t163);
t168 = t98.*u4_2;
t99 = t165-t166-t168+t97.*u3_2;
t170 = t4.*t54.*(t34+t35-t37-t38-x1_2+x2_2+x3_2-x4_2).*2.0;
t172 = t101.*u1_2;
t173 = u2_2.*(t59-t170);
t175 = t104.*u4_2;
t105 = t172-t173-t175+t103.*u3_2;
t176 = t66-t106;
t180 = t176.*u1_2;
t181 = t108.*u2_2;
t182 = t109.*u3_2;
t183 = t110.*u4_2;
t111 = t180+t181-t182+t183;
t184 = t69-t112;
t188 = t184.*u1_2;
t189 = t114.*u2_2;
t190 = t115.*u3_2;
t191 = t116.*u4_2;
t117 = t188+t189-t190+t191;
t119 = t2.*(1.0./3.0);
t118 = t119+1.0;
t120 = t119-1.0;
t122 = t118.^2;
t123 = t120.^2;
t125 = t118.*t120.*u2_1.*(1.0./4.0);
t126 = t118.*t120.*u4_1.*(1.0./4.0);
t121 = -t125-t126+t122.*u1_1.*(1.0./4.0)+t123.*u3_1.*(1.0./4.0);
t128 = t118.*t120.*u1_1.*(1.0./4.0);
t129 = t118.*t120.*u3_1.*(1.0./4.0);
t124 = -t128-t129+t122.*u2_1.*(1.0./4.0)+t123.*u4_1.*(1.0./4.0);
t127 = t125+t126-t123.*u1_1.*(1.0./4.0)-t122.*u3_1.*(1.0./4.0);
t130 = t128+t129-t123.*u2_1.*(1.0./4.0)-t122.*u4_1.*(1.0./4.0);
t167 = t57-t163;
t174 = t61-t170;
out1 = reshape([t33.^2.*(-1.0./2.0)-t41.^2.*(1.0./2.0),t75.^2.*(-1.0./2.0)-t80.^2.*(1.0./2.0),t121.^2.*1.0./(t122.*x1_1.*(1.0./4.0)+t123.*x3_1.*(1.0./4.0)-t118.*t120.*x2_1.*(1.0./4.0)-t118.*t120.*x4_1.*(1.0./4.0)).^2.*(-1.0./2.0),(t133-t134+t131.*u1_1+t132.*u2_1).*(t135+t136+t137-t138)+(t141-t142+t139.*u1_1+t140.*u2_1).*(t143+t144+t145-t146),0.0,0.0,t49.^2.*(-1.0./2.0)-t52.^2.*(1.0./2.0),t87.^2.*(-1.0./2.0)-t93.^2.*(1.0./2.0),t124.^2.*1.0./(t122.*x2_1.*(1.0./4.0)+t123.*x4_1.*(1.0./4.0)-t118.*t120.*x1_1.*(1.0./4.0)-t118.*t120.*x3_1.*(1.0./4.0)).^2.*(-1.0./2.0),t93.*(t156-t157-t155.*u1_1+u2_1.*(t51-t89))+t49.*t87,0.0,0.0,t58.^2.*(-1.0./2.0)-t62.^2.*(1.0./2.0),t99.^2.*(-1.0./2.0)-t105.^2.*(1.0./2.0),t127.^2.*1.0./(t123.*x1_1.*(1.0./4.0)+t122.*x3_1.*(1.0./4.0)-t118.*t120.*x2_1.*(1.0./4.0)-t118.*t120.*x4_1.*(1.0./4.0)).^2.*(-1.0./2.0),(t162-t164+t167.*u3_1-u2_1.*(t55-t163)).*(t165-t166-t168+t167.*u3_2)+(t169-t171+t174.*u3_1-u2_1.*(t59-t170)).*(t172-t173-t175+t174.*u3_2),0.0,0.0,t67.^2.*(-1.0./2.0)-t70.^2.*(1.0./2.0),t111.^2.*(-1.0./2.0)-t117.^2.*(1.0./2.0),t130.^2.*1.0./(t123.*x2_1.*(1.0./4.0)+t122.*x4_1.*(1.0./4.0)-t118.*t120.*x1_1.*(1.0./4.0)-t118.*t120.*x3_1.*(1.0./4.0)).^2.*(-1.0./2.0),(t177-t178+t179+t176.*u1_1).*(t180+t181-t182+t183)+(t185-t186+t187+t184.*u1_1).*(t188+t189-t190+t191),0.0,0.0],[6,4]);