MATLAB的答案

帮助优化代码回到1号线如果条件不满足。

1视图(30天)
詹姆斯·库克
詹姆斯·库克 2023年4月4日
回答: 沃尔特·罗伯森 2023年4月4日
我下面的代码需要迭代如果不满足一个条件。
如果任何价值1 x3x8 Tsai_Wu_Total矩阵>我需要迭代LT(代码)的顶部和重新运行整个代码iterively直到Tsai_Wu_Total矩阵< 1。
所以只有改变在Tsai_Wu_Total LT,直到所有值< 1。完整代码如下。谢谢。
清晰的所有,关闭所有,clc
格式短g
% %用户输入
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%的厚度是首先列出所需的方向
%层压板- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
LT = 0.180;%板厚度
TLT = LT * 8;%层厚度
层压板= [TLT / 8 TLT / 8 TLT / 8 TLT / 8 TLT / 8 TLT / 8 TLT / 8 TLT / 8;
45 0 45 -45 90 90 -45 0);% 8层厚度相等
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% %的距离计算- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
电脑=大小(层压板,2);%厚度计算
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%材料特性- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
核心= 0;% < = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =输入核心厚度
E_1 = 3 * 10 ^ 7;
E_2 = 7.5 * 10 ^ 5;
G_12 = 3.75 * 10 ^ 5;
v_12 = 0.25;
v_21 = v_12 * (E_2 / E_1);
XT = 1.5 * 10 ^ 5;
XC = -1.0 * 10 ^ 5;
欧美= 6.0 * 10 ^ 3;
YC = -1.7 * 10 ^ 4;
S = 1.0 * 10 ^ 4;
SIG_B = 0;
Density_Laminate = 0.0625;
Mass_Laminate = Density_Laminate *总和(强化(1:));
Density_core = 0;
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%热,湿热输入- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
a_1 = -1.7 * 10 ^ 7;
a₂= 1.6 * 10 ^ 5;
b_1 = 0;
b_2 = 0;
Dt = -200;
Dc = 0;
为i = 1:电脑
a_x (i) = a_1 * cosd(层压板(我)2)^ 2 + a₂*信德(层压板(我)2)^ 2;
得+ (i) = a_1 *信德(层压板(我)2)^ 2 + a₂* cosd(层压板(我)2)^ 2;
a_xy (i) = 2 * cosd(层压板(我))*信德(强化(2,i)) * (a_1-a_2);
提出(i) = b_1 * cosd(层压板(我)2)^ 2 + b_2 *信德(层压板(我)2)^ 2;
的话(我)= b_1 *信德(层压板(我)2)^ 2 + b_2 * cosd(层压板(我)2)^ 2;
b_xy (i) = 2 * cosd(层压板(我))*信德(强化(2,i)) * (b_1-b_2);
α(:,i) = (a_x(我),得+(我),a_xy (i)];
结束
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%力量和时间- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
F_M = [2000、3000、0 0 5000 0];
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% % Q矩阵计算
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Q11 = E_1 / (1-v_12 * v_21);
如= E_2 / (1-v_12 * v_21);
12 = (E_1 * v_21) / (1-v_12 * v_21);
温度系数= (E_1 * v_21) / (1-v_12 * v_21);
Q66 = G_12;
Q_M = [Q11 Q12 0;如对方篮里0;0 0 Q66];
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%这个For循环计算层的总高度- - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
t (1,1) = 0;
为i = 1:电脑
t = t +层压板(1,);
sum_t (i + 1) = t;
结束
为k = 1:尺寸(sum_t, 2)
如果sum_t (k) = = 1000
sum_t (k) = [];
其他的sum_t (k) = sum_t (k);
结束
结束
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%这For循环索引mid-surf“0”需要- - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
t_2 = t / 2;
为j = 1: PC + 1
如果sum_t (j) < t_2
Zero_check (j) = 1;
其他的
Zero_check (j) = 0;
结束
Top_Plies =总和(Zero_check);
结束
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%这个For循环检查厚度层之间“边界——谎言
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
为m = 1: PC + 1
如果sum_t(1米)= = t_2
Layer_check(1米)= 1;
其他的
Layer_check(1米)= 0;
结束
Ply_layer =总和(Layer_check);
结束
如果核心= = 0
Mid_Surf =“层间”;
其他的
Mid_Surf =“中间的核心”;
结束
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%这一步计算层的距离mid-surf下降
%的厚度
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
如果Mid_Surf = =“中间的核心”
Ply_D = 0 (1、PC + 1);
sum_t (1、PC + 2) = 1000;
为k = 1: PC + 2
p = k - 1;
如果sum_t (k) < t_2
Ply_D (k) = sum_t (k) t_2;
其他的
如果k = = Top_Plies + 1%这个语句添加一个“0”
% mid-surf
Ply_D (k) = 0;
其他的
Ply_D (k) = sum_t (1, p) t_2;
结束
结束
结束
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%这一步计算层的距离mid-surf下降
%的两层之间
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
其他的
为n = 1: PC + 1
如果sum_t < t_2
Ply_D (1, n) = t_2-sum_t (1, n);
其他的
Ply_D (1, n) = sum_t t_2 (1, n);
结束
结束
结束
如果核心= = 0
Ply_D = Ply_D;
其他的
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%计算距离复合层的核心结构- - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
为i = 1: PC + 2
如果Ply_D (i) < 0
Ply_D1 (i) = Ply_D(我)生水/ 2;
elseifPly_D (i) = = 0
Ply_D1 (i) = Ply_D(我)+核心/ 2;
elseifPly_D (i) > 0
Ply_D1 (i) = Ply_D(我)+核心/ 2;
结束
结束
Ply_D = [Ply_D1(1, - 1:电脑/ 2)0 Ply_D1 (1、PC / 2 + 3:结束))
结束
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%的变换矩阵- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
R = [1 0 0
0 1 0
0 0 2];
为j = 1:电脑
c (j) = cosd(层压板(j));
s (j) =信德(层压板(j));
结束
为i = 1:电脑
T (:,:, i) = (c (i) ^ 2 s (i) ^ 2 c(我)* 2 * s (i); s(我)^ 2 c (i) ^ 2 c(我)* 2 * s(我);…
c - s(我)*(我)(我)* (i) c (i) ^ 2 s (i) ^ 2);
结束
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% Q_BAR矩阵- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
为z = 1:电脑
Q_BAR (:,:, z) =发票(T (:,:, z)) * Q_M * R * T (:,:, z) *发票(R);
结束
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%一个矩阵- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
一个= 0 (3、3);
为k = 1:电脑
A_1 = Q_BAR (1 1 k) * (Ply_D (k + 1) -Ply_D (k));
(1)= (1)+ A_1;
a₂= Q_BAR (1、2、k) * (Ply_D (k + 1) -Ply_D (k));
(1、2)= (1、2)+ A₂;
A_3 = Q_BAR (1、3 k) * (Ply_D (k + 1) -Ply_D (k));
(1、3)= (1、3)+ A_3;
A_4 = Q_BAR (2, 1, k) * (Ply_D (k + 1) -Ply_D (k));
(2,1)= (2,1)+ A_4;
A_5 = Q_BAR (2 2 k) * (Ply_D (k + 1) -Ply_D (k));
(2,2)= (2,2)+ A_5;
A_6 = Q_BAR (2 3 k) * (Ply_D (k + 1) -Ply_D (k));
(2、3)= (2、3)+ A_6;
A_7 = Q_BAR (3 1 k) * (Ply_D (k + 1) -Ply_D (k));
一个(3,1)= (1)+ A_7;
A_8 = Q_BAR (3 2 k) * (Ply_D (k + 1) -Ply_D (k));
(2)= (2)+ A_8;
A_9 = Q_BAR (3 3 k) * (Ply_D (k + 1) -Ply_D (k));
(3)= (3、3)+ A_9;
结束
B = 0 (3、3);
为k = 1:电脑
B_1 = Q_BAR (1 1 k) * ((Ply_D (k + 1) ^ 2-Ply_D (k) ^ 2) / 2);
B (1, - 1) = (1) + B_1;
B_2 = Q_BAR (1、2、k) * ((Ply_D (k + 1) ^ 2-Ply_D (k) ^ 2) / 2);
(1、2)= B (1、2) + B_2;
B_3 = Q_BAR (1、3 k) * ((Ply_D (k + 1) ^ 2-Ply_D (k) ^ 2) / 2);
(1、3)= B (1、3) + B_3;
B_4 = Q_BAR (2, 1, k) * ((Ply_D (k + 1) ^ 2-Ply_D (k) ^ 2) / 2);
(2,1)= B (2, 1) + B_4;
B_5 = Q_BAR (2 2 k) * ((Ply_D (k + 1) ^ 2-Ply_D (k) ^ 2) / 2);
(2,2)= B (2, 2) + B_5;
B_6 = Q_BAR (2 3 k) * ((Ply_D (k + 1) ^ 2-Ply_D (k) ^ 2) / 2);
(2、3)= B (2、3) + B_6;
B_7 = Q_BAR (3 1 k) * ((Ply_D (k + 1) ^ 2-Ply_D (k) ^ 2) / 2);
B (3,1) = (1) + B_7;
B_8 = Q_BAR (3 2 k) * ((Ply_D (k + 1) ^ 2-Ply_D (k) ^ 2) / 2);
(2)= B (3 2) + B_8;
B_9 = Q_BAR (3 3 k) * ((Ply_D (k + 1) ^ 2-Ply_D (k) ^ 2) / 2);
(3,3)= B (3,3) + B_9;
结束
D = 0 (3,3);
为k = 1:电脑
D_1 = Q_BAR (1 1 k) * ((Ply_D (k + 1) ^ 3-Ply_D (k) ^ 3) / 3);
(1)= D (1,1) + D_1;
D_2 = Q_BAR (1、2、k) * ((Ply_D (k + 1) ^ 3-Ply_D (k) ^ 3) / 3);
(1、2)= D (1、2) + D_2;
D_3 = Q_BAR (1、3 k) * ((Ply_D (k + 1) ^ 3-Ply_D (k) ^ 3) / 3);
(1、3)= D (1、3) + D_3;
D_4 = Q_BAR (2, 1, k) * ((Ply_D (k + 1) ^ 3-Ply_D (k) ^ 3) / 3);
D (2, 1) = D (2, 1) + D_4;
D_5 = Q_BAR (2 2 k) * ((Ply_D (k + 1) ^ 3-Ply_D (k) ^ 3) / 3);
D (2, 2) = D (2, 2) + D_5;
D_6 = Q_BAR (2 3 k) * ((Ply_D (k + 1) ^ 3-Ply_D (k) ^ 3) / 3);
(2、3)= D (2、3) + D_6;
D_7 = Q_BAR (3 1 k) * ((Ply_D (k + 1) ^ 3-Ply_D (k) ^ 3) / 3);
D (3,1) = D (3,1) + D_7;
D_8 = Q_BAR (3 2 k) * ((Ply_D (k + 1) ^ 3-Ply_D (k) ^ 3) / 3);
(2)= D (3 2) + D_8;
D_9 = Q_BAR (3 3 k) * ((Ply_D (k + 1) ^ 3-Ply_D (k) ^ 3) / 3);
(3)= D (3、3) + D_9;
结束
ABD = (A, B, B, D);
ABD_I =发票(ABD);
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%热力和时刻- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
为i = 1:电脑
N_xT (i) = (Q_BAR (1, 1, 1) * a_x(我)+ Q_BAR(1、2) *得+(我)+ Q_BAR(1, 3,我)* a_xy(我))*…
((Ply_D -Ply_D (i + 1)(我)));
N_yT (i) = (Q_BAR (1、2) * a_x(我)+ Q_BAR(2, 2,我)*得+(我)+ Q_BAR(2、3、我)* a_xy(我))*…
((Ply_D -Ply_D (i + 1)(我)));
N_xy (i) = (Q_BAR(1, 3,我)* a_x(我)+ Q_BAR(2、3、我)*得+(我)+ Q_BAR (3 3 i) * a_xy(我))*…
((Ply_D -Ply_D (i + 1)(我)));
M_xT (i) = (Q_BAR (1, 1, 1) * a_x(我)+ Q_BAR(1、2) *得+(我)+ Q_BAR(1, 3,我)* a_xy(我))*…
((Ply_D (i + 1) ^ 2-Ply_D(我)^ 2));
M_yT (i) = (Q_BAR (1、2) * a_x(我)+ Q_BAR(2, 2,我)*得+(我)+ Q_BAR(2、3、我)* a_xy(我))*…
((Ply_D (i + 1) ^ 2-Ply_D(我)^ 2));
M_xyT (i) = (Q_BAR(1, 3,我)* a_x(我)+ Q_BAR(2、3、我)*得+(我)+ Q_BAR (3 3 i) * a_xy(我))*…
((Ply_D (i + 1) ^ 2-Ply_D(我)^ 2));
结束
N_xxT = Dt *总和(N_xT);
N_yyT = Dt *总和(N_yT);
N_xyT = Dt *总和(N_xy);
M_xxT = Dt / 2 *总和(M_xT);
M_yyT = Dt / 2 *总和(M_yT);
M_xyT = Dt / 2 *总和(M_xyT);
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%湿热力和时刻- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
为i = 1:电脑
N_xH (i) = (Q_BAR提出(1,1,1)*(我)+ Q_BAR的话(1、2)*(我)+ Q_BAR(1, 3,我)* b_xy(我))*…
((Ply_D -Ply_D (i + 1)(我)));
N_yH (i) = (Q_BAR提出(1、2)*(我)+ Q_BAR的话(2,2,我)*(我)+ Q_BAR(2、3、我)* b_xy(我))*…
((Ply_D -Ply_D (i + 1)(我)));
N_xyH (i) = (Q_BAR提出(1,3,我)*(我)+ Q_BAR的话(2、3、我)*(我)+ Q_BAR (3 3 i) * b_xy(我))*…
((Ply_D -Ply_D (i + 1)(我)));
M_xH (i) = (Q_BAR提出(1,1,1)*(我)+ Q_BAR的话(1、2)*(我)+ Q_BAR(1, 3,我)* b_xy(我))*…
((Ply_D (i + 1) ^ 2-Ply_D(我)^ 2));
M_yH (i) = (Q_BAR提出(1、2)*(我)+ Q_BAR的话(2,2,我)*(我)+ Q_BAR(2、3、我)* b_xy(我))*…
((Ply_D (i + 1) ^ 2-Ply_D(我)^ 2));
M_xyH (i) = (Q_BAR提出(1,3,我)*(我)+ Q_BAR的话(2、3、我)*(我)+ Q_BAR (3 3 i) * b_xy(我))*…
((Ply_D (i + 1) ^ 2-Ply_D(我)^ 2));
结束
N_xxH =直流*总和(N_xH);
N_yyH =直流*总和(N_yH);
N_xyH =直流*总和(N_xyH);
M_xxH =直流/ 2 *总和(M_xH);
M_yyH =直流/ 2 *总和(M_yH);
M_xyH =直流/ 2 *总和(M_xyH);
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%的总力量和时间- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
N_xxTotal = F_M (1,1) + N_xxT + N_xxH;
N_yyTotal = F_M (1、2) + N_yyT + N_yyH;
N_xyTotal = F_M (1、3) + N_xyT + N_xyH;
M_xxTotal = F_M (1、4) + M_xxT + M_xxH;
M_yyTotal = F_M (1、5) + M_yyT + M_yyH;
M_xyTotal = F_M (1,6) + M_xyT + M_xyH;
F_M_T = [N_xxTotal; N_yyTotal; N_xyTotal M_xxTotal; M_yyTotal; M_xyTotal);
% %
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%解决强化应变和曲率- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
Ep_k = ABD_I * F_M_T;
% %
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%解决个人的全球和当地的应变和应力- - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
为n = 1: PC + 1
如果sum_t < t_2
Ply_Dis (1, n) = t_2-sum_t (1, n);
其他的
Ply_Dis (1, n) = sum_t t_2 (1, n);
结束
结束
如果核心= = 0
Ply_Dis = Ply_D;
其他的
为i = 1: PC + 1
如果Ply_Dis (i) < 0
Ply_Dis (i) = Ply_Dis(我)生水/ 2;
elseifPly_Dis (i) > 0
Ply_Dis (i) = Ply_Dis(我)+核心/ 2;
结束
结束
结束
为i = 1:长度(Ply_Dis) 1
P(:,我)= (Ply_Dis(我),(Ply_Dis(我)+ Ply_Dis (i + 1) / 2, Ply_Dis (i + 1)];
结束
为i = 1:长度(P)
Global_Top_Strains_Ply (:, i) = (Ep_k (1:3) + P(我)。* Ep_k (4:6));
Global_Top_Strains_Ply(:,我)= Global_Top_Strains_Ply(:,我)——(Dt *α(:,i));
Global_Middle_Strains_Ply (:, i) = (Ep_k (1:3) + P(我)。* Ep_k (4:6));
Global_Middle_Strains_Ply(:,我)= Global_Middle_Strains_Ply(:,我)——(Dt *α(:,i));
Global_Bottom_Strains_Ply (:, i) = (Ep_k (1:3) + P(我)。* Ep_k (4:6));
Global_Bottom_Strains_Ply(:,我)= Global_Bottom_Strains_Ply(:,我)——(Dt *α(:,i));
Total_Global_Strains_Ply(:,:我)= (Global_Top_Strains_Ply(:,我)…
Global_Middle_Strains_Ply(:,我),Global_Bottom_Strains_Ply (:, i)];
Total_Local_Strains (:,:, i) = T(:,:我)* Total_Global_Strains_Ply(:,:我);
Total_Global_Stresses(:,:我)= Q_BAR (:,:, i) *…
Total_Global_Strains_Ply(:,:我);
Total_Local_Stresses (:,:, i) = T(:,:我)* Total_Global_Stresses(:,:我);
结束
% %
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%解决如果失败发生在任何位置厚度- - - - - - - - - - - - - - - - - - - - - - -
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
F1 = (1 / XT) + (1 / XC);
季= 1 / (XT * XC);
F2 =(1 /次)+ (1 / YC);
——F22 = 1 /(欧美* YC);
F6 = 0;
F66 = 1 / S ^ 2;
如果SIG_B = = 0;
F12 = 0;
其他的
F12 = 1 / (2 * SIG_B ^ 2) * (1 - ((1 / XT) + (1 / XC) +(1 /次)+ (1 / YC)) * SIG_B + ((1 / XT * XC) +…
(1 /次* YC) * SIG_B ^ 2));
结束
为i = 1:长度(Total_Local_Stresses)
Tsai_Wu_Top (:, i) = F1 * Total_Local_Stresses (1, 1, 1) + F2 *…
Total_Local_Stresses (2, 1, i) + F11 * Total_Local_Stresses(1, - 1,我)^ 2 +…
——F22 * Total_Local_Stresses(2, 1,我)^ 2 + F66 * Total_Local_Stresses (3、1, i) ^ 2…
+ 2 * F12 * Total_Local_Stresses (1, 1, 1) * Total_Local_Stresses(2, 1,我);
Tsai_Wu_Middle (:, i) = F1 * Total_Local_Stresses (1、2) + F2 *…
Total_Local_Stresses (2, 2, i) + F11 * Total_Local_Stresses (1、2) ^ 2 +…
——F22 * Total_Local_Stresses(2, 2,我)^ 2 + F66 * Total_Local_Stresses (3 2 i) ^ 2…
+ 2 * F12 * Total_Local_Stresses (1、2) * Total_Local_Stresses(2, 2,我);
Tsai_Wu_Bottom (:, i) = F1 * Total_Local_Stresses(1、3、我)+ F2 *…
Total_Local_Stresses (2、3、i) + F11 * Total_Local_Stresses(1、3、我)^ 2 +…
——F22 * Total_Local_Stresses(2、3、我)^ 2 + F66 * Total_Local_Stresses (3 3 i) ^ 2…
+ 2 * F12 * Total_Local_Stresses(1, 3,我)* Total_Local_Stresses(2、3、我);
Tsai_Wu_Total(:,:我)= (Tsai_Wu_Top(:,我),Tsai_Wu_Middle(:,我)…
Tsai_Wu_Bottom (:, i)];
结束

在回答这个问题。

答案(1)

沃尔特·罗伯森
沃尔特·罗伯森 2023年4月4日
有一些您可以使用不同的逻辑结构。
  • 你可以开始而真正的。在某种程度上测试是否都在范围的变量;如果是这样,那么打破循环
  • 你可以开始而真正的。在某种程度上测试是否变量的范围;如果是这样,那么继续回到顶部的循环
  • 你可以初始化变量的宽容,然后使用而的条件,测试他们的宽容。不需要测试循环内的条件。
  • 你可以设置一个标志等坚持= true;和而坚持。在循环中,测试的变量,如果他们都在范围内,集坚持= false;。你不需要打破或继续如果你使用这个

类别

找到更多的在日期和时间在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!