後藤資料

柴田の修論日誌 - *ハニカム完成(一応形だけは)(12/3/4)

目次

ハニカム完成(一応形だけは)(12/3/4)

これはx方向，y方向分割数を増やしても（偶数であれば）大丈夫なプログラム。右の画像はx方向に16分割，y方向に10分割したときの結果。左の画像はx方向に10分割，y方向に20分割した結果。x方向の分割数よりy方向の分割数の方が大きくないとダメかと思ったが，逆でもイケル。まあとりあえず形だけはなんとか卒業までにできたので自分としては満足。 でもこれは，s3要素で作ったためobjファイル形式にして節点と節点の間の点を作ったほうがいいと思う。一発でハニカムが作れるプログラムだと楽だがs6要素で作るとなると点と点の間も自分で作らないといけないので結構大変だと思う。結局ダイヤカットや蛇腹折りのようにするのが楽かも・・でもそうなると節点番号を出したのも意味なくなるのだが・・。作ったのはいいがこれを使うこともなくいなくなるので少し寂しい気もする。まあ仕方がないが。

・・・・もし万が一使う人がいれば，この日誌の一番下を見てほしい。

暇なのでハニカム作ってみた(12/2/21)

1方向だけ（x方向だけかな？）の分割数を増やしても大丈夫なプログラムを作成． これを基礎としてy方向の分割数を増やせたらいいなぁ．もうちょっと時間があれば作りたいけどあんまりやるとM1に怒られそうなのでここらへんでやめておこうかなと思う（笑）あと正直飽きてきたのでここらが潮時だと思う．

暇なので(引越しとか掃除とかしたほうがいいのだろうけど・・・あと東北支部のスライド作ったり)ハニカムっぽいのを作ってみた．まだ1列しかないけど・・・．この先増やせていけたらカッコイイ感じになるのでは？自分はいないだろうけど．これはs4のシェル要素で構成されている．s4のシェル要素でもできることが分かったが，s8にした方がいいので考えなおす必要がある．objファイル形式にした方がいいのかもしれない．

シェル要素片持ち梁(プログラム)について(一応完成？)(12/2/19)

http://www.str.ce.akita-u.ac.jp/~gotouhan/shibata/kadai2.f

4年生の時につくったプログラムを移植して改造した．下に書いてある箇所を直そうと思って取り組む．もうちょっと時間がかかるかと思ったが，一応できた． このプログラムは拘束条件も載荷条件もちゃんとしてる・・と思う． 使わない節点まで出てきていたが，そこも修正済み．すぐに計算に使えると思う．だが，ヤング率(弾性係数)は自分で変える必要あり． 実際に使う人が確かめた後に使ったほうがいいと思う．必要であれば改造した方がいいかも・・・．上の画像はx方向分割数4のz方向分割数2の余計な節点がない片持ち梁(一応成功か？)

シェル要素片持ち梁(プログラム)について(12/2/18)

片持ち梁のような長方形断面のものはs6要素よりもs8要素の方が適しているような気がするのでs8要素で片持ち梁ができるようなものを作った．

形しかできていないので拘束条件，載荷条件はこれから足していくつもり． あと，片持ち梁の一番最後の列は使わないのに節点を出しているので消す必要あり．ホントはプログラムの中でしっかりできてればよかったんだけど，自分が4年の時に作ったものを移植しているのでそこら辺に不備が・・・ でも，x方向分割数とz方向分割数を同じにしなくても形をつくることができるのでどちらかの分割数を増やすことで傾向がみたいときはいいかも？ 上の画像はx方向分割数が4でz方向分割数が2のときだが，右端の節点は余計なので消す必要がある．(失敗作)

http://www.str.ce.akita-u.ac.jp/~gotouhan/shibata/kadai.f

このプログラムであれば片持ち梁をどれだけ細長くして分割を増やしても計算可能． 自分がやったのは軸長1m，桁幅2mmで分割数x方向40，z方向40でも大丈夫だった．心配な人はやってみたほうがいいかも． ただx方向及びz方向の分割数は同じにしないとダメなところは相変わらずなので注意が必要．x方向,z方向どちらかの分割数を増やして傾向をみたいときなんかは困るのでつくり変える必要があるかも？． 新しくしたプログラムの方はb=の部分をなくした形になっている．

