#author("2020-01-20T12:29:32+09:00","default:kouzouken","kouzouken")
#contents
*卒論日誌 [#d9858f5b]
,4/11,14:30~16:00,1.5,ブラインドタッチ
,4/15,14:30~16:00,1.5,ブラインドタッチ
,4/16,12:00~13:30,1.5,ブラインドタッチ
,4/17,13:00~14:00,1.0,ブラインドタッチ、VI
,5/13,15:00~16:00,1.0,VI
,5/16,14:00~15:00,1.0,VI
,5/27,11:30~13:30,2.0,VI
,6/6,15:00~16:00,1.0,VI
,6/7,10:00~16:30,6.5,VI
,6/9,19:00~22:00,3,VI
,6/16,18:00~2:00,8,salome
,6/17,11:00~13:00,2,salome
,6/18,13:00~15:00,2,VI
,6/20,17:00~18:30,1.5,VI
,7/12,13:00~7:00,4.0,tex
,7/15,14:00~18:00,4.0,tex
,7/16,14:00~16:00,2.0,tex
,7/18,17:00~21:00,4.0,tex
,7/19,10:00~12:00,2.0,tex
,8/6,12:00~16:00,4.0,tex
,9/2,15:00~18:00,3.0,tex
,9/6,13:00~16:00,3.0,tex
,9/9,13:00~17:00,4.0,tex
,9/20,13:00~16:00,3.0,tex
,9/25,13:00~16:00,3.0,tex
,10/3,9:00~13:00,4,tex
,10/4,13:00~16:00,3,tex
,10/8,13:00~16:00,3,tex
,10/9,13:00~16:00,3,tex、VI
,10/10,13:00~16:00,3,tex、VI
,10/15,13:00~16:00,3,tex、VI
,10/16,13:00~16:00,3,tex、VI
,10/17,13:00~16:00,3,tex、VI
,10/23,13:00~17:00,4,tex、VI
,10/24,13:00~17:00,4,tex、VI
,10/25,13:00~17:00,4,tex、VI
,10/28,13:00~17:00,4,tex、VI
,10/30,13:00~17:00,4,tex、VI
,10/31,13:00~17:00,4,tex、VI
,11/1,13:00~18:00,5,tex、VI
,11/5,13:00~18:00,5,tex、VI
,11/6,13:30~18:30,5,計算、VI
,11/7,10:00~17:00,7,計算
,11/8,12:30~18:30,6,計算、VI
,11/11,13:00~17:00,4,tex、VI
,11/12,13:00~17:00,4,tex、VI
,11/13,10:00~16:00,6,tex
,11/14,13:00~16:00,3,計算、VI
,11/14,14:30~16:30,2:30,3dプリンター関連
,11/15,15:00~17:30,2:30,3dプリンター関連
,11/19,13:00~16:30,3:30,3dプリンター関連
,11/20,13:00~16:30,3:30,3dプリンター関連
,11/21,11:00~15:30,3:30,3dプリンター関連
,11/22,11:00~18:30,6:00,3dプリンター関連
,11/25,14:00~18:30,4:30,3dプリンター関連
,11/27,13:00~17:00,4:00,3dプリンター関連
,12/2,13:00~17:00,4:00,3dプリンター関連
,12/3,13:00~18:00,5:00,3dプリンター関連
,12/4,13:00~17:00,4:00,3dプリンター関連
,12/8,13:00~17:00,4:00,3dプリンター関連
,12/9,12:00~16:00,4:00,比較用グラフづくり
,12/10,12:00~18:00,6:00,比較用グラフづくり
,12/11,12:00~19:00,7:00,有限要素によるばね定数計算
,12/12,12:00~19:00,7:00,有限要素によるばね定数計算
,12/13,12:00~20:00,8:00,有限要素によるばね定数計算
,12/14,12:00~18:00,6:00,有限要素によるばね定数計算
,12/15,13:00~19:00,6:00,有限要素によるばね定数計算
,12/16,11:00~20:00,9:00,ばね定数計算&実験準備
,12/17,10:00~18:00,8:00,ばね定数計算&実験
,12/18,12:00~18:00,6:00,実験装置製作
,12/19,10:00~20:00,10:00,実験装置製作
,12/21,12:00~20:00,8:00,実験装置製作
,12/22,14:00~20:00,6:00,実験装置製作
,12/23,12:00~18:00,6:00,実験装置製作
,12/24,9:30~19:00,9:30,実験
,12/25,9:30~22:00,12:30,中間発表スライド製作
,1/6,12:30~18:00,5:30,蛇腹折りstl製作、卒論概要
,1/7,10:00~17:00,7:00,蛇腹折り製作、座屈モード図製作
,1/8,10:00~18:00,8:00,座屈モード図製作
,1/9,13:00~20:00,7:00,座屈モード図製作
,1/10,9:30~18:00,4:00,実験
,1/12,14:00~18:00,4:00,実験したモデルの座屈モード図を製作
,1/14,14:00~21:00,7:00,パターン数、変位の3Dグラフ製作
,1/15,14:00~21:00,7:00,破壊試験 パターン数、変位の3Dグラフ製作
,1/16,14:00~21:00,7:00,破壊試験 パターン数、変位の3Dグラフ製作
,1/20,14:00~19:00,5:00,卒論概要、実験まとめ
,1/21,14:00~18:00,4:00,卒論概要、実験まとめ
,1/22,13:00~20:00,7:00,卒論概要、実験まとめ
,1/23,13:00~20:00,7:00,卒論概要
,1/24,14:00~20:00,6:00,卒論概要
,1/25,12:00~24:00,12:00,卒論概要
,1/26,14:00~18:00,4:00,卒論概要
,1/27,11:00~20:00,9:00,卒論概要
,1/28,9:00~15:00,6:00,卒論概要
,1/29,11:00~24:00,13:00,卒論概要
,1/30,24:00~15:00,15:00,卒論概要
,1/31,12:00~19:00,5:00,発表スライド
,2/1,14:00~20:00,6:00,発表スライド
,2/2,14:00~19:00,5:00,発表練習
,2/3,10:00~15:00,5:00,発表練習
,合計時間,486h

*今後の計画 [#g4a724b2]
卒論のパワーポイント製作、発表練習
*使用プログラム [#jfc24bda]

&link(jabaram4n6.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabaram4n6.inp)

&link(jabaram6n8.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabaram6n8.inp)

&link(zakutum6n10a7.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/zakutum6n10a7.inp)

&link(zakutum6n10a9.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/zakutum6n10a9.inp)

&link(zakutum6n10a99.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/zakutum6n10a99.inp)

*モデル製作 [#x3011528]
orita4ziku.f90をviで開きパターン数や高さ円周を決定し、コンパイル実行。

そうするとobjファイルができる。
それをble6oritata.f90で載荷条件や材料の条件を決定。
以下のようにコンパイルするとinpファイルができる。
このinpファイルはヤング率や荷重の数値が桁数が大きくなりすぎてい
るのでccxで計算すると変になる。なので、有効数値何桁かに決めて
調整したほうが良い。

例えばble6のファイルでヤング率4Gpaとすると
inpファイル内では4.00000000134887992のように
なっているので13~を消して保存する。

~$ gfortran -o orita orita4ziku.f90 

~$ ./orita   <-これによりorita4ziku.objが生成される

~$ gfortran -o ble ble6oritata.f90

~$ ./ble <orita4ziku.obj >jabaraMNA.inp

*破壊試験 [#hc78e19b]
CalculiXより座屈荷重は以下のようになる

,初期高さ,0.7,0.9,0.99
,座屈荷重[N],296,402,368

実験結果

,初期高さ,0.7,0.9,0.99
,座屈荷重[N],118,225,472
-破壊後
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/DSCN3297.JPG

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/DSCN3286.JPG

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/DSCN3391.JPG
-破壊直前
--周方向6高さ方向10初期高さ0.9
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabaram6n10a09-hakai.png
--周方向6高さ方向10初期高さ0.99
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabaram6n10a099-hakai.png

-横軸が変位(mm)で縦軸が荷重(N)
--周方向6高さ方向10初期高さ0.7
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/hakaim6n10a07.png

&link(hakaim6n10a07.dat,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/hakaim6n10a07.dat)

--周方向6高さ方向10初期高さ0.9
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/hakaim6n10a09.png

&link(hakaim6n10a09.dat,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/hakaim6n10a09.dat)

--周方向6高さ方向10初期高さ0.99
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabarazikken3.png

&link(hakaim6n10a099.dat,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/hakaim6n10a099.dat)

*CalculiXによる計算結果 [#ya6f5936]
実験と同じように10kgを載荷した際の高さ方向変位をパターン毎に表したもの。1列目が周方向分割数で2列目が高さ方向分割数、3列目が変位

&link(jabaraheni.dat,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabaraheni.dat)

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabaraheni.png
-座屈荷重3dplot

&link(jabarazakutukaju.dat,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara-zakutukaju.dat)

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara-zakutukaju.png

変位0の地点は製作出来なかった。

グラフ四隅
&link(ori4-4.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/ori4-4.inp)

&link(ori4-20.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/ori4-20.inp)

&link(ori14-18.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/ori14-18.inp)

&link(ori14-20.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/ori14-20.inp)
*実験したパターン数のinpファイル [#p27e017d]
&link(ori6-10-7.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/ori6-10-7.inp)

&link(ori6-10-9.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/ori6-10-9.inp)

&link(ori6-10-99.inp,http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/ori6-10-99.inp)

*実験結果 [#jb627875]
**12/17 [#f956a684]
 500g載荷
,経過時間,0分,5分,10分,15分,20分,25分,30分
,変位,4.65,4.51,4.48,4.46,4.45,4.44,4.43

**12/24 [#x6417584]
 10kg載荷 初期高さ0.9
,載荷,1kg,2kg,3kg,4kg,5kg,6kg,7kg,8kg,9kg,10kg
,変位,0.25,0.50,0.69,0.91,1.11,1.33,1.56,1.78,2.07,2.32

,経過時間,0分,5分,10分,15分,20分,25分,30分
,変位,2.32,2.86,3.10,3.28,3.44,3.58,3.70

,除荷,0分,5分,10分,15分,20分,25分,30分
,変位,3.70,0.96,0.81,0.72,0.67,0.62,0.60

バネ定数[除荷手前30分時]26486
k/k円=0.031621
ただの円筒よりk/k円=0.062161

**12/24 [#yf3741e9]
 10kg載荷 初期高さ0.7
,載荷,1kg,2kg,3kg,4kg,5kg,6kg,7kg,8kg,9kg,10kg
,変位,0.02,0.12,0.63,1.22,1.69,2.25,2.64,3.39,4.10,4.56

,経過時間,0分,5分,10分,15分,20分,25分,30分
,変位,4.56,7.04,8.01,8.74,9.60,10.25,10.66

,除荷,0分,5分,10分,15分,20分,25分,30分
,変位,10.66,1.81,1.40,1.16,1.01,0.90,0.82

バネ定数[除荷手前30分時]9193
k/k円=0.010975
ただの円筒よりk/k円=0.021575

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabarasiken.png
試験中に蛇腹の端の部分が歪んだ
**1/10 [#ta6328af]
 10kg載荷 初期高さ0.99
,載荷,1kg,2kg,3kg,4kg,5kg,6kg,7kg,8kg,9kg,10kg
,変位,0.13,0.23,0.32,0.44,0.51,0.63,0.71,0.75,0.83,0.99

,経過時間,0分,5分,10分,15分,20分,25分,30分
,変位,0.99,1.27,1.39,1.49,1.57,1.69,1.70

,除荷,0分,5分,10分,15分,20分,25分,30分
,変位,1.70,0.52,0.46,0.42,0.40,0.38,0.36

バネ定数[除荷手前30分時]57647
k/k円=0.068824
ただの円筒よりk/k円=0.13529

*座屈モード [#pe1d1723]
**3Dプリンターで製作、実際試験したものと同じ条件 [#sf0d393e]
上から周方向6高さ方向10初期高さ0.7荷重98N座屈モード倍率0.001

周方向6高さ方向10初期高さ0.9荷重98N座屈モード倍率0.001

周方向6高さ方向10初期高さ0.99荷重98N座屈モード倍率0.001
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara6107.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara6109.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara61099.png
**周方向4高さ方向6~20[2刻み]荷重98N座屈モード倍率0.001 [#sc35fcd4]
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara467.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara487.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara4107.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara4127.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara4147.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara4167.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara4187.png

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara4207.png

**周方向6高さ方向6~20[2刻み]荷重98N座屈モード倍率0.001 [#vecacbd3]
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara667.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara687.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara6107.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara6127.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara6147.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara6167.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara6187.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara6207.png
**周方向8高さ方向10~20[2刻み]荷重98N座屈モード倍率0.001 [#odf31c34]

高さ方向分割数6と8は節点の位置的に製作不可
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara8107.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara8127.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara8147.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara8167.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara8187.png
http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/jabara8207.png

*3Dプリンター [#u868c175]
 昭和アルミ缶の190ml缶の寸法に合わせて蛇腹円筒を作成
 予備試験をして試験内容を決定する。
 蛇腹円筒を製作。厚さを0.8mmにしたが斜め部分では厚さが0.4~0.5mmと半分程度になってしまっている
 蛇腹の初期高さを0.7から0.9へと変更
 それによる斜め部分の厚さの変化が無いように厚さを変更0.8から0.622…

 初期高さ0.9 厚さ0.4 ヤング率1.3GPa にて製作

 http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/hcpy_1.png
**養生2日目の5分ごとのヤング率 [#c8a1d6fc]
開始0秒のひずみで初期化したもの
,,5分,10分,15分,20分,25分,30分,平均(GPa)
,たてヤング率,1.468,1.429,1.404,1.429,1.455,1.496,1.447
,よこヤング率,-2.222,-1.839,-1.601,-1.391,-1.312,-1.212,-1.596

開始30秒のひずみで初期化したもの
,,5分,10分,15分,20分,25分,30分,平均(GPa)
,たてヤング率,1.468,1.429,1.404,1.429,1.455,1.496,1.447
,よこヤング率,-2.222,-1.839,-1.600,-1.391,-1.312,-1.212,-1.596
**バネ定数 [#e0ded3c9]
ヤング率1.3GPa 厚さ0.4mm 載荷荷重0.5kg 周方向分割数6 高さ方向分割数10 形状は昭和アルミニウム缶の190ml缶を基準とした蛇腹折り円筒をcalculixで計算しバネ定数を求める。
,初期高さ,ばね定数
,1,777777
,0.95,105150
,0.9,62420
,0.85,43750
,0.8,33108
,0.75,26344
,0.7,21585
,0.65,18148
,0.6,15457
,0.55,13315
,0.5,11556
,0.45,10082
,0.4,8844
,0.35,7852
,0.3,6911
,0.25,6102
,0.2,5378
,0.15,4666
,0.1,3920
,0.5,2882

*ポアソン比 [#b24a9be1]
ポアソン比=横ひずみ/縦ひずみ
,,1分,5分,10分,30分,1時間,2時間,6時間,24時間-30秒
,造形直後,0.48076,0.55635,0.68184,0.84513,0.97130,1.0905,1.4565,1.1130
,養生2日,0.53579,0.66064,0.77690,1.2337,1.5962,1.9718,2.6667,3.7284
,養生7日,0.51336,0.57151,0.62119,0.83037,0.097879,0.24472,1.3779,0.90340

*6/17課題 [#u8182d2b]
E=206GP ν=0.3 -b=0.2m h=0.2m
理論値12.515mm

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/kekka4.png
-横軸:lengh
-縦軸:たわみ

*課題 [#va76bcbd]
E=206GP ν=0.3 片持ち梁に1MNの荷重をかけた際のたわみを分割数を増やしながらccx_2.3で解く。またl=1mとし、断面が縦長、正方形、横長とした場合についても同様の計算をする。今回はnx=10,ny=20,としnzの分割数を増やしていった。

**縦長断面 [#x9640da3]
-b=0.1m h=0.2m
理論値25.029mm

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/hoge.png

-横軸:要素数
-縦軸:たわみ

**正方形断面 [#p2d388ed]
-b=0.2m h=0.2m
理論値12.515mm

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/hoge2.png

-横軸:要素数
-縦軸:たわみ

**横長断面 [#u3746213]
-b=0.4m h=0.1m
理論値48.922mm

http://www.str.ce.akita-u.ac.jp/~gotouhan/tsuruta/hoge3.png

-横軸:要素数
-縦軸:たわみ



トップ   編集 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS