数理解析研究所講究録
第 744 巻 1991 年 48-62

48
平均曲率一定の曲面の安定性0.

電気通信大学

小林 昭彦 (Akihiko Kobayashi)

電気通信大学渡辺 二郎 (Jiro Watanabe)

はじめに

図 0.1 のような容器に気体が詰められているとする.

図 $0.1$

容器の壁は固定境界の部分と薄い膜からできている自由境界の部分からなるものとする. 薄膜 $S$ の表面張力が一様であるとすれば,$S$

は容器の口の枠$\Gamma$に張られる平均曲率一定の曲面である. また, この平均曲率は容器の内部と外部の圧力差 $p_{\infty}-p$ に比例する.

容器の外部の圧力$p_{\infty}$は一定であるとし, 容器の内部の圧力 $P$ をゆっくり変化させる.

このとき, 境界面は安定な状態をたどって変化するものとする.

容器の口の枠が適当な形であるとき, ある圧力差のところで境界面 $S$ は不連続的に変化することが数値計算の結果として示される (図 6.3-4 と図 8.3-4 参照). このことにっいて報告する.
注意 本稿では平均曲率一定の曲面の安定性を論ずるが, 曲面が安定であるとか不安定であるとかいうのは, どのような物理系で考えるかに依存する.

したがって, 同一の曲面でも, 考える自由エネルギーを変えれば安定にも不安定にもなりうる. すなわ

ち, 曲面の安定性は曲面が本来もっている性質ではない.

-1-

49
1. 与えられた 「体積」 をもっ平均曲率一定の曲面
3 次元空間内にひとつの Jordan 閉曲線 P が与えられている. F に張られる曲面 S があり, そのパラメータ表現のひとつを $x$ とする. すなわち
$x:\overline{B}=\{u^{2}+v^{2}\leqq 1\}arrow R^{3}$.

$x$はなめらかで 1 対 1 であり,$x_{u}\cross x_{v}\neq 0_{*}$$x(\overline{B})=S$$x|_{\text{\^{a}} B}$

は.

$\Gamma$のパラメータ表現であるとする.
便宜上, 曲面 $S$ というべきところを $x$ ということがある. 本稿ではとくにことわらな

いかぎり, 曲面といえば\Gamma に張られているものとする. さらに記号として$\alpha(x)=S$の面積,$v(x)= \frac{1}{3}\int_{S}x\cdot ndS$

( $S$ (または ) の体積 ).

(1.1)$x$を定義する. ここで は $S$ の単位法線ベクトルである.
$n$

図 1.1

H.C.Wente [12] と K.Steffen [8] は. 任意の $V\in R$ に対して,$\Gamma$に張られる曲面 $x$ で$v(x)=V$ をみたすもののなかに面積 $a(x)$ を最小にするものがあることを示した. この曲面を $x(V)$ とかく:

a$(x(V))= \min_{vtx=V}\mathfrak{a}(x)$.

(1.2)

この曲面 $x(V)$ の平均曲率 $H(V)$ は一定である.
各 $V\in R$ に対して曲面 $x(V)$ とその平均曲率 $H(V)$ を数値計算で求めることを考える.

注意

上記のように Wente [12] と Steffen [8] は, 与えられた体積をもつ曲面のう$-2^{-}-$

50ちで面積最小のものを求めることにより平均曲率一定の曲面を求めた.

しかし他の多くの存在定理は, 与えられた平均曲率をもっ曲面を求める ([1], [4], [9], [10], [11] など). このとき平均曲率の絶対値が十分小さくなければ, これを平均曲率とする曲面は存在しない ([3], [4], [11] など).

2. 数値計算の方法
条件 $v(x)=V$ のもとで面積 $a(x)$ を最小とするような曲面 $x$ を求める. 、この極値問題に対する Lagrangian は$\varphi(x.,H)=a(x)+2H[v(x)-.V]$である.
$\varphi$の停留点 $(x.H)$ を求めたい.

そのために Lagrangian$\varphi$の代わりに augm ented Lagrangian$\psi(x,H)=\varphi(x.H)+r\cdot[v(x)-V]^{2}$を考え,$\psi$$\psi$.

