24小时热门版块排行榜    

查看: 850  |  回复: 7

bingmei527

铁虫 (初入文坛)

[求助] 偏微分求解已有2人参与

问题在附件图片上,希望哪位高手帮帮忙!

偏微分求解
IMG_20140414_230851.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

巅子~

金虫 (小有名气)

癫子

【答案】应助回帖

感谢参与,应助指数 +1
做出来了不知道怎么回复给你

[ 发自小木虫客户端 ]
2楼2014-04-15 00:48:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

感谢参与,应助指数 +1
这应该是弹性力学中的无限长圆柱杆件的轴对称问题的力的平衡方程吧。由广义胡可定律,εθ=1/E*(σθ-ν*σr)
依题意εθ=k*σr*(1+ν)/E*r^2/R^2
由此得到σθ=σr*[ν+k*(1+ν)*(r/R)^2]                           (1)
代入力的平衡方程中得到:
Pσr/Pr+[1-ν-k*(1+ν)*(r/R)^2]/r*σr+[ν+k*(1+ν)*(r/R)^2]/r*Pσr/Pθ=0
令σr=f(r)*g(θ)                                                           (2)
g(θ)*df(r)/dr+[1-ν-k*(1+ν)*(r/R)^2]/r*f(r)*g(θ)+[ν+k*(1+ν)*(r/R)^2]/r*f(r)*dg(θ)/dθ=0
两边同除以 f(r)*g(θ),得到 :
{[df(r)/dr]/f(r)+[1-ν-k*(1+ν)*(r/R)^2]/r}/{[ν+k*(1+ν)*(r/R)^2]/r}
= - [dg(θ)/dθ]/g(θ)
由于左边为r的函数,右边为θ的函数,两者相等说明均为某一个常数A。由此得到:
- [dg(θ)/dθ]/g(θ)=A                                        (3)
{[df(r)/dr]/f(r)+[1-ν-k*(1+ν)*(r/R)^2]/r}/{[ν+k*(1+ν)*(r/R)^2]/r}
   =A                                                              (4)   
由(3):
   g(θ)=B*e^(-A*θ)                                         (5)
由(4):
    f(r)=C*e^{(A+1)*r-R/Sqrt[ν*(1+ν)*k]*Arctg{sqrt[k*(1+ν)/ν]/R}}                                                             (6)
将(5)(6)代入(2)得到σr,再将σr代入(1)即得到σθ。
3楼2014-04-15 10:17:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

补遗:计算中出现的A、B、C均为积分常数,要由方程定解问题的边界条件来确定。
4楼2014-04-15 10:21:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmei527

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by peterflyer at 2014-04-15 10:17:24
这应该是弹性力学中的无限长圆柱杆件的轴对称问题的力的平衡方程吧。由广义胡可定律,εθ=1/E*(σθ-ν*σr)
依题意εθ=k*σr*(1+ν)/E*r^2/R^2
由此得到σθ=σr*                           (1)
代入力的平衡 ...

谢谢~
5楼2014-04-15 15:57:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


引用回帖:
5楼: Originally posted by bingmei527 at 2014-04-15 15:57:44
谢谢~...

更正:前面f(r)的积分计算有问题,现更正如下:
f(r)=C*r^[(A+1)*ν-1]*e^[(A+1)*k*(1+ν)*r^2/(2*R^2)]                                                            (6)
6楼2014-04-15 16:25:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmei527

铁虫 (初入文坛)

引用回帖:
6楼: Originally posted by peterflyer at 2014-04-15 16:25:28
更正:前面f(r)的积分计算有问题,现更正如下:
f(r)=C*r^*e^                                                            (6)...

谢谢~
7楼2014-04-16 15:36:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nagami

木虫 (正式写手)

这字写得不错,比我蚂蚁爬的好多了
女靠衣装;男靠金装
8楼2014-04-16 16:53:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bingmei527 的主题更新
信息提示
请填处理意见