【求助】如何求解这样的方程组
mathematica吧
全部回复
仅看楼主
level 11
代码:
Solve[{Integrate[1/
Sqrt[(2 P)/(
E0 Iy) (Cos[\[Theta]A] -
Cos[\[Theta]])], {\[Theta], \[Theta]A, \[Theta]B}] == L/2(*AB*),
Integrate[-(1/
Sqrt[(2 P)/(E0 Iy) (1 - Cos[\[Theta]]) +
d\[Theta]C^2]), {\[Theta], \[Theta]B, 0}] == \[Pi] r1/2(*BC*),
Sqrt[(2 P)/(E0 Iy) (Cos[\[Theta]A] - Cos[\[Theta]B])] +
Sqrt[(2 P)/(E0 Iy) (1 - Cos[\[Theta]B]) + d\[Theta]C^2] ==
1/r1(*B*),
Integrate[Cos[\[Theta]]/
Sqrt[(2 P)/(
E0 Iy) (Cos[\[Theta]A] -
Cos[\[Theta]])], {\[Theta], \[Theta]A, \[Theta]B}] +
Integrate[-(Cos[\[Theta]]/
Sqrt[(2 P)/(E0 Iy) (1 - Cos[\[Theta]]) +
d\[Theta]C^2]), {\[Theta], \[Theta]B, 0}] ==
1/2 \[Epsilon] Lb}, {P, \[Theta]A, \[Theta]B, d\[Theta]C}]
2019年04月21日 06点04分 1
level 11
附上常量的值:
R = 3; r = 0.5; r1 = 0.075; r2 = 0.045; E0 = 2500; L = 0.09; Lb = \
0.43; Iy = 3.125`*^-10; \[Epsilon] = 0.00031288629055851536`;
2019年04月21日 06点04分 2
吧务
level 15
1. 先设法弄清楚你的积分能不能积出来,积不出来就别往Solve里塞了,因为Solve不能处理积分方程。
2. 就算积出来了,如果是超越方程也请做好算不出来的准备,因为只有很少的一部分超越方程有符号解。
3. 如果积不出来则需要考虑数值解,FindRoot慢慢试吧。(找合适初值这步目测又是够呛。)
2019年05月04日 07点05分 3
1