稲荷の卒論日誌 - *座屈荷重-折れ角

2枚の板を重ねたやつを上に引っ張る(高さ0の状態から)
2枚の板を重ねたやつを上に引っ張る(初期不整を与えた状態から)
10/3 15時
ccxの増分繰り返しの制御
10/3
10/1
9/30 13時
9/26 13時
円筒を引っ張った
8段のやつを引っ張った
厚さ0にしたらこうなった
gnuplot 要素数と変位
拘束条件、両端固定梁
拘束条件、単純梁
08/6/27
リンク
LaTeX
fortranメモ
稲荷の日誌

座屈荷重-折れ角

折れ角

ψ

は下の図の

ψ

。

初期高さを $a$ とすると、この $ψ$ は、

$ψ = a r c sin (a)$

になる。ちなみに $ψ = a r c sin (a)$ の時は、 $a = sin (ψ)$ 。

$ψ = a r c sin (a)$ から、初期高さを折れ角になおした↓

初期高さ	折れ角 $ψ$
0	0
0.01	0.57297
0.1	5.73917
0.2	11.5370
0.3	17.4576
0.4	23.5782
0.5	30
0.6	36.8699
0.7	44.4270
0.8	53.1301
0.9	64.1581
0.95	71.8051
0.98	78.5217
0.99	81.8904
1	90

ばね定数-折れ角

座屈荷重-折れ角

tubu4sub_2.f。初期高さを倍率で指定できるようにした。実行すると、最小の高さ方向分割数も求めることができる。

文字コード

支承について

Ctrl+f→耐用年数 http://www.j-rba.com/php/faq/faq.php

角柱の角数-圧縮による座屈荷重

28角柱の圧縮による座屈荷重 $9.66977 \times 10^{- 3}$ MN 30角柱の圧縮による座屈荷重 $9.98460 \times 10^{- 3}$ MN

最小の高さ方向分割数を求める計算式

×高さ方向分割数nを求める式

$n > \frac{h m}{s} (\frac{1}{tan \frac{2 π}{m}} + \frac{1}{sin \frac{2 π}{m}}) \times \sqrt{1 - a^{2}}$

a:初期高さ
s:円周の長さ
h:角柱の高さ
m:周方向分割数

右辺に出た値よりnを大きくしないと、折り畳み円筒が作れない。初期高さaには、0を代入して計算したほうがいい。h～0.0hまで解析するなら

ble6ccxfについて

inpファイルの座標値について

ble6ccxf.fを使うと有効数字９ケタまで出てくる。今まで使っていたble6ori3.fを使うと有効数字５ケタまで出てくる。

ble6ccxfとble6ori3を、引張の座屈荷重で比較

4角形20段

初期高さ	座屈荷重MN(ble6ccxf)	座屈荷重MN(ble6ori3)
0.9h	$9.051698 \times 10^{- 3}$	$9.050731 \times 10^{- 3}$
0.8h	$1.018282 \times 10^{- 2}$	$1.018381 \times 10^{- 2}$
0.7h	$1.105013 \times 10^{- 2}$	$1.105171 \times 10^{- 2}$
0.0h	$2.943180 \times 10^{- 3}$	$2.943783 \times 10^{- 3}$

6角形20段

初期高さ	座屈荷重MN(ble6ccxf)	座屈荷重MN(ble6ori3)
0.9h	$1.054768 \times 10^{- 2}$	$1.054118 \times 10^{- 2}$
0.8h	$1.020788 \times 10^{- 2}$	$1.020751 \times 10^{- 2}$
0.0h	$3.385171 \times 10^{- 3}$	$3.386029 \times 10^{- 3}$

8角形20段

初期高さ	座屈荷重MN(ble6ccxf)	座屈荷重MN(ble6ori3)
0.9h	$1.238677 \times 10^{- 2}$	$1.239338 \times 10^{- 2}$
0.8h	$7.300154 \times 10^{- 3}$	$7.304122 \times 10^{- 3}$
0.0h	$4.406272 \times 10^{- 3}$	$4.406603 \times 10^{- 3}$

円管を座屈

円管(周方向分割数26、高さ方向分割数36)を圧縮で座屈させた。

角数とバネ定数の関係

4角柱から、(30角柱)に増やしていったときのバネ定数の変化。

28角柱のばね定数= $2.13711 \times 10^{1}$ MN/m

30角柱のばね定数= $2.13940 \times 10^{1}$ MN/m