ちょっとした疑問？

4年生が使っていたプログラムの中で，細長くして分割数を増やすとおかしくなる現象が起きていたので調査を開始．その中で気になった点を書く． 幅と軸長と厚さ？桁高？を最初の方に入力しているが，シェル要素で桁高は厚さと同じような意味だし，プログラムを見てもその後の計算で桁高を使っている感じもないので入れる必要があるのか？とちょっと疑問に思った． その後厚さは入力するので二重で入力していることにはならないのか？ 上のプログラムの中ではb=を入力しているし， 
*SHELL SECTION,ELSET=Eall,MATERIAL=EL
の後にも入れているのでb=の値はいるのか？

交番載荷について(12/1/30)

*RESTART
*RESTART，WRITE，FREQUENCY=1

を使っていく．

Calculix上でのリスタート(再計算)のやり方(詳細)

発表練習だけやってても暇なのでリスタートについて書いておく．これがあればすこしは研究の幅が広がるかもしれない，そうでもないかもしれないけど．

1，計算するreiwrite.inpファイル(名前は何でもいい)の中に

　　　*RESTART

　　　*RESTART，WRITE，FREQUENCY=1

　　を入れて計算を行う.その時にreiwrite.routファイルが出力される.

2，リスタートファイルreiread.inp(これも名前は何でもいい)を作っておく，その中身は

　　　*RESTART，READ，STEP=1

　　　*RESTART，WRITE，FREQUENCY=1 を書いて，後は荷重の設定をしておく．

　　　ここで節点情報や要素情報はいらず(rinファイルから読み込んでいるからか) 荷重条件だけでいい．ついでに拘束条件もいらない．

3，リスタートファイルはrinファイルを読み込んでその情報を基にリスタートを行っているので，reiwrite.routをreiread.rinにコピーしておく．

4，コピーした後，reiread.inpの計算を行う．そのときreiread.routが出力される.

5，その後またリスタートをする場合，reiread.routをreiread.rinにコピーする． そしてreiread.inpの中身を

　　　*RESTART，READ，STEP=2

　　　*RESTART，WRITE，FREQUENCY=1

　　　にしてSTEP=の数字を増やす．

6，何回か繰り返してリスタートをする場合には，5の作業を繰り返してSTEP=の数字を 増やしていく．

交番載荷の荷重制御(11/11/29)

4角形

6角形

8角形

文献等

交番荷重の制御など

Calculix2.3による解析結果との比較(11/11/20)

リスタートの計算結果

リスタートまでは1.7と2.3で違いは出なかった。しかし、その後のリスタートの計算が1.7ではできたが、2.3はできなかった。動解析は難しいのかもしれない．他の人は動解析ではない結果も違っていたとかいなかったとかの報告をちらほら聞くので正直よく分からないが自分は単純な計算結果は1.7と2.3で同じ結果となった．
ccx_1.7ccx_2.3
折りたたみ円筒○○
折りたたみ円筒(リスタート)○×
片持ち梁○○
片持ち梁(リスタート)○×

蓋?板付き折りたたみ円筒(11/11/7)

載荷の向きを変えた交番載荷(6角)(11/11/4)

切れ目なし円筒

切れ目あり円筒

初期高さを変えた交番載荷(切れ目あり)(11/10/13)

4角形の場合(切れ目あり)

6角形の場合(切れ目あり)

8角形の場合(切れ目あり)

初期高さを変えた交番載荷(11/10/12)

4角形の場合

6角形の場合

8角形の場合

xfigのエラーについて(11/9/16)

 グラフがxfigで加工できなかったが、viで加工したいfigファイルを開いて

# Begin plot #1
# End plot #1

 の前と後にある

6 1612 1351 6974 4550
-6

 を消すか、コメント分にするとxfigのEditで加工できる。

切れ目のある円筒の厚さを変化(4角)(11/9/9)

切れ目のある円筒の厚さを変化(6角)(11/9/9)

切れ目のある円筒の厚さを変化(8角)(11/10/11)

厚さ1mmのリスタートができない(11/9/5)

切れ目のありなしの比較(11/8/31)

切れ目のある円筒のリスタート結果(11/8/31)

切れ目のある円筒(4角)のリスタート(11/8/30)

厚さを変えて比較

厚さを変えて比較(8角)(11/9/5)

