Math Pro 數學補給站

標題: 多圓柱相交的體積 [打印本頁]

作者: bugmens 時間: 2018-3-17 12:37 標題: 多圓柱相交的體積

之前做兩圓柱相交的體積
https://math.pro/db/viewthread.php?tid=162&page=4#pid17835

最近在做三圓柱相交的體積，有問題想請各位網友提供意見
我參考這裡的資料
http://calculus.nctu.edu.tw/uplo ... ns_porj_cylind.html
在18.bmp的\( z=\sqrt{1-x^2} \)轉換成19.bmp的\( z=\sqrt{1-(r cos \theta)^2} \)
要怎麼用動畫直觀表達出\( x=r cos \theta \)

參考資料的圓柱坐標用\( (r,\theta,z) \)，wiki用的是\( (\rho,\phi,z) \)
https://zh.wikipedia.org/wiki/圓柱坐標系
大家覺得這個有需要改嗎？

107.4.6
目前沒有網友回應，先做成gif動畫
　\(\displaystyle 16 \int_0^{\frac{\pi}{4}} \int_0^1 \int_0^{\sqrt{1-(r cos \theta)^2}}r dz dr d \theta \)
\(\displaystyle =16 \int_0^{\frac{\pi}{4}} \int_0^1 rz \Bigg\vert\;_{z=0}^{z=\sqrt{1-(r cos \theta)^2}} dr d \theta\)
\( \displaystyle =16 \int_0^{\frac{\pi}{4}} \int_0^1 r\sqrt{1-(rcos\theta)^2} dr d\theta \)
\(\displaystyle =16 \int_0^{\frac{\pi}{4}} -\frac{(1-(rcos\theta)^2)^{\frac{3}{2}}}{cos \theta} \Bigg\vert\;_{r=0}^{r=1}d \theta \)
\(\displaystyle =16 \int_0^{\frac{\pi}{4}} \frac{1-sin^3 \theta}{3cos^2 \theta}d\theta\)
\(\displaystyle =\frac{16}{3}(tan\theta-sec\theta-cos\theta)\Bigg\vert\;_{\theta=0}^{\theta=\frac{\pi}{4}}\)
\(\displaystyle =\frac{16}{3}(3-\frac{3\sqrt{2}}{2})\)
\(=16-8\sqrt{2}\)

附件: 三圓柱相交的體積SketchUp檔.zip (2018-4-6 07:34, 398.17 KB) / 該附件被下載次數 11183
https://math.pro/db/attachment.php?aid=4335&k=385dc437313a41866f3d514b6c2fc576&t=1732197383

圖片附件: 三圓柱相交的體積.gif (2018-4-6 08:07, 152.43 KB) / 該附件被下載次數 5911
https://math.pro/db/attachment.php?aid=4345&k=ede3c83c3446467e912e84ab4e24cfa0&t=1732197383

作者: bugmens 時間: 2018-9-23 14:21

Moreton Moore在1974年發表一篇關於多圓柱相交的體積的論文，他設定圓柱的直徑為2\((d=2)\)，\(S_n\)代表\(n\)對稱圓柱相交的體積。
二圓柱相交的體積\( \displaystyle S_2=\frac{2}{3}d^3=\frac{16}{3}\)

三圓柱相交的體積\( S_3=(2-\sqrt{2})d^3=8(2-\sqrt{2}) \)
四圓柱相交的體積\( \displaystyle S_4=\frac{3}{2}\sqrt{2}(2-\sqrt{3})d^3=12(2\sqrt{2}-\sqrt{6}) \)

六圓柱相交的體積\( \displaystyle S_6=\frac{2}{3}(3+2\sqrt{3}-4\sqrt{2})d^3=\frac{16}{3}(3+2\sqrt{3}-4\sqrt{2})\)

徒手畫四圓柱和六圓柱相交的圖形非常麻煩，要借助SketchUp的Ruby程式畫出圖形。

1.設定立方體8個頂點坐標為
\((1,-1,1),(1,1,1),(-1,1,1),(-1,-1,1)\)
\((1,-1,-1),(1,1,-1),(-1,1,-1),(-1,-1,-1)\)
立方體邊長2，中心為\((0,0,0)\)

2.從立方體8個頂點對角兩兩連線形成四圓柱相交的模型
圓柱1：\((-1,-1,-1)\)和\((1,1,1)\)連線
圓柱2：\((1,-1,-1,)\)和\((-1,1,1)\)連線
圓柱3：\((1,1,-1)\)和\((-1,-1,1)\)連線
圓柱4：\((-1,1,-1)\)和\((1,-1,1)\)連線

3.轉換成\(SketcuUp\)程式碼
ent=Sketchup.active_model.entities
center=[[-1,-1,-1],[1,-1,-1],[1,1,-1],[-1,1,-1]]
normal=[[2,2,2],[-2,2,2],[-2,-2,2],[2,-2,2]]
for i in 0..3
  circle=ent.add_circle center[ i ],normal[ i ],1,60
  circle_face=ent.add_face circle
  circle_face.pushpull 4
  ent.add_group circle_face.all_connected
end
程式碼解釋可以看底下參考資料的"當SketchUp遇見Ruby－邁向程式化建模之路"

4.將程式碼儲存成4intersection.rb，將檔案放到SketchUp目錄下
SketchUp2017
C:\Program Files\SketchUp\SketchUp 2017\ShippedExtensions
SketchUp 8
C:\Program Files (x86)\Google\Google SketchUp 8\Plugins

