24小时热门版块排行榜    

查看: 1493  |  回复: 12

zhengjx

金虫 (正式写手)

[求助] 曲面积分计算 已有2人参与

已知:r、R、a、b>0;R>r。
1、求球体x^2+y^2+z^2=r^2被圆柱体(x-a)^2+z^2=R^2和平面y=0所截下部分的曲面面积。
2、求球体x^2+y^2+z^2=r^2被圆柱体(x-a)^2+z^2=R^2和平面y=0、y=-b所截下部分的曲面面积。
请高手给出详细的求解步骤!书上只有一个基本公式:
曲面积分计算
Snap2.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

科研仅是生活的一小部分
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

hank612

至尊木虫 (著名写手)

【答案】应助回帖

人力有时而穷
曲面积分计算-1
emuch2.jpg

We_must_know. We_will_know.
2楼2014-01-07 13:58:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
zhengjx: 金币+20, ★★★很有帮助, 非常感谢!交点坐标好像不对。 2014-01-08 11:17:24
(1)第一问   
此处由于圆柱体为(x-a)^2+z^2=R^2,积分区域S在XOZ平面上,是圆(x-a)^2+z^2=R^2和圆x^2+z^2=r^2的交集,设它们的交点分别为P和Q,
    故xP=(r^2-R^2)/(2*a) ; xQ=xP
         zP=sqrt{r^2-xP^2}   ; zQ=-zP
  故而用下列公式求面积而不用楼主写的公式。
A=DoubleIntegral{sqrt[1+(Py/Px)^2+(Py/Pz)^2]*dS,S}
   其中,y=sqrt{r^2-x^2-z^2}, dS=dx*dz
所以:
A=DoubleIntegral{r/sqrt[r^2-x^2-z^2]*dx*dz , S}
  =Integral{dz*Integral{r/sqrt[r^2-x^2-z^2]*dx, a-sqrt[R^2-z^2], sqrt[r^2-z^2]}, zQ, zP}
   =Integral{π/2-ArcSin{[a-sqrt(R^2-z^2)]/sqrt(r^2-z^2)}, zQ, zP}
   =π/2*(zP-zQ)-Integral{ArcSin{[a-sqrt(R^2-z^2)]/sqrt(r^2-z^2)}, zQ, zP}
  =π*zP-2*Integral{ArcSin{[a-sqrt(R^2-z^2)]/sqrt(r^2-z^2)},0, zP}
   后面的积分利用了被积函数的偶函数性质以及积分区间的对称性,但它是个椭圆积分,积不出解析式出来,只能求数值积分。
3楼2014-01-07 22:06:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zhengjx

金虫 (正式写手)