初期高さと座屈荷重の関係(20段)

座屈荷重ー初期高さグラフ

引張による座屈荷重が正で、圧縮による座屈荷重が負

0.4h～0.8hの時、引張に対して座屈しにくいのは4角形。圧縮に対しては、8角形>6>4の順で座屈しにくい。

圧縮の座屈で、hの時に座屈荷重が落ちるのは、座屈のモードが違うから。0.1～0.9のときは局部座屈で、hは柱の座屈

引張による座屈荷重

初期高さ	座屈荷重MN(10角形)	座屈荷重MN(八角形)	座屈荷重MN(六角形)	座屈荷重MN(四角形)
h	座屈しない	座屈しない	座屈しない	座屈しない
0.99h		$1.569016 \times 10^{- 2}$	$1.158171 \times 10^{- 2}$	$1.006837 \times 10^{- 2}$
0.98h	$2.062936 \times 10^{- 2}$	$1.280695 \times 10^{- 2}$	$9.755915 \times 10^{- 3}$	$8.727733 \times 10^{- 3}$
0.95h		$1.297384 \times 10^{- 2}$	$9.618086 \times 10^{- 3}$	$8.437849 \times 10^{- 3}$
0.9h		$1.239338 \times 10^{- 2}$	$1.054117 \times 10^{- 2}$	$9.050730 \times 10^{- 3}$
0.8h		$7.304123 \times 10^{- 3}$	$1.020751 \times 10^{- 2}$	$1.018368 \times 10^{- 2}$
0.7h		$5.843784 \times 10^{- 3}$	$6.865188 \times 10^{- 3}$	$1.105169 \times 10^{- 2}$
0.6h		$5.227489 \times 10^{- 3}$	$5.177106 \times 10^{- 3}$	$9.394533 \times 10^{- 3}$
0.5h		$4.886570 \times 10^{- 3}$	$4.145332 \times 10^{- 3}$	$7.280426 \times 10^{- 3}$
0.4h		$4.738051 \times 10^{- 3}$	$3.439166 \times 10^{- 3}$	$5.270162 \times 10^{- 3}$
0.3h		$4.473502 \times 10^{- 3}$	$2.964714 \times 10^{- 3}$	$3.514957 \times 10^{- 3}$
0.2h		$3.247805 \times 10^{- 3}$	$2.581031 \times 10^{- 3}$	$1.945069 \times 10^{- 3}$
0.1h		$9.910846 \times 10^{- 4}$	$1.076251 \times 10^{- 3}$	$9.301209 \times 10^{- 4}$
0.01h		$2.220658 \times 10^{- 4}$	$1.876647 \times 10^{- 4}$	$1.729322 \times 10^{- 4}$
0.0h		$4.406609 \times 10^{- 3}$	$3.386023 \times 10^{- 3}$	$2.943783 \times 10^{- 3}$

圧縮による座屈荷重

初期高さ	座屈荷重MN(10角形)	座屈荷重MN(八角形)	座屈荷重MN(六角形)	座屈荷重MN(四角形)
h		$9.227003 \times 10^{- 4}$	$5.492898 \times 10^{- 4}$	$2.7766868 \times 10^{- 4}$
0.99h		$6.546836 \times 10^{- 3}$	$4.717111 \times 10^{- 3}$	$1.933918 \times 10^{- 3}$
0.98h	$8.419665 \times 10^{- 3}$	$7.403405 \times 10^{- 3}$	$5.826513 \times 10^{- 3}$	$3.204657 \times 10^{- 3}$
0.95h		$8.048174 \times 10^{- 3}$	$6.838534 \times 10^{- 3}$	$4.604031 \times 10^{- 3}$
0.9h		$7.971688 \times 10^{- 3}$	$6.896642 \times 10^{- 3}$	$4.982872 \times 10^{- 3}$
0.8h		$7.263279 \times 10^{- 3}$	$6.337780 \times 10^{- 3}$	$4.710710 \times 10^{- 3}$
0.7h		$6.640925 \times 10^{- 3}$	$5.559279 \times 10^{- 3}$	$4.187311 \times 10^{- 3}$
0.6h		$6.113434 \times 10^{- 3}$	$4.903180 \times 10^{- 3}$	$3.556752 \times 10^{- 3}$
0.5h		$5.597126 \times 10^{- 3}$	$4.205012 \times 10^{- 3}$	$2.809136 \times 10^{- 3}$
0.4h		$3.850159 \times 10^{- 3}$	$3.266090 \times 10^{- 3}$	$1.960684 \times 10^{- 3}$
0.3h		$2.136936 \times 10^{- 3}$	$2.046746 \times 10^{- 3}$	$1.160451 \times 10^{- 3}$
0.2h		$8.826734 \times 10^{- 4}$	$9.318319 \times 10^{- 4}$	$5.764478 \times 10^{- 4}$
0.1h		$2.166468 \times 10^{- 4}$	$2.907486 \times 10^{- 4}$	$2.288958 \times 10^{- 4}$
0.01h			$5.742953 \times 10^{- 5}$	$5.518308 \times 10^{- 5}$