執行SketchUp可以看到四圓柱的模型，但要得到四圓柱交集後的結果，要另外安裝OScoolean plugins。

6.安裝OScoolean plugins
到http://www.plotyuki.sakura.ne.jp ... =28&category=01下載OSCoolean_2013.zip
將解壓縮的檔案放到SketchUp目錄下
SketchUp2017
C:\Program Files\SketchUp\SketchUp 2017\ShippedExtensions
SketchUp 8
C:\Program Files (x86)\Google\Google SketchUp 8\Plugins

7.重新執行SketchUp，全選四個圓柱，按下滑鼠右鍵選"實心工具/交集"，得到四圓柱交集後的模型。

===========================================================

取立方體12個邊中點對角兩兩連線形成六圓柱相交的圖形
圓柱1：\((1,-1,0)\)和\((-1,1,0)\)連線
圓柱2：\((-1,-1,0)\)和\((1,1,0)\)連線
圓柱3：\((1,0,1)\)和\((-1,0,-1)\)連線
圓柱4：\((-1,0,1)\)和\((1,0,-1)\)連線
圓柱5：\((0,-1,1)\)和\((0,1,-1)\)連線
圓柱6：\((0,-1,-1)\)和\((0,1,1)\)連線

轉換成\(SketcuUp\)程式碼
ent=Sketchup.active_model.entities
center=[[1,-1,0],[-1,-1,0],[1,0,1],[-1,0,1],[0,-1,1],[0,-1,-1]]
normal=[[-2,2,0],[2,2,0],[-2,0,-2],[2,0,-2],[0,2,-2],[0,2,2]]
for i in 0..5
  circle=ent.add_circle center[ i ],normal[ i ],1,60
  circle_face=ent.add_face circle
  circle_face.pushpull 3
  ent.add_group circle_face.all_connected
end
其餘部分就和前面相同
但交集後的模型在6條線的交點處會有折線或斷線，這是比較可惜的地方

參考資料
Moore, M. "Symmetrical Intersections of Right Circular Cylinders." Math. Gaz. 58, 181-185, 1974.
http://mathworld.wolfram.com/SteinmetzSolid.html
當SketchUp遇見Ruby－邁向程式化建模之路，https://books.google.com.tw/book ... e&q&f=false

附件: 多圓柱相交的體積moore1974.pdf (2018-9-23 14:21, 253.67 KB) / 該附件被下載次數 11127
https://math.pro/db/attachment.php?aid=4701&k=c8866f7c24c61dc8cbfead18fc55debb&t=1732197383

附件: OSCoolean_2013.zip (2018-9-23 14:21, 25.87 KB) / 該附件被下載次數 8231
https://math.pro/db/attachment.php?aid=4702&k=e0b679f6386387507176e660e3844d2e&t=1732197383

附件: 四圓柱相交的體積SketcuUp檔.zip (2018-9-23 14:21, 49.72 KB) / 該附件被下載次數 8273
https://math.pro/db/attachment.php?aid=4703&k=c6cc82dba32494317aa7b9141821b080&t=1732197383

圖片附件: 四圓柱相交的體積.gif (2018-9-23 14:21, 17.53 KB) / 該附件被下載次數 5463
https://math.pro/db/attachment.php?aid=4704&k=151549630e4056c8e6c155dac10178fd&t=1732197383

附件: 六圓柱相交的體積SketcuUp檔.zip (2018-9-23 14:21, 67.31 KB) / 該附件被下載次數 8092
https://math.pro/db/attachment.php?aid=4705&k=48086c00c87d316543bad7e8e6a9d0b4&t=1732197383

圖片附件: 六圓柱相交的體積.gif (2018-9-23 14:21, 20.11 KB) / 該附件被下載次數 5683
https://math.pro/db/attachment.php?aid=4706&k=cd477dd5f962e85da36e1c633894b2eb&t=1732197383

圖片附件: 二圓柱相交的體積moore.gif (2018-9-23 14:21, 23.28 KB) / 該附件被下載次數 6477
https://math.pro/db/attachment.php?aid=4707&k=6bc3238eda8e418f9f4ca04e60dcc212&t=1732197383

圖片附件: 四圓柱相交的體積moore.gif (2018-9-23 14:21, 39.77 KB) / 該附件被下載次數 6372
https://math.pro/db/attachment.php?aid=4708&k=7b2474783428aa7952389262895edfa3&t=1732197383

圖片附件: 六圓柱相交的體積moore.gif (2018-9-23 14:21, 48.13 KB) / 該附件被下載次數 6536
https://math.pro/db/attachment.php?aid=4709&k=e461649a1a73b593ecd7873a7391f978&t=1732197383

作者: Joy091 時間: 2022-5-6 16:49 標題: 利用牟合方蓋的結果，計算三圓柱相交的體積

另一種想法供參考：
利用對稱性，可將三圓柱相交的體積拆解為6塊局部的牟合方蓋，
再一一貼於某個正立方體的六個面上。

附件: [利用牟合方蓋計算三圓柱相交的體積] 牟合方蓋與牟合圓塊體積推導.pdf (2022-5-6 16:51, 350.22 KB) / 該附件被下載次數 4543
https://math.pro/db/attachment.php?aid=6347&k=8725760daf7eaae2a938a037bfc0d651&t=1732197383

歡迎光臨 Math Pro 數學補給站 (https://math.pro/db/)

論壇程式使用 Discuz! 6.1.0