関原の卒論日誌

再発表†

結果はまだ出せていませんが、大まかなスライドの流れと疑問?を作りました。もしこれでだいたいの方向性がいいならテストが終わってからすぐにでも結果もスライドも作りたいです。

ダイヤカット円筒†

シェル要素†

• http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/sekihara/daiyat6c.f90
  • このファイルで自分の作りたいダイヤカット円筒のサイズやパターン数を決めて実行ファイルにして./daiyastl>???.inpと書きだすとシェル要素でメッシュが切られてるダイヤカット円筒のinpファイルが出来る。

ｓａｒｏｍｅで要素を決める†

• http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/sekihara/daiyastl.f90
  • このファイルで自分の作りたいダイヤカット円筒のサイズやパターン数を決めて実行ファイルにして./daiyastl>???.stlのように書き出すとダイヤカット円筒のstlファイルが出来る。
• http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/sekihara/c3d4unv19a3.f90
  • その後stlファイルをfreecadでbrepファイルにエクスポートしsaromeでメッシュ分割し、unvファイルにエクスポートし、メッシュの要素数をこのファイルに書き込み実行するとダイヤカット円筒のinpファイルが出来る。詳しい方法は、下に書いてある。

頂点を吐き出すプログラム†

• http://www.str.ce.akita-u.ac.jp/~gotouhan/j2015/sekihara/tyouten.f90
  • このプログラムである程度の頂点は吐き出すことができたが、座標の値が０や０に限りなく近い数値だけは吐き出すことが出来なかった。ある程度は吐き出せているから残りは目視でやっても大幅に作業は短縮できた。
• 
  • この画像はダイヤカット円筒を真上からみた画像で、頂点を吐き出して目視で頂点を全部書き込み載荷したが変な挙動になった。色々手をつけてみてがわからない。
• 
  • 目視での頂点が間違ってたみたいで修正するときれいに載荷できた形になった。

ccx,cgx†

• 
• 
  • 上が坪井さんのダイヤカット円筒で、それと比べると均一に載荷できてないように見えるから何が原因なのか調べたい。
  • 今回はcgxで節点セットを作成してその節点セットを解析ファイルのcload入力に書いていたがこのやり方自体厳しい気がする。c3d4unv.f90コンパイルして出てくる実行ファイルを実行してできるinpファイルに最初から拘束や荷重が書きこまれているようにしたいができない。c3d4unv.f90から出来るinpファイルをccxで計算すると、エラーが出て多分そのせいで拘束や荷重が書き込まれずにできているんだと思う。でも一応cgxで図は見れる。
  • 
  • 最初からinpファイルに書きこまれているようにできたが、画像のようになってしまった。何が原因なのかわからない。

卒論テーマ†

• ダイヤカットや蛇腹折り円筒を今までより厳密なFEMで解析し、従来法と比較
  • 今まではシェル要素で解いているが、四面体要素で厚みや折れ目を考慮して解析してみたい。
  • 折りたたみ構造の折り目を弾塑性要素でFEM解析
  • これまで、ダイヤカットや蛇腹折り円筒は、シェル要素で解析しているが、 実際には、金属類をプレスで折り曲げるとすれば、折り目は塑性変形しており、 折り目部分を弾塑性要素で解析すべきではないか（初期状態を塑性硬化の状態にする？）

卒論日誌†

バネ定数†

• 坪井さんの卒論で使われていたdaiyat6a.f90を少し変えてコンパイルし実行すると、大量のinpファイルとjobdaiya1が出てきた。
• それらをgFTPからk2に移動してjobdaiya1を実行すれば、すべてのパターンのdatファイルができた。
• datファイルの中身は様々な点での変位が計算されていた。
• 荷重（２００Ｎ）÷鉛直変位＝バネ定数より下の表にまとめた。
• 53π×10^-3を25π×10^-3にして高さを104,5×10^-3を50×10^-3にした。

バネ定数 表 グラフ†

３次元プロット グラフの画像†

• ( set pm3d splot "hoge" w l
• ) で３次元の色ありのグラフが出来る。

その後に

• ( set terminal png
• ) と打ち
• ( set out "image1.png" replot exit
• ) でpngファイルとして出力される。

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で解ける。

これまでの目標†

後藤ちゃちゃ(16/2/1)†

作業記録なので、これまでの経緯も削除しないで記録しておいていただけますか。

12/18までの目標†

• k2でcalculixで計算して結果を持ってきて見てみる。

12/4までの目標†

• calculixを使えたので座屈モードなど計算する。

11/6までの目標†

• 後藤資料のCAELinuxのFEM解法のページにある「Salomeでメッシュを切ったunvファイルからCalculiXのinpに変換」の方法について試してみる。

10/23までの目標†

• まずcalculixで計算できるとこまでやる。

10/16までの目標†

• ダイヤカット円筒のstlファイルを作り、それをsalomeに持ってきて四面体要素で解いてみる。
• 弾塑性要素で解析する方法についても探ってみる。

10/9までの目標†

• 坪井の卒論日誌 坪井さんのdaiyat6a.f90を去年の坪井さんの卒論と同じ 円筒（同じ材料定数、違うサイズ）にして解いてみる。
  • 大量のinpファイルが生成されるので、daiyat6aスクリプトを使ってそれをCalculiXに自動入力し計算
  • 計算されたdatファイルから、上端部の鉛直変位を読み取る。
  • 荷重/鉛直変位でばね定数を求める
  • それらをデータにして、gnuplotで3次元プロットしてみる（坪井さんと同じグラフになるか）
  • ここまでの結果を東北支部概要スタイルで\( \TeX \)でまとめる。

最新の100件