$r>0$の停留点を宇沢アルゴリズムで求める [2].の停留点 $(x_{0},H_{0})$ は(P)$\{\nabla_{H}^{x}\psi(x.H)=2[v(x_{0}^{0})\nabla\psi(x_{0_{0}}.H_{0_{0}})=\nabla_{x}a(x);_{V]=0^{+r[v(x_{0})-V]\}\nabla_{x}v(x_{0})=0}}2\{H_{0}$をみたすものであるから, これから容易に$\psi$$\varphi$と$\psi$.の停留点が同一であることがわかる.

は $x$ に関して下に有界であり, 次の定理がなりたっ.

任意の H\in R と正の r に対して \Gamma に張られる曲面 xo で定理(2.1),$\psi(x_{0},H)=\min_{x}\psi(x,H)$

をみたすものが存在する. このどき,$x_{0}$の平均曲率は$H+r[v(x_{0})-V]$である.
この最小化問題を数値的に解くには, 例えば共役勾配法 (Fletcher-Reeves 法) や準Newton 法 ( $Davidon-Fletcher$ -Powell 法) といった反復解法を用いればよい [7].
$\psi$の停留点を求める問題 (P) を解くために宇沢アルゴリズムを用いる.

宇沢アルゴリズム(1)$H_{0}\in R$を求める.

-3–

51
(2)

$H_{\mathfrak{n}}\in R$が与えられたとき,$\{\begin{array}{l}\psi(x_{n}.H_{\mathfrak{n}})=\min\psi(x.H_{\mathfrak{n}})xH_{n*1}=H_{n}+r[v(x_{n})-V]\end{array}$

によって$x_{\mathfrak{n}}$と$H_{\mathfrak{n}+1}$を求める.

このようにして得られる列$v(x_{\infty})=V$$(x_{\mathfrak{n}}.H_{\mathfrak{n}})$が$narrow\infty$のとき収束すれば, その極限$(x_{\infty},H_{\infty})$

は,をみたしていれば, 問題 (P) の解であることを証明することができる.

3. 数値計算結果の概要
Jordan 閉曲線 として図 3.1 の枠を考える.
$\Gamma$

この枠に張られる極小曲面は, $1.40<\ell<1.53$ のとき 3 種類あるが, [5] ではこれらのうちの 2 種類を数値計算で求めている.

この枠に近い形のものが図 3.2 の枠である.
たものである.

$\Gamma$の線分の部分を catenary でおきかえ

Wiener はこの枠 ( $1.0554<l<1.3254$ の範囲で, 適当な開き角の場合)に少なくとも 2っの極小曲面が張られることを注意した [6]. このことに因んで図 3.1 の枠を,問題 (P) を数値的に解くことにより, 問題 (1.2) の解である曲面 $x(V)$ とその平均曲率$H(V)$ が求められた.

関数 $Varrow>a(x(V))$ と$V\mapsto H(V)$ の計算結果のグラフが図 3.3と図 3.4 である.
これから $da(x(V))/dV=-2H(V)$ がなりたっことが読み取れるが, このことは, 微分可能性があれば, 証明することができる.
以後, 本稿では $x(V),H(V)$ は $V$ に関して微分可能であると仮定する.

求められた曲面 $x(V)$ が汎関数-4-52$\varphi(x.H(V))=a(x)+2H(V)\cdot v(x)$の極小を与える点かどうかを数値実験で確かめてみた. すなわち $x(V)$ を出発点として勾配法によって $\varphi(x,H(V))$ がより小さい値をとる点をたどらせたが, その軌跡は
1) $H’(V)<0$ のときは. 出発点 $x(V)$ から動かなかった. したがって $x(V)$ は$\varphi(x.H(V))$の極小点であることが数値的に確かめられた.

2) $H’(V)>0$ のときは. 出発点 $x(V)$ からしだいにはずれていった. したがっては $\varphi(x.H(V))$ の鞍点であることが数値的に確かめられた.

ao$1.\emptyset$

図 3.3

9.01.0

図 3.4$a\emptyset$$3.\emptyset$$4.\emptyset$

a=a(x(V)) のグラフ

$2.\emptyset$$3.\emptyset$

H=H(V) のグラフ

-5-$4.\emptyset$$\vee$$x(\overline{V})$53
4. 境界面の安定性
いま考えている物理系において温度は一様かっ一定であるとする. 系のエネルギー変化は, 自由境界面 $S$ の面積と位置が変化することによって引き起こされるものとする.

