石坂の卒論日記

ねじれ座屈荷重の比較†

• (

グリーンヒル問題の座屈モーメント理論値＝865.281(Nm)

                               解析値＝868.232(Nm)   近い値が解析できた。

• )

グリーンヒル問題†

後藤ちゃちゃ(15/1/9)

• http://www.google.co.jp/url?sa=t&rct=j&q=&esrc=s&source=web&cd=2&ved=0CC4QFjAB&url=http%3A%2F%2Flink.springer.com%2Fcontent%2Fpdf%2F10.1007%252FBF02426915.pdf&ei=soGvVM7XE8LMmAWD2YHACQ&usg=AFQjCNEFGoPXYmIJadw7-R8uSNiA7DXOPA&bvm=bv.83339334,d.dGY

• ( Yoshiaki GOTO, Xiao-Song LI, Toshihiro KASUGAI and Makoto OBATA: Analysis of Greenhill Problem by a Co-Rotational Method, 構造工学論文集, Vol. 41A, 1995.
• ) に、ねじりを受ける円断面のグリーンヒルによる座屈モーメントが書かれています。 \( M_{cr}=2.861\frac{\pi}{\ell}EI \) ただし、\( I=\frac{\pi r^{4}}{4} \)

解析プログラムを調査†

• (

  x軸回りの曲げ試験で測定されるヤング率 ex=206.d9
  y軸回りの曲げ試験で測定されるヤング率：ey=ex
  ねじり試験で測定されるせん断弾性係数：g=1.d0
  ell=1.d0

  c 断面積：a

         a=1.5d-4

  c x軸（強軸）回りの断面二次モーメント：xi

         yi=1.25d-9

  c y軸（弱軸）回りの断面二次モーメント：yi

         xi=2.8125d-9

  c ねじり定数：ej

         ej=1.d0

  c そりねじり定数：wi

         wi=1.d0

  c write(*,*)' 一発目の弧長を入れて下さい' c read(*,*) r

          r=eul*0.999d0

  c write(*,*)' 増分弧長を入れて下さい' c read(*,*) rn

         rn=eul*0.1d-5

• )

片持ち梁鋼材†

• 分割数64
  • (

    ----------------------------------------------------
    N-R反復数= 2
    det=      2.865d1940  ,  hu koyuuti:   0
     p=  6.35497E+02   PL/EA=  2.05662E-05    Euler=  6.35356E+02
    64右  0.000E+00  0.000E+00 -2.057E-05  0.000E+00  0.000E+00  0.000E+00  0.00E+00
    ----------------------------------------------------

  • )

• 分割数256
  • (

    ----------------------------------------------------
    N-R反復数= 2
    det=     -1.886d8092  ,  hu koyuuti:   1
     p=  6.35377E+02   PL/EA=  2.05624E-05    Euler=  6.35356E+02

    256右 0.000E+00 0.000E+00 -2.056E-05 0.000E+00 0.000E+00 0.000E+00 0.00E+00

    座屈したみでだよ................
    ----------------------------------------------------

• )

後藤メモ†

k2のgakusei/bunken/g02j97oiraa1.pdf がグリーンヒルの問題を梁要素で解いている例題ですが、つりあい状態がどんどん変化する座屈後変形の問題なので、CalculiXでは、解けるとしても簡単には解けないでしょう。（座屈したら、座屈モードを出力して、それを\( 10^{-3} \)倍とか小さくした値を初期不正としてinpの初期データに加えて、静的問題として解いて座屈後変形を起こして....みたいなやり方）

• http://www.str.ce.akita-u.ac.jp/~gotou/programoj/#hari ここのhari.fのプログラム（上記の論文で用いられているプログラム）で、上記の例題を同じように解けるかどうかをまず確認して、引張を加えたらどうなるかというアプローチでもいいかもしれません。 まあ、一緒に考えますが、まずは、文献を読んで、どんな例題なのかイメージできるようにしといてください。

卒論日誌†

四面体メッシュ†

• グラフ

直方体メッシュ†

• 初期値15以下でセグメントを設定すると誤差が大きくなることがわかった。
• 直方体分割の場合セグメント数を大きくする事によって理論値の値に近づいていくという事がわかった。

直方異性体†

• グラフ

2つの材料でできた部材†

• 試験結果 x.12mm y.50mm z.50mm たわみはy軸方向とする

  length

  節点数

  要素数

  たわみ

  18

  464

  1728

  1.247

• hari.fの解析
• modo.dの値が反転して出ている可能性があるので、今までうまく行かなかったオイラー座屈の座屈後の挙動にいけるようなmodoの組み合わせを見つける。

最新の100件