| 查看: 785 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
wangzgenyu新虫 (小有名气)
|
[求助]
mathematica不能化简得问题 已有1人参与
|
|
|
程序编了有一段时间了，只是一个简单的二元二次方程组求解问题，但就是不能很好的化简，现在代码拷贝在下面，希望有人能指点一些。另外将代码的规则写出以后，还带有TextCell[""],不知道是为什么。 Clear["Global`*"] Di = 1.98; tc = 0.01; Ri = 1.99; hs = 0.01; ts = 0.01; p = 1.6 10^6; nu = 0.3; r = 0.5*(Di + tc); R = Ri + 0.5*hs; nfy = ArcSin[r/R]; kc = Power[3 (1 - nu^2), (4)^-1]/Sqrt[r*tc]; nlam = Power[3 (1 - nu^2), (4)^-1] Sqrt[R/hs]; k21 = 1 - (1 - 2 nu)/(2 Power[3 (1 - nu^2), (4)^-1]) Sqrt[hs/R] Cot[nfy]; k22 = 1 - (1 + 2 nu)/(2 Power[3 (1 - nu^2), (4)^-1]) Sqrt[hs/R] Cot[nfy]; k2d = 1 - nu/Power[3 (1 - nu^2), (4)^-1] Sqrt[hs/R] Cot[nfy] + (1 - 2 nu)^2/( 8 Sqrt[3 (1 - nu^2)]) hs/R (Cot[nfy])^2; k23 = 1 + (2 - nu)/(2 nu Power[3 (1 - nu^2), (4)^-1]) hs/R Cot[nfy]; a11 = Sqrt[r/tc]*(1 + Sqrt[Sin[nfy]]*k2d); a12 = (1 - tc^2/(k21*hs^2)) ; a21 = (1 - tc^2/(k21*hs^2)); a22 = 2 ((3 (1 - nu^2))^0.75)*Sqrt[tc/ r] (1 + 1/(k21 Sqrt[Sin[nfy]]) (tc/hs)^2.5); a10 = 1/4 ((1 - nu)/Sin[nfy] tc/hs - 2 Power[3 (1 - nu^2), (4)^-1] Sqrt[R/hs] tc/hs k2d* Cos[nfy] - (2 - nu)); a20 = (p*r)/(2.0*k21)*Cot[nfy]; Solve[{a11 x + a12 y == a10, a21 x + a22 y == a20}, {x, y}] |
回复此楼» 猜你喜欢
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有3人回复
-
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
wangzgenyu
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 581
- 帖子: 178
- 在线: 20小时
- 虫号: 1606005
- 注册: 2012-02-08
- 性别: GG
- 专业: 控制论中的数学方法
3楼2016-05-18 16:29:00
wangzgenyu
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 581
- 帖子: 178
- 在线: 20小时
- 虫号: 1606005
- 注册: 2012-02-08
- 性别: GG
- 专业: 控制论中的数学方法
2楼2016-05-18 16:18:57
kill_speed
铁虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 306.8
- 红花: 1
- 帖子: 343
- 在线: 29.4小时
- 虫号: 4280372
- 注册: 2015-12-10
- 性别: GG
- 专业: 信号理论与信号处理
4楼2016-05-19 12:36:43