近藤の卒論日誌

あつこ弁当 696枚×5円=3400円分のポイント†

　構造研の先輩・同期の皆さん、ありがとうございます。質費に充てました。

Salome-Meca2015.2のポストプロセッサ"ParaViS"†

　Salome-Meca2015.2をインストールしてみたところ、PostProがなくなっていた。(Code_AsterのMEDファイルを読めるようにParaviewを作り替えた?)ParaViSを使え、ということらしいので、FEMで使えそうな機能をメモする。

読み込む†

• ParaViSに行く(赤・緑・青 のボタン)
• ファイル→ParaViewファイルを開く→Code_Asterで計算した*.rmedファイルを選択して開く

変形図に応力のコンター図を重ね合わせる†

• パイプブラウザーのrmedファイルを選択したまま、フィルター→Common→Warp By Vector
• PropertiesのVectorsを「RESU____DEPL」(変位)に設定→Apply
• PropertiesのScale Factorは変形の倍率なので、1倍がそのまま。適度に大きくしたり小さくすると見やすくなる。→Apply
• Coloringを見たい応力の種類によって変える。
  • 例えば、ミーゼス応力なら、「RESU____SIEQ_NOEU」の「VMIS」と、いうふうに。

グラフ化（軸方向距離-変位or応力or...）†

• パイプブラウザーのrmedファイルを選択したまま、フィルター→Data Analysis→Plot Over Line
• PropertiesのPoint1とPoint2は、線分の起点と終点なので、適宜座標値を入力。これでグラフの横軸が決まる。
  • 例えば、梁の中立軸の両端の座標値とかを入力。
• DisplayのSeries Parametersを見たい変位や応力(例えば、荷重方向の変位なら"RESU____DEPL_DY"だけ)に選択。
• ViewのLeft Axis Lange(横軸)やBottom Axis Lange(縦軸)の値を変えて、グラフの範囲を適宜修正。
  • チェックを外すと、自動で調度良い範囲になる。

グラフを数値化†

• グラフの右上の「日」みたいなボタン→SpreadSheet View を押すと表が出てくる

csvファイルに出力†

• 表を選択したまま、ファイル→エクスポート→名前をつけて（csvしかない）ok
  • 形式を見てみると、gnuplotの入力データがfortranで簡単に作れそう。

ある断面の軸方向応力をベクトル表示†

• 「伊藤の卒論日誌」あたりにもあったような。やるべし。

断面の変位を平均する†

　そんな機能があってもいいではないかと、模索中。以下、暫定版。

• パイプブラウザーのrmedファイルを選択したまま、フィルター→Common→Slice
• 見たい断面の、原点座標(origin)を入力し、面の法線(nomal)を選択。
• パイプブラウザーのSliceを選択したまま、フィルター→Data Analysis→Plot Data
• 欲しい変位だけを選択。
• グラフ右上の「日」→SpreadSheet View→ファイル→エクスポート→csv
• LibreOfficeのCarcで open→"f(x)"(InsertのFunction Wizardでもいい)のAverageで平均化
  • fortranで計算してもいいけど

3D地図をstlファイルにダウンロード†