cgxで座屈モードを見た(八角形)

h(八角柱)

引張では座屈しなかった
圧縮

0.9h

0.8h

0.7h

0.6h

0.5h

0.4h

0.3h

0.2h

0.1h

0.0h

引張

引張の時は谷折り部が変形、する傾向がある
圧縮の時は山折り部が変形(同時に谷折り部も少し変形)、する傾向がある
初期高さが大きい場合は固定端側、小さい場合は載荷側に変形、する傾向がある

cgxで座屈モードを見た(六角形)

h(六角柱)

引張では座屈しなかった
圧縮

0.9h

0.8h

0.7h

0.6h

0.5h

0.4h

0.3h

0.2h

0.1h

0.0h

引張

引張の時は谷折り部が変形、する傾向がある
圧縮の時は山折り部が変形(同時に谷折り部も少し変形)、する傾向がある
初期高さが大きい場合は固定端側、小さい場合は載荷側に変形、する傾向がある

cgxで座屈モードを見た(四角形)

h(四角柱)

引張では座屈しなかった
圧縮

0.9h

0.8h

0.7h

0.6h

0.5h

0.4h

0.3h

0.2h

0.1h

0.0h

引張

初期高さ・圧縮・引張に関係なく、すべての座屈モードが、山折り部が変形
初期高さが低くても高くても載荷側に変形、する傾向がある

初期高さ-ばね定数関係

縦軸がばね定数(MN/m)、横軸が初期高さ。
3つとも段数は20段。
圧縮も引張もバネ定数は一緒

初期高さ	ばね定数MN/m(4角形)	ばね定数MN/m(6角形)	ばね定数MN/m(8角形)
0.0h	$2.850 \times 10^{- 2}$	$5.162 \times 10^{- 2}$	$9.037 \times 10^{- 2}$
0.1h	$8.414 \times 10^{- 2}$	$1.009 \times 10^{- 1}$	$1.050 \times 10^{- 1}$
0.2h	$1.175 \times 10^{- 1}$	$1.322 \times 10^{- 1}$	$1.243 \times 10^{- 1}$
0.3h	$1.514 \times 10^{- 1}$	$1.647 \times 10^{- 1}$	$1.471 \times 10^{- 1}$
0.4h	$1.876 \times 10^{- 1}$	$2.067 \times 10^{- 1}$	$1.786 \times 10^{- 1}$
0.5h	$2.317 \times 10^{- 1}$	$2.632 \times 10^{- 1}$	$2.238 \times 10^{- 1}$
0.6h	$2.914 \times 10^{- 1}$	$3.411 \times 10^{- 1}$	$2.950 \times 10^{- 1}$
0.7h	$3.837 \times 10^{- 1}$	$4.667 \times 10^{- 1}$	$4.140 \times 10^{- 1}$
0.8h	$5.521 \times 10^{- 1}$	$6.879 \times 10^{- 1}$	$6.427 \times 10^{- 1}$
0.9h	$9.527 \times 10^{- 1}$	$1.228$	$1.230$
0.95h	$1.449$	$1.926$	$2.076$
0.98h	$2.234$	$2.956$	$3.416$
0.99h	$2.992$	$3.813$	$4.445$
1.0h	$1.989 \times 10^{1}$	$2.044 \times 10^{1}$	$2.072 \times 10^{1}$

四段四角形の、初期高さと座屈荷重

座屈荷重