したがって系のエネルギーの変化する部分は$\alpha\cdot a(x)+(p_{\infty}-p)v(x)[=\alpha\cdot\varphi(x.H_{*})$.

$H_{*}= \frac{p_{\infty}-1p}{2\alpha}$(一定)$]$である. 第 1 項の増加は, 表面張力に逆らう仕事によって引き起こされ, 第 2 項の増加は圧力差に逆らう仕事によって引き起こされる.
この系の平衡状態は(4. 1)$\nabla_{x}\varphi(x_{0}.H_{*})=0$をみたすような$x_{0}$であり, この状態

$x_{0}$が安定か不安定かは,$x=x_{0}$が $\varphi(x,H_{*})$ の極小点であるか鞍点であるかである.
曲面 $x(V)$ とその平均曲率 $H(V)$ は (4.1) をみたすが,$x(V_{0})$において(又は $>0$ )$\frac{dH}{dV}(V_{0})<0$

ならば,$V=V_{0}$は安定 (又は不安定) である. すなわち $x(V_{0})$ は $\varphi(x.H(V_{0}))$ の極小点 (又は鞍点) である.
証明
$f(V)=\varphi(x(V).H(V_{0}))=a(x(V))+2H(V_{0})\cdot V$

とおく. このとき$f’(V_{0})=-2H(V_{0})+2H(V_{0})=0$.

$f”(V_{0})=-2H’(V_{0})$がなりたっ. したがって$V_{0}$

ならば,$V$が$V_{0}$に十分近いとき.

$f(V)\geqq f(V_{0})$

したがって $V$ が$H’(V_{0})<0$に十分近く, かつ $v(x)=V$ であるとき,$\varphi(x,H(V_{0}))\geqq\varphi(x(V),H(V_{0}))\geqq\varphi(x(V_{0}),H(V_{0}))$がなりたち,$x(V_{0})$.の近傍として, 十分小さい $\delta>0$ に対して, $\{x:|V(x)-V_{0}|<\delta\}$をとれば, ここで $\varphi(x,H(V_{0}))$ は $x=x(V_{0})$ で極小となる.

逆に$H’(V_{0})>0$$f(V)<f(V_{0})$であるから, V が$V_{0}$

ならば.

.

$V$が$V_{0}$に十分近いとき$V\neq V_{0}$に十分近く. かっ $v(x)=V_{0}$ であれば,$\varphi(x(V),H(V_{0}))\leqq\varphi(x(V_{0}),H(V_{0}))\leqq\varphi(x_{*}H(V_{0}))$-6-

54がなりたっ. これは,

注童

$x=x(V_{0})$ が鞍点であることを意味する.

上記の結論が数値計算の結果と一致することは前節で述べた通りである.

5. 最小固有値の符号
いま, 曲面の平均曲率 H。が一定であるとすると$S_{0}$$\nabla_{x}\varphi(x_{0*}H_{0})=0$(5. 1)$\nabla_{x^{2}}\varphi(x_{0}.H_{0})(\zeta.\zeta)=\int_{S_{0}}\zeta\cdot L\zeta dS$(5.2)

がなりたっ. ここで$x_{0}$は$S_{0}$のパラメータ表現であり,$L\zeta=-\triangle\zeta-2(2H_{0}^{2}-K_{0})\zeta$$\triangle$

は$S_{0}$.

上の L. aplace-Beltram 作用素,$i$

$K_{0}$

は$S_{0}$

方向外向きの変位を表す滑らかな実数値関数で$\Gamma$のガウス曲率である.

$\zeta$は$S_{0}$上の法線上では $0$ であるとする.

