| 查看: 954 | 回复: 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}] |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
面上合作单位盖章
已经有5人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
申请删除本帖
已经有7人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
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
投票:
10