初期高さ	座屈荷重MN(引張)	座屈荷重MN(圧縮)
h(四角柱)	座屈しない	$3.159421 \times 10^{- 4}$
0.9h	$1.836944 \times 10^{- 3}$	$6.389281 \times 10^{- 4}$
0.8h	$1.123588 \times 10^{- 3}$	$4.688874 \times 10^{- 4}$
0.6h	$7.472452 \times 10^{- 4}$	$3.162731 \times 10^{- 4}$
0.4h	$5.415855 \times 10^{- 4}$	$2.106308 \times 10^{- 4}$
0.2h	$3.424033 \times 10^{- 4}$	$1.116131 \times 10^{- 4}$
0.0h	$3.780499 \times 10^{- 4}$	$3.780503 \times 10^{- 4}$

座屈モードをcgxで見る

h(四角柱)

引張では座屈しなかった
圧縮

0.9h

0.8h

0.6h

0.4h

0.2h

0.0h

どの初期高さでも、引張の時は谷折り部が変形

どの初期高さでも、圧縮の時は山折り部が変形(同時に谷折り部も少し変形)

0.9hの引張座屈を、横から見た画像。

0.8h、0.6h、0.4h、0.2hの引張座屈はこういうふうに左右対象になる

ハードディスクの容量をチェックしたいときは

df -m

ディレクトリごと圧縮して.tar.gzアーカイブファイルを作るには、

初期高さとばね定数の関係(四段の四角形)その3

hが四角柱の場合

1個の台形につき3要素と、12要素を比較

hが円管の場合

縦軸がバネ定数(MN/m)、横軸が初期高さ。

赤線が今までの要素数で、緑線が要素数を増やしたやつ。

グラフ

初期高さとばね定数の関係(四段の四角形)その2

縦軸がバネ定数(MN/m)、横軸が初期高さ。四角柱がすごく硬い。

グラフ

4角柱の画像

初期高さとばね定数の関係(四段の四角形)

縦軸がバネ定数(MN/m)、横軸が初期高さ。

グラフ

0.9hをcgxで見ると、円管というより四角柱に近づいてるように見える。だから0.9hのばね定数が大きかったのかな。

初期不整による、荷重-変位の関係の変化(四段の四角形)

円管の高さをhとして、0.8h、0.6h…というふうに折り畳んでいく。そして、それらに圧縮力を加えた。円管の材料諸元↓

高さ	12cm
半径	3cm
厚さ	0.2mm
ヤング率	69GPa

拘束条件は、固定側x,y,z、引張側x,yの並進変位の固定。

赤線がh(円管)、緑線が0.8h、濃い青線が0.6h、紫線が0.4h、薄い青線が0.2h、茶線が0h。

縦軸が荷重(N)で、横軸が変位(mm)。

赤線が緑線と重なっていて見にくいが、数値でみると、0.8hよりh(円管)の方が変位が大きい。
円管の載荷節点数と、折り畳み円筒の載荷節点数が違うので、CLOADは(1/載荷節点数)で指定している。
この図は圧縮による変位のグラフだけど、引張では数値の＋－が変わるだけだった。

後藤ちゃちゃ08/11/13

グラフは、だいたい予想通りだけど、圧縮も引張もバネ定数は同じなんだ(まるでバネみたい)。縦軸をバネ定数(つまり荷重/変位)、横軸を初期高さにしてみたら、予想通りの曲線になりそうかな。

inpのCLOADと、datの座屈荷重との関係

いつもは、CLOADで(1/載荷節点数)と入力したものを、座屈荷重とみなしている。

もし10/載荷節点数とかにしたら、座屈荷重はどう変わるのか見てみた↓。

CLOAD	0.001/n	0.01/n	0.1/n	1/n	10/n	100/n
理論値との誤差	座屈せず	3%	3%	3%	3%	座屈せず

n:載荷節点数

10/nの座屈荷重は、1/nの座屈荷重の0.1倍になる。0.1/nにすれば10倍、0.01/nにすれば100倍の座屈荷重になる。

だから、inpのCLOADで0.1/nで入力したら、datファイルの座屈荷重に0.1かければよい、みたい

後藤ちゃちゃ08/11/14

ということは、座屈荷重が計算されるぶんには、どの倍率で計算しても精度的には大丈夫そうだということでいいかな。前に、梁の横ねじれ座屈をやったときは、倍率によって精度が変わったりしたような気もするけど、それはccxのバージョンのせいかな。 ccx_1.5とか古いバージョンでも同じ結果になるかな。まあ、これを信じるなら、刈屋さんのダイヤカットの座屈も信じてもいいのかも....

板の座屈(2次)

2次の座屈モード。162要素の時、誤差は2.7%になった。収束傾向にある。