厚さを変えて比較(6角)(11/9/4)

厚さを変えて比較(4角)(11/8/28)

円筒との比較

円筒との比較3(11/8/25)

円筒との比較2(11/8/2)

円筒との比較(11/7/27)

20角形の円筒(円筒との比較)(11/7/26)

ダイヤカットのプログラムを変えて出力

伸縮できる円筒のプログラムにて出力

8角の円筒のリスタート計算結果(11/7/25)

6角の円筒のリスタート計算結果(11/7/20)

4角の円筒のリスタート計算結果(11/7/15)

4角と6角と8角の円筒(11/7/9)

RESTART(伸縮できる円筒ver)(11/7/6)

RESTARTについて2(11/6/28)

RESTARTについて(11/6/9)

*RESTART,WRITE,FREQUENCY=1
を入れて計算すると~.routが出力される。 
*RESTART,READ,STEP=1
を入れて計算すると~.rinが出力される。

RESTARTのエラーについて(11/6/9)

WRITE→READの順に計算すると思う。 しかしそのまま計算すると、 
*ERROR in restartshort: requested step
        is not in the restart file
のエラーが出てしまうので、 ~.routの中を~.rinにコピーしたが、 
*ERROR in nodes: *NODE should be placed
   before all step definitions
のエラーが出てしまって計算できない。

荷重と変位(初期高さ0.9 円周が40cm)(11/5/26)

切れ目なし

切れ目あり

切れ目に対して垂直に載荷(11/5/30)

切れ目に沿って載荷(11/5/31)

初期高さと$\frac{k}{k_円}$の関係

切れ目なし(高さ方向分割数10)(11/5/25)

切れ目あり(高さ方向分割数10)(11/5/25)

高さ方向分割数を少なくするとうまく出力できない(11/5/17)

切れ目に垂直な方向に載荷したとき(11/5/13)

切れ目に沿う方向に載荷したとき(11/5/13)

座標変換プログラム

4,6,8角形の初期高さが0.9の円筒(高さ方向分割数が10ver)

切れ目なしの荷重と変位の関係(11/4/25)

切れ目ありの荷重と変位の関係(11/4/25)

4,6,8角形の切れ目のある円筒(4/3)

切れ目のある円筒(再計算)(3/21)

とりあえず4角形を再計算。(3/17)ここからM2の作業です．

tex

高さ方向分割数が少ないときちんと描けない

高さ方向分割数と高さの関係(11/01/13)

周方向分割数と円周の関係(11/01/13)

木橋関係

節点を共有させるプログラム2(10/12/13)

節点を共有させるプログラム(10/12/12)

ダイヤカットの座屈関係

材料諸元は、

半径(cm)高さ(cm)ヤング率(アルミ)(MPa)ポアソン比厚さ[mm]
1650700000.3000000121

座屈に関しては、

座屈モード精度倍率
10.10.003

inpのCLOADと、datの座屈荷重(10/11/30)

CLOAD10/n1/n0.1/n0.05/n0.01/n0.001/n0.0001/n
座屈荷重座屈せず5.8705220.5870469×$10^{2}$0.1174076×$10^{3}$0.5870234×$10^{3}$0.5872686×$10^{4}$座屈せず

座屈荷重については、ここにかいてあるので参考に。

精度について

荷重が0.01/nのとき

精度0.010.0010.00010.00001
座屈荷重0.5870234×$10^{3}$0.5870234×$10^{3}$0.5870234×$10^{3}$0.5870234×$10^{3}$

モードの倍率 0.003で描いてみた。(10/11/24)

返答と分からないこと(10/11/24)

後藤ちゃちゃ(10/11/24)

さっき、江村さんとこで、ccx_1.7+tasudatw.fで書かせようとしたときは、 no buckling occursエラーでしたが、どうやったらモード出せましたか。

なるほど、そうすると、ここにあるように、倍率を 調整すれば、ccx_1.7+tasudatw.fでモード図をかけるということですね。

伸縮できる円筒の塑性関係

*DEFORMATION PLASTICITY
210000.,0.3,235,100,0.5

後藤ちゃちゃ(10/11/15)

アクセス制限下のccxマニュアル中の画像を引用するとパスワードを求められるので、 ここに式を書けばいいのでは?

