- モデルの回転や移動は、blenderでやるとかなり動作が速いのでおすすめ。
- 例えば、国土地理院の3Dの富士山↓は元々、面の数が"133110"もあるので、細かすぎて使いにくい。
- 難しい
- 楽に減らせる。面の数が1/10でも↓結構きれい。
- 適宜、「その1」と「その2」を組み合わせる。
- blenderでは、右上の「面:」に面の数が書かれている。
- blenderでは、"z"でエッジが見れる。
| 日付 | 用途 | 出金(円) | 残高(円) |
| 3/28 | カードのポイントに換金 | 0 | 3400 |
| 3/29 | ティッシュ×2 | 740 | 2660 |
| Scalars | radius | |||
| Glyph Type | Sphere | |||
| Radius | 1.0 | |||
| Theta Resolution | 16 | 多いほうが球体の角が取れる | ||
| Scalar mode | scalar | |||
| Set scale factor | 1.0 |
| vectors | v | |||
| Glyph Type | Arrow | |||
| Scalar mode | vector | |||
| Set scale factor | 0.02 | 矢印の大きさ。適度に小さくして、調整。 |
| title | text | タイトルの名前 | |
| Font | times | 好きな書体に | |
| B | 太字 | ||
| I | 斜め | ||
| S | 影 | ||
| size | 20 | 文字サイズ | |
| Color | 黒とかに | ||
| label | 上記と同じ | 帯の方 |
| Animation Duration(sec) | アニメのコマ数 |
| Resolution(pixels) | 横幅・縦幅のピクセル値 |
| 1/29まで | 発表練習とスライドの修正。 |
| 達成度 | した。時間内に終わらないので、あとは、各スライド10秒ぐらい縮めて喋れるようにする。 |
| 1/22まで | 支部TeX・卒論TeXを後藤さんに提出して修正する。 |
| 達成度 | 完成。支部概要は提出。 |
| 1/15まで | 支部TeXを「1.はじめに 2.実習環境 3.実習内容 4.実習課題」ぐらいまで、大まかに書き進める。 |
| 達成度 | あんまり進んでいない。 |
| 1/8まで | libreoffice→powerpointで読める形式を見つける。 |
| 達成度 | libreofficeで、pptファイルとして保存すると、powerpointで普通に開ける。 |
| 12/25まで | odpファイルの画像を小さくする。powerpointでodpファイルを読み込む。 |
| 達成度 | odpファイルのサイズが半分ぐらいになった。odpファイルを開くと、画像が暗くなる。 |
| 12/22まで | 中間発表のスライドとTeXをつくる。 |
| 達成度 | できた。 |
| 12/18まで | 動画をつくる。その２。 |
| 達成度 | 動画を増やした。 |
| 12/11まで | paraviewやblenderで何か面白そうな動画をつくる。横ねじれ座屈、エラスティカ問題......? |
| 達成度 | code-asterとblenderの石積み動画をつくった。 |
| 12/4まで | blenderで何らかの動画ファイルをつくれるようにする。出来ればliggghts用にpraviewの動画ファイルも作りたい。 |
| 達成度 | paraviewからjpgとかを吐き出して、gif化。libbreofficeへの挿入もうまくいった。 |
| 11/27まで | blenderでラーメンの柱のすり抜けをなくす。 |
| 達成度 | 床のコリジョンの設定で「内側」の値を0.2から1.0にすると、すり抜けを防止できた。FEMよりも計算が遅い。 |
| 11/20まで | blenderで片持ち梁をつくる。固定端をmodifierのhookという機能でやってみる。 |
| 達成度 | 壁と梁をくっつけるとどちらかの設定に依存してしまうので、ラーメンにした。曲げを見るには、辺を内部にたくさん生成しないといけないし、そのせいで梁が重すぎて、柱部分がすり抜けるなどの問題がある。 |
| 11/13まで | 2D円筒要素のDEMをつくりたい。(上手くいったら2D多角形要素も考える) |
| 達成度 | 円筒同士の衝突計算、円筒と床の衝突計算をした。それをparaviewでの映像化して観た。 |
| 11/6まで | femとの比較用にblenderで簡単な梁をつくりたい |
| 達成度 | 片持ち梁を作ろうとするも、片面の固定がわからない。 |
| 10/30まで | science directあたりで、physics engineとfemをどう連携しているか調べる。 |
| 達成度 | physicsとfemのハイブリッドしたようなものはなさそう。 |
| 10/23まで | ゲームとFEMをやっている人の文献を調べる。 |
| 達成度 | 日本語のページを調べた。 |
| 10/16まで | blenderのfluidを試す。粘性を上げて、どこまで固体に近づくのか。 |
| 達成度 | 粘性を最も高い値にしてみたが、流体は使えそうにないことがわかった。 |
| 10/9まで | 背面を球体の自由落下でpackingしてみる。それをエクスポートして、pythonでrigidの特性を与えてみる。 |
| 達成度 | 球体の自由落下によるpackingはできたが、その後の石をactiveに変更する作業は自動化できていない。 |
| 9/28まで | 夏休みの課題達成度のTeX・発表用スライドを作る。 |
| 達成度 | とりあえず完成。 |
| 9/21まで | blender --pythonで直方体のソフトボディを無数に作り、裏込め＋背面土を再現する。 |
| 達成度 | 要素数が3桁とかになると、生成を自動化しても計算が終わらない上、どうしてもpythonで制御できない条件がある。そのため、剛体で土っぽくなる条件を探した。 |
| 9/14まで | blenderのソフトボディをsalomeのメッシュを使って先週の続きをやる。 |
| 達成度 | salomeのメッシュを色々な形式に変換してblenderにインポートしたが、ソフトボディの要素が増えるわでないとうことが解った。ソフトボディの要素を細かくする場合、blenderのナイフ機能を使って切り分けなければならない。 |
| 9/7まで | blenderのソフトボディによる裏込めをより詳細に再現してみる。出来ればasterの塑性解析を試す。 |
| 達成度 | どちらもできていない。ソフトボディは物体内部のメッシュを手動で制御(生成)しなければいけないことがわかった。 |
| 8/31まで | blenderでソフトボディで裏込めを再現してみる。 |
| 達成度 | 形としては出来たものの、asterに近いとはいえない。 |
| 8/24まで | 先週のに石を想定した物体を1・2個埋め込み、2種材料として自重で沈む問題をasterで解く。 |
| 達成度 | 石の1･2個どころか、沢山の場合でもヤング率の異なる複合材料として解けた。 |
| 8/17まで | 石積みなし(背面土+裏込だけ)で自重で沈む問題をasterで解く。出来れば過渡応答もつくる。 |
| 達成度 | お盆気分が抜けず、しばらくサボっていたため8/20に達成。石を繰り抜いた部分が細かい場合には解けない。 |
| 8/10まで | salome-mecaで自重でたわむ単純梁をつくって手計算と比較。ついでに、同問題を動的に出来るのか調べる(あったら試す) |
| 達成度 | 自重は単位体積あたりに密度分の荷重をかけると出来る。過渡応答解析も重くなるけどできた。自重を10倍とかにしたのを一瞬だけかけた後、除荷することで、自重方向のインパルス応答もみれた。 |
| 8/3まで | オープンキャンパス用に何か(自由落下->坂->ドミノとか?)つくり、blenderの操作方法の勉強する。 |
| 達成度 | できた |
| 7/27まで | ドミノ倒しの論文を探す、物理エンジンについて調べる、球を繋げる(継続中) |
| 達成度 | blender2.7以上で剛体シミュレーション(弾性、摩擦等を変更できる)ができることを確認。 |
| 7/13まで | 球を2つ繋げてみる、将棋倒しの計算方法を探す |
| 達成度 | 球をつなげる方法が分かるも、formatの問題で出来ず。ドミノ倒しはBlenderが多いことがわかる。 |
| 7/6まで | 球以外のもの(円柱とか)をつくって、坂を転がす。 |
| 達成度 | 球以外できず。Catherine O'Sullivan『粒子個別要素法』をパラパラめくり、球をくっつけて物体をつくるのではないかと想像する。 |
| 6/29まで | 粒子の投入コマンドを調べる。複数の物体の経時変化ををグラフ化する。 |
| 達成度 | 投入コマンドは分からず。stlファイルで坂を入れられた。gnuplotで複数の粒子の経時変化をグラフ化できた。 |
| 6/22まで | 球をピンポイントで配置するためのコマンドを探して、それをできるようにする。1個の物体の経時変化をグラフ化する。 |
| 達成度 | 入力ファイルの編集は進歩なし。gnuplotでグラフ化は成功。 |
$ mogrify -format jpg moto*.png
$ mogrify -crop 400x300+500+400 moto*.jpg
$ mogrify -resize 400 moto*.jpg
$ convert -delay 10 moto*.jpg filename.gif
$ paraview kekka..vtk&
| Glyph Type | Sphere-Radius | Scale Mode | Set Scale Factor |
| Sphere | 1 | scalar | 1(Editにチェックを入れて) |
| 日付 | 時間帯 | 作業時間(hr) | 内容 | 立会 |
| 4/13(月) | 13:30-17:30 | 1.5 | 戸川隼人『有限要素法へのガイド』の1章を読む | |
| 4/16(木) | 14:30-17:30 | 3.0 | 戸川隼人『有限要素法へのガイド』の2章を読む | |
| 4/17(金) | 11:00-12:00 | 1.0 | 戸川隼人『有限要素法へのガイド』の3章を読む | |
| 4/17(金) | 15:30-17:30 | 2.0 | 戸川隼人『有限要素法へのガイド』の3章を読む | |
| 4/21(火) | 9:00-10:00 | 1.0 | 戸川隼人『有限要素法へのガイド』のプログラム2.1を修正しつつ模写 | |
| 4/21(火) | 14:00-17:30 | 3.0 | 戸川隼人『有限要素法へのガイド』のプログラム2.2を修正しつつ模写 | |
| 4/24(金) | 13:30-17:30 | 4.0 | 戸川隼人『FORTRANによる有限要素法入門』の1章を読む | |
| 4/27(月) | 9:00-10:00 | 1.0 | 戸川隼人『FORTRANによる有限要素法入門』のプログラム1.1〜1.2を模写 | |
| 4/27(月) | 13:30-17:00 | 1.0 | 戸川隼人『FORTRANによる有限要素法入門』のプログラム1.3〜1.4を模写 | |
| 4/28(火) | 9:00-12:00 | 3.0 | LIGGGHTSの入手 | |
| 4/28(火) | 13:30-17:30 | 4.0 | LIGGGHTSのコンパイルの準備およびMakefileの修正 | |
| 4/29(水) | 10:00-12:00 | 2.0 | LIGGGHTSのコンパイルの準備から(/opt に)やり直し | |
| 4/29(水) | 13:30-17:30 | 4.0 | LIGGGHTSのコンパイルを色々と試す | |
| 4/30(木) | 13:30-18:30 | 5.0 | LIGGGHTSのコンパイル(後藤ちゃちゃ)をしたり、ライブラリをつくったり、パスを追加したり、サンプルを実行したり、paraviewを入手したり | |
| 4月計 | 35.5 | |||
| 5/1(金) | 9:00-12:00 | 3.0 | pizza.pyのインストール、およびサンプルの実行を試みる | |
| 5/1(金) | 13:30-18:00 | 4.5 | サンプルを計算しparaviewで再生するも失敗(vtkファイルは開けるし、Timeも動いているが、物体？が現れない) | |
| 5/6(水) | 14:00-18:00 | 3.5 | LIGGGHTSの読み方を調べたり、linuxのリンクについて調べたり | |
| 5/7(木) | 16:00-18:00 | 2.0 | LIGGGHTSのサンプルの実行(が出来ず。dumpファイルが見つからない？) | |
| 5/8(金) | 13:00-18:00 | 5.0 | LIGGGHTSをインストールからやり直す。examplesの計算、結果ファイルdumpをvtkへの形式変換、paraviewで(赤い点々がうようよ動く)の確認 | |
| 5/9(土) | 10:00-12:00 | 2.0 | LIGGGHTSのサンプル「conveyor」のFiletersのGlyphをいじる | |
| 5/9(土) | 16:00-18:00 | 2.0 | PENGUINITISさんのサンプル「box」を試すも、実行出来なかったのでin.boxをいじる | |
| 5/12(火) | 11:30-19:00 | 7.0 | LIGGGHTSのサンプルの中身を簡単な問題に書き換えてみる(物性値、壁の形状stl、球の投入速度？、等ができた)。salome-mecaをいじる | |
| 5/13(水) | 13:30-18:00 | 4.0 | 柴田良一『「salome-meca」ではじめる構造解析』のp.121までやる | |
| 5/14(木) | 9:30-12:00 | 2.0 | 柴田良一『「salome-meca」ではじめる構造解析』の付録以外をやる | |
| 5/14(木) | 14:30-17:30 | 3.0 | 榎戸正一『開発現場で活用できるオープンソースCAE�◆戮�case1をやる | |
| 5/15(金) | 10:30-12:00 | 1.5 | 榎戸正一『開発現場で活用できるオープンソースCAE�◆戮�case2-1をやる | |
| 5/15(金) | 13:30-17:30 | 2.0 | 榎戸正一『開発現場で活用できるオープンソースCAE�◆戮�case2-2・2-3・2-4をやる | |
| 5/18(月) | 10:00-12:00 | 2.0 | 榎戸正一『開発現場で活用できるオープンソースCAE� 戮�p.90まで (静的)をやる | |
| 5/19(火) | 13:30-18:00 | 4.0 | 榎戸正一『開発現場で活用できるオープンソースCAE� 戮鬚笋�(固有値解析も) | |
| 5/20(水) | 15:00-18:00 | 3.0 | 鈴木さんとgimagereaderをインストールし、日本語化してみる。lammpsのコマンド(liggghtsとコマンドが同じで情報量が多い？)を調べる | |
| 5/21(木) | 10:30-12:00 | 1.5 | lammpsのコマンドを眺める | |
| 5/21(木) | 14:30-18:00 | 3.5 | liggghtsのコマンドを見つけて眺める | |
| 5/22(金) | 10:30-12:00 | 1.5 | liggghtsの入力ファイルをいじる。チュートリアルなるものを見つける | |
| 5/22(金) | 13:30-18:00 | 4.0 | Ram Chand Raguelさんのチュートリアルをやったり、入力ファイルをいじったり | |
| 5/25(月) | 8:30-17:30 | 4.0 | チュートリアルの計算が出来ないのはおかしいと思い、LIGGGHTS-PUBLICとCFDEMcoupling-PUBLICをインストール | |
| 5/27(水) | 14:00-18:00 | 4.0 | lppをインストール、examplesで動作確認 | 河原さん |
| 5/28(木) | 8:00-12:00 | 3.5 | salomeで壁条件をつくり、そのstlファイルを使ってexamplesを計算 | |
| 5/28(木) | 15:30-17:30 | 0.5 | tutorial(序文)でliggghtsのコマンドを勉強する | |
| 5/29(金) | 8:00-17:30 | 7.5 | tutorial(bin flow)でliggghtsのコマンドを勉強する | |
| 5/30(土) | 10:00-12:00 | 2.0 | tutorial(bin flow)を計算・修正 | |
| 5月計 | 82.5 | |||
| 6/1(月) | 8:30-17:00 | 4.0 | tutorial(bin flow)を修正・計算、マニュアルのcontact model exampleを計算 | |
| 6/3(水) | 13:30-18:00 | 4.0 | マニュアルのcontact model exampleを改造して、一個の球を自由落下させた。目視と手計算を比較(反発係数OK、速度?) | |
| 6/4(木) | 8:00-12:00 | 4.0 | paraviewのグラフ化を試みる(ダメ) | |
| 6/4(木) | 14:30-17:30 | 3.0 | paraviewのチュートリアルを見つける。29/106までやってみる | |
| 6/8(月) | 13:30-17:30 | 3.0 | liggghtsの計算がparaviewで自由落下する球の時間とか反発係数が手計算と合っているか、の確認 | |
| 6/9(火) | 8:00-10:00 | 2.0 | liggghtsの計算が出力ファイルにどう書かれているかの確認(速度と力) | |
| 6/9(火) | 14:00-17:30 | 3.0 | liggghtsの出力ファイルより、粒子が変形しているときに働いている力を確認(途中) | |
| 6/10(水) | 14:00-17:30 | 3.5 | liggghtsの出力ファイルより、粒子が変形しているときに働いている力を確認(勉強中) | |
| 6/11(木) | 14:30-17:30 | 3.0 | liggghtsの出力ファイル(たくさん出てくる)を一つにまとめるプログラムをつくる(openfileの勉強中) | |
| 6/15(月) | 8:30-18:00 | 6.0 | liggghtsの出力ファイルを途中まで一つにまとめるプログラムをつくる(grepで解決と思いきや順番が変になるので) | 後藤さん |
| 6/16(火) | 13:30-18:00 | 4.0 | gnuplotについて調べる。gnuplot用にプログラム(6/15)を改造して、gnuplotで自由落下する球をグラフ化する | |
| 6/17(水) | 13:30-17:30 | 4.0 | gnuplotをいじる。liggghtsマニュアルから座標指定のコマンドを探す | |
| 6/18(木) | 9:00-12:00 | 3.0 | liggghtsについてきたtutorials publicの中で、座標指定しているものがないかを計算してみて探す | |
| 6/18(木) | 13:30-17:30 | 4.0 | liggghtsのマニュアルから座標指定のコマンドを探す | |
| 6/19(金) | 13:30-19:00 | 5.0 | stl→freecad×salomeの組み合わせからメッシュを切る方法を探す | |
| 6/20(土) | 11:00-18:00 | 4.0 | daiyastl.f90をfreecad経由でsalome-mecaにいれ、solid要素として厚みを持たせて解いてみる | |
| 6/22(月) | 8:30-18:00 | 5.0 | liggghtsのマニュアルを眺めたり、fortranのwrite文のformatを調べたり | |
| 6/23(火) | 13:30-17:30 | 4.0 | 複数の粒子の挙動をgnuplot用にまとめるプログラムをつくる | |
| 6/24(水) | 14:00-18:00 | 4.0 | 入力ファイルに坂(stlファイルでつくった)を読み込ませ、一個の球を転がす。坂の方向の移動量とか速度をグラフ化 | |
| 6/25(木) | 9:30-17:30 | 5.0 | 摩擦係数を変化させ、坂で球を転がす。摩擦係数をある程度大きくすると剛体の運動に従う(小さい時は質点に同じ)ことを確認 | |
| 6/26(金) | 9:00-18:30 | 6.5 | 球の自由落下と坂を転がした結果を日誌に書いたり、日誌を整理したり | |
| 6/29(月) | 10:30-17:30 | 4.0 | wikiの数式を勉強し、日誌を書き直す。liggghtsで物体をつくる方法を調べる | |
| 6/30(火) | 13:30-17:30 | 4.0 | 図書館に入れてもらった、林真『はじめてのParaView?』を見ながらparaviewをいじったり、鈴木さんとcode-asterの直交異方性をやってみたり | |
| 6月計 | 82.0 | |||
| 7/2(木) | 10:00-18:00 | 5.0 | code-asterの直交異方性の軸LTNがどうなっているのか試す。liggghtsの塊を入れるコマンドを調べる | |
| 7/3(金) | 13:00-17:30 | 3.0 | liggghtsの座標指定したり、塊を入れたりするコマンドを探す | |
| 7/6(月) | 8:00-18:00 | 6.0 | commファイルの出力をいじったり、熱応力の例題をやってみたり | |
| 7/7(火) | 8:00-9:00 | 1.0 | 熱応力の例題をやってみるも出来ず | |
| 7/8(水) | 14:00-18:00 | 4.0 | ドミノ倒しのシミュレーションを調べる(blenderとか自作とか)、商用(liggghts-premium)でなくても(liggghts-publicでも)コマンドrigidで塊を作れるかも | |
| 7/9(木) | 10:00-18:00 | 4.0 | liggghtsでread_data→group→rigidの順で塊を作るために、read_dataの書式を調べる | 7 |
| 7/10(金) | 13:00-18:00 | 4.0 | liggghtsのinsert/packとかfix growとかをいじったり、read_dumpについて調べたり | |
| 7/13(月) | 9:00-18:00 | 4.0 | Bullet Physics(Blenderの物理演算エンジン)、ODE、Newtonとかについて調べる、石積の資料(付箋:なくさないように ごとう)を見つけて見る | |
| 7/14(火) | 13:00-18:00 | 4.0 | bullet physicsとかphysXとかについて調べたり、blenderをいじってみたり | |
| 7/15(水) | 13:00-17:00 | 3.0 | liggghtsのread_dumpの使い方を調べて色々と試す | |
| 7/16(木) | 10:00-16:00 | 4.0 | liggghtsのread_dataの使い方を調べて色々と試す | |
| 7/17(金) | 10:00-17:00 | 6.0 | 3Dプリンターwikiを確認したり、LaTeXを少しやってみたり | |
| 7/18(土) | 14:00-17:00 | 3.0 | githubにliggghtsのexamplesがあるのを見つける | |
| 7/19(日) | 13:00-17:00 | 4.0 | githubのliggghtsのexamplesとかbenchmarkとかを色々と見たり、試したり | |
| 7/21(火) | 14:00-17:30 | 3.0 | githubのliggghtsのexamplesを試す | |
| 7/22(水) | 13:00-22:00 | 6.0 | 個人のPCにUbuntu15.04を入れてみる。liggghtsやblenderは使えるようになったものの、salomeやcode-asterのインストールができない | |
| 7/23(木) | 14:00-17:30 | 3.0 | blenderの物理エンジンをaddon(?)しなくても、物理計算ができることを自由落下問題で確認。salomeでドミノをつくりstlファイルをblenderにインポートして物理計算してみるが、動きが変(操作や設定に問題あり?) | |
| 7/24(金) | 13:00-18:00 | 3.0 | オープンキャンパスの振り子をプログラムを書き換えて、計算する | |
| 7/25(土) | 15:00-22:00 | 4.0 | 個人のPCにCAELinuxをいれようとする。Ubuntu15.04ベースだと無理そうだと確認 | |
| 7/26(日) | 8:30-14:00 | 3.0 | 個人のPCにUbuntu12.04.04とCAELinux2013を入れてみる | |
| 7/27(月) | 10:00-17:30 | 6.0 | 個人のPCにubuntu12.04の公式パッケージではない新しいバージョンのblenderやliggghts、paraviewをいれて動作確認する。振り子の調整 | |
| 7/28(火) | 14:00-23:30 | 6.0 | 振り子の調整、pythonでblenderの制御する方法を試す。blenderで、円柱と球を坂で転がしてみる | |
| 7/29(水) | 11:30-25:00 | 7.0 | 落下する球を作り直す(減衰機能を0にしないと空気抵抗のようなものが働き手計算と合わない)。salomeで円筒・球殻をつくり、blenderにインポートして転がす | |
| 7/31(金) | 15:00-18:00 | 3.0 | ドミノ倒しを作る、asterの座屈に挑戦(エラー)、bullet physicsの原理について書かれたページを探す | 斉藤さん |
| 7月計 | 99.0 | |||
| 8/4(火) | 11:00-17:30 | 5.0 | code-asterで自重でたわむ梁を静的に求めた | |
| 8/5(水) | 13:30-19:00 | 5.0 | code-asterで自重でたわむ梁の過渡応答解析をした | |
| 8/6(木) | 14:00-19:00 | 5.0 | salomeで石積みのモデルを作り、静的な解析をするも失敗。多分、石同士が線で接触(要素の線を共有)しているのが原因？ | |
| 8/7(金) | 14:00-18:00 | 2.0 | salomeで石積みのモデル(石同士が少し離れた)を作り、自重による沈下を静的に解析をするも失敗。多分、隙間のメッシュが粗すぎるから？ | |
| 8/9(日) | 13:00-18:00 | 4.0 | salomeで石積みのモデル(石の周りだけsubmeshで細かくする)を作り、自重による沈下を静的に解析をするも失敗。 | 斎藤さん |
| 8/20(木) | 19:00-23:00 | 3.0 | 拘束等の確認のために台形柱のみで解析をする。これは出来たので、拘束条件はOK。しかし、台形柱の石の部分を繰り抜いた複雑な形だとダメ。 | |
| 8/21(金) | 13:00-16:00 | 3.0 | 石の部分を粗めに繰り抜いた台形柱を解析する。これは出来たが、石をはめ込むとダメ。ヤング率を同じにしてもダメ。 | 斉藤さん |
| 8/24(月) | 14:00-24:00 | 8.0 | 石積み+裏込め(+背面土)で互いのヤング率を変えて静的と動的に解析する。繰り抜いたところをパーティションしないとダメ。 | |
| 8/25(火) | 14:00-23:30 | 7.0 | salomeからblenderに石積みをインポートして、soft bodyで裏込めを再現してみる。asterで変形したようになる設定をsoft bodyで色々試すも、思うような柔らかい物体にならない。多分、メッシュが原因。 | |
| 8/26(水) | 10:00-17:30 | 6.0 | 裏込めの応力歪曲線を設定して、asterで石積みの塑性変形を解いてみる。インパルス応答と静的自重で比較をしてみるも、結果は変わらず。多分、荷重が弾性域だからか設定のミスか？ | |
| 8/27(木) | 8:00-16:00 | 2.0 | 木ごころ橋(山形県金山町)やら、いくつかの橋を見学。 | |
| 8/28(金) | 10:00-20:00 | 7.0 | 斉藤さんのプログラミング講座(ccx用の矩形メッシュ分割)を受ける。asterの塑性解析について調べる。 | 斉藤さん |
| 8/29(土) | 11:30-15:00 | 3.0 | 斉藤さんのプログラミング講座(ccxの拘束条件のカード)を受ける。ccxの拘束条件以外のカードとかを調べる。 | 斉藤さん |
| 8/30(日) | 21:00-23:00 | 2.0 | 斉藤さんのプログラムを眺めたり、いじったり。 | |
| 8/31(月) | 10:00-19:00 | 7.0 | 斉藤さんのプログラミング講座(ccx用のL字・箱型メッシュ分割)を受ける。 | 斉藤さん |
| 8月計 | 69.0 | |||
| 9/1(火) | 13:00-18:00 | 5.0 | ccx用の箱桁をつって片持ち梁として計算。 | |
| 9/2(水) | 13:00-16:00 | 3.0 | blenderでsoft bodyの計算間隔を短くしてみるも、時間だけかかって精度はあんまり変わらない。 | |
| 9/3(木) | 13:00-16:00 | 3.0 | blenderで内部に細かなメッシュを生成するのが面倒なので、soft bodyのself colligionを小さくしてみれば、メッシュを細かくしたようになるのではとやってみるも、自分自身との衝突判定が正確になるだけ。 | |
| 9/7(月) | 10:00-16:00 | 4.0 | blenderのsoft bodyで内部にメッシュを生成するために、salomeのmesh→freecad→blenderを繰り返す | |
| 9/8(火) | 10:00-16:00 | 4.0 | Gmshをいじっていたらasterのmedファイルがインポートできることが分かり、rmedファイル(計算結果)をいろいろなテキストファイルにエクスポートして、眺める。 | |
| 9/9(水) | 9:00-16:00 | 4.0 | 様々なファイル形式に対応しており、ccxとの連携にも便利(inpファイルも作れる)だと思い、freecadをいじる。freecadの画面では解けてはいない。asterのBarre要素を試す。 | |
| 9/10(木) | 9:00-19:00 | 6.0 | blenderのsoft bodyで裏込め＋背面土内部にメッシュを作るもこれといって変化なし。 | |
| 9/11(金) | 11:00-16:00 | 3.0 | blenderのfluidを試す | |
| 9/16(水) | 8:00-19:30 | 8.0 | blender --pythonで剛体を積み上げる。ソフトボディも積み上げる。 | |
| 9/17(木) | 15:00-21:00 | 5.0 | blender --pythonでソフトボディを積み上げる。collisionでソフトボディ同士がぶつかるようにして、soft bodyで衝突判定をuse_face_collisionにすると、いい感じに柔らかい土っぽくなる。 | |
| 9/18(金) | 13:00-19:00 | 5.0 | blender --pythonでuse_face_collisionを入れようとするが、出来そうにない。 | |
| 9/22(火) | 20:00-24:00 | 3.0 | blender --pythonで(use_face_collisionを入れずに)、300*400*400の立方体要素をソフトボディで作ってみる。計算が終わらない。 | |
| 9/23(水) | 10:00-18:00 | 4.0 | blenderのrigidで直方体要素を使って土を表現できるようにpythonプログラムをつくり、色々な条件を試す。 | |
| 9/24(木) | 13:00-18:00 | 5.0 | 夏休み明け発表用のスライドを作る。9/23の計算が途中だが諦める。 | |
| 9/25(金) | 13:00-18:00 | 5.0 | 夏休み明け発表用のスライドを作る。とりあえず完成。 | |
| 9/26(土) | 13:30-25:00 | 7.0 | 夏休み明け発表用のTeXを作る。Xfigで図を作ったり、epsファイルを作ったり。 | |
| 9/27(日) | 9:00-22:00 | 9.0 | 夏休み明け発表用のTeXを作る。図表と文章は完成。あとは書式を整えるだけ。 | |
| 9/28(月) | 10:30-16:00 | 4.0 | 夏休み明け発表用のTeXを作る。とりあえず完成。発表練習をして、スライドとかを修正する。 | |
| 9/29(火) | 10:30-12:00 | 1.0 | 中間発表の練習 | |
| 9/30(水) | 13:00-26:00 | 8.0 | pythonで直方体の石を段違いに積み上げたり、壁を生成したりするプログラムを作る。 | |
| 9月計 | 96.0 | |||
| 10/1(木) | 14:00-22:00 | 6.0 | pythonで裏込め(剛体球体)を単純立方格子状に充填するプログラムを作る。 | |
| 10/2(金) | 16:00-18:00 | 2.0 | pythonで球体packingプログラムを作り、試しにobjファイルにエクスポート。 | |
| 10/4(日) | 11:00-16:00 | 4.0 | objファイルにrigidの設定を付け加えるプログラムを作ろうとするも、インポートまでしかできず。 | |
| 10/5(月) | 10:00-16:30 | 5.0 | 10/4と同じ作業をするも、進展なし。 | |
| 10/6(火) | 9:00-18:00 | 4.0 | blenderのsoftbodyをいじってみる。salomeで切ったメッシュをfreecad経由で色んな形式にしてblenderにエクスポートしてみるも、内部にメッシュが生成されないせいで動きが変。 | |
| 10/7(水) | 13:30-17:30 | 4.0 | 背面土を詰め込んだ後、それまで固定していた石を手動で動く設定に変更した。お試し用なのでモデルは雑。保存のやり方がわからず。 | |
| 10/8(木) | 13:00-16:00 | 3.0 | salomeからblenderにメッシュをインポートさせても内部にメッシュが生成されないので、手動で内部にメッシュを作る方法を試しみた(soft body用に) | |
| 10/9(金) | 10:00-18:30 | 3.0 | Elmerを使ってみるが、解析は出来ず。Blenderの保存の仕方を色々と試すが、物理演算結果がリセットされてしまう。 | |
| 10/10(土) | 11:00-17:30 | 3.0 | Gmshを使用してみる。code-asterの結果ファイルrmedを開いて、色々な形式に保存してみた。txtやposファイルとかは一応テキストファイルっぽい。 | |
| 10/13(火) | 10:30-18:30 | 5.0 | blenderの流体機能を使ってみる。動粘性係数を一番高い10×10$^10$m/s$^2$に設定してみるも、固体っぽくない。uv_sphereをico_sphereに変えてみる。 | |
| 10/14(水) | 10:00-23:00 | 7.0 | 大久保『石積の秘法とその解説』を読む。本を参考に新しく石のモデルを作る(石＋裏込)。ジオメトリーのみ。 | |
| 10/15(木) | 13:30-23:30 | 6.0 | 新しく作ったモデルをsalomeでメッシュに切り、asterで自重を載荷して静的解析。blenderでは、pythonでプログラムを作ってみた。日誌整理。 | |
| 10/16(金) | 10:00-17:00 | 3.0 | blenderで新しいモデルをpackingする。裏込は楕円球に変えてみた。計算中。動的解析用commファイルを作って解析。 | |
| 10/19(月) | 13:00-17:00 | 3.0 | 「fem」、「有限要素」、「ゲームエンジン」、「物理エンジン」とかのワードで色々調べる。 | |
| 10/20(火) | 13:00-23:00 | 7.0 | blenderで石積の計算してみるも上手く行かないので、asterの接触の例題をやってみるも計算失敗。 | |
| 10/21(水) | 11:00-17:00 | 3.0 | caelinuxのインストール。blender用のmesher "swiftsnap"を試そうとして失敗。asterの接触の例題その1をやってみるが出来ない。 | |
| 10/22(木) | 17:30-19:30 | 1.0 | 10/21の続きをやる。 | |
| 10/23(金) | 8:30-18:30 | 6.0 | asterで接触の問題をやってみる。寒風の文献をめくる。 | 後藤さん |
| 10/24(土) | 13:30-17:30 | 3.0 | blender --pythonで寒風の資料みたいなのを作って計算してみる。球体の半径を乱数で生成するプログラムを作る。 | |
| 10/25(日) | 13:30-17:30 | 3.0 | 10/24のを半径乱数verでも計算してみる。ccx_2.5/test/contact*.inpを計算してみるも、よくわからない。 | |
| 10/26(月) | 11:30-17:30 | 4.0 | 半径乱数verの半径とかを調整(側面の摩擦係数の影響大)。asterでウェブ上にある接触解析(計算途中)と塑性解析をやってみる。 | |
| 10/27(火) | 13:00-17:30 | 4.0 | 試しにやった接触解析が途中で終わっていたのでやり直す。science directを漁る | |
| 10/28(水) | 9:30-19:00 | 5.0 | 先日解けなかったweb上にあったasterでの解析例を、バージョンをSTA10から11に変えて計算してみるが、ダメだった。2D接触の物性を塑性に書き換えて、計算してみる。science directを漁る。山内さんと折り目をつくる。 | |
| 10/29(木) | 10:30-16:30 | 5.0 | science directでphysics、fem、hybrid、multi、fusion...とかを調べる。山内さんと折り目をつくる。 | |
| 10/30(金) | 16:00-18:30 | 2.0 | 滝田さんのunvファイルの体積求めるプログラムの復習ついでに、実はメッシュファイルの体積を求める機能がsalomeかどこかにあるのではないかと、色々探してみる。 | |
| 10月計 | 101.0 | |||
| 11/2(月) | 13:00-17:00 | 4.0 | blenderで片持ち梁を作るために、頂点を固定する方法を探すが、見つからない。 | |
| 11/4(水) | 13:00-16:00 | 3.0 | 昨日につづき、blenderで片持ち梁を作るために、頂点を固定する方法を探すが、見つからない。 | |
| 11/5(木) | 10:00-17:30 | 6.0 | 棚にあった桝谷浩『落石覆工の落石による衝撃力に関する研究』をめくっていると、自作DEMについての記述を発見。fortranで2D円筒要素の衝突の組み合わせと初期条件の読み込みプログラムを作ってみる。 | |
| 11/6(金) | 16:00-18:30 | 2.0 | 創造工房の来週までの課題をasterと手計算で解いて、メッシュとか単位とかを確認する。初期設定のメッシュ長さだと明らかに粗い。 | |
| 11/8(日) | 15:00-23:00 | 6.0 | 桝谷先生の論文をわりと真面目に読む。円筒要素DEMプログラムをつくる。加速度、速度、変位増分の積分とか、についてかく。 | |
| 11/9(月) | 11:00-17:30 | 5.0 | 円筒要素DEMプログラムをつくる。要素間の接触力についてかく。dt後にどうなっているのかの確認中。 | |
| 11/10(火) | 9:00-17:30 | 7.0 | 円筒要素DEMプログラムをつくる。各ステップでの円筒中心座標をgnuplot用ファイルに書き込むようにする。同時におかしな所を確認+修正中。円筒同士の衝突判定と自由落下はOK。接触力が変。 | |
| 11/11(水) | 9:00-17:30 | 7.0 | 円筒要素DEMプログラムをつくる。出力をparaview用のvtkフォーマットにし、動画として見ることで、計算結果の確認をする。 | |
| 11/12(木) | 13:30-17:30 | 4.0 | 円筒要素DEMプログラムをつくる。円筒要素と床および壁との接触についてをかく。接触力が変(そのうち直す)。 | |
| 11/13(金) | 9:30-17:30 | 3.0 | 円筒要素DEMプログラムをつくる。床と要素の接触を修正。blenderで壁と梁を結合してみるが、どちらかの属性にひとまとめになってしまう。 | |
| 11/16(月) | 16:00-17:30 | 1.0 | blenderのbone機能を調べるが、物理演算をしていない模様。固定端をつくるには、hookという機能が使えそう。 | |
| 11/17(火) | 10:00-17:30 | 5.0 | blenderのhookという機能を試してみるが、なぜか固定できない。 | |
| 11/18(水) | 13:30-17:30 | 4.0 | blenderでラーメン構造を作ってみるが、1つ1つの辺は曲がらず、伸び縮みのみする模様。柱部のみが上下に伸び縮みする。かといって、梁部分の辺を細かく切り分けると、表面が波打ってしまう。 | |
| 11/19(木) | 13:30-17:30 | 4.0 | blenderでラーメン構造の梁部分が波打たないように、中に辺をたくさんつくってみるが、そこばかり重くなってしまい、柱部分がひっくり返る。また、辺同士が重複するなど、思うように辺を作れない。 | |
| 11/24(火) | 9:30-16:30 | 6.0 | DEMプログラムの各減衰係数、逓減率を計算する(理解していない)。佐藤佑さんの実験の手伝い。 | |
| 11/25(水) | 10:00-17:00 | 6.0 | DEMプログラムの法線方向のバネ定数を工夫する。ヘルツの接触から導けないか、とか。 | |
| 11/26(木) | 14:00-17:00 | 3.0 | blenderで作ったラーメンの柱のすり抜けを解消。梁の内部に辺を生成するも、辺が重複したりと、やっかい。柱の内部にも辺を生成してみる。 | |
| 11/27(金) | 9:30-17:00 | 3.0 | blenderで作ったラーメンは自重しかかけていないのに、変形が大きすぎるので、曲げ強さを大きくしてみた。 | |
| 11/28(土) | 14:00-17:00 | 2.0 | asterの結果をテキストファイルにしたいので、commファイルをいじって、format='GMSH'というのを試してみた。エラーはでないが、結果ファイルが出来上がらない。 | |
| 11/30(月) | 13:00-17:00 | 4.0 | commファイルをいじって、format='IDEAS'にしてみたら、asterの結果ファイルがunv形式になった。libreofficeやadobe readerで使えそうな動画形式やら動画作成ソフトやらを探す。 | |
| 11月計 | 84.0 | |||
| 12/1(火) | 9:30-17:00 | 4.0 | imagemagickのconvertコマンドを色々と試す。mpeg形式ならimagemagick→libreofficeに出来そう。スクリーンショットを連続で画像に保存する方法を探す。 | |
| 12/2(水) | 10:00-20:30 | 5.0 | fortranでimagemagickコマンドを書き込みスクリーンショットをたくさん撮る。それをgif化してlibreofficeで観た。 | |
| 12/3(木) | 9:30-17:30 | 5.0 | 12/2のfortranをスクリーンショットをウィンドウ指定にし、gif化する範囲を指定した。liggghts・自作DEMの計算結果を、直接paraviewから画像にたくさんはき出してgif化。 | |
| 12/4(金) | 14:30-17:30 | 2.0 | 心臓シミュレーションと仕分けの話について調べる。 | |
| 12/6(日) | 16:00-17:00 | 1.0 | 横ねじれ座屈について調べる。 | |
| 12/7(月) | 10:00-18:00 | 6.0 | 3Dプリンターの印刷について行くだけの人をやる。面白そうな動画ということで、石積のFEM動的解析のgifとかをつくる。これは、commファイルを見なおして、再計算した。 | |
| 12/8(火) | 10:30-17:00 | 5.0 | blenderの動画をつくるための計算をする。並行して、画像保存機能について、いろいろ試す。カメラやランプの向き配置をためしたり、壁をやや透明にしたり、影を消したりしてみる。 | |
| 12/9(水) | 10:00-17:00 | 5.0 | blenderの計算継続中。伊藤さんの曲げ試験(練習)を手伝う。 | |
| 12/10(木) | 10:30-17;00 | 5.0 | blenderの計算中PCが動かなくなるので辞めて、昨日一昨日とは別のblenderを計算(途中)。伊藤さんの曲げ試験(本番)を手伝う。 | |
| 12/11(金) | 10:30-17:00 | 3.0 | 昨日計算をしたblenderをgif化。表面の線を細かくして計算。 | |
| 12/12(土) | 10:30-17:00 | 5.0 | asterで石積のモーダル解析の計算しなおし。昨日計算したblenderをgif化。それだと、柔らかすぎたのでパラメータを変えて計算。 | |
| 12/14(月) | 10:30-17:00 | 3.0 | 支部発表申し込み。blenderのパラメータを変えて計算(途中)。 | |
| 12/15(火) | 10:00-17:30 | 5.0 | blenderのパラメータを変えて計算。動画化してみるもsoftbodyだと、うまく見えない。 | |
| 12/16(水) | 10:30-17:30 | 4.0 | 3Dプリンターの造形を見に行くだけの人をやる。blenderの剛体を半径とかを調整して計算(途中)。 | |
| 12/17(木) | 9:00-17:00 | 6.0 | blenderの剛体の計算(途中)。 | |
| 12/18(金) | 9:30-17:30 | 4.0 | 12/17のblenderをgif化。asterの石積をモーダル解析の密度を修正して、過渡応答解析の粘性減衰を変えて計算しなおし。中間発表のスライドを考える。 | |
| 12/19(土) | 10:30-17:30 | 5.0 | 中間発表のスライドに使うbullet physicsの検証用動画等をつくる。 | |
| 12/20(日) | 9:00-16:30 | 6.0 | 中間発表のスライドをつくる。完成。texを考える。 | |
| 12/21(月) | 9:00-18:00 | 8.0 | texをつくる。発表練習をする。時間オーバーする。texつくる。完成。 | |
| 12/22(火) | 8:30-20:00 | 5.0 | 発表練習をする。佐藤さんの引張試験を手伝う。スライドの画像を小さくするとファイルサイズが半分ぐらいになった。powerpointでodpファイルを開くと画像が暗くなった。しかし、pptxファイルに保存したものを読み込むよりはかなりまし。河原さんの支部TeXを少し書き換える。 | 後藤さん・斉藤さん・河原さん |
| 12月計 | 92.0 | |||
| 1/7(木) | 14:00-17:00 | 3.0 | blenderのキーフレームという機能を発見。使えるかもしれないので、調べて試す。 | |
| 1/8(金) | 10:30-17:00 | 4.0 | impress→powerpointで読める形式(.ppt)を確認。blenderのキーフレームをいじる。 | |
| 1/9(土) | 13:00-17:00 | 3.0 | 支部のTeXを考える。「ccxとasterの違い」にまつわる質問に備えて、というのもあるが、純粋に気になったのでasterの結果ファイル(unv)の中身を調べてまとめる。 | |
| 1/10(日) | 13:00-17:00 | 4.0 | 昨日、unvファイル上でグループ情報がどう書かれているのかを確認しいなかったので、引き続き調べた。シェル要素を試しに解いてみるも、失敗。 | |
| 1/11(月) | 13:00-17:00 | 4.0 | astkの勉強。asterでシェル要素を解く。指定した線で回転しないようにはできるが、ヒンジにするやり方が分からない。asterでは、荷重の代わりに初めから変位を与えられる模様。モーメント力も入れられる模様 | |
| 1/12(火) | 10:30-17:30 | 5.0 | asterで高橋さんの木組の接触解析を試す。commファイルは問題なく完成。収束しないのは、ニュートン法とやらの設定の問題だと思う。 | |
| 1/13(水) | 10:30-17:00 | 5.0 | asterの接触に摩擦係数を入れてみる。commファイルは問題なく完成。メッシュを細かくしたり、時間の間隔を見直すも、結果が出ない。他にもペナルティ法(収束しやすいらしい)とかラグランジェ法とかがあるらしいので、そのうちに試す。 | |
| 1/15(金) | 10:00-17:00 | 5.0 | 創造工房のスクリーンショットをとる、藤田さんのを計算する。卒論概要・スライドを修正。 | |
| 1/16(土) | 21:00−23:00 | 2.0 | CAE懇話会の会誌第11号を読む。 | |
| 1/17(日) | 16:00−23:00 | 4.0 | 一昨年の卒論発表を途中まで見る。CAE懇話会の会誌第19号を読む。 | |
| 1/18(月) | 21:00-24:00 | 3.0 | 支部概要を大体手書きで書く。 | |
| 1/19(火) | 10:30-18:00 | 6.0 | 支部概要をTeXで書く。DEXCS-Salomeとかについて調べる。 | |
| 1/20(水) | 10:30-19:00 | 7.0 | 後藤さんによる支部及び卒論概要の添削。その後、丁度2ページにならなかったのを無理やり合わせた。佐藤さんの実験手伝い。 | 後藤さん・斉藤さん・河原さん |
| 1/21(木) | 10:30-18:30 | 6.0 | 後藤さんによる支部概要の添削。文中の式とか変なところを直す。DEXCS-SalomeをDVD起動して使ってみる。でも、エラーはフランス語のまま。asterのmessファイルの英語化方法を調べると、STA11.3以降ならできるとか(いつもはSTA10)。 | 後藤さん |
| 1/22(金) | 10:30-18:00 | 5.0 | 支部概要の直しと提出。支部の宿の予約など。スライドの修正。 | |
| 1/23(土) | 13:00-17:00 | 3.0 | 支部の質問に備えて、ccxの理解を深めようと斉藤さんのccxプログラムを実行してみる。 | |
| 1/25(月) | 11:00-17:30 | 5.0 | 高橋さん・伊藤さんと発表練習をする。スライドの修正。gifが暗かったので動画の作り直し。 | |
| 1/26(火) | 9:00-18:00 | 7.0 | 後藤さんによる発表指導を受け、伊藤さんのも見る。スライドの指摘されたところを直す。主にgimpで画像を見やすくする。 | 後藤さん |
| 1/27(水) | 9:00-18:00 | 7.0 | 鈴木さんの発表を見る。ImpressでTeXを書けるようにする。スライドの修正つき、動画を見て、salomeで物体の色を変えたり、線を太くしたり、スカラーバーの設定をいじってスクリーンショットを撮る。 | |
| 1/28(木) | 9:00-18:30 | 7.0 | 概要の連絡先を消して、濃く印刷。スライドの修正につき、動画を見て、gifを作りなおしたり、femの画像を変えたり。発表練習をする。 | |
| 1/29(金) | 9:00-18:30 | 7.0 | 奥村さん、山内さんらの発表を見る。斉藤さん、鈴木さんと発表練習をする。日誌の整理。伊藤さん、高橋さん、菊地さんと発表練習。 | 斉藤さん |
| 1/30(土) | ||||
| 1/31(日) | ||||
| 1月計 | ||||
| 2/1(月) | ||||
| 総計 |
| スクロール | 拡大縮小 |
| スクロールでドラッグ | 見る角度を変える |
| shft + スクロールでドラッグ | 見る位置を移動 |
| 1 or 3 or 7 | xz・yz・xyとかの正面に視点を移動 |
| 5 or 9 | 選択した物体に接近 or 選択した物体の裏側に視点を移動 |
| 2 or 8 | 下or上に視点を移動 |
| 4 or 6 | 左or右に視点を移動 |
| aa | 全物体を選択 |
| z | 物体のedgeのみの表示にする。 |
| g (+ x or y or z + 移動量) | 物体の移動。座標と移動量を指定できる。 |
| shift + d | 物体のコピー |
| s (+ 倍率) | 物体の大きさ変更。数値を入力する。 |
| shift + a | 物体の追加 |
| shift + c | 原点を選択 |
| del | 物体の削除 |
| esc | キャンセル |
| alt + a | 動画の再生 |
| esc | 動画の終了 |
| Rigid Body | Type | Active | 動く物体 |
| Passive | 動かない物体 | ||
| Dynamic | on | 動く物体 | |
| Mass | 重量(kg)。左のPhysics→Calculate Mass で、物体を選択すると自動で計算できる。あるいは、その下の欄に密度の入力を入力してもいい。 | ||
| Rigid Body Colligions | Shape | Mesh | 衝突判定が三角形ポリゴンのメッシュ。困ったら、なんでもこれ。 |
| Convex Hull | 衝突判定が凸包。真っ平らな床とか壁とかはMeshでなくてもこれでいい。 |
| Surface Response | Friction | 摩擦係数(0.000~1.000) | |
| Bounciness | 反発係数(0.000~1.000) | ||
| Sensitivity(Colligion Margin) | Margin | 衝突判定の距離(0.001~1.000) | |
| Rigid Body Dynamics | Damping | Tlanslation | 移動の減衰係数(0.000~1.000) |
| Rotation | 回転の減衰係数(0.000~1.000) |
$ blender --python filename.py&
$ blender filename1.blend --python filename2.py&
$ cd /opt $ sudo mkdir LIGGGHTS $ cd LIGGGHTS $ git clone git://github.com/CFDEMproject/LIGGGHTS-PUBLIC.git
$ cd /opt/LIGGGHTS/LIGGGHTS-PUBLIC/src/MAKE $ cp Makefile.linux Makefile.linux.org $ vi Makefile.linux
| 項目 | 変更前 | 変更後 | 補足説明 |
| CC= | icc | gcc | コンパイラ |
| LINK= | icc | gcc | リンカー |
| MPI_INC= | -DMPICH_IGNORE_CXX_SEEK | -DMPICH_IGNORE_CXX_SEEK -l/usr/lib/openmpi/include | MPIのincludeの設定 |
| MPI_LIB= | -lmpich -lpthread | -L/usr/lib -lmpi -L/usr/lib/openmpi/lib -lmpi_cxx_lpthread | MPIのlibの設定 |
$ cd /opt/LIGGGHTS/LIGGGHTS-PUBLIC/src $ make fedora_fpic $ make makelib $ make -f Makefile.lib fedora_fpicliblmp_fedora_fpic.aという実行ファイルが出来ればOK。
$ ln -s lmp_fedora_fpic liggghts $ ln -s lmp_fedora_fpic liblmp.a
$ vi ~/.bashrcのどこかしらに、以下の文を追加。人によってはパスは違うかも。
alias lpp="python $CFDEM_LPP_DIR/lpp.py" alias pizza="python $CFDEM_PIZZA_DIR/pizza.py"
export PATH=/opt/LIGGGHTS/LIGGGHTS-PUBLIC/src:$PATH export PATH=/opt/LIGGGHTS/LPP/src:$PATH export CFDEM_LPP_DIR=/opt/LIGGGHTS/LPP/src export CFDEM_PIZZA_DIR=$CFDEM_LPP_DIRLPPを使わない人は、export PATH〜LIGGGHTS-PUBLIC/srcの行だけでいいかも。それと、lpp.pyしか使っていないので、pizza.pyの行(alias pizza〜 と 〜CFDEM_PIZZA〜)もいりませんでした。
sudo apt-get install liggghtsとかで入んないですかね。足りないものはあるかもしれませんが。
$ sudo mkdir /opt/LIGGGHTS $ sudo chmod 777 /opt/LIGGGHTSすべてのユーザーが実行可能になるらしい。上記参考サイトを真似てやったけど、いらなかったかも。
$ mkdir /home/kouzou/2015/dem $ cp /opt/LIGGGHTS/LIGGGHTS-PUBLIC/examples/LIGGGHTS/Tutorials_public /home/kouzou/2015/demkouzou/のどこかにDEMの作業用のディレクトリを作って、そこにTutorials_public(サンプルが色々入っている)をコピーしてから使うことにした。
$ liggghts < in.filename
$ cd /opt/LIGGGHTS $ git clone git://github.com/CFDEMproject/LPP.git
$ lpp.py dump*.filename
| command | 可 | 不可 | 備考 | |
| fix | wall/gran model hertz tangential history | wall/gran/hertz/history | ||
| pair_style | gran model hertz tangential history | gran/hertz/history | ||
| fix | mesh/surface/stress | mesh/gran/stressanalysis | ??? |
上の図のように、球体を自由落下させた。球に働いている力は重力のみである。反発係数だけを変えたのが下の図である。
横軸が経過時間、縦軸が鉛直座標である。反発係数$e$を$\sqrt{\frac{1}{1}}=1.000,\sqrt{\frac{3}{4}}=0.866,\sqrt{\frac{1}{2}}=0.707,\sqrt{\frac{1}{4}}=0.500$に変化させて計算した。
上の図のように、球体を坂に置いた。球に働いている力は重力のみである。摩擦係数だけを変えたのが下の図である。
横軸が経過時間、縦軸が鉛直座標である。摩擦係数$\mu$のみを$0.00,0.05,0.10,0.30,0.50,1.00$に変化させて計算した。
| 摩擦片 | 摩擦面 | 摩擦係数 |
| 石 | 石 | 0.6~0.7 |
| 石 | 土 | 0.5(乾)~0.3(湿) |
| 木 | 石 | 0.4 |
$ cat dump* |grep 1\ 1> tatoeba
$ xwininfo |grep '^xwininfo: Window id:' | awk '{print $4}'
で、ウィンドウのidを調べて、実行。
$ convert moto*.jpg -crop 400x300+500+400 kiritori.jpg
$ convert -delay 10 kiritori*.jpg filename.gif
$ mogrify -resize 400 *jpg *gif
| 条件 | 4面体要素の種類 | 要素数 | たわみ(m) | 相対誤差(%) |
| 理論値 | --- | --- | 0.206 | --- |
| 点載荷 | 1次 | 8908 | 0.148 | 28.2 |
| 面載荷 | 1次 | 8908 | 0.142 | 31.1 |
| 点載荷 | 1次 | 75562 | 0.194 | 5.8 |
| 面載荷 | 1次 | 75562 | 0.193 | 6.3 |
| 点載荷 | 2次 | 8436 | 0.206 | 0.0 |
| 面載荷 | 2次 | 8436 | 0.206 | 0.0 |
| 拘束条件 | 載荷条件 | 要素の種類 | 要素数 | 節点数 | たわみ(mm) | 体積(mm$^3$) |
| 線 | 面(1cmのひもを想定) | 4面体2次 | 32705 | 55184 | 1.16 | 11276 |
| 初等梁理論 | 静的解析 | 過渡応答解析 | |
| 中央のたわみ(mm) | 4.0578 | 4.1960 | 4.1958 |