画像

板の座屈

モード1の座屈。162要素の時、誤差は3%になった。収束傾向にある。真横から見た画像

斜め上から

画像

2枚の板がはじめから開いている状態から、引っ張った。

1枚あたり2要素の場合と、1枚あたり16要素の場合をくらべた。

16要素のほうが滑らかに開いている。

08/10/29

圧縮による座屈。→画像

NLGEOMの有無に関わらず、ccx1.7で計算すれば、frdファイルの無駄な節点はできてしまう。

08/10/28 WARNINGについて

6角形の、初期不整を与えた折り畳み円筒を引っ張った。

NLGEOMを付けた方を、ccxが計算しているときにでるWARNING。小さい荷重でもでる。
NLGEOMを付けない場合は、でない。

*WARNING in statics: a nonlinear geometric analysis is requested
          but no time increment nor step is specified
          the defaults (1,1) are used

6角形の、高さ0の折り畳み円筒を引っ張った。

NLGEOMを付けた方を、ccxが計算しているときにでるWARNINGとINFO。小さい荷重でもでる。
NLGEOMを付けない場合は、INFOだけでる。

 *INFO in gen3dnor: in some  nodes opposite normals are defined

*WARNING in statics: a nonlinear geometric analysis is requested
          but no time increment nor step is specified
          the defaults (1,1) are used

08/10/27 15時

6角形の、初期不整を与えた折り畳み円筒を引っ張った。

板厚1mm、半径5cm、高さ2cm。
拘束条件は、固定側x,y,z、引張側x,yの並進変位の固定。

横軸が変位(mm)、縦軸が力(N)。赤線がNLGEOMを付けた場合で、緑線がNLGEOMをはずした場合。

NLGEOMを付けた場合、40000N以上は、too many cutbacks。

6角形の、高さ0の折り畳み円筒を引っ張った。

板厚1mm、半径5cm。
拘束条件は、固定側x,y,z、引張側x,yの並進変位の固定。

横軸が変位(mm)、縦軸が力(N)。赤線がNLGEOMを付けた場合で、緑線がNLGEOMをはずした場合。

NLGEOMを付けた時、80000N以上でエラー。エラーは以下のようになる。

 *ERROR: increment size smaller than minimum
 best solution and residuals are in the frd file

08/10/27

2枚の板を重ねて高さ0の状態から思いっきり引っ張った。→アニメーション

さらに大きな荷重で引っ張るとこうなる

2枚の板がはじめから開いている状態から、引っ張った。アニメーション

さらに大きな荷重で引っ張るとこうなる

08/10/22

frdファイルについて

ccx1.7で作ったfrdファイルは、NLGEOMがついているいないにかかわらず、節点がよけいにでてくる。378個の節点のやつをccx1.7で計算し、frdファイルを見たら節点が4000くらいまで増えていた。

引張側も拘束した場合について

高さ0の円筒を固定側x,y,z、引張側x,yを固定して引っ張っても、うまく伸びない。

ccx1.5でも1.7でも、非線型を入れても入れなくても、14角形でも6角形でも、どんな組合せでも、うまくいかない。

かなり小さい初期不整を入れた場合も同じである。

08/10/21

固定側の全節点のx,y,zの並進変位を拘束したもの→アニメーション

高さ0の状態から引っ張った

6角形、4段、半径5cm、210GPa、板厚1mmの円管に、高さ0から0.5MNでひっぱった。→add-displacement

これは、ccx1.7で、NLGEOMをいれないで計算した。これに、すごく小さい初期不整をいれても、変形状態はほとんど変わらない。

08/10/20 15時

全節点数は378で382という節点はなかった。

08/10/20

「1点でx,y,θz拘束」という境界条件では、ccx1.7でNLGEOMをいれていないと、

*INFO in cascade: linear MPCs and
       nonlinear MPCs depend on each other
       common node: 382 in direction 2

 Factoring the system of equations using spooles

のようにでて計算できないが、NLGEOMを入れると、計算できた。

計算できたといっても、何度も上のINFOがでできて、時間がかかって計算が終わる。

これのアニメーションを見ると、回転もしないし上にも動かない。

でもadd-displacementを見ると、上に伸びているのが確認できた。

10/17 その1

高さ3cm、半径5cm、板厚1mm、ヤング率69GPaで、ccx1.7でのアニメーション。上に引っ張られないで横に動く。

同じ物をccx1.5でやると、しっかり上に動く。