$E\epsilon=\sigma+\alpha(\frac{|\sigma |}{\sigma_{0}})^{n-1}\sigma$

指数っていうのは、この(n-1)の部分を指定できるということかな。

そうだと思います。 nが1だと$(\frac{|\sigma |}{\sigma_{0}})^{n-1}$の部分が常に1になって、 一番下のグラフみたいに直線になります。 

*DEFORMATION PLASTICITY
210000.,0.3,800.,12,0.4

salome関係

salomeのプログラム(11/01/14)

cgxでsalomeの円柱を出力

3 102 4 1
1 102 1 3
9 203 8 6 2
19 203 9 4 3

2つのdatファイルのデータを1つのinpファイルにするプログラム

変位が大きくでてしまう理由について

折り紙構造関係

波形鋼板(両端まっすぐver)(10/11/2)

open(7,file='yousosuu.txt') !yousosuu.fで数えた要素数ファイル
read(7,*) nyou !要素数の読み込み
の下に 
read(7,*) m,n,l
いれておけば、ちょっとは楽になる。(ble6ccxをいちいちコンパイルしなくて済む)

波がいくつかある波形鋼板(10/9/17)

まっすぐ切ったダイヤカット円筒 (10/9/8)

波形鋼板ウェブ(1/4ver) (10/9/8)(10/9/10)

波形鋼板ウェブ(1/2ver)とまっすぐ切ったダイヤカット（10/7/30）

応力-ひずみ曲線2

s6要素の片持ち梁を引っ張ったときの応力-ひずみ曲線 

*DEFORMATION PLASTICITY
210000.,0.3,800.,12.,0.4
を入れて計算すると上図のような結果になった。

ダイヤカット円筒（半分）

ダイヤカットを切って半分にしてみた。が少し違う?

応力-ひずみ曲線

*PLASTIC
235.,0.
245.,0.5
をいれて計算してみたところ、応力が200MPaを超えると下記のエラーが表示される。 
*ERROR: increment size smaller than minimum
best solution and residuals are in the frd file

2010後期分

時間月火水木金
1.2時限8:50-10:20
3.4時限10:30-12:00環境材料学(ゼミ)
5.6時限12:50-14:20物理学実験(TA)技術者倫理特論Ⅱ木構造工学
7.8時限14:30-16:00物理学実験(TA)構造設計学特論

2010前期分

時間月火水木金
1.2時限8:50-10:20建設構造工学
3.4時限10:30-12:00構造力学特論材料設計学交通システム計画特論
5.6時限12:50-14:20土質工学特論物理学実験(TA)ゼミ
7.8時限14:30-16:00物理学実験(TA)土木環境工学実験(TA)

構造実験(アクリル片持ち梁のたわみ試験)

http://www.str.ce.akita-u.ac.jp/~gotou/zikken/

hanikamu.f

ここにくるなんて物好きッスね。それか相当焦っているか。一応、ここにハニカム(柴田)を置いていくので利用したい人は利用してみて。しかし、このプログラムでできているのは節点番号と要素番号だけで拘束条件、載荷条件は正しくないので作る必要あり。

と、この最後の文は旅立つ朝7時頃に書いております。ねむ・・くはあまりない。 この3年間(大学の半分も)研究室にお世話になりました。色んなことがここであったような気がします。プログラムができなくてあ〜ってなったりtexの図が思うようにでなくてあ〜ってなったり4年の頃卒論をまとめる時期深夜に研究室で変にテンション上がったり最初人数いなくて最近は増えたなぁ〜とか思ったりあんまり書くと涙腺が・・・(嘘です)。 M1とは2年一緒にいたけど、いい子達でしたしもっとしゃべってみたかったかな。 4年生とは1年だったけどおとなしい子達で(自分にはそう見えたが)やっぱりもっとしゃべってみたかったかな。もっと教えてあげられるとことか手伝ってあげられるとこがあったように思って後悔してることとかもあります。そんなことないよって言ってもらえたら泣いちゃいますけど。あんまり書くと恥ずかしいこととか書きそうになるのでこんな最後にこっそり書いてます。ここを見る人が何人いることやら。

え〜と、とにかく全部ひっくるめて楽しかったです。

ありがとうございました。

Powered by YukiWiki 2.1.3
Modified by GOTOU Humihiko.

*