level 3
朱老四__
楼主
代码如下,是一个力平衡和力矩平衡的方程组。 力和力矩都是迎角的函数
clWing[a_] := 0.25704 + 0.11369*a + 0.00106*a*a (*机翼升力系数,a为迎角*)
cmWing[a_] := -0.04954 - 0.00124*a - (1.6305*10^-4)*a*a(*机翼升力力矩系数,a为迎角*)
clTail[b_] := 0.25704 + 0.11269*b*-+0.00106*b*b(*尾翼升力系数,b为尾翼迎角*)
cmTail[b_] := (5.037*10^-5) + 0.00215*b(*尾翼翼升力力矩系数,b为迎角*)
p := 0.5*1.225*30*30(*气动系数,动压*)
wingLift[a_] := 8.7*p*clWing[a] (*机翼升力,a为迎角*)
tailLift[b_] := 1.5*p*clTail[b] (*尾翼升力,a为迎角*)
mwing[a_] := p*cmWing[a]*8.7 (*机翼升力力矩 ,a为迎角*)
mTail[b_] := p*cmTail[b]*1.5 (*尾翼升力力矩 ,b为迎角*)
In[5]:= SolveValues[
wingLift[a] + tailLift[b] == 9.8*400 &&
mwing[a] + mTail[b] + 0.05*Winglift[a] - 2.5*tailLift[b] == 0, {a,
b}, Reals]
(*力相等,力矩相等*)
\:6B63\:5728\:8BA1\:7B97In[5]:= SolveValues::inex: SolveValues 无法求解具有不精确系数的系统或通过将系统中出现的不精确数直接合理化处理得到的系统. 由于 Solve 所用的许多方法需要精确输入,为 Solve 提供一个精确版本的系统可能会有帮助.
Out[5]= SolveValues[
tailLift[b] + wingLift[a] == 3920. &&
mTail[b] + mwing[a] - 2.5 tailLift[b] + 0.05 Winglift[a] == 0, {a,
b}, Reals]
2022年10月22日 02点10分
1
clWing[a_] := 0.25704 + 0.11369*a + 0.00106*a*a (*机翼升力系数,a为迎角*)
cmWing[a_] := -0.04954 - 0.00124*a - (1.6305*10^-4)*a*a(*机翼升力力矩系数,a为迎角*)
clTail[b_] := 0.25704 + 0.11269*b*-+0.00106*b*b(*尾翼升力系数,b为尾翼迎角*)
cmTail[b_] := (5.037*10^-5) + 0.00215*b(*尾翼翼升力力矩系数,b为迎角*)
p := 0.5*1.225*30*30(*气动系数,动压*)
wingLift[a_] := 8.7*p*clWing[a] (*机翼升力,a为迎角*)
tailLift[b_] := 1.5*p*clTail[b] (*尾翼升力,a为迎角*)
mwing[a_] := p*cmWing[a]*8.7 (*机翼升力力矩 ,a为迎角*)
mTail[b_] := p*cmTail[b]*1.5 (*尾翼升力力矩 ,b为迎角*)
In[5]:= SolveValues[
wingLift[a] + tailLift[b] == 9.8*400 &&
mwing[a] + mTail[b] + 0.05*Winglift[a] - 2.5*tailLift[b] == 0, {a,
b}, Reals]
(*力相等,力矩相等*)
\:6B63\:5728\:8BA1\:7B97In[5]:= SolveValues::inex: SolveValues 无法求解具有不精确系数的系统或通过将系统中出现的不精确数直接合理化处理得到的系统. 由于 Solve 所用的许多方法需要精确输入,为 Solve 提供一个精确版本的系统可能会有帮助.
Out[5]= SolveValues[
tailLift[b] + wingLift[a] == 3920. &&
mTail[b] + mwing[a] - 2.5 tailLift[b] + 0.05 Winglift[a] == 0, {a,
b}, Reals]