ふたつとも幾何学非線型をいれていない。

ccx1.7に幾何学非線型をいれると、ccx1.5のように上に動くようになる。

後藤ちゃちゃ08/10/17

上に伸びている方も、固定端部の1点の回りに回転しているように見える。境界条件は、

1)「1点でx,y,θz拘束」ではなくて、
2)「直径を挟んで向かい合う2点で、それぞれx方向とy方向を拘束」でしたっけ。

だとすると、直径を挟む2点が、一方は半径方向にだけ動けて、他方は接線方向にだけ動けるので、半径方向に動ける方を回転中心にして、向かい側の点が接線方向に動くとやや回転しているように見えるということだろうか。 1)だとどうなるだろう。

08/10/17 14時

2)でやっていました。なので1)の拘束条件をccx1.7でやってみました。すると、

*INFO in cascade: linear MPCs and
       nonlinear MPCs depend on each other
       common node: 382 in direction 2

 Factoring the system of equations using spooles

のようにでて計算されませんでした。

10/17 その2

その1にあるように、69GPaで、ccx1.7の幾何学非線型を入れない状態でやると、おかしいので210GPaでやった。

高さ3cm、半径5cm、板厚1mm、ヤング率210GPa、ccx1.7でNLGEOMを入れないで、荷重0.1MNで引っ張った。

上と同じ物を荷重0.05MNで。

片持ち梁を変形(STEP,NLGEOMで)

後藤ちゃちゃ(08/10/10)

バージョン依存のエラーは、まとめると以下のような感じでしょうか。折れ角を持って接合されているシェル要素を引っ張ったり押したりした場合、

バージョン	幾何学非線形	出力結果
ccx_1.5	入れない	too many cutbacks
ccx_1.5	入れる	too many cutbacks
ccx_1.7	入れない	出る
ccx_1.7	入れる	too many cutbacks

幅5cm、軸長10cm、ヤング率69GPa、1要素で、板厚0.1mmの時。力は0.06Nで押している。0.6Nにすると、too many cutbacks。

幅5cm、軸長10cm、ヤング率69GPa、1要素で、板厚1mmの時。力は600Nで押している。6000Nにすると、too many cutbacks。

1段あたり10cmの初期不整をあたえた板

幅5cm、軸長14cmで、板厚0.1mmの時。力は0.6Nで引っ張っている。6Nにすると、too many cutbacksがでる。

後藤ちゃちゃ(08/10/9)

そうすると、折れ角が小さすぎることが原因ではないですね。変位(特に回転)が大きくなりすぎるせいかどうかは、片持ち梁を極端に変形させて確認したいですが、 STEPのところを

*STEP,NLGEOM

として幾何学非線型解析にして解いてみますか(そうすれば、変位や回転の大きさには制約がなくなりますので)。

幅5cm、軸長14cmで、板厚1mmの時。力は6Nで引っ張っている。60Nにすると、too many cutbacksがでる。

2枚の板を重ねたやつを上に引っ張る(高さ0の状態から)

幅5cm、軸長10cmで、板厚0.1mmの時。力は0.6Nで引っ張っている。6Nにすると、too many cutbacksがでる。

幅5cm、軸長10cmで、板厚1mmの時。力は6Nで引っ張っている。60Nにすると、too many cutbacksがでる。

2枚の板を重ねたやつを上に引っ張る(初期不整を与えた状態から)

幅5cm、軸長10cmで、板厚0.1mmの時。力は0.6Nで引っ張っている。6Nにすると、too many cutbacksがでる。

幅5cm、軸長10cmで、板厚1mmの時。力は6Nで引っ張っている。60Nにすると、too many cutbacksがでる。

10/3 15時

高さ0の折り畳み円筒を、大きい荷重で引っ張るときにでるエラーの文。

ERROR: too many cutbacks
best solution and residuals are in the frd file

ccxの増分繰り返しの制御

*CONTROLS,PARAMETERS=TIME INCREMENTATION
4,8,9,16,10,

10/3

アニメ

分布

z方向の変位の分布。

1点はx軸方向の並進変位を拘束し、その点に対象な点はy軸方向の並進変位を拘束した。

10/1

分布

z方向の変位の分布。

1点はx軸方向の並進変位を拘束し、その点に対象な点はy軸方向の並進変位を拘束した。

tasu.fの改造したやつ

tasu2.f

9/30 13時

アニメ