次の固有値問題を考える :$\{\Gamma^{0}S$

上上でで.

$L\zeta=\lambda\zeta\zeta=0$

この問題のスペクトル集合は実固有値の列$\lambda_{1}<\lambda_{2}\leqq\lambda_{3}\leqq\cdotsarrow\infty$

であり

(5.3)

$\lambda_{1}=\min\{\int_{s_{0}}\zeta\cdot L\zeta dS;\int_{S_{0}}\zeta^{2}dS=1\}$である. 曲面 $x(V)$ 上の$L$の最小固有値

H’(V)<O(又は>0) ならば$\lambda_{1}$にっいて$\lambda_{1}\geqq 0($又は $<0)$がなりたつ.
証明$H’(V)<0$ のとき, 前節の結果から $v(x)$ が V に十分近いとき$\varphi(x,H(V))\geqq\varphi(x(V),H(V))$

であるから, (5. 1-2) により$\int_{s_{0}}\zeta\cdot L\zeta dS=\nabla_{x^{2}}\varphi(x(V).H(V))(\zeta.\zeta)\geqq 0$

したがって, (5.3) から$\lambda_{1}\geqq 0$.

である.

次に$\zeta(V)=\frac{dx(V)}{dV}\cdot n$( は $x(V)$ の単位法線ベクトル)$n$-7-(5.4)

55とおけば,$\nabla_{x^{2}}\varphi(x(V).H(V))(\zeta(V),\zeta(V))=-2H’(V)$(5.5)

であるから, $H’(V)>0$ のときは上の等式と (5.2) から$\int_{s(\gamma)}\zeta(V)\cdot L\zeta(V)dS<0$

したがって (5.3) から$\lambda_{1}<0$.

である.

6. Wiener の枠に張られる境界面の安定性
第 4 節の結果を枠が Wiener 枠である場合に適用する.
いま考えている物理系 (図 0.1) の温度は一様で一定であるとする. 容器外部の圧力は一定であるとする. 容器内部の圧力$P$$p_{\infty}$をゆっくり変化させるときの自由境界 $S$ の形の変化を調べる. 容器の体積は十分大きくて,$S$の形の変化に伴う容器内部の圧力 $p$ の変化はないものとする.
容器内部の圧力が $p$ であるとき, 境界面 $S$ は平均曲率一定の曲面であり, その平均曲率は$H=(p_{\infty}-p)/2\alpha$である. この曲面には曲線 $H=H(V)$ 上の 1 点 (V, $H(V)$ ) が対応する.

$\sqrt{}$

図 6.1

もし図 6.1 のように 3 点 A.B.

$C$で条件$H(V)=(p_{\infty}-p)/2\alpha$をみたすならば, 点 $B$ では $H’(V)>0$ であり, 対応する曲面は不安定であるので,ような曲面はいま考えている系において現れることはない.

-8-

この

56
点 $A_{*}C$ では $H’(V)<0$ であり, 対応する曲面は安定であるので, このような曲面は現れることになる.

いま容器内部の圧力 $p$ を$A_{0}$から$A_{2}$$p_{0}$から$p_{2}$まで徐々に増加させると, 境界面 $S$ に対応する点はまで進み, さらに $P$ を増加させると $S$ に対応する点は$A_{2}$から$C_{2}$に移り. そこから C3 の方向に進む (図 6.2. 図 6.4).

逆に容器内部の圧力 p を p3 から p1 まで減少させれば, 境界面に対応する点は C3 から$C_{1}$まで進み, さらに $p$ を減少させれば,$S$に対応する点は$C_{1}$から$A_{1}$に移り. そこから$A_{0}$の方向に進む (図 6.3).

図 6.2$H=\frac{P_{\Phi}rightarrow P}{2\wedge}\{$

図 6.3

2

図 6.4-9-

57
7. 容器の体積が比較的小さい場合の境界面の安定性
図 0.1 のような容器に気体が入っているとする. 容器の壁は固定境界の部分と自由境
界 $S$ の部分からなり,$S$$V.\nu$の境界$\Gamma$$S$は薄膜でできていて, その表面張力は一様であるとする.

は変化しないとする. 容器内部の気体の圧力, 体積, モル数をそれぞれとし. 容器外の気体の圧力$p_{\infty}$$P$,はっねに一定とする. 容器内外の気体の温度は一様で一定とする.
モル数\mbox{\boldmath $\nu$} が一定であるとして,$S$を法線方向外向きに$\zeta$だけ変化させたときの自由エネルギー $F$ の第 1 変分は(7.1)$\nabla_{x}F(\zeta)=-2\alpha\int_{S}H\zeta dS+(p_{\infty}-p)\int_{S}\zeta dS$である. ここで$\alpha$は境界面 $S$ の表面張力係数であり, $H$ は $S$ の平均曲率 (曲率の中心が容器の外側にあるとき曲率は正) とする.

(7.1) の右辺の第 1 項は,$S$の面積変化に伴う仕事であり, 第 2 項は体積 $V$ の変化に伴うものである. 平衡状態は$\nabla_{x}F=0$で与えられる. (7. 1) から$\nabla_{x}F=0$.

$p_{\infty}-p=2\alpha H$$\Leftrightarrow$

いま平衡状態にある曲面を$S_{0}$

そのパラメータ表現を

この節では, 容器内の気体のモル数積$1^{\gamma}$

.

$\nu$

$x_{0}$とする.

は一定であるとする. 容器内の圧力 $p$ はその体だけの関数であるとし,$p=\nu RT/V$(7.2)(R は気体定数)がなりたっものとする. すなわち容器内の気体は理想気体で温度は変化しないとする.
いま$\overline{H}_{\nu}(V)=\frac{p_{\infty}-\nu RT/V}{2\alpha}$

とおく. 条件 (7.2) のもとでは$p_{\infty}-p=2\alpha\overline{H}_{\nu}(V)$である. 条件 (7.2) のもとで

$x_{0}$

における $F$ の第 2 変分は$\nabla_{x^{2}}F(\zeta.\zeta)=\alpha\int_{s_{0}}(\zeta L\zeta+4H_{0}^{2}\zeta^{2})dS+2\alpha\overline{H}_{\nu}’(\dot{V}_{0})\cdot(\int_{S_{0}}\zeta dS)^{2}$$+2 \alpha\overline{H}_{\nu}’(V_{0})(-2\int_{s_{0}}H_{0}\zeta^{2}dS)$

-10-

58
(7.3)

$= \alpha\int_{S_{0}}\zeta\cdot L\zeta dS+2\alpha\overline{H}_{\nu}’\cdot(\int_{S_{0}}\zeta dS)^{2}$

ここで $V_{0}=v(x_{0}),$

いま,

とくに

$\zeta$$p_{0}=\nu RT/V_{0}$.

$\overline{H}_{\nu}(V_{0})=(p_{\infty}-p)/2\alpha=H_{0}$として (5.4) の $\zeta(V)$ をとる.

である.

(7.3) において $x_{0},S_{0}.V_{0*}H_{0}$ はそれぞれ$x(V).S(V).V.H(V)$ である. (5.2) と (5.5) に注意すれば, (7.3) から$\nabla_{x^{2}}F(\zeta(V).\zeta(V))=2\alpha[\overline{H}_{\nu}’(V)-H’(V)]$がわかる.

またさらに条件 (7.2) のもとでは, (7.1) から$F(x)=\alpha\cdot\alpha(x)+p_{\infty}v(x)-\nu RT\cdot\log v(x)+const$であることがわかる.
以上のことから$\frac{dF(x(V))}{dV}=\alpha\cdot\frac{da(x(V))}{dV}+p_{\infty}-\frac{\nu RT}{V}=-2\alpha H(V)+2\alpha\overline{H}_{\nu}(V)=0$,$\frac{d^{2}F(x(V))}{dV^{2}}=\alpha\cdot\frac{d^{2}a(x(V))}{dV^{2}}+\frac{\nu RT}{V^{2}}=-2\alpha H’(V)+2\alpha\overline{H}_{\nu}’(V)$

がわかる.

したがって

$\overline{H}_{\nu}’(V)-H’(V)<0$であるとき. 条件 (7.2) のもとで $F$ は $x=x(V)$ において極小とはならない. 逆に$\overline{H}_{\nu}’(V)-H’(V)>0$であるときは,$V_{1}$が $V$ に十分近いとき$F(x(V_{1}))\geqq F(x(V))$

であり,また $v(x)=V_{1}$ であるとき$F(x)\geqq F(x(V_{1}))$であるから, $v(x)$ が $V$ に十分近いときっねに$F(x)\geqq F(x(V))$

がなりたっ. すなわち条件 (7.2) のもとで $F$ は $x=x(V)$ において極小である.

8. Wiener 枠に張られる境界面 (容器の体積が比較的小さい場合)
容器の口の枠が Wiener 枠である場合に前節の結果を適用する.

図 8.1 の点 A. では$C$$\overline{H}_{\nu}’(V)-H’(V)>0$であるから,A.

$C$に対応する曲面は安定であり, 点 $B$ では$-11-$

59$\overline{H}_{\nu}’(V)-H’(V)<0$であるから,$B$に対応する曲面は不安定である. ただし条件 (7.2) がっねになりたっよ

うな変化に関する安定性である.
いま容器内部の気体のモル数をゆっくり変化させる過程を考える. 温度は一様で一定であるとし, 条件 (7.2) が成り立っものとする.

図 8.1
図 8.2 において
$A_{2}$$\nu$が\mbox{\boldmath $\nu$}0 からに徐々に変化していく.に移る (図 8.3). さらに
逆に$\nu$を$\nu_{3}$$\nu$$\nu_{2}$$\nu$をまで徐々に増加していけば, 曲面に対応する点は$\nu_{2}$からさらに増加させれば. $S$ に対応する点はを増加させれば,から減少させるときは,$S$$S$に対応する点はに対応する点は

図 8.2

-12-

$C_{2}$

から$C_{3}$$C_{3}arrow C_{1}arrow A_{1}arrow A_{0}$$A_{0}$$A_{2}$

からから$C_{2}$の方向に進む.と進む (図 8.4).

60

$O$
図 8.3

.

$C_{2}$

図 8.4$\ovalbox{\tt\small REJECT}$A1$C_{t}$

9. シャボン玉の安定性
枠$\Gamma$は円周とする. 原点を$\Gamma$の中心にとって体積 $V$ を定義すれば, 曲線 $H=H(V)$ は

図 9.1 のようになる.

図 9.1

図 9.2

石鹸膜とストローからなる図 9.2 のような系を考える. この系においては石鹸膜が自由境界であり, ストローとその底面が固定境界である. 温度 $T$ は一様かっ一定とする.

内部の体積は $V+V_{0}$.

気体のモル数は$\nu$とする.

-13-61関数$H= \overline{H}_{\nu}(V)=\frac{1}{2\alpha}(p_{\infty}-\frac{\nu RT}{V_{0}+V})$のグラフを図 9.1 の中にかきこめば, 図 9.3 のようになる.

図 9.3

図 9.3 からわかるように 2っの曲線 $H=H(V)$ と $H=\overline{H}_{\nu}(V)$ の交点において$\overline{H}_{\nu}’(V)>H’(V)$であるから, この交点に対応する境界面は第 7 節の結果により安定である.

したがって円形枠に張られる石鹸膜は常に安定である.
次に球形のシャボン玉の安定性について示す. シャボン玉の半径を r, 体積を V, 平
均曲率は $H=H(V)=-1/r$ であるとする. シャボン玉内部の圧力は $p=\nu RT/V$ である.
$\overline{H}_{\nu}(V)=\frac{1}{2\alpha}(p_{\infty}-\frac{\nu RT}{V})$を $V$ で微分して$\overline{H}_{\nu}’(V)=\frac{\nu RT}{2\alpha V^{2}}=\frac{p}{2\alpha V}$を得る. 一方$H’(V)= \frac{1}{4\pi r^{4}}=-\frac{H}{3V}$であり,$H=(p_{\infty}-p)/2\alpha$$H’(V)= \frac{p-p_{\infty}}{6\alpha V}$

であるから

.

したがって

-14-

62

$\overline{H}_{\nu}’(V)-H’(V)=\frac{2p+p_{\infty}}{6\alpha V}>0$となりシャボン玉は安定である.

文献[11H.Br\’ezis, J.M.Coron: Multiple solutions of H-systems and Rellich’sconjecture, Comm. Pure Appl. Math. 37(1984), 149-187.

[2] M.Fortin,

R.Glowinski: Augmented Lagrangian methods, applications tothe numerical solution of boundary-value problems, Amsterdam, NorthHolland, 1983.
[31E.Heinz: On the nonexistence of a surface of constant mean curvatureinite area and prescribed rectifiable boundary, Arch. Rat.
Mech. Anal. 35 (1969), 249-252.

with[41$f$

S.Hildebrandt: On the Plateau problem for surfaces of constant mean
curvature, Comm. Pure Appl. Math. 23 (1970), 97-114.

[5] M. Hinata, M. Shimasaki, T. Kiyono:Numerical solution of Plateau’$s$problem by a finite element method, Math. Comp. , 28(1974) , 45-60.
[61 J. C. C. Nitsche: Lectures

on minimal surfaces, Vol. 1, Cambridge Univ.

Press, 1989.
[7] 坂和 正敏 :非線形システムの最適化, 森北出版, 1986.

[81 K. Steffen: Fl\"achen konstanter mittlerer$Kr\ddot{u}I\ddagger mlung$mit vorgegebenem

Volumen oder Fl\"achenirAalt, Arch. Rat. Mech. Anal. 49 (1972), 99-128.
[9] K.Steffen: On the nonuniqueness of surfaces with prescribed

constant

mean curvature spanning a given contour, Arch. Rat. Mech. Anal.
94 (1986), 101-122.
[101 M. Struwe: Large H-surfaces via the $mountain-pass-le\mathfrak{n}ma$ , Math. Ann.

270 (1985), 441-459.
[11] M.Struwe: Plateau’s problem and the calculus of variations, Math.
Notes 35, Princeton Univ. Press, 1988.
[121 H.C.Wente: A general existence theorem for surfaces of constant meancurvature, Math. Z. 120(1971), 277-288.

-15-