引用回帖:
3楼: Originally posted by peterflyer at 2014-01-07 22:06:53
(1)第一问   
此处由于圆柱体为(x-a)^2+z^2=R^2,积分区域S在XOZ平面上,是圆(x-a)^2+z^2=R^2和圆x^2+z^2=r^2的交集,设它们的交点分别为P和Q,
    故xP=(r^2-R^2)/(2*a) ; xQ=xP
         zP=sqrt{r^2-xP^ ...

参考您的思路
第一问:
1、交点应为Xp=(a^2-R^2+r^2)/(2*a) 。
2、A=DoubleIntegral{sqrt[1+(Py/Px)^2+(Py/Pz)^2]*dS,S}
   其中,y=sqrt{r^2-x^2-z^2}, dS=dx*dz
考虑对称性,
A=DoubleIntegral{r/sqrt[r^2-x^2-z^2]*dx*dz , S}
  =2*Integral{dz*Integral{r/sqrt[r^2-x^2-z^2]*dx}
积分上限x1=r,下限x2=xp,上限z1=zp,下限z2=0。

下一步怎么求,请继续指导!非常感谢!
科研仅是生活的一小部分
4楼2014-01-08 11:15:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengjx

金虫 (正式写手)

可否采用柱面坐标求解?高数知识都忘得差不多了,重新看书太难了。
科研仅是生活的一小部分
5楼2014-01-08 11:19:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


引用回帖:
4楼: Originally posted by zhengjx at 2014-01-08 11:15:56
参考您的思路
第一问:
1、交点应为Xp=(a^2-R^2+r^2)/(2*a) 。
2、A=DoubleIntegral{sqrt*dS,S}
   其中,y=sqrt{r^2-x^2-z^2}, dS=dx*dz
考虑对称性,
A=DoubleIntegral{r/sqrt*dx*dz , S}
  =2*Integral{ ...

楼主是个细心人,Xp应该是:   Xp=[r^2+a^2-R^2]/(2*a), Zp=Sqrt[r^2-Xp^2].  谢谢指正。关于积分,实在是积不出来。可考虑将积分上下限通过代换变为-1和+1,然后用高斯-勒让德多项式的零点为插值点进行所谓的高斯数值积分运算,这在数值分析教材中都已讲过的,楼主查阅一下教材吧。
6楼2014-01-08 11:49:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengjx

金虫 (正式写手)

参考您的思路
第一问:
1、交点应为Xp=(a^2-R^2+r^2)/(2*a) 。
2、A=DoubleIntegral{sqrt[1+(Py/Px)^2+(Py/Pz)^2]*dS,S}
   其中,y=sqrt{r^2-x^2-z^2}, dS=dx*dz
考虑对称性,
A=DoubleIntegral{r/sqrt[r^2-x^2-z^2]*dx*dz , S}
  =2*Integral{dz*Integral{r/sqrt[r^2-x^2-z^2]*dx}
积分上限x1=xp,下限x2=-r,上限z1=zp,下限z2=0。

积分上下限搞错了。
附图片:
曲面积分计算-2
Snap3.jpg

科研仅是生活的一小部分
7楼2014-01-08 11:49:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

引用回帖:
4楼: Originally posted by zhengjx at 2014-01-08 11:15:56
参考您的思路
第一问:
1、交点应为Xp=(a^2-R^2+r^2)/(2*a) 。
2、A=DoubleIntegral{sqrt*dS,S}
   其中,y=sqrt{r^2-x^2-z^2}, dS=dx*dz
考虑对称性,
A=DoubleIntegral{r/sqrt*dx*dz , S}
  =2*Integral{ ...

根据积分区域的形状,应该先对x求积分,将被积函数中的x用z表示。积分下限是柱面与y=0的交线:x=a-sqrt(R^2-z^2),为何根式取负号,楼主画一下图就清楚了;积分上限为球面与y=0的交线:x=sqrt(r^2-z^2),被积函数为r/sqrt[r^2-x^2-z^2],求出后的一元被积函数变为z的函数,积分上下限分别为ZQ=-Zp和Zp。
8楼2014-01-08 12:00:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

引用回帖:
5楼: Originally posted by zhengjx at 2014-01-08 11:19:44
可否采用柱面坐标求解?高数知识都忘得差不多了,重新看书太难了。

用极坐标也得不到简化,因为被积分区域的表达式复杂,不像圆心在坐标原点时那样能大幅减少计算量。
9楼2014-01-08 12:02:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

引用回帖:
4楼: Originally posted by zhengjx at 2014-01-08 11:15:56
参考您的思路
第一问:
1、交点应为Xp=(a^2-R^2+r^2)/(2*a) 。
2、A=DoubleIntegral{sqrt*dS,S}
   其中,y=sqrt{r^2-x^2-z^2}, dS=dx*dz
考虑对称性,
A=DoubleIntegral{r/sqrt*dx*dz , S}
  =2*Integral{ ...

内积分,将Z看作常数,只对x积分,楼主应该知道1/sqrt[a^2-u^2]的不定积分为arcSin(u/a)+C吧?
10楼2014-01-08 12:05:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhengjx 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 17:35 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 16:15 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 14:58 by CgyNCDVNhGVg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:26 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
信息提示
请填处理意见