主題
Search

最小旋轉曲面

下載 Mathematica 筆記本下載 Wolfram 筆記本

變分法 可以用來找到從點 (x_1,y_1) 到點 (x_2,y_2) 的曲線，當這條曲線繞 x 軸 旋轉時，會產生具有最小 表面積 A 的曲面（即 最小曲面）。這等價於找到穿過兩個圓形線框的 最小曲面。面積 元素是

dA=2piyds=2piysqrt(1+y^('2))dx,
(1)

因此 表面積 是

A=2piintysqrt(1+y^('2))dx,
(2)

而我們正在最小化的量是

f=ysqrt(1+y^'^2).
(3)

這個方程有 f_x=0，所以我們可以使用 貝爾特拉米恆等式

f-y_x(partialf)/(partialy_x)=a
(4)

得到

ysqrt(1+y^('2))-y^'(yy^')/(sqrt(1+y^('2)))=a
(5)
y(1+y^('2))-yy^('2)=asqrt(1+y^('2))
(6)
y=asqrt(1+y^('2))
(7)
y/(sqrt(1+y^('2)))=a
(8)
(y^2)/(a^2)-1=y^('2)
(9)
(dx)/(dy)=1/(y^')=a/(sqrt(y^2-a^2))
(10)
x=aint(dy)/(sqrt(y^2-a^2))=acosh^(-1)(y/a)+b
(11)
y=acosh((x-b)/a),
(12)

這被稱為 懸鏈線，透過旋轉它生成的曲面被稱為 懸鏈面。兩個常數 a 和 b 由以下兩個隱式方程確定

y_1=acosh((x_1-b)/a)
(13)
y_2=acosh((x_2-b)/a),
(14)

這兩個方程無法解析求解。

CatenoidXAxis

一般情況比這個解所暗示的要複雜一些。為了理解這一點，考慮兩個半徑相等 半徑 y_0 的環之間的 最小曲面。不失一般性，取兩個環的中點為原點。那麼兩個端點位於 (-x_0,y_0) 和 (x_0,y_0)，並且

y_0=acosh((-x_0-b)/a)=acosh((x_0-b)/a).
(15)

但是 cosh(-x)=cosh(x)，所以

cosh((-x_0-b)/a)=cosh((-x_0+b)/a).
(16)

反轉每一邊

-x_0-b=-x_0+b,
(17)

因此 b=0 （由於對稱性，必須如此，因為我們選擇了兩個環之間的原點），並且 最小曲面 的方程簡化為

y=acosh(x/a).
(18)

在端點處

y_0=acosh((x_0)/a),
(19)

但是對於某些 x_0 和 y_0 的值，這個方程沒有解。這個事實的物理意義是，曲面斷裂並形成每個環中的圓盤，以最小化 面積。變分法 不能用於找到這種不連續的解（在這種情況下稱為 戈爾德施密特解）。上面顯示了幾個端點選擇的最小曲面。前兩種情況是 懸鏈面，而第三種情況是 戈爾德施密特解。

為了找到可以獲得 懸鏈線 解的 x_0/y_0 的最大值，令 p=1/a。然後 (17) 給出

y_0p=cosh(px_0).
(20)

現在，將 x_0 的最大值記為 x_0^*。那麼 dx_0/dp=0 將成立。對 (20) 取 d/dp，

y_0=sinh(px_0)(x_0+p(dx_0)/(dp)).
(21)

現在設定 dx_0/dp=0

y_0=x_0sinh(px_0^*).
(22)

從 (20)，

py_0^*=cosh(px_0^*).
(23)

取 (23) ÷ (22)，

px_0^*=coth(px_0^*).
(24)

定義 u=px_0^*,

u=cothu.
(25)

這有解 u=1.1996789403...。從 (22)，y_0p=coshu。將此除以 (25) 得到 y_0/x_0=sinhu，因此 x_0/y_0 的最大可能值是

(x_0)/(y_0)=cschu=0.6627434193....
(26)

因此，僅當 x_0/y_0>0.6627... 時，才存在戈爾德施密特環解。

最小 懸鏈面 曲面的 表面積 由下式給出

A=2(2pi)int_0^(x_0)ysqrt(1+y^('2))dx,
(27)

但是由於

y=sqrt(1+y^('2))a
(28)
=acosh(x/a),
(29)
A=(4pi)/aint_0^(x_0)y^2dx
(30)
=4piaint_0^(x_0)cosh^2(x/a)dx
(31)
=4piaint_0^(x_0)1/2[cosh((2x)/a)+1]dx
(32)
=2pia[int_0^(x_0)cosh((2x)/a)dx+int_0^(x_0)dx]
(33)
=2pia[a/2sinh((2x)/a)+x]_0^(x_0)
(34)
=pia^2[sinh((2x)/a)+(2x)/a]_0^(x_0)
(35)
=pia^2[sinh((2x_0)/a)+(2x_0)/a].
(36)
CatenoidSolution

在求解 (◇) 中的 a 時需要謹慎。如果我們取 x_0=1/2 和 y_0=1，則 (◇) 變為

1=acosh(1/(2a)),
(37)

這有兩個解：a_1=0.2350...（“深”），和 a_2=0.8483...（“平坦”）。但是，將這些代入 (◇) 和 x_0=1/2，我們發現 A_1=6.8456... 和 A_2=5.9917...。所以 A_1 實際上不是區域性最小值，而 A_2 是唯一真正的最小解。

懸鏈面 解的 表面積 等於 戈爾德施密特解 的表面積，當 (◇) 等於兩個圓盤的 面積 時，

pia^2[sinh((2x_0)/a)+(2x_0)/a]=2piy_0^2
(38)
a^2[2sinh((x_0)/a)cosh((x_0)/a)+(2x_0)/a]-2y_0^2=0
(39)
a^2[cosh((x_0)/a)sqrt(cosh^2((x_0)/a)-1)+(x_0)/a]-y_0^2=0.
(40)

代入

(y_0)/a=cosh((x_0)/a)
(41)
(y_0)/asqrt(((y_0)/a)^2-1)+cosh^(-1)((y_0)/a)-((y_0)/a)^2=0.
(42)

定義

u=(y_0)/a
(43)

給出

usqrt(u^2-1)+cosh^(-1)u-u^2=0.
(44)

這有一個解 u=1.2113614259。對於

A_(catenary)=A_(2 disks)
(45)

的 x_0/y_0 值是

(x_0)/(y_0)=((x_0)/a)/((y_0)/a)=(cosh^(-1)((y_0)/a))/((y_0)/a)=(cosh^(-1)u)/u=0.5276973967.
(46)

因此 x_0/y_0 in (0.52770,0.6627)，懸鏈線 解的 面積 大於兩個圓盤的面積，因此它僅作為 區域性最小值 存在。

也存在具有一個圓盤（半徑為 r）的解，該圓盤位於由兩個旋轉 懸鏈面 支撐的環之間。面積 大於簡單 懸鏈面 的面積，但它是 區域性最小值。這條曲線的 正 半部分的方程是

y=c_1cosh(x/(c_1)+c_3).
(47)

在 (0,r) 處，

r=c_1cosh(c_3).
(48)

在 (x_0,y_0) 處，

y_0=c_1cosh((x_0)/(c_1)+c_3).
(49)

兩個 懸鏈面 的 面積 是

A_(catenoids)=2(2pi)int_0^(x_0)ysqrt(1+y^('2))dx
(50)
=(4pi)/(c_1)int_0^(x_0)y^2dx
(51)
=4pic_1int_0^(x_0)cosh^2(x/(c_1)+c_3)dx.
(52)

現在令 u=x/c_1+c_3，所以 du=dx/c_1

A=4pic_1^2int_(c_3)^(x_0/x_1+c_3)cosh^2udu
(53)
=4pic_1^21/2int_(c_3)^(x_0/x_1+c_3)[cosh(2u)+1]du
(54)
=2pic_1^2[1/2sinh(2u)+u]_(c_3)^(x_0/x_1+c_3)
(55)
=2pic_1^2{1/2sinh[2((x_0)/(c_1)+c_3)]-1/2sinh(2c_3)+(x_0)/(c_1)}
(56)
=pic_1^2{sinh[2((x_0)/(c_1)+c_3)]-sinh(2c_3)+(2x_0)/(c_1)}.
(57)

中心 圓盤 的 面積 是

A_(disk)=pir^2=pic_1^2cosh^2c_3,
(58)

因此總 面積 是

A=pic_1^2{sinh[2((x_0)/(c_1)+c_3)]+[cosh^2c_3-sinh(2c_3)]+(2x_0)/(c_1)}.
(59)

根據 普拉託定律，懸鏈面 以 120 degrees 的 角度 相遇，所以

tan30 degrees=[(dy)/(dx)]_(x=0)
(60)
=[sinh(x/(c_1)+c_3)]_(x=0)
(61)
=sinhc_3=1/(sqrt(3))
(62)

並且

c_3=sinh^(-1)(1/(sqrt(3))).
(63)

這意味著

cosh^2c_3-sinh(2c_3)=[1+sinh^2c_3]-2sinhc_3sqrt(1+sinh^2c_3)
(64)
=(1+1/3)-2(1/(sqrt(3)))sqrt(1+1/3)
(65)
=4/3-2/(sqrt(3))2/(sqrt(3))=0,
(66)

所以

A=pic_1^2{sinh[2((x_0)/(c_1)+c_3)]+(2x_0)/(c_1)}.
(67)

現在檢查 x_0/y_0,

(x_0)/(y_0)=((x_0)/(c_1))/((y_0)/(c_1))=((x_0)/(c_1))/(cosh((x_0)/(c_1)+c_3))=usech(u+c_3),
(68)

其中 u=x_0/c_1。找到 x_0/y_0 的最大比率給出

d/(du)((x_0)/(y_0))=sech(u+c_3)-utanh(u+c_3)sech(u+c_3)=0
(69)
utanh(u+c_3)=1,
(70)

其中 c_3=sinh^(-1)(1/sqrt(3)) 如上所述。解是 u=1.0799632187，因此帶有中心圓盤的兩個 懸鏈面 的 x_0/y_0 的最大值是 y_0=0.4078241702。

如果我們有興趣找到從點 (x_1,y_1) 到點 (x_2,y_2) 的曲線，當這條曲線繞 y 軸（而不是 x 軸）旋轉時，會產生具有最小 表面積 A 的曲面，我們按上述步驟進行。請注意，該解在物理上等效於繞 x 軸 旋轉的情況，但採用不同的數學形式。面積 元素是

dA=2pixds=2pixsqrt(1+y^('2))dx
(71)
A=2piintxsqrt(1+y^('2))dx,
(72)

而我們正在最小化的量是

f=xsqrt(1+y^'^2).
(73)

求導得到

(partialf)/(partialy)=0
(74)
d/(dx)(partialf)/(partialy^')=d/(dx)((xy^')/(sqrt(1+y^('2)))),
(75)

因此 尤拉-拉格朗日微分方程 變為

(partialf)/(partialy)-d/(dx)(partialf)/(partialy^')=d/(dx)((xy^')/(sqrt(1+y^('2))))=0.
(76)
(xy^')/(sqrt(1+y^('2)))=a
(77)
x^2y^('2)=a^2(1+y^('2))
(78)
y^('2)(x^2-a^2)=a^2
(79)
(dy)/(dx)=a/(sqrt(x^2-a^2))
(80)
y=aint(dx)/(sqrt(x^2-a^2))+b=acosh^(-1)(x/a)+b.
(81)

求解 x 得到

x=acosh((y-b)/a),
(82)

這是 懸鏈線 的方程。懸鏈面 旋轉產生的 表面積 是

A=2piintxsqrt(1+y^('2))dx=2piintxsqrt(1+(a^2)/(x^2-a^2))dx
(83)
=2piintx/(sqrt(x^2-a^2))sqrt((x^2-a^2)+a^2)dx
(84)
=2piint(x^2dx)/(sqrt(x^2-a^2))
(85)
=[x/2sqrt(x^2-a^2)+(a^2)/2ln(x+sqrt(x^2-a^2))]_(x_1)^(x_2)
(86)
=1/2[x_2sqrt(x_2^2-a^2)-x_1sqrt(x_1^2-a^2)+a^2ln((x_2+sqrt(x_2^2-a^2))/(x_1+sqrt(x_1^2-a^2)))].
(87)

Isenberg (1992, p. 80) 討論了找到穿過軸線彼此偏移的兩個環的 最小曲面。

參見

最小曲面, 旋轉曲面

在 中探索

參考文獻

Arfken, G. Mathematical Methods for Physicists, 3rd ed. Orlando, FL: Academic Press, pp. 931-937, 1985.Goldstein, H. Classical Mechanics, 2nd ed. Reading, MA: Addison-Wesley, p. 42, 1980.Isenberg, C. The Science of Soap Films and Soap Bubbles. New York: Dover, pp. 79-80 and Appendix III, 1992.

在 上引用

最小旋轉曲面

請引用為

Weisstein, Eric W. "最小旋轉曲面。" 來自 Web 資源。 https://mathworld.tw/MinimalSurfaceofRevolution.html

主題分類