• &link(国土地理院のサイト,http://cyberjapandata.gsi.go.jp/3d/index.html)から、3D地図のstlファイルがダウンロードできる。これを3Dプリンターで造形するのが流行っているらしい。
• 下記の方法 blender → freecad → salome で遊ぶと楽しそう。
  • 面の数を1/10に減らして、ようやくsalomeで動いたが、クリエイトグループが大変かも？

細かすぎるstlファイルblenderで粗目にする†

• モデルの回転や移動は、blenderでやるとかなり動作が速いのでおすすめ。

• blender → ファイル → インポート"stl"
• モディファイ → 追加 → ポリゴン数削減
  • 例えば、国土地理院の3Dの富士山↓は元々、面の数が"133110"もあるので、細かすぎて使いにくい。
• その1：平面 → 角度制限 → 適度な角度に → 適用
  • 難しい
• その2：束ねる → 「三角面化」にチェック → 「比率」を適度に小さく → 適用
  • 楽に減らせる。面の数が1/10でも↓結構きれい。
  • 適宜、「その1」と「その2」を組み合わせる。
    • blenderでは、右上の「面:」に面の数が書かれている。
    • blenderでは、"z"でエッジが見れる。
• ファイル → エクスポート"stl"

stlファイルをsalomeにいれて使う†

freecad†

• stlファイルを開く
• partにて
  • 部品→メッシュから形状を作成する
    • 「縫い合わせのトレランス」は、何だかよくわかりません。
  • 部品→シェイプビルダー→シェルからソリッドを作成
• ソリッドを選択したまま……ファイル→エクスポート→BREP format(*.brep *.brp)でsave

salome†

• geometoryにて
  • ファイル→インポート→先ほどのbrepファイルをインポート
  • measures→basic propertiesで、volumeを確認して、値が正になっていればOK
  • volumeが負の値の時は、もう一度freecadに戻って、
    • ソリッドを選択したまま……部品→図形を反転
    • 反転したソリッドを選択したまま……ファイル→エクスポート→BREP format(*.brep *.brp)でsave
  • これをインポートすると、多分volumeが正の値になっている。
• あとは、普通にcreate groupして、meshに切れば、asterで解ける。

paraview の使い方†

インストール(CAELinuxには元々入っていますが)†

• ダウンロード
  • http://www.paraview.org/download/
• 解凍
  • tar xvzf ParaView-4.3.1-Linux-64bit-glibc-2.3.6.tar.gz
    • バージョンによって、新しいのだと5.0.0とか。
• 端末から連続したvtkファイルを開くには、
  • paraview filename..vtk&
    • paraviewをアイコンから開いて、openからファイルを指定するよりも楽。
  • 再生Modeは、Real Timeが速くておすすめ。

各種フィルター†

　フィルターなるものでデータを色々と見やすくできる。グラフとかもかけるらしい。 vtkフォーマットは分かりやすくて便利。

• 基本的には、
  • Filters→Alphabetical

球体をつくる†

• vtkファイルを押したまま、Glyph .

  Scalars	radius
  Glyph Type	Sphere
  		Radius	1.0
  		Theta Resolution	16	多いほうが球体の角が取れる
  Scalar mode	scalar
  	Set scale factor	1.0

力をスカラーで色付けしたい†

• 上記 球体の設定に加え、Solid Color を f に変更。
  • 右隣のMagunitudeは合成。X,Y,Zで成分別にできる。

速度ベクトルを描きたい†

• 新しく Glyph を追加

  vectors	v
  Glyph Type	Arrow
  Scalar mode	vector
  	Set scale factor	0.02		矢印の大きさ。適度に小さくして、調整。

• Solid Color を v に変更。
  • 右隣のMagunitudeは合成。X,Y,Zで成分別にできる。

色・文字の設定†

数値と色の帯(スカラーバー)†

• Toggle Color Legend visibility(左上のカラフルな帯のアイコン)
  • 位置とか大きさは、マウスの左クリックで調整できる。
  • Color Scale Editor(その右隣)→Color Lengend

    title	text	タイトルの名前
    	Font	times	好きな書体に
    		B	太字
    		I	斜め
    		S	影
    	size	20	文字サイズ
    	Color	黒とかに
    label	上記と同じ	帯の方

画面の背景†

• View Setting(Layout#1の下の電卓みたいなアイコン)
• Choose Color→白とかに

たくさんの静止画を保存†

• Start Time(開始時間)とEnd Time(終了時間)を設定して、南京錠マークをチェックしておくと、撮りたい範囲を指定できる。
• File→Save Animation

  Animation Duration(sec)	アニメのコマ数
  Resolution(pixels)	横幅・縦幅のピクセル値

gif形式に動画化†

• imagemagickより、
  • convert filename*.jpg filename.gif
• 画像ビューアーで見る。
  • eog filename.gif&

#########ここより上はただのメモ###########†

卒論のテーマ†

• 近藤：&link(LIGGGHTS,http://www.geocities.jp/penguinitis2002/study/LIGGGHTS/LIGGGHTS-example/LIGGGHTS-example.html)等のDEMツールを使って石積み護岸の安定性解析等できないかの基礎研究
  • DEM（個別要素法）を使ったことがないので、ためしに使ってみたい。
  • 石積み護岸に対して、なんらかの数値的アプローチをしてみたい。

• 石積みを例にしたFEMとDEMの比較(当初のテーマ)
  • FEM:積んだだけの石同士の間に柔らかい弾塑性要素を挟み、インパルスをいれる。
    • salome-mecaで出来そう。
  • DEM:積んだだけの石にインパルスをいれる。
    • liggghtsでは(機能的に)多角形要素が出来ないし、(近藤の能力的に?)正確な座標指定が出来ない。
    • では、よくあるドミノとかのアニメはどうやってる?　→　物理エンジン

テーマの修正とか目的とか方向性みたいなもの†

物理エンジンの特徴†

　物理エンジン(blenderのbullet)はDEMとかFEMとかと比べると、本当にそうなるのか、ちゃんと計算しているのか、というと、ちょっとあやしい。 それっぽく見せているだけなような気もする。 剛体の動きは結構いい線を行っているように見えるが、軟体(? soft body)などの内側に辺が無い物体の変形は、とくにあやしい。 しかし、物理エンジンは、比較的計算は速い(簡単なものならほとんど対話型)し、アニメーションをつくるための色んな機能が(ほとんど使いこなせていないぐらい)ある。

• 曲げを表現しようとすると、むしろ有限要素法よりもメッシュ分割・計算ともに遅い?→あんまり、メリットがない?

何かに生かせないか?†

　たとえ厳密に物理現象を再現していなくても、それっぽい災害とかのアニメーションをつくって見せることは、防災意識を促すといった観点から、それなりに意義がある。 アニメーションは専門知識とかがなくても、見るだけで直感的に何が起きているのかわかるので、そこがいいところ。

方向性としては……†

　まず、code-asterとblender、liggghtsとblenderで簡単な問題の比較検証をして、どれくらい厳密なのかを一応確認してみる。 次に、DEMとかFEMの計算結果と同じになるように、物理エンジンのパラメータを探す。 一度ちょうどいいパラメータを探しておけば、ちょっと違う構造で映像を見せたい時に、すぐにそれっぽい動きを再現することが出来る……のでは?

修正の修正?†

• 目的 　動画の作成に際して、3DCGソフトを使うものと思われるが、その1つのblenderではbulletという物理エンジンが組み込まれている。bulletでは、剛体の衝突判定や物体の変形を計算してくれるが、どれだけ厳密なのかは明らかでない。 構造物の主要な部分は従来どおりFEMやDEMで計算する必要があるが、それ以外のところを補うツールとして物理エンジンを利用できるのでは?

修正の修正の修正?†

• 目的とか背景とか 　シミュレーション計算結果の可視化方法の模索。 　計算結果は複雑かつ数が多いため、可視化しないと評価が難しいし、人に伝わらない。しかし、シミュレーション(模擬実験)を専門にしている人とビジュアライゼーション(可視化)を専門にしている人はそれぞれ違い、お互いがそれだけでも難しいことをやっている。そのため、この2つの連携は必ずしもうまくいかず、現状ではフィードバックを何度も繰り返しそれらしい動画をつくっている。 　こういった課題に対して、シミュレーションをやっている我々はどうしたらいいかを、色々と模索する。
  • こういった課題は&link(心臓シミュレーションと仕分けの例,http://www.sciement.com/blog/3dcg/ut-heart-2/1250/)に見られるように、なかなか理解されない。　
• 方針 　とりあえず、salomemecaやらparaviewやらblenderやらで面白そうな動画をつくる。
• シナリオ
  • まずは正攻法! 　asterとかで真面目に計算したのを、paraviewに入れてパラパラ漫画にする。計算結果に忠実だが、パラパラのコマ分の全てを計算するのは大変。
  • 足りないところは補う 　なので、ある程度はちゃんと計算をして、それを真似て3DCGソフトblenderの物理エンジンで補間できないか、という提案。
    • ただし、多分blenderにfemの結果とかの座標データはいれられないので、なんとなくで動画を作ることになるが、そもそもblender自体の操作も大変......。

目標†

一週間の目標とその達成度†

• 卒論本体を書く。
• 支部のスライドをつくる。
• ちょっとした案
  • blenderのキーフレームという機能が使えないか?
    • 物体の位置・向き・スケールを2コマ以上設定すると、途中のコマをblenderが補間してくれるという便利機能。
    • 柔らか部分についてasterの計算結果(座標)をいくつか取り出して、キーフレームを設定できないだろうか?
    • 物体をインポートして、消して、インポートして、消して、の繰り返しは出来ないかも。
• そのうちやりたいこと
  • asterで弾塑性解析をしたい。web上にあった例は解いてみたが、他のモデルに適用していない。
  • asterで自重の他に、横揺れのインパルスを入れたい。追従荷重というやつ?
  • 多角形DEMプログラムをつくりたい。2D円verは作ってみたが、色々とあやしいのでまだまだ修正の必要あり。

    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でグラフ化は成功。

物理エンジンの検証†

剛体運動†

反発係数†

• 反発係数を左から順に、\( \sqrt{0.00},\sqrt{0.25},\sqrt{0.50},\sqrt{0.75},\sqrt{1.00} \)とした。
• 手計算通りになる。自由落下だと元の高さの反発係数2乗倍の高さまで跳ね上がる。

摩擦係数†

• 摩擦係数を左から順に、\( 0.00,0.25,0.50,0.75,1.00 \)とした。
• 摩擦係数が低いほど、よく滑る。回転運動も再現している。

弾性変形†

• ラーメン構造を自重で変形させる。
  • 下記のようにパラメータの設定がとても曖昧。ヤング率やポアソン比ではない。
  • そこで、edgeという辺の数や入れ方で調整してみる。

edge:初期†

edge:表面たくさん†

• blenderの機能で自動細分化できる。

edge:内部も格子状にたくさん†

• 重なった点や線を除去しないと、物体が重なっておかしなことになるので注意！
• salomeで作ったメッシュはstlファイルにエクスポートしてblenderにインポートできるが、そもそもstlファイルは表面形状しか持っていない。なので、ずっと手作業で点とか線をつくっていたけど、プログラムから作ると分割数も変えられて、楽。
  • &link(naibu.py,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/blender/naibu.py)
    • "aa"で全選択。
    • 左側の「ツール」→「編集」→「統合」
    • 下の「オブジェクト」→「トランスフォーム」→「原点をジオメトリに移動」
    • 下の「オブジェクトモード」を「編集モード」にかえて、
    • 左側の「ツール」→「メッシュツール」→「重複頂点を削除」
    • あとは、好きなようにラーメンのソフトボディ、床のコリジョンを設定すべし。
• プログラムで作ると、動きが断然滑らか。手作業で作った上のは、どこかを間違えていたらしい。

石積†

salomeでモデルをつくる†

• geometoryにて、
  • 0.次の1～4までは、別のsalomeで作ってbrepファイルにエクスポートしたものをインポートすると動作が軽い。
    • ついでに、一部をstlファイルにエクスポートしておくと、blenderにインポート出来るので便利。
  • 1.vertex→edge→face→solid→translation等々で石段をつくり、compoundする。
  • 2.face→roatation→extrusion等々で裏込部分の平行四辺形をつくる。
  • 3.2の平行四辺形を石をツールにcutする。
  • 4.1の石をfaceにexplodeして、裏込との境界部分のみを選択しcompoundする。
  • 5.1の石と3の裏込をcompoundし、石・裏込部・下面・後面・左面・右面にそれぞれcreate groupする。
    • 下の画像は、境界条件を設定する4つの面
  • 6.5でまとめた物体を、4の境界面をツールにpartitonする。
  • 7.6でpartitionした物体を、5でcreate groupしたグループを使用して、再びcreate groupする。
  • partitionした物体を選択して、いつも通り四面体要素にcreate meshする。
    • これも、medファイルにエクスポートしておくと、asterで計算するとき、別のsalome-mecaで開けるので、動作が軽いかも。

asterで計算†

静的解析:自重で沈む†

• &link(sei.comm,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/salome-meca/sei/sei.comm) 　背面が2mmぐらい沈んでいる。
• 変位図……単位:m、スケール:100倍

動的弾性解析:途中で上下に自重の10倍のインパルスを入れる†

• &link(dou.comm,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/salome-meca/dou/dou.comm)
  • 静的解析とあんまり変わらない。裏込を弾塑性体にする必要あり? インパルスを横揺れに変える必要あり?
  • 変位図……単位:m,sec、スケール:100倍↑,粘性減衰:土0.04,石0.00002

動的弾性解析:途中で上下に自重の100倍のインパルスを入れる。↑より、粘性減衰を小さくした。†

• 変位図……単位:m,sec、スケール:10倍↑,粘性減衰:土0.00004,石0.00002
  • もうちょっと、先まで見たい。

モーダル解析†

• &link(sei.comm,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/salome-meca/sei/sei.comm)。上の動的解析と一緒にやった。
• 13.4Hz.1次モード的な
• 22.4Hz.2次モード的な
• 29.3Hz.ヤング率が異なる積石と裏込(本当は石だけど変形するように土にした)が食い込んでいる。
• 30.4Hz.ねじれ。しかし、実物の奥行き方向はもっと長いので、信用できないと思う。
• 34.6Hz.3次モード的な

参考†

• 積石のサイズ……1尺x1尺5寸x目測1尺5寸→0.30m x 0.45m x 0.45m　の4角錐
  • 大久保森造「石積みの秘法とその解説」理工図書,p.17
• 裏米の入れ方……「自然土質の場合における石積の裏込めの入れ方は積石の長さだけ(倍)」、図が斜め
  • 大久保森造「石積みの秘法とその解説」理工図書,p.102
• FEM解析での土の物性値……E=16.8MPa,\( \nu \)=0.33,\( \gamma \)=14.7kN/m\( ^3 \)とか?
  • ttp://www.forum8.co.jp/product/uc1/jiban/pdf/geo200707.pdf
• 土の応力ひずみ曲線……(\( \epsilon \),E):(0,0)→(0.01,150)and(0.01,150)→(0.15,150)とか?
  • ttp://shear.nagaokaut.ac.jp/Archive/PDF/Soilv2.pdf
• 男鹿石(安山岩)のヤング率……特に無いので石の平均ということで、36.75GPaぐらい?
  • http://www.minami-m.co.jp/reference/tokusei.html
  • http://www.green-site.com/manual/rock.htm
  • ttp://www.do-geo.com/cgi-bin/admin/survey/06rock_engineering.pdf
  • http://www.minami-m.co.jp/reference/tokusei.html
• 男鹿石(安山岩)の密度……2.5g/cm\( ^3 \)ぐらい?
  • http://www.monotsukuri.net/wbt/wbt_soseki/s0329/s0329.htm
  • ttp://www.do-geo.com/cgi-bin/admin/survey/06rock_engineering.pdf
  • http://www.minami-m.co.jp/reference/tokusei.html

blenderにモデルを作成†

• 1.石のstlファイルをインポート。スケールはsalomeに同じ。
• 2.edit mode→"p"→「構造的に分離したパーツで」にして分離する
• 3.運動の原点を物体の重心に持っていく(必須):"aa"で全選択→object→transform→origin to center of mass
• 4.床、壁をadd meshで手動で入れる。
• 5.以上は物理演算のrigidのpassiveにしておく。
• 6.適宜、反発係数・摩擦係数・質量等を変える。

pythonで栗石の自由落下†

• &link(ransuupack.py,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/blender/ransuupack.py)
  • 半径を0.075~0.030mの乱数にした。回転減衰:1.0,移動減衰:0.9
• 7.再生して、栗石を敷き詰める。

参考†

• 石の摩擦係数……blenderだと0.8ぐらいにすると、0.8×0.8＝0.64で丁度いい?
  • http://www.minami-m.co.jp/reference/tokusei.html
  • ttp://library.jsce.or.jp/jsce/open/00127/2001/47-A3-1605.pdf
• 石の反発係数……blenderだと0.2ぐらいにすると、0.2×0.2＝0.04で丁度いい?
  • ttp://www.skr.mlit.go.jp/yongi/duties/research/rakuseki/gaiyo2.pdf

blenderで計算†

• 8.最終コマで停止し、再生コマ(Dimensions/Frame Range)と計算コマ(rigid body cache)をそこからに変更する。
• 9.石を物理演算のactiveにする。
• 10.再生。
• 11.計算が終わったら、ベイクしてから保存する。
  • あとで開いてみてみると、なぜだか1個ぐらい石がずれている。

裏込にsoftbodyを使う場合†

• 7'.裏込のstlファイルをインポート。
• 8'.裏込についても、運動の原点を物体の重心に持っていく(必須):左クリックで選択→object→transform→origin to center of mass
• 9'.裏込をsoftbodyに、積石をactiveに設定。
• 10'.裏込および積石を全部collisionに設定。

blenderから動画保存†

• softbody:ほぼ初期状態：曲げだけ0.000→0.500に変えてみた
  • stlファイルの線の配置が荒いためか、すり抜ける。しかも、線が表面にしかないので左上の背面がペラペラしている。
• softbody:表面だけ細かくしてみた
  • 表面がただれるようだ。もうちょっと硬くてもいい。ちなみに、これだけでもすごく重い。
• softbody:表面だけ細かくして、曲げだけを0.500→1.000に変えてみた。
  • 徐々に曲げだけを大きくしていった。
  • どうやら曲げの強さパラメータを最大にしても、表面がだらーとするし、背面側に積石がすり抜ける
• softbody+rigid:背面にもrigidの属性を与えた
  • 積石の背面側へのすり抜けはなくなったが、飛び出してしまう。
• softbody+rigid:softbodyのgoalを1.000にした。
  • ほとんど、動かなくなってしまう。
• softbody+rigid:softbodyのgoalを0.900にした。
  • 今度は、動き過ぎる。
• softbody+rigid:softbodyのgoalを0.990にした。
  • goalを1.000に近づけても、まだまだ動き過ぎる。

• rigid:上のpackingのあと、積石をactive状態にした。
  • 石の回転減衰and移動減衰:1.000
  • softbodyよりいい。背面に、裏込なのに高く積みすぎ? もうちょっと先まで見たい。
    • 計算時間が約20時間ぐらい?

imagemagickでgif動画をつくる†

元画像の加工†

• jpg化
  • ( $ mogrify -format jpg moto*.png
  • )
  • jpgにしたほうが軽い。
• 画像の切り取り
  • ( $ mogrify -crop 400x300+500+400 moto*.jpg
  • )
  • 横幅x縦幅+横位置+縦位置 (ピクセル値、位置は左上から)
• 画像サイズを変える
  • ( $ mogrify -resize 400 moto*.jpg
  • )
  • 横が400ピクセルになる
    • このくらい小さくしないと、odpファイルの動作が遅い。

画像のgif化†

• ( $ convert -delay 10 moto*.jpg filename.gif
• )
• "-delay X"は再生コマ間隔がX/100秒となる。
  • 画像と同じで"$ eog filename.gif"で観られる。

officeで使う†

• gifはlibreofficeのimpressに画像として貼れる。
• impressから"microsoft powerpoint 97/2000/XP (.ppt)"で保存すると、そのままpowerpointでも読める。

LIGGGHTSでも動画を作ってみる†

　下記paraview+上記imagemagickの方法でgif化

• 例えば、こんなのとか?

2次元DEM(円筒要素)†

計算†

• &link(en.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/DEM/en.f90)
• &link(syoki.dat,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/DEM/syoki.dat)
  • おかしな動きをするので、そのうちに修正。

描画†

• paraview上で「ファイルを開く」:file→open→kekka..vtk→ok
• あるいは端末上で、
  • ( $ paraview kekka..vtk&
  • )
• 球体を表示させる:fileters→alphabetical→glyph→apply
  • 円筒もあるようだが、2Dなので真横から見れば一緒。
• 表示の設定:propertiesを下の表のように変更する。
  • paraviewのバージョン(新しいもの)によっては、Sphere-Radiusの項目が無いので、代わりにScale Modeを2にする。

    Glyph Type	Sphere-Radius	Scale Mode	Set Scale Factor
    Sphere	1	scalar	1(Editにチェックを入れて)

• アニメーション再生:play
  • modeをreal timeにすると動きが速くなる。

動画保存†

• file→save animation→animation durution(コマ数):100 (仮)→save animation
  • ↓お試し計算をimagemagickでgif化してみる。

卒論日誌†

------------------------------おまけ------------------------------†

Blenderの入手†

• https://www.blender.org/download/
• 解凍:tar xvfi blender-2.75a-linux-glibc211-x86_64.tar.bz2
  • ↑の「2.75a」はバージョン

日本語化†

• file(左上)
• preferences
• system(右上)
• international fonts(右下)にチェック
• languageを日本語に
• インターフェース・ツールチップ・新規データ(右下)にチェック
• ユーザー設定の保存(左下)

物理演算エンジン「bullet physics」の原理†

• http://www.slis.tsukuba.ac.jp/~fujisawa.makoto.fu/lecture/iml/index.html

Blenderの使い方†

参考URL†

• http://blender.jp/
• http://wiki.blender.org/index.php/Doc:JA/2.6/Manual
• http://mm-tech.co.jp/Blender.html

基本的な操作方法†

覚えないと使いこなせないキーボード操作がたくさんある。よく使うものを随時書き足すべし。

• 視点の移動

  スクロール	拡大縮小
  スクロールでドラッグ	見る角度を変える
  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	動画の終了

Object†

• Transform:物体の位置とか向きとか大きさを変更できる。

Physicsの設定†

Rigid Body(剛体)†

soft body(柔らかい物体)†

• colligion:soft bodyに衝突を入れたい物体には、collisionを設定。
  • それでもすり抜ける場合は、innerの値を大きくする。
• 曲げを見るには、edit mode で辺をたくさん生成しなければならない。

Materials†

• Diffuse:物体の色を変えられる。
• Transparency→alpha:値を小さくすると、透明になる。壁とかに使うといい。

Sceneの設定†

• Gravity:重力の大きさと方向を指定できる。
• Rigid Body World/Step Per Second:1秒あたりに計算する回数。精度に影響する。
• Rigid Body Cache:物理演算の開始(Start)と終了(End)コマを設定。

動画保存†

• レンダー/レンダー:カメラで撮る枠が見れる。
  • 表示:別ウィンドウ　がおすすめ
• 寸法
  • フレーム範囲:画像保存のステップと開始終了コマを設定。
• シェーディング/影:全ての物体の影を消す。
  • 影を消すには、他に、各物体ごとの設定(Materials)・ランプの設定(Lamp)等でも出来るらしい。
• 出力:画像を保存するディレクトリを指定
  • 形式、品質とうを指定。
• アニメーション:画像を保存

Renderの設定†

• Dimensions/Frame Range:動画の開始(Start Frame)と終了(End Frame)コマを設定。

ミスした時に戻る方法†

• 左側の「tool」→下の方の「history」→「undo」

キーフレーム†

　物体の位置・向き・スケールを2コマ以上設定しておくと、間の動作を補間してくれる便利機能。

• 1.物体を選択し、移動させたい位置や変化させたい向き、スケールを設定する。
• 2.物体を選択した状態で"i"を押す。
• 3.指定したいコマをクリック(下の帯)する。

キーフレーム＋物理エンジン†

　キーフレームを設定した物体と、何もしていない物体で、衝突とかの物理エンジンができる。

• キーフレームを設定した物体(例、剛体)
  • 剛体のタイプの「アニメ」にチェックを入れる。
  • あとは、通常通りに剛体の質量・摩擦・反発等を設定する。
• 何もしていない物体。
  • 普通に剛体の設定をする。

細かな設定†

• 1.左下の時計マークをクリック→グラフエディター
• 2.各動きの速度がグラフになっているので、手動で変更
  • 右クリックで選択、左クリックで決定
• 3.編集が終了次第、左下のグラフエディターをタイムラインに戻して、再生。

pythonでblenderを制御†

• 実行
  • ( $ blender --python filename.py&
  • )
• blendファイルを開いて、そのファイル上でpythonスクリプトを実行
  • ( $ blender filename1.blend --python filename2.py&
  • )

LIGGGHTS(ライツ?)等のインストール†

• 参考
  • http://www.geocities.co.jp/penguinitis2002/study/LIGGGHTS/LIGGGHTS/LIGGGHTS.html
  • http://blog.goo.ne.jp/dalab-gnct/e/e7883c4285c436a1b147688d4ce97b16
  • http://blog.goo.ne.jp/dalab-gnct/e/a4eabbb27eadb0af269a6eb9e9f1319e

準備†

• システム→システム管理→Synaptic
• FFTW:fftw-dev,fftw-docs,fftw2
• OpenMPI:libopenmpi-dev,libopenmpi1.3,openmpi-common,libopenmpi-dbg,openmpi-bin,openmpi-checkpoint,openmpi-doc 　をインストール指定し、適用。

LIGGGHTS-PUBLICの入手†

• ( $ cd /opt $ sudo mkdir LIGGGHTS $ cd LIGGGHTS $ git clone git://github.com/CFDEMproject/LIGGGHTS-PUBLIC.git
• )

Makefileの修正†

• ( $ 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_fpic
• ) 　liblmp_fedora_fpic.aという実行ファイルが出来ればOK。

リンク(?)をつくる†

• ( $ ln -s lmp_fedora_fpic liggghts $ ln -s lmp_fedora_fpic liblmp.a
• )

.bashrcに追加(下記LPPをインストールした後)†

• ( $ 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_DIR
• ) 　LPPを使わない人は、export PATH～LIGGGHTS-PUBLIC/srcの行だけでいいかも。それと、lpp.pyしか使っていないので、pizza.pyの行(alias pizza～ と ～CFDEM_PIZZA～)もいりませんでした。

後藤ちゃちゃ(15/4/30)†

• ( sudo apt-get install liggghts
• ) とかで入んないですかね。足りないものはあるかもしれませんが。

• ubuntu12.04のパッケージにあるliggghts1.5-1が手に入る模様。

chmod 777　?†

• ( $ sudo mkdir /opt/LIGGGHTS $ sudo chmod 777 /opt/LIGGGHTS
• ) 　すべてのユーザーが実行可能になるらしい。上記参考サイトを真似てやったけど、いらなかったかも。

LIGGGHTSの計算†

• 準備(一例)
  • ( $ mkdir /home/kouzou/2015/dem $ cp /opt/LIGGGHTS/LIGGGHTS-PUBLIC/examples/LIGGGHTS/Tutorials_public /home/kouzou/2015/dem
  • ) 　kouzou/のどこかにDEMの作業用のディレクトリを作って、そこにTutorials_public(サンプルが色々入っている)をコピーしてから使うことにした。
• 計算方法(in.filenameをvi等で編集した後、同ディレクトリにて)
  • ( $ liggghts < in.filename
  • )
• 計算過程は、log.liggghtsに書かれ、結果はpost/dump*.filename(変更可)に出力される。

関連ソフト(可視化)†

lpp(vtkファイルに変換)†

• インストール
  • ( $ cd /opt/LIGGGHTS $ git clone git://github.com/CFDEMproject/LPP.git
  • )
  • LPP/srcにあるlpp.pyが実行ファイル。LPPごと/opt/LIGGGHTSに移動させた。
• 計算方法(dump*.filenameと同ディレクトリにて)
  • ( $ lpp.py dump*.filename
  • )

paraview(vtkファイルを描出)†

• 準備
  • システム→システム管理→Synaptic
  • python-pmw,python-pmw-doc
  • python-imaging-tk,python-imaging-tk-dbg
  • python-opengl 　をインストール指定し、適用。
• インストール
  • http://www.paraview.org/download/
• 解凍
  • tar xvzf ParaView-4.3.1-Linux-64bit-glibc-2.3.6.tar.gz
• 使い方 　以下、おすすめの設定。詳しくは林真『はじめてのParaView』@図書館を参照。
  • paraviewを開く(アプリケーションのグラフィックス にある）
  • File → Open → dump..vtk → OK
  • Apply ……物体が見えるようになる
  • Filters → Common → Glyph
    • Glyph type : sphere ……点が球体になる
    • Scalars : radius
    • Scale Mode : scalar
    • Scale Factor : 1 ……球体の半径が計算結果と同じ大きさになる
    • Apply ……変更が反映される
  • View → Animation View → Mode:real time ……早く動くので見やすくなる
  • Play：再生

LIGGGHTS-PUBLIC 3.2.0　のコマンドのこと†

参考†

• Manual:http://www.cfdem.com/media/DEM/docu/Manual.html
• GitHub(benchmarkとかexamplesとかに計算例がある):
  • https://github.com/CFDEMproject/LIGGGHTS-3-beta-PUBLIC
  • https://github.com/metrdia/LIGGGHTS_B
  • https://github.com/ligg-tst/liggghts-examples/tree/master/USER
• 粒子DEMの解説：Catherine O'Sullivan『粒子個別要素法』@図書館

対応表ver3.2.0†

　ネット上にあるのを参考にすると、バージョンの問題で「そのコマンドは定義されていません」というエラーが出ることが結構ある。なので、下の表ように書きなおさなければならない。

LIGGGHTSと手計算の比較†

自由落下†

• 入力ファイル:&link(in.rakka,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/LIGGGHTS/nyuuryoku/in.rakka)

　上の図のように、球体を自由落下させた。球に働いている力は重力のみである。反発係数だけを変えたのが下の図である。

　横軸が経過時間、縦軸が鉛直座標である。反発係数\( 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 \)に変化させて計算した。

• 自由落下させたとき、跳ね返る高さは元の高さの\( e^2 \)倍となる。図より、LIGGGHTSの計算でもそうなっている(\( \frac{1}{1},\frac{3}{4},\frac{1}{2},\frac{1}{4} \))ことが分かる。
• \( t=\sqrt{\frac{2z}{g}} \)　に\( z=0.14m,g=9.81m/s^2 \)を代入すると、落下までの時間は\( t=0.17sec \)……LIGGGHTSの計算と一致

坂と球†

• 入力ファイル：&link(in.saka,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/LIGGGHTS/nyuuryoku/in.saka)
• 壁：&link(saka.stl,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/LIGGGHTS/kabe/saka.stl)

　上の図のように、球体を坂に置いた。球に働いている力は重力のみである。摩擦係数だけを変えたのが下の図である。

　横軸が経過時間、縦軸が鉛直座標である。摩擦係数\( \mu \)のみを\( 0.00,0.05,0.10,0.30,0.50,1.00 \)に変化させて計算した。

• 形状を考慮しない質点の力学では、時間tと鉛直変位の関係は
• \( t=\sqrt{\frac{2z}{g(sin\theta-\mu cos\theta)sin\theta}} \)
  • 上の式に\( z=0.28m,\theta=25.8^\circ,\mu=0.00 \)を代入すると、落下までの時間は\( t=0.55sec \)……\( \mu=0.00 \)の結果に一致
  • 上の式に\( z=0.28m,\theta=25.8^\circ,\mu=0.05 \)を代入すると、落下までの時間は\( t=0.58sec \)……\( \mu=0.05 \)の結果に一致
• 形状による回転を考慮した剛体の力学では、時間tと鉛直変位の関係は
• \( t=\sqrt{\frac{14z}{5g sin^2\theta}} \)
  • 上の式に\( z=0.28m,\theta=25.8^\circ \)を代入すると、落下までの時間は\( t=0.65sec \)……\( \mu=0.30,0.50,1.00 \)の結果に一致 　傾きがより急な場合はどうなるか分からないが、傾斜角25.8度では摩擦係数が0.3以上では、ほぼ剛体の運動に収束していることが分かった。また、摩擦係数にも色々あるような気がするが、下の表とLIGGGHTSの計算結果と比較すると石積みの場合は転がりのみで、ほとんど滑りが発生しないのではないかと思う。

    摩擦片	摩擦面	摩擦係数
    石	石	0.6~0.7
    石	土	0.5(乾)~0.3(湿)
    木	石	0.4

• ttp://www.labnotes.jp/pdf2/friction.pdf

LIGGGHTS関連に使用したFORTRAN90のプログラム†

項目別経時変化をgnuplot用にまとめる†

• &link(t.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/LIGGGHTS/fortran/t.f90)
  • 上手いことやればgrepでも出来るかも。

grep†

• 出力ファイルに書かれたデータをまとめる
• 例えば、
  • ( $ cat dump* |grep 1\ 1> tatoeba
  • )
• で、「dump〇〇」というファイルで、「１　１」つづきになる行を「tatoeba」という名前のファイルにまとめてくれる。

imagemagickのコマンドで、動画を作ってみる†

• スクリーンをjpgに連写
  • ( \( xwininfo |grep '^xwininfo: Window id:' | awk '{print \)4}'
  • ) 　で、ウィンドウのidを調べて、実行。
  • &link(rensya.f90,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/rensya.f90)
    • imagemagickのスクリーンショットに元々タイムラグがあるため、連写間隔が広くてよくない方法!
• 画像の切り取り
  • ( $ convert moto*.jpg -crop 400x300+500+400 kiritori.jpg
  • )
  • 横幅x縦幅+横位置+縦位置 (ピクセル値、位置は左上から)
• 画像のgif化
  • ( $ convert -delay 10 kiritori*.jpg filename.gif
  • )
  • "-delay X"は再生コマ間隔がX/100秒となる。
    • 画像と同じで"$ eog filename.gif"で観られる。
    • gifはlibreofficeのwriterに画像として貼れる。

libreofficeのimpressに画像を貼り付ける†

• サイズを小さくしないとodpファイルが動かなくて、発表がスムーズに出来ない。なので、
  • 400ピクセルぐらいにしたほうがいい。
  • pngよりjpgの方が軽い。
  • ( $ mogrify -resize 400 *jpg *gif
  • )

salome-meca課題メモ†

code-asterを使ってみる†

• d=0.02m,l=1.0mの円柱(鋼材：E=206GPa)の片持ち梁先端に1000N(≒100kg)の集中載荷をした。

  条件	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

• まとめ
  • 面載荷よりも点載荷の方が少しだけ理論値に近い
  • 要素数が多いほうが理論値に近い
  • メッシュが粗くても、4面体2次要素は理論値に極めて近い

6/22までの課題†

　昨年度の創造工房でFEM解析(鈴木さんの)に使用した値

• ヤング率：2840MPa
• ポアソン比：0.313
• 荷重：1000N を用いた。
  • モデルの結果

    拘束条件	載荷条件	要素の種類	要素数	節点数	たわみ(mm)	体積(mm\( ^3 \))
    線	面(1cmのひもを想定)	4面体2次	32705	55184	1.16	11276

振り子†

　&link(Chipmunk BASICのプログラム,http://www.str.ce.akita-u.ac.jp/~gotou/programoj/bas/hurikoc.bas)を書き換えて、作れそうな長さの振り子でもっとも最短で振り子が揃うものを探してみた。

• 条件
  • 長い方：45～70cm
  • 短い方：20～35cm
  • 変化させた長さ0.1cmずつ
• &link(印刷用の振り子の長さ,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/opencampus/huriko)
  • 長さの幅を広げれば、もっと短い時間で済むのがあるかも？
• &link(12個の振り子プログラム,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/kondo/opencampus/huriko12c.bas)

夏休みの課題†

• 1.blenderで石積み(裏込めあり、杭なし、背面土なし、水なし)をつくる
  • 裏込めをヤング率の低いsoftbodyにするか、粘性係数の高いfluidにするか。
    • softbodyでやったが、動きが変だし、計算が終わらなかったりする。
  • まずは自重のみ。出来ればインパルスを入れたい。
  • 変位や応力とかの確認が出来るのかを調べる。
    • ポストプロセッサーは内蔵していないと思う。
• 2.salome-mecaで自重で変形する石積み(同条件)をつくり静的な解析をする……OK
  • まずは、密度を入れて単純梁で手計算と比較する。
  • 裏込めはヤング率の低くする。
• 3.salome-mecaで動的な解析が出来るかの確認……OK
  • まずは自重のみ。出来ればインパルスを入れたい。
    • 過渡応答解析は出来たが、弾性解析だと結局元に戻ってしまうので、塑性解析が必要？
  • 単純梁でこれができたら、いずれは石積みを解きたい。

code-asterで動的な解析†

自重でたわむ単純梁†

• 10mm角のスパン100mm(創造工房に同じ)
• ヤング率:1MPa(山内さんの折り目を参考に適当にやわらかくした)
• ポアソン比:0.3(一般的な土)
• 参考:ttp://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin=attach&refer=SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2&openfile=12-00.pdf
• 密度:2.65g/cm\( ^3 \)(一般的な土)
• 粘性減衰:0.01(適当)
• タイムステップ0.001秒(適当)で0.1秒まで(適当)計算

  	初等梁理論	静的解析	過渡応答解析
  中央のたわみ(mm)	4.0578	4.1960	4.1958

• 一瞬だけ荷重をかけたり(自重をかける時間を0.001とかにする)、一瞬だけ荷重のスケールを大きく(or小さく)すると、ちゃんと梁が上下に揺れるので、垂直方向のインパルスはいれられる。
  • 自重をかけた状態で横揺れを入れるには、もう少しasterの勉強が必要。

最新の100件