線形な計算ではこういうメリコミを防げない。BFGSでやってみるか。
全エネルギー＝バルクの線形弾性＋界面の開きに対するLJポテンシャル、を最小化。LJの根拠は無し。計算が簡単でヘシアンが連続で負で無限で正で０になるから。界面の横ズレはEnergy cost を考えない。いや、本当は粒界面は平面でなく、ボロノイ的にできた後に若干緩和して曲面となる。そのコストは出てくるか。
しかしこれだけでも線形弾性では出ない極端な応力集中が出そうだ。
いやいやいや、この自由度数で普通のBFGSは無理だがな。自由度数の二乗のメモリがいる。
http://en.wikipedia.org/wiki/BFGS_method
