1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
| ...
j = 1;
for n = trianglesInLay
Kc = E/(3*(1-2*nu));
G = E/(2*(1+nu));
A = Kc + 4/3*G;
B = Kc - 2/3*G;
C = 2*G;
X2 = 0;
syms X1;
F1(1, n) = eval( int( -(1/2)*delta*(A*(2*a(1, n)*X1+c(1, n)*X2+d(1, n))*u(1, 1, n)+2*B*(a(1, n)*X1^2+b(1, n)*X2^2+c(1, n)*X1*X2+d(1, n)*X1+h(1, n)*X2+g(1, n))*u(2, 1, n)/delta+A*(2*a(2, n)*X1+c(2, n)*X2+d(2, n))*u(1, 2, n)+2*B*(a(2, n)*X1^2+b(2, n)*X2^2+c(2, n)*X1*X2+d(2, n)*X1+h(2, n)*X2+g(2, n))*u(2, 2, n)/delta+A*(2*a(3, n)*X1+c(3, n)*X2+d(3, n))*u(1, 3, n)+2*B*(a(3, n)*X1^2+b(3, n)*X2^2+c(3, n)*X1*X2+d(3, n)*X1+h(3, n)*X2+g(3, n))*u(2, 3, n)/delta+A*(2*a(4, n)*X1+c(4, n)*X2+d(4, n))*u(1, 4, n)+2*B*(a(4, n)*X1^2+b(4, n)*X2^2+c(4, n)*X1*X2+d(4, n)*X1+h(4, n)*X2+g(4, n))*u(2, 4, n)/delta+A*(2*a(5, n)*X1+c(5, n)*X2+d(5, n))*u(1, 5, n)+2*B*(a(5, n)*X1^2+b(5, n)*X2^2+c(5, n)*X1*X2+d(5, n)*X1+h(5, n)*X2+g(5, n))*u(2, 5, n)/delta+A*(2*a(6, n)*X1+c(6, n)*X2+d(6, n))*u(1, 6, n)+2*B*(a(6, n)*X1^2+b(6, n)*X2^2+c(6, n)*X1*X2+d(6, n)*X1+h(6, n)*X2+g(6, n))*u(2, 6, n)/delta)*(2*a(1, n)*X1+c(1, n)*X2+d(1, n)) ,X1, p(1, ed2(j)), p(1, ed2(j+1)) ) );
F1(2, n) = eval( int( -(1/2)*delta*(A*(2*a(1, n)*X1+c(1, n)*X2+d(1, n))*u(1, 1, n)+2*B*(a(1, n)*X1^2+b(1, n)*X2^2+c(1, n)*X1*X2+d(1, n)*X1+h(1, n)*X2+g(1, n))*u(2, 1, n)/delta+A*(2*a(2, n)*X1+c(2, n)*X2+d(2, n))*u(1, 2, n)+2*B*(a(2, n)*X1^2+b(2, n)*X2^2+c(2, n)*X1*X2+d(2, n)*X1+h(2, n)*X2+g(2, n))*u(2, 2, n)/delta+A*(2*a(3, n)*X1+c(3, n)*X2+d(3, n))*u(1, 3, n)+2*B*(a(3, n)*X1^2+b(3, n)*X2^2+c(3, n)*X1*X2+d(3, n)*X1+h(3, n)*X2+g(3, n))*u(2, 3, n)/delta+A*(2*a(4, n)*X1+c(4, n)*X2+d(4, n))*u(1, 4, n)+2*B*(a(4, n)*X1^2+b(4, n)*X2^2+c(4, n)*X1*X2+d(4, n)*X1+h(4, n)*X2+g(4, n))*u(2, 4, n)/delta+A*(2*a(5, n)*X1+c(5, n)*X2+d(5, n))*u(1, 5, n)+2*B*(a(5, n)*X1^2+b(5, n)*X2^2+c(5, n)*X1*X2+d(5, n)*X1+h(5, n)*X2+g(5, n))*u(2, 5, n)/delta+A*(2*a(6, n)*X1+c(6, n)*X2+d(6, n))*u(1, 6, n)+2*B*(a(6, n)*X1^2+b(6, n)*X2^2+c(6, n)*X1*X2+d(6, n)*X1+h(6, n)*X2+g(6, n))*u(2, 6, n)/delta)*(2*a(2, n)*X1+c(2, n)*X2+d(2, n)) ,X1, p(1, ed2(j)), p(1, ed2(j+1)) ) );
F1(3, n) = eval( int( -(1/2)*delta*(A*(2*a(1, n)*X1+c(1, n)*X2+d(1, n))*u(1, 1, n)+2*B*(a(1, n)*X1^2+b(1, n)*X2^2+c(1, n)*X1*X2+d(1, n)*X1+h(1, n)*X2+g(1, n))*u(2, 1, n)/delta+A*(2*a(2, n)*X1+c(2, n)*X2+d(2, n))*u(1, 2, n)+2*B*(a(2, n)*X1^2+b(2, n)*X2^2+c(2, n)*X1*X2+d(2, n)*X1+h(2, n)*X2+g(2, n))*u(2, 2, n)/delta+A*(2*a(3, n)*X1+c(3, n)*X2+d(3, n))*u(1, 3, n)+2*B*(a(3, n)*X1^2+b(3, n)*X2^2+c(3, n)*X1*X2+d(3, n)*X1+h(3, n)*X2+g(3, n))*u(2, 3, n)/delta+A*(2*a(4, n)*X1+c(4, n)*X2+d(4, n))*u(1, 4, n)+2*B*(a(4, n)*X1^2+b(4, n)*X2^2+c(4, n)*X1*X2+d(4, n)*X1+h(4, n)*X2+g(4, n))*u(2, 4, n)/delta+A*(2*a(5, n)*X1+c(5, n)*X2+d(5, n))*u(1, 5, n)+2*B*(a(5, n)*X1^2+b(5, n)*X2^2+c(5, n)*X1*X2+d(5, n)*X1+h(5, n)*X2+g(5, n))*u(2, 5, n)/delta+A*(2*a(6, n)*X1+c(6, n)*X2+d(6, n))*u(1, 6, n)+2*B*(a(6, n)*X1^2+b(6, n)*X2^2+c(6, n)*X1*X2+d(6, n)*X1+h(6, n)*X2+g(6, n))*u(2, 6, n)/delta)*(2*a(3, n)*X1+c(3, n)*X2+d(3, n)) ,X1, p(1, ed2(j)), p(1, ed2(j+1)) ) );
F2(1, n) = eval( int( -(2*A*(a(1, n)*X1^2+b(1, n)*X2^2+c(1, n)*X1*X2+d(1, n)*X1+h(1, n)*X2+g(1, n))*u(2, 1, n)/delta+B*(2*a(1, n)*X1+c(1, n)*X2+d(1, n))*u(1, 1, n)+2*A*(a(2, n)*X1^2+b(2, n)*X2^2+c(2, n)*X1*X2+d(2, n)*X1+h(2, n)*X2+g(2, n))*u(2, 2, n)/delta+B*(2*a(2, n)*X1+c(2, n)*X2+d(2, n))*u(1, 2, n)+2*A*(a(3, n)*X1^2+b(3, n)*X2^2+c(3, n)*X1*X2+d(3, n)*X1+h(3, n)*X2+g(3, n))*u(2, 3, n)/delta+B*(2*a(3, n)*X1+c(3, n)*X2+d(3, n))*u(1, 3, n)+2*A*(a(4, n)*X1^2+b(4, n)*X2^2+c(4, n)*X1*X2+d(4, n)*X1+h(4, n)*X2+g(4, n))*u(2, 4, n)/delta+B*(2*a(4, n)*X1+c(4, n)*X2+d(4, n))*u(1, 4, n)+2*A*(a(5, n)*X1^2+b(5, n)*X2^2+c(5, n)*X1*X2+d(5, n)*X1+h(5, n)*X2+g(5, n))*u(2, 5, n)/delta+B*(2*a(5, n)*X1+c(5, n)*X2+d(5, n))*u(1, 5, n)+2*A*(a(6, n)*X1^2+b(6, n)*X2^2+c(6, n)*X1*X2+d(6, n)*X1+h(6, n)*X2+g(6, n))*u(2, 6, n)/delta+B*(2*a(6, n)*X1+c(6, n)*X2+d(6, n))*u(1, 6, n))*(a(1, n)*X1^2+b(1, n)*X2^2+c(1, n)*X1*X2+d(1, n)*X1+h(1, n)*X2+g(1, n)) ,X1, p(1, ed2(j)), p(1, ed2(j+1)) ) );
F2(2, n) = eval( int( -(2*A*(a(1, n)*X1^2+b(1, n)*X2^2+c(1, n)*X1*X2+d(1, n)*X1+h(1, n)*X2+g(1, n))*u(2, 1, n)/delta+B*(2*a(1, n)*X1+c(1, n)*X2+d(1, n))*u(1, 1, n)+2*A*(a(2, n)*X1^2+b(2, n)*X2^2+c(2, n)*X1*X2+d(2, n)*X1+h(2, n)*X2+g(2, n))*u(2, 2, n)/delta+B*(2*a(2, n)*X1+c(2, n)*X2+d(2, n))*u(1, 2, n)+2*A*(a(3, n)*X1^2+b(3, n)*X2^2+c(3, n)*X1*X2+d(3, n)*X1+h(3, n)*X2+g(3, n))*u(2, 3, n)/delta+B*(2*a(3, n)*X1+c(3, n)*X2+d(3, n))*u(1, 3, n)+2*A*(a(4, n)*X1^2+b(4, n)*X2^2+c(4, n)*X1*X2+d(4, n)*X1+h(4, n)*X2+g(4, n))*u(2, 4, n)/delta+B*(2*a(4, n)*X1+c(4, n)*X2+d(4, n))*u(1, 4, n)+2*A*(a(5, n)*X1^2+b(5, n)*X2^2+c(5, n)*X1*X2+d(5, n)*X1+h(5, n)*X2+g(5, n))*u(2, 5, n)/delta+B*(2*a(5, n)*X1+c(5, n)*X2+d(5, n))*u(1, 5, n)+2*A*(a(6, n)*X1^2+b(6, n)*X2^2+c(6, n)*X1*X2+d(6, n)*X1+h(6, n)*X2+g(6, n))*u(2, 6, n)/delta+B*(2*a(6, n)*X1+c(6, n)*X2+d(6, n))*u(1, 6, n))*(a(2, n)*X1^2+b(2, n)*X2^2+c(2, n)*X1*X2+d(2, n)*X1+h(2, n)*X2+g(2, n)) ,X1, p(1, ed2(j)), p(1, ed2(j+1)) ) );
F2(3, n) = eval( int( -(2*A*(a(1, n)*X1^2+b(1, n)*X2^2+c(1, n)*X1*X2+d(1, n)*X1+h(1, n)*X2+g(1, n))*u(2, 1, n)/delta+B*(2*a(1, n)*X1+c(1, n)*X2+d(1, n))*u(1, 1, n)+2*A*(a(2, n)*X1^2+b(2, n)*X2^2+c(2, n)*X1*X2+d(2, n)*X1+h(2, n)*X2+g(2, n))*u(2, 2, n)/delta+B*(2*a(2, n)*X1+c(2, n)*X2+d(2, n))*u(1, 2, n)+2*A*(a(3, n)*X1^2+b(3, n)*X2^2+c(3, n)*X1*X2+d(3, n)*X1+h(3, n)*X2+g(3, n))*u(2, 3, n)/delta+B*(2*a(3, n)*X1+c(3, n)*X2+d(3, n))*u(1, 3, n)+2*A*(a(4, n)*X1^2+b(4, n)*X2^2+c(4, n)*X1*X2+d(4, n)*X1+h(4, n)*X2+g(4, n))*u(2, 4, n)/delta+B*(2*a(4, n)*X1+c(4, n)*X2+d(4, n))*u(1, 4, n)+2*A*(a(5, n)*X1^2+b(5, n)*X2^2+c(5, n)*X1*X2+d(5, n)*X1+h(5, n)*X2+g(5, n))*u(2, 5, n)/delta+B*(2*a(5, n)*X1+c(5, n)*X2+d(5, n))*u(1, 5, n)+2*A*(a(6, n)*X1^2+b(6, n)*X2^2+c(6, n)*X1*X2+d(6, n)*X1+h(6, n)*X2+g(6, n))*u(2, 6, n)/delta+B*(2*a(6, n)*X1+c(6, n)*X2+d(6, n))*u(1, 6, n))*(a(3, n)*X1^2+b(3, n)*X2^2+c(3, n)*X1*X2+d(3, n)*X1+h(3, n)*X2+g(3, n)) ,X1, p(1, ed2(j)), p(1, ed2(j+1)) ) );
j = j + 1;
end
... |
Причём в коде я ничего не менял.
(