いままでの座屈した円管で、Animateのstartをおしてから、Viewingのadd-displacementを押したらこうなった。

9/26 13時

アニメ

円筒を引っ張った

cgxのmovie

8段のやつを引っ張った

BUCKLEができない。アニメ

厚さ0にしたらこうなった

cgxのmovie

gnuplot 要素数と変位

幅0.5m、桁高0.5m、軸長10mの片持ち梁の自由端に集中荷重が作用している。横軸を要素数、縦軸を自由端の鉛直変位[m]とする。緑の直線はティモシェンコ梁の式の変位で、青の直線はせん断無視の梁理論の変位。 32要素の変位が1.7609で、64要素の変位が1.7641になり、収束値を1.76にした。

要素数を32で、l/hをかえた。縦軸を誤差(%)、横軸をl/hとする。赤線が梁理論との誤差で、緑線がティモシェンコとの誤差である。

拘束条件、両端固定梁

軸方向がx、幅方向がy。固定支点の端の節点(Nfixx)は

Nfixx,1
Nfixx,4,6
8,1,6

cgxのmovie

拘束条件、単純梁

軸方向がx、幅方向がy。ヒンジ支点の端の節点(Nfixx)は、

Nfixx,1
Nfixx,4
Nfixx,6
8,1,4
8,6
Nfixx2,4
Nfixx2,6
102,2,4
102,6

cgxのmovie

08/6/27

幅0.5m、桁高0.5m、軸長10mの片持ちばりの曲げが、理論値との誤差-4%になった。幅0.5m、桁高0.01m、軸長10mの引張りは、誤差-0.9%になった。

稲荷の卒論日誌 - *座屈荷重-折れ角

目次

座屈荷重-折れ角

文字コード

支承について

角柱の角数-圧縮による座屈荷重

最小の高さ方向分割数を求める計算式

ble6ccxfについて

inpファイルの座標値について

ble6ccxfとble6ori3を、引張の座屈荷重で比較

円管を座屈

角数とバネ定数の関係

初期高さと座屈荷重の関係(20段)

座屈荷重ー初期高さ グラフ

引張による座屈荷重

圧縮による座屈荷重

cgxで座屈モードを見た(八角形)

cgxで座屈モードを見た(六角形)

cgxで座屈モードを見た(四角形)

初期高さ-ばね定数 関係

四段四角形の、初期高さと座屈荷重

座屈荷重

座屈モードをcgxで見る

ハードディスクの容量をチェックしたいときは

初期高さとばね定数の関係(四段の四角形)その3

hが四角柱の場合

hが円管の場合

初期高さとばね定数の関係(四段の四角形)その2

初期高さとばね定数の関係(四段の四角形)

初期不整による、荷重-変位の関係の変化(四段の四角形)

後藤ちゃちゃ08/11/13

inpのCLOADと、datの座屈荷重との関係

後藤ちゃちゃ08/11/14

板の座屈(2次)

板の座屈

2枚の板がはじめから開いている状態から、引っ張った。

08/10/29

08/10/28 WARNINGについて

6角形の、初期不整を与えた折り畳み円筒を引っ張った。

6角形の、高さ0の折り畳み円筒を引っ張った。

08/10/27 15時

6角形の、初期不整を与えた折り畳み円筒を引っ張った。

6角形の、高さ0の折り畳み円筒を引っ張った。

08/10/27

08/10/22

frdファイルについて

引張側も拘束した場合について

08/10/21

高さ0の状態から引っ張った

08/10/20 15時

08/10/20

10/17 その1

後藤ちゃちゃ08/10/17

08/10/17 14時

10/17 その2

片持ち梁を変形(STEP,NLGEOMで)

後藤ちゃちゃ(08/10/10)

1段あたり10cmの初期不整をあたえた板

後藤ちゃちゃ(08/10/9)

2枚の板を重ねたやつを上に引っ張る(高さ0の状態から)

2枚の板を重ねたやつを上に引っ張る(初期不整を与えた状態から)

10/3 15時

ccxの増分繰り返しの制御

10/3

10/1

9/30 13時

9/26 13時

円筒を引っ張った

8段のやつを引っ張った

厚さ0にしたらこうなった

gnuplot 要素数と変位

拘束条件、両端固定梁

拘束条件、単純梁

08/6/27

リンク

LaTeX

fortranメモ

稲荷の日誌

座屈荷重ー初期高さグラフ

初期高さ-ばね定数関係