“我”总是等于4因为for循环没有条件,将阻止它。因为“我”永远是4月底循环,QbarA, QbarB, QbarC永远不会设置,因此,当你得到“A”的行被设置,你会得到一个错误。
下面的代码可能是更接近你想要的,但你的输出函数将返回一个数组失配误差,因为,B和D是3 x3数组,N和M是1 x3数组。需要背景对你的算法进一步帮助。
clc,关闭所有、清晰变量
E1 = 40;
E2 = 1.8;
V12 = 0.25;
V21 = 0.01125;
G12 = 0.7;
t = 0.005;
z0 = - t * 2;
z1 = - t;
z2 = 0;
z3 = t;
z4 = t * 2;
θ= [0,90,-30];
Q11 = E1 / (1-V12 * V21);
12 = (V21 * E2) / (1-V12 * V21);
如= E2 / (1-V12 * V21);
Q66 = G12;
Nx = 1200 * 4 * t;
纽约= -600 * 4 * t;
Nxy = 600 * 4 * t;
我= 36/12;
Mx = -18/18;
Mxy = 0;
N = [Nx;纽约;Nxy);
M = [Mx;我;Mxy];
为i = 1:4
m = cosd(θ(i));
n =信德(θ(i));
Qbar11 = Q11。* m ^ 4 + 2。* (12 + 2 * Q66)。* m ^ 2。* n ^ 2 +的时候* n ^ 4;
Qbar12 = (Q11 + Q22-4 * Q66) * m ^ 2 * n ^ 2 + 12 * (m ^ 4 + n ^ 4);
Qbar16 = (Q11-Q12-2 * Q66) * * n + m ^ 3 (Q12-Q22 + 2 * Q66) * m * n ^ 3;
Qbar22 = Q11。* n ^ 4 + 2。* (12 + 2 * Q66)。* m ^ 2。* n ^ 2 +的时候* m ^ 4;
Qbar26 = (Q11-Q12-2 * Q66) * m * n ^ 3 + (Q12-Q22 + 2 * Q66) * m ^ 3 * n;
Qbar66 = (Q11 + Q22-2 * Q12-2 * Q66) * m ^ 2 * n ^ 2 + Q66 * (m ^ 4 + n ^ 4);
一个= [Qbar11 Qbar12 Qbar16
Qbar12 Qbar22 Qbar26
Qbar16 Qbar26 Qbar66];
=一个* t +一个* t t +一个+一个* *;
B = 1/2 *(一个* (z1 ^ 2-z0 ^ 2) +一个* (z2 ^ 2-z1 ^ 2) +一个* (z3 ^ 2-z2 ^ 2) +一个* (z4 ^ 2-z3 ^ 2));
D = 1/3 *(一个* (z1 ^ 3-z0 ^ 3) +一个* (z2 ^ 3-z1 ^ 3) +一个* (z3 ^ 3-z2 ^ 3) +一个* (z4 ^ 3-z3 ^ 3));
输出= (A, B, B, D) ^ 1。* [N M]
结束
