我试图从为每个力(V
)生成的表中提取F
的最低值,以便建立最佳结果,例如:对于F=5000
,V=0.06575
是最佳结果
我希望能够为F
的每个实例显示此内容。我不太确定我哪里出错了。
%===========================================
%This Script calculates
%===========================================
%Parameter Values
%F:force (N)
%L1: Length of beam 1 (m)
%L2: Length of beam 2 (m)
%H: Height of beam 1 (m)
%h: Height of beam 2 (m)
%b: Cross section width (m)
%S: Yield strength (Pa)
%N: Factor of safety
%V:Volume of the beam (m^2)
%=========================================
clc %clear command window
clear all %clear all variables previously generated and stored in Matlab
format short g
L1=2;
L2=3;
b=0.05;
S=248;
K=1.42;
results=[];
optimalresults=[];
nd=3.6; %desired factor of safety calculated by group number
for F=5000:100:6000
for h=0.2:0.005:0.3 %loop of h
for H=0.2:0.005:0.3 % loop of H
M1=F*(L1+L2);
M2=F*L2;
I1=b*H^3/12;
I2=b*h^3/12;
Y1=H/2;
Y2=h/2;
Max_stress_1=M1*Y1/I1;
Max_stress_2=K*(M2*Y2/I2);
Max_total_stress=max(Max_stress_1,Max_stress_2);
Max_total_stress=Max_total_stress/1e6;
N=0.577*S/Max_total_stress;
%if N is greater than desired factor of safety then display results
if N >= nd
V=(b*L1*H)+(b*L2*h);
results=[results; F,V,N];sort(F);
end
end
end
results_sorted=sortrows(results,1); %sorts factor of safety from min to max
disp(' F V N');
disp(results_sorted);
optimum=results_sorted(1,:)
end
我需要更改哪些内容才能显示每个V
F
的最小值?
答案 0 :(得分:0)
我的猜测是你可能希望这个位在3 for
个循环之外:
results_sorted=sortrows(results,1); %sorts factor of safety from min to max
disp(' F V N');
disp(results_sorted);
optimum=results_sorted(1,:)
另外,我不确定你在循环中使用sort(F);
取得了什么。
答案 1 :(得分:0)
我认为您的代码需要进行两项更改才能获得每个F的最小V值:
在使用sortrows的地方,您应该使用第二列(V)对表进行排序,以便将最低的V放置在顶行。 results_sorted=sortrows(results,2);
您应该只对当前F的结果进行排序,因此要么清除变量"结果"对于每个F迭代results=[];
或您使用新变量。一定要检查"结果"在排序行之前不是空的。
变量"最佳"现在每个F迭代的V值最低。