| 查看: 801 | 回复: 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}] |
回复此楼» 猜你喜欢
有院领导为了换新车，用横向课题经费买了俩车
已经有9人回复
-
酰胺脱乙酰基
已经有13人回复
-
博士延得我，科研能力直往上蹿
已经有8人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有4人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有5人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有6人回复
-
面上项目申报
已经有3人回复
-
遇见不省心的家人很难过
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
wangzgenyu
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 581
- 帖子: 178
- 在线: 20小时
- 虫号: 1606005
- 注册: 2012-02-08
- 性别: GG
- 专业: 控制论中的数学方法
2楼2016-05-18 16:18:57
wangzgenyu
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 581
- 帖子: 178
- 在线: 20小时
- 虫号: 1606005
- 注册: 2012-02-08
- 性别: GG
- 专业: 控制论中的数学方法
3楼2016-05-18 16:29:00
kill_speed
铁虫 (正式写手)
- 应助: 11 (小学生)
- 金币: 306.8
- 红花: 1
- 帖子: 343
- 在线: 29.4小时
- 虫号: 4280372
- 注册: 2015-12-10
- 性别: GG
- 专业: 信号理论与信号处理
4楼2016-05-19 12:36:43