0

我有一个看起来像这样的方程：

w1   = (1/EI)*(heaviside(x - a1).*(((x - a1).^3)*F1(1))/6 + (C(1, 1)*x.^3)/6 + (C(2, 1)*x.^2)/2 + C(3, 1)*x + C(4, 1)) + (1/GAs).*(- heaviside(x - a1).*(x - a1)*F1(1) - C(1, 1).*x) - C(1, 1)/c_wa;

我使用来自测量的七对值来求解这个方程。我设法用 cftool-toolbox 解决了这个问题。但我并不是真的想知道 Cs，而是 c1、c2、k1 和 k2。我有四个 Cs 方程。现在我想为 c1、c2、k1 和 k2 求解这些方程。到目前为止，我不知道该怎么做，因为我要求解的变量在每个方程中出现多次。

我试图将这四个方程插入我的第一个方程并使用 cftool 来求解它，但这不能正常工作，因为它变得太复杂了。

 C(1, 1) = ((c1*F*(-2*a^3*c2*(EI*(k1 + k2) + k1*k2*l) + 3*a^2*c2*k1*l*(2*EI + k2*l) + 6*a*c2*EI*l*(2*EI + k2*l) - ...
     12*EI^2*(c2*l^2 + k1 + k2) - 4*c2*EI*l^3*(k1 + k2) - c2*k1*k2*l^4 - 12*EI*k1*k2*l))/ ...
     (12*EI^2*(c1 + c2)*(k1 + k2) + 12*c1*c2*EI^2*l^2 + 4*c1*c2*EI*l^3*(k1 + k2) + ...
     12*EI*k1*k2*l*(c1 + c2) + c1*c2*k1*k2*l^4));

C(2, 1) = ((F*k1*(2*EI*l*(a^3*c1*c2 + 6*a*k2*(c1 + c2) - 6*c2*EI) + c2*l^2*(a^3*c1*k2 - 6*a^2*c1*EI - 6*EI*k2) - ...
     2*a*c1*c2*l^3*(a*k2 - 2*EI) - 6*a*EI*(c1 + c2)*(a*k2 - 2*EI) + a*c1*c2*k2*l^4))/ ...
     (12*EI^2*(c1 + c2)*(k1 + k2) + 12*c1*c2*EI^2*l^2 + 4*c1*c2*EI*l^3*(k1 + k2) + ...
     12*EI*k1*k2*l*(c1 + c2) + c1*c2*k1*k2*l^4));

C(3, 1) = ((EI*F*(2*EI*l*(a^3*c1*c2 + 6*a*k2*(c1 + c2) - 6*c2*EI) + c2*l^2*(a^3*c1*k2 - 6*a^2*c1*EI - 6*EI*k2) - ...
     2*a*c1*c2*l^3*(a*k2 - 2*EI) - 6*a*EI*(c1 + c2)*(a*k2 - 2*EI) + a*c1*c2*k2*l^4))/ ...
     (12*EI^2*(c1 + c2)*(k1 + k2) + 12*c1*c2*EI^2*l^2 + 4*c1*c2*EI*l^3*(k1 + k2) + ...
     12*EI*k1*k2*l*(c1 + c2) + c1*c2*k1*k2*l^4));

C(4, 1) = ((EI*F*(2*a^3*c2*(EI*(k1 + k2) + k1*k2*l) - 3*a^2*c2*k1*l*(2*EI + k2*l) - 6*a*c2*EI*l*(2*EI + k2*l) + ...
     12*EI^2*(c2*l^2 + k1 + k2) + 4*EI*l*(c2*l^2*(k1 + k2) + 3*k1*k2) + c2*k1*k2*l^4))/ ...
     (12*EI^2*(c1 + c2)*(k1 + k2) + 12*c1*c2*EI^2*l^2 + 4*c1*c2*EI*l^3*(k1 + k2) + ...
     12*EI*k1*k2*l*(c1 + c2) + c1*c2*k1*k2*l^4));

我将非常感谢一些提示！如果此类问题只能通过 Mathematica 或 Derive 解决，我也很乐意获得有关这些问题的解决方案！

编辑（给比尔 8 月 7 日 23:08）：

感谢您的回复！

  • “...”表示换行
  • EI, F1(1), GAs, F, c_wa 给定变量（你在哪里看到 F？）

所以：

C(1, 1) = ((c1*F*(-2*a^3*c2*(EI*(k1 + k2) + k1*k2*l) + 3*a^2*c2*k1*l*(2*EI + k2*l) + 6*a*c2*EI*l*(2*EI + k2*l) - 12*EI^2*(c2*l^2 + k1 + k2) - 4*c2*EI*l^3*(k1 + k2) - c2*k1*k2*l^4 - 12*EI*k1*k2*l))/ (12*EI^2*(c1 + c2)*(k1 + k2) + 12*c1*c2*EI^2*l^2 + 4*c1*c2*EI*l^3*(k1 + k2) + 12*EI*k1*k2*l*(c1 + c2) + c1*c2*k1*k2*l^4));

C(2, 1) = ((F*k1*(2*EI*l*(a^3*c1*c2 + 6*a*k2*(c1 + c2) - 6*c2*EI) + c2*l^2*(a^3*c1*k2 - 6*a^2*c1*EI - 6*EI*k2) - ...
     2*a*c1*c2*l^3*(a*k2 - 2*EI) - 6*a*EI*(c1 + c2)*(a*k2 - 2*EI) + a*c1*c2*k2*l^4))/ (12*EI^2*(c1 + c2)*(k1 + k2) + 12*c1*c2*EI^2*l^2 + 4*c1*c2*EI*l^3*(k1 + k2) + 12*EI*k1*k2*l*(c1 + c2) + c1*c2*k1*k2*l^4));

C(3, 1) = ((EI*F*(2*EI*l*(a^3*c1*c2 + 6*a*k2*(c1 + c2) - 6*c2*EI) + c2*l^2*(a^3*c1*k2 - 6*a^2*c1*EI - 6*EI*k2) - 2*a*c1*c2*l^3*(a*k2 - 2*EI) - 6*a*EI*(c1 + c2)*(a*k2 - 2*EI) + a*c1*c2*k2*l^4))/ (12*EI^2*(c1 + c2)*(k1 + k2) + 12*c1*c2*EI^2*l^2 + 4*c1*c2*EI*l^3*(k1 + k2) + 12*EI*k1*k2*l*(c1 + c2) + c1*c2*k1*k2*l^4));

C(4, 1) = ((EI*F*(2*a^3*c2*(EI*(k1 + k2) + k1*k2*l) - 3*a^2*c2*k1*l*(2*EI + k2*l) - 6*a*c2*EI*l*(2*EI + k2*l) + 12*EI^2*(c2*l^2 + k1 + k2) + 4*EI*l*(c2*l^2*(k1 + k2) + 3*k1*k2) + c2*k1*k2*l^4))/ (12*EI^2*(c1 + c2)*(k1 + k2) + 12*c1*c2*EI^2*l^2 + 4*c1*c2*EI*l^3*(k1 + k2) + 12*EI*k1*k2*l*(c1 + c2) + c1*c2*k1*k2*l^4));

即使您不需要他们来解决问题以便理解：

EI = 8.0934e+11
F1 = [10000 345] 
=> F = F1(1) = 10000 
=> a = F1(2) = 345
GAs = 7.0852e+08 
c_wa = 2.115e+05;

亲切的问候，

AMK

4

0 回答 0