-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathStrainsAlmansi.m
More file actions
222 lines (220 loc) · 7.89 KB
/
Copy pathStrainsAlmansi.m
File metadata and controls
222 lines (220 loc) · 7.89 KB
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
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
function out1 = StrainsAlmansi(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)
%STRAINSALMANSI
% OUT1 = STRAINSALMANSI(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:17
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_1.*(1.0./3.0);
t34 = t2.*x3_1.*(1.0./3.0);
t36 = t2.*x2_1.*(1.0./3.0);
t37 = t2.*x4_1.*(1.0./3.0);
t35 = t33+t34-t36-t37+x1_1-x2_1-x3_1+x4_1;
t38 = t33+t34-t36-t37+x1_1+x2_1-x3_1-x4_1;
t39 = t4.*t17.*t35.*2.0;
t88 = t4.*t17.*t38.*2.0;
t89 = t17.*t30.*t35.*2.0;
t90 = t17.*t30.*t38.*2.0;
t91 = t89-t90;
t162 = t39-t90;
t163 = t162.*u4_1;
t166 = t91.*u3_1;
t40 = -t163+t166+u2_1.*(t88-t17.*t30.*t35.*2.0)+u1_1.*(t39-t4.*t17.*(t33+t34+x1_1+x2_1-x3_1-x4_1-t2.*x2_1.*(1.0./3.0)-t2.*x4_1.*(1.0./3.0)).*2.0);
t41 = t4.*t17.*t31.*2.0;
t42 = t17.*t18.*t30.*2.0;
t45 = t17.*t30.*t31.*2.0;
t43 = t42-t45;
t44 = t43.*u3_1;
t96 = t32-t41;
t97 = t41-t42;
t170 = t96.*u1_1;
t171 = t97.*u2_1;
t172 = u4_1.*(t32-t45);
t46 = t44+t170+t171-t172;
t47 = t2.*x1_2.*x3_1.*(2.0./3.0);
t50 = t2.*x1_1.*x3_2.*(2.0./3.0);
t48 = t7+t8+t9+t10+t11-t12-t13+t15-t19-t20-t21-t22-t23+t24+t25-t27+t47-t50;
t49 = 1.0./t48;
t51 = t5+t6-t28-t29-x1_2-x2_2+x3_2+x4_2;
t52 = t4.*t18.*t49.*2.0;
t53 = t18.*t30.*t49.*2.0;
t54 = t33+t34-t36-t37-x1_1-x2_1+x3_1+x4_1;
t55 = t4.*t49.*t54.*2.0;
t99 = t30.*t35.*t49.*2.0;
t100 = t4.*t35.*t49.*2.0;
t101 = t30.*t49.*t54.*2.0;
t102 = t100-t101;
t103 = t99-t101;
t175 = t102.*u3_1;
t177 = t103.*u4_1;
t56 = t175-t177+u2_1.*(t55-t4.*t35.*t49.*2.0)-u1_1.*(t55-t30.*t35.*t49.*2.0);
t58 = t4.*t49.*t51.*2.0;
t57 = t53-t58;
t60 = t30.*t49.*t51.*2.0;
t59 = t52-t60;
t108 = t52-t58;
t109 = t53-t60;
t182 = t57.*u1_1;
t183 = t108.*u2_1;
t184 = t59.*u3_1;
t185 = t109.*u4_1;
t61 = t182-t183+t184-t185;
t62 = t7+t8+t9+t10-t11-t12-t13-t14-t15-t19-t20-t21-t22+t23+t24+t25+t26+t27;
t63 = 1.0./t62;
t64 = t30.*t51.*t63.*2.0;
t65 = t5+t6-t28-t29-x1_2+x2_2+x3_2-x4_2;
t66 = t4.*t51.*t63.*2.0;
t67 = t30.*t54.*t63.*2.0;
t68 = t33+t34-t36-t37-x1_1+x2_1+x3_1-x4_1;
t69 = t4.*t54.*t63.*2.0;
t111 = t4.*t63.*t68.*2.0;
t112 = t67-t111;
t113 = t30.*t63.*t68.*2.0;
t114 = t67-t113;
t115 = t69-t111;
t116 = t69-t113;
t186 = t114.*u1_1;
t189 = t116.*u4_1;
t70 = -t186+t189+t112.*u2_1-t115.*u3_1;
t73 = t30.*t63.*t65.*2.0;
t71 = t64-t73;
t72 = t66-t4.*t63.*t65.*2.0;
t118 = t4.*t63.*(t5+t6-t28-t29-x1_2+x2_2+x3_2-x4_2).*2.0;
t119 = t66-t73;
t188 = t64-t118;
t196 = t71.*u1_1;
t198 = t119.*u4_1;
t202 = t188.*u2_1;
t74 = t196-t198-t202+t72.*u3_1;
t75 = t7+t8+t9+t10-t11+t12+t13-t15-t19-t20-t21-t22+t23-t24-t25+t27-t47+t50;
t76 = 1.0./t75;
t77 = t30.*t31.*t76.*2.0;
t78 = t4.*t31.*t76.*2.0;
t79 = t30.*t38.*t76.*2.0;
t80 = t4.*t38.*t76.*2.0;
t121 = t4.*t68.*t76.*2.0;
t122 = t79-t121;
t123 = t30.*t68.*t76.*2.0;
t124 = t79-t123;
t125 = t80-t123;
t126 = t80-t121;
t203 = t122.*u1_1;
t204 = t124.*u2_1;
t206 = t125.*u3_1;
t208 = t126.*u4_1;
t81 = t203-t204+t206-t208;
t86 = t4.*t65.*t76.*2.0;
t82 = t77-t86;
t83 = t82.*u1_1;
t85 = t30.*t65.*t76.*2.0;
t84 = t77-t85;
t130 = t78-t85;
t131 = t78-t86;
t213 = t84.*u2_1;
t214 = t130.*u3_1;
t215 = t131.*u4_1;
t87 = t83-t213+t214-t215;
t92 = t39-t88;
t93 = t88-t89;
t94 = u4_2.*(t39-t90);
t167 = t92.*u1_2;
t168 = t93.*u2_2;
t169 = t91.*u3_2;
t95 = -t94+t167+t168+t169;
t160 = t96.*u1_2;
t161 = t97.*u2_2;
t164 = t43.*u3_2;
t165 = u4_2.*(t32-t45);
t98 = t160+t161+t164-t165;
t104 = t55-t99;
t173 = t55-t100;
t105 = t173.*u2_2;
t106 = t102.*u3_2;
t180 = t104.*u1_2;
t181 = t103.*u4_2;
t107 = t105+t106-t180-t181;
t174 = t57.*u1_2;
t176 = t108.*u2_2;
t178 = t59.*u3_2;
t179 = t109.*u4_2;
t110 = t174-t176+t178-t179;
t187 = t4.*t63.*(t33+t34-t36-t37-x1_1+x2_1+x3_1-x4_1).*2.0;
t193 = t114.*u1_2;
t194 = u2_2.*(t67-t187);
t195 = t116.*u4_2;
t117 = t193-t194-t195+t115.*u3_2;
t197 = t66-t118;
t190 = t197.*u3_2;
t199 = t71.*u1_2;
t200 = u2_2.*(t64-t118);
t201 = t119.*u4_2;
t120 = t190+t199-t200-t201;
t127 = t124.*u2_2;
t128 = t126.*u4_2;
t211 = t122.*u1_2;
t212 = t125.*u3_2;
t129 = t127+t128-t211-t212;
t205 = t82.*u1_2;
t207 = t84.*u2_2;
t209 = t130.*u3_2;
t210 = t131.*u4_2;
t132 = t205-t207+t209-t210;
t134 = t2.*(1.0./3.0);
t133 = t134+1.0;
t135 = t134-1.0;
t136 = t133.^2;
t137 = t135.^2;
t138 = t136.*u1_1.*(1.0./4.0);
t139 = t137.*u3_1.*(1.0./4.0);
t150 = t133.*t135.*u2_1.*(1.0./4.0);
t151 = t133.*t135.*u4_1.*(1.0./4.0);
t140 = t138+t139-t150-t151;
t141 = t136.*x1_1.*(1.0./4.0);
t142 = t137.*x3_1.*(1.0./4.0);
t152 = t133.*t135.*x2_1.*(1.0./4.0);
t153 = t133.*t135.*x4_1.*(1.0./4.0);
t143 = t141+t142-t152-t153;
t144 = t136.*u2_1.*(1.0./4.0);
t145 = t137.*u4_1.*(1.0./4.0);
t155 = t133.*t135.*u1_1.*(1.0./4.0);
t156 = t133.*t135.*u3_1.*(1.0./4.0);
t146 = t144+t145-t155-t156;
t147 = t136.*x2_1.*(1.0./4.0);
t148 = t137.*x4_1.*(1.0./4.0);
t157 = t133.*t135.*x1_1.*(1.0./4.0);
t158 = t133.*t135.*x3_1.*(1.0./4.0);
t149 = t147+t148-t157-t158;
t154 = t150+t151-t137.*u1_1.*(1.0./4.0)-t136.*u3_1.*(1.0./4.0);
t159 = t155+t156-t137.*u2_1.*(1.0./4.0)-t136.*u4_1.*(1.0./4.0);
t191 = u2_1.*(t67-t187);
t192 = t69-t187;
out1 = reshape([t44+u2_1.*(t41-t17.*t18.*t30.*2.0)-u4_1.*(t32-t17.*t30.*t31.*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)+t40.^2.*(1.0./2.0)+t46.^2.*(1.0./2.0),t94-t92.*u1_2-t93.*u2_2-t91.*u3_2+t95.^2.*(1.0./2.0)+t98.^2.*(1.0./2.0),t140.^2.*1.0./t143.^2.*(1.0./2.0)-t140./t143,-t160-t161-t163-t164+t165-t95.*(-t163+t166+t92.*u1_1+t93.*u2_1)-t46.*t98+u1_1.*(t39-t88)+u2_1.*(t88-t89)+u3_1.*(t89-t90),0.0,0.0,-t57.*u1_1-t59.*u3_1+u2_1.*(t52-t4.*t49.*t51.*2.0)+u4_1.*(t53-t30.*t49.*t51.*2.0)+t56.^2.*(1.0./2.0)+t61.^2.*(1.0./2.0),t105+t106-t104.*u1_2-t103.*u4_2+t107.^2.*(1.0./2.0)+t110.^2.*(1.0./2.0),t146.^2.*1.0./t149.^2.*(1.0./2.0)-t146./t149,t174-t175-t176+t177+t178-t179-t107.*(t175-t177-t104.*u1_1+t173.*u2_1)-t61.*t110-t173.*u2_1+u1_1.*(t55-t99),0.0,0.0,-t71.*u1_1-t72.*u3_1+u2_1.*(t64-t4.*t63.*t65.*2.0)+u4_1.*(t66-t30.*t63.*t65.*2.0)+t70.^2.*(1.0./2.0)+t74.^2.*(1.0./2.0),-t112.*u2_2+u3_2.*(t69-t4.*t63.*(t33+t34-t36-t37-x1_1+x2_1+x3_1-x4_1).*2.0)+u1_2.*(t67-t30.*t63.*(t33+t34-t36-t37-x1_1+x2_1+x3_1-x4_1).*2.0)+u4_2.*(t113-t4.*t63.*(t33+t34-t36-t37-x1_1-x2_1+x3_1+x4_1).*2.0)+t117.^2.*(1.0./2.0)+t120.^2.*(1.0./2.0),t154.^2.*1.0./(t152+t153-t137.*x1_1.*(1.0./4.0)-t136.*x3_1.*(1.0./4.0)).^2.*(1.0./2.0)-(t150+t151-t137.*u1_1.*(1.0./4.0)-t136.*u3_1.*(1.0./4.0))./(t152+t153-t137.*x1_1.*(1.0./4.0)-t136.*x3_1.*(1.0./4.0)),-t186+t189+t190+t191-t188.*u2_2-t192.*u3_1+u1_2.*(t64-t30.*t63.*(t5+t6-t28-t29-x1_2+x2_2+x3_2-x4_2).*2.0)+u4_2.*(t73-t4.*t63.*(t5+t6-t28-t29-x1_2-x2_2+x3_2+x4_2).*2.0)-t120.*(t196-t198-t202+t197.*u3_1)-(t186-t189-t191+t192.*u3_1).*(t193-t194-t195+t192.*u3_2),0.0,0.0,t83-t84.*u2_1-u4_1.*(t78-t4.*t65.*t76.*2.0)+u3_1.*(t78-t30.*t65.*t76.*2.0)+t81.^2.*(1.0./2.0)+t87.^2.*(1.0./2.0),t127+t128-t122.*u1_2-t125.*u3_2+t129.^2.*(1.0./2.0)+t132.^2.*(1.0./2.0),t159.^2.*1.0./(t157+t158-t137.*x2_1.*(1.0./4.0)-t136.*x4_1.*(1.0./4.0)).^2.*(1.0./2.0)-(t155+t156-t137.*u2_1.*(1.0./4.0)-t136.*u4_1.*(1.0./4.0))./(t157+t158-t137.*x2_1.*(1.0./4.0)-t136.*x4_1.*(1.0./4.0)),t203-t204-t205+t206+t207-t208-t209+t210-t87.*t132+t81.*(t127+t128-t211-t212),0.0,0.0],[6,4]);