![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
#author("2020-01-20T12:31:16+09:00","default:kouzouken","kouzouken") #contents *卒論日誌 ※忘れないうちに気づいたことは卒論日誌にメモする！ [#p57c9c32] ,日付,時間帯,作業時間,内容,立会 ,4/10,14:30-16:00,1.5,パソコン練習, ,4/11,14:00-15:00,1.0,パソコン練習, ,4/14,14:30-15:30,1.0,パソコン練習, ,4/15,13:30-15:00,1.5,パソコン練習, ,4/16,14:00-16:00,2.0,パソコン練習, ,4/17,15:00-17:00,2.0,パソコン練習, ,4/18,14:00-17:00,3.0,パソコン練習, ,4/19,14:00-16:30,2.5,パソコン練習, ,4/20,14:30-17:00,2.5,パソコン練習,斉藤 ,4/21,9:30-11:00,1.5,パソコン練習, ,4/21,14:30-16:30,2.5,プログラミング練習, ,4/21,17:00-19:30,2.5,構造力学, ,4/22,14:30-16:30,2.5,プログラミング練習, ,4/23,14:30-16:30,2.0,パソコン練習, ,5/5,17:00-23:00,6.0,構造力学, ,5/7,14:00-15:30,1.5,パソコン練習, ,5/11,20:30-22:30,2.0,構造力学, ,5/14,14:00-17:00,3.0,Salome練習＋課題（途中まで）, ,5/15,14:30-16:30,2.0,Salome練習（解析は終了）, ,5/15,20:00-23:00,3.0,構造力学, ,5/16,13:30-16:00,2.5,Salome練習（まとめ1）, ,5/17,14:00-15:00,1.0,Salome練習（まとめ2）, ,5/19,16:30-18:00,1.5,グラフ作成、貼付け, ,5/19,23:30-24:30,1.0,構造力学, ,5/20,13:30-15:00,1.5,Salome課題, ,5/21,13:30-16:00,2.5,Salome課題, ,5/22,14:30-16:00,1.5,Salome課題, ,5/22,23:00-24:00,1.0,構造力学, ,5/23,13:30-16:00,2.5,Salome課題, ,5/24,18:00-24:00,6.0,構造力学, ,5/26,14:00-17:00,3.0,Salome課題, ,5/27,13:30-16:30,3.0,Salome課題, ,6/2,16:00-18:00,2.0,Tex課題, ,6/2,22:00-25:00,3.0,構造力学, ,6/3,14:00-16:30,2.5,Tex課題, ,6/4,14:00-16:30,2.5,Tex課題, ,6/5,15:30-17:00,1.5,Tex課題, ,6/13,16:00-17:00,1.0,Salome課題, ,7/6,14:00-17:00,3.0,Tex練習, ,7/8,14:00-17:00,3.0,箱型断面の解析, ,7/13,18:30-19:30,1.0,箱型断面の解析, ,7/17,16:00-22:00,6.0,3点曲げと半解析, ,8/28,15:00-17:00,2.0,Mentat練習,尾山 ,9/2,12:30-17:00,4.5,Mentatで接触解析,尾山 ,9/3,11:00-15:30,4.5,スライド作成, ,9/4,13:00-17:30,4.5,Mentatで接触解析&スライド作成,尾山 ,9/16,12:00-14:30,2.5,摩擦試験方法の詳細を調べる, ,9/18,13:00-16:00,3.0,Mentatで接触解析,尾山 ,9/22,10:30-12:00,1.5,Mentatで接触解析, ,9/25,10:30-12:00,1.5,Mentatで接触解析, ,9/25,14:00-16:00,2.0,Mentatで接触解析, ,9/29,10:30-12:30,2.0,Mentatで接触解析, ,9/29,15:00-19:00,4.0,Mentatで接触解析,尾山 ,9/30,10:30-12:30,2.0,Mentatで接触解析, ,9/30,15:00-19:30,4.5,Mentatで接触解析,尾山 ,10/1,10:30-12:30,2.0,Mentatで接触解析, ,10/1,14:00-17:00,3.0,Mentatで接触解析,尾山 ,10/2,10:30-12:30,2.0,スライド作成&後藤さんに質問, ,10/2,14:00-24:30,10.0,Mentatで接触解析&スライド作成,尾山 ,10/6,14:00-15:00,1.0,摩擦試験の文献を読む, ,10/7,14:00-18:00,4.0,摩擦試験の文献を読む&Mentatで接触解析&荷重と水平変位の関係を調べる,尾山 ,10/8,14:00-20:00,6.0,摩擦係数について調べる→静止摩擦係数,尾山・斉藤 ,10/9,12:00-20:00,8.0,摩擦試験の水平変位についてMentatで接触解析&グラフの作成,尾山・斉藤 ,10/15,13:30-18:30,5.0,摩擦試験をMentatで解析→力をプレストレスにすれば上手くいく感じ・・・,尾山 ,10/16,10:00-12:00 13:30-15:00,3.5,摩擦試験をMentatで解析&結果をグラフにまとめる,尾山 ,10/20,13:30-18:00,4.5,Mentatで鉄の引張試験 ※午前中に研究室の模様替えをする,尾山・斉藤・坪井 ,10/21,10:00-11:30,1.5,Mentatで鉄の引張試験 ※その後、3年生が使うパソコンの準備,斉藤・坪井 ,10/22,14:00-19:00,5.0,Mentatで鉄の引張試験 ※やりながら坪井に引張試験のやり方を教わる,尾山・坪井 ,10/23,10:00-12:00 13:00-17:00,6.0,Memtatで鉄の引張試験 弾塑性解析(失敗?) 坪井から引張試験の方法を教わる,坪井 ,10/27,10:30-12:00 14:00-17:00,5.5,Mentatで弾塑性解析→とりあえずいったかな,尾山 ,10/28,10:00-12:30 14:00-18:30,7.0,Mentatで弾塑性解析,尾山 ,10/29,10:00-11:00 14:30-17:00,3.5,Mentatで弾塑性解析,尾山・斉藤 ,10/31,13:00-16:00,3.0,森吉で木橋見学,尾山・斉藤・石坂 ,11/4,13:00-16:00,3.0,Mentatで弾塑性解析, ,11/5,14:00-16:00,2.0,Mentatで弾塑性解析, ,11/6,10:00-12:00 13:30-17:00,5.5,Mentatで弾塑性解析・グラフの作成,尾山 ,11/10,10:00-12:00 13:00-17:00,6.0,Mentatで弾塑性解析 降伏後の傾きを調べる グラフ作成,尾山 ,11/11,10:00-12:00 13:30-16:00,4.5,Mentatで弾塑性解析 降伏後の傾きを調べる グラフの作成,尾山 ,11/12,10:30-12:30 13:30-15:30,4.0,Mentatで弾塑性解析 降伏後の傾きを調べる グラフの作成, ,11/13,10:30-12:30 14:00-18:00,6.0,Mentatで弾塑性解析 降伏後の傾きを調べる グラフの作成 後藤さんと今後の方針の確認,尾山 ,11/17,13:30-18:00,4.5,Mentatで摩擦係数を与えたときの水平変位の再確認 Mentatでの応力の表示方法を調べる,尾山 ,11/18,10:30-12:00 13:30-18:00,6.0,Mentatで摩擦係数を与えたときの水平変位の再確認 Mentatで応力の解析,尾山 ,11/19,10:00-12:00 14:00-22:30,10.5,Mentatで摩擦係数を与えたときの水平変位の再確認 Mentatで応力の解析 応力と変位の関係のグラフの作成,尾山・斉藤 ,11/20,10:30-12:30 14:30-18:30,6.0,Mentatで応力の解析 応力と変位の関係のグラフの作成,尾山 ,11/25,13:30-19:00,6.5,挟んだモデルの作成,尾山 ,11/26,9:30-12:00 13:30-22:00,12.0,挟んだモデルの解析と接触解析,尾山 ,11/27,10:30-12:00 13:30-19:00,7.0,挟んだモデルの解析 グラフの作成,尾山 ,11/28,10:30-20:00,9.5,挟んだモデルの解析 接触解析 グラフの作成,尾山 ,12/1,10:00-12:00 13:30-20:00,8.5,挟んだモデルの解析,尾山・斉藤 ,12/2,10:30-12:30 13:30-19:00,7.5,挟んだモデルの解析,尾山 ,12/3,10:30-23:00,12.5,挟んだモデルの解析 条件を変えての弾塑性解析 坪井の実験の手伝い,尾山・斉藤・坪井 ,12/4,11:00-23:00,12.0,弾塑性解析,尾山 ,12/5,14:30-19:00,4.5,弾塑性解析,尾山 ,12/8,10:30-19:00,8.5,弾塑性解析 グラフ作成,尾山 ,12/9,10:30-11:30,13.0,弾塑性解析 グラフ作成 挟んだモデルの解析,尾山 ,12/10,11:00-21:00,10.0,挟んだモデルの解析 後藤さんに質問,尾山 ,12/11,10:00-20:00,10.0,挟まないでの解析 後藤さんに質問 土木学会・支部発表の申込,尾山・斉藤・河原 ,12/12,10:00-20:30,10.5,挟んだモデルの解析 支部発表の手続き,尾山・斉藤・河原・大竹 ,12/15,10:00-19:00,9.0,挟んだモデルの解析 スライド作成,尾山 ,12/16,10:00-20:30,10.5,荷重をかける位置を変えたときの弾塑性解析 支部発表の手続き 後藤さんに質問,尾山・河原・大竹 ,12/17,11:00-20:00,9.0,Mises応力、相当応力の解析 後藤さんに質問,尾山 ,12/18,11:00-25:00,14.0,Mises応力の解析 相当応力の解析 スライド作成,尾山 ,12/20,14:00-17:00,3.0,スライド作成, ,1/7,11:00-15:30,4.5,硬い材料を挟んだ場合の解析, ,1/9,14:30-16:30,2.0,木材の引張試験,尾山 ,1/13,13:00-19:30,6.5,概要作成,尾山 ,1/14,11:00-18:30,7.5,概要作成,尾山 ,1/15,13:00-24:00,11.0,ヤング率を小さくした場合の解析 概要作成 スライド作成,尾山 ,1/19,14:00-21:30,7.5,概要作成 接触解析 オンサイト木橋に降伏点を与えた材料を挟んだモデルの解析,尾山・斉藤 ,1/20,9:30-19:30,10.0,概要作成 接触解析 オンサイト木橋に降伏点を与えた材料を挟んだ場合の解析,尾山 ,1/21,10:00-20:00,10.0,概要作成 スライド作成,尾山・斉藤 ,1/22,10:00-23:30,13.0,概要作成 スライド作成 3D組の実験の手伝い,尾山・坪井 ,1/23,11:00-17:00,4.0,概要作成,尾山 ,1/26,13:00-20:00,7.0,概要作成 東北支部に概要提出 スライド作成,尾山 ,1/27,11:00-19:30,8.5,卒論本体作成 スライド作成 発表練習,尾山・後藤さん ,1/28,14:00-19:30,5.5,卒論本体作成,尾山 ,1/29,11:00-20:00,9.0,卒論本体作成,尾山 ,1/30,11:00-20:00,9.0,発表練習 卒論本体作成,尾山・坪井 ,2/1,14:00-22:00,8.0,発表練習, ,2/2,14:00-22:00,8.0,卒論本体作成 発表練習, ,合計作業時間(h) ,583.5 *東北支部発表 [#qb55ce14] -申込方法等 http://www.jsce.or.jp/branch/tohoku/info/giken/H26/index.html -東北学院大学多賀城キャンパスアクセス方法 http://www.tohoku-gakuin.ac.jp/access/tagajo.html -東北学院大学 http://www.tohoku-gakuin.ac.jp -発表申請書 下の例を見て書いてね --&link(happyousien.doc,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/happyousien.doc) -発表申請書(例) --&link(happyousien(sample).pdf,http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/happyousien(sample).pdf) *卒業論文発表会・東北支部技術発表会メモ [#i949c303] -1月15日(木)9:30まで研究室集合 --支部発表の作戦を練る -1月23日(金)13:00まで 論文タイトルの提出 --提出内容 ---学籍番号・氏名・論文タイトル・指導教員名 --後藤さんに提出 -1月26日(月) 東北支部の原稿締切 --この日まで卒論発表の概要もでかす -1月第5週に発表練習 --スライドもちまちま作る -1月29日(木)13:00まで 論文概要の提出 --提出内容 ---片面印刷した原稿・pdf化したファイル --ファイル名 ---学籍番号「75○○○○○」の「75」を抜いた番号 ---おそらく11746になる --後藤さんに提出 -2月2日(月)9:00～15:00まで 発表用ファイルの提出 --工学資源学部1号館4F学科図書館(土木の図書館)に提出 --学科の図書館にあるパソコンに落とす ---デスクトップの中にある「卒業論文」フォルダに保存 --ファイル名 ---学籍番号「75○○○○○」の「75」を抜いた番号 ---おそらく11746になる -2月3日(火)～4日(水)9:00～17:00 修士論文・卒業論文発表会 --会場:総合研究棟2F大セミナー室 --発表時間 ---発表6分・質問4分 --卒論日誌を持参 ---事前に該当ページを印刷し、後藤さんからハンコをもらう -3月7日(月) 平成26年度土木学会東北支部技術研究発表会 --会場:東北学院大学工学部多賀城キャンパス *今後の方針 [#f7aa7d18] -Marc Mentatで接触解析 -木材の間にゴムみたいなものを挟めて接触解析の代わりになる方法の提案 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/masatsu.png ←こういう感じのもの *次のゼミまでやること(ゼミは毎週金曜日) [#vb0290a1] -JIS K 7125 http://kikakurui.com/k7/K7125-1999-01.html -Marc Mentatで弾塑性解析→できた -降伏後に傾きが与えられるかcheck(bilinearモデルに想定するような・・・)→できた -Mentatの弾塑性解析においてbi-linearモデルと完全弾塑性モデルに変えられる方法を調べる→変えられた -弾塑性材料の側面をローラーで固定、底面を固定した際、弾塑性材料は塑性するかcheck http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/sok.png -概要を書く *荷重と変位の関係 [#mfe6600e] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/JIS.png -ファイル名:jis2 -節点番号 --上:127 --下:1104 ,荷重(N),変位(m) ,0.1N,7.737714416E-13 ,0.2N,1.46830883E-12 ,0.3N,2.20246325E-12 ,0.4N,2.93661766E-12 ,0.5N,3.67077208E-12 ,0.6N,4.4049265E-12 ,0.7N,5.13908091E-12 ,0.8N,5.87323533E-12 ,0.9N,6.60738974E-12 ,1.0N,2.182469255E-10 ,1.1N,0.0850635 →1.1Nで動き始める http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/JISg.png *鉄の引張試験(とりあえず練習として弾塑性とかの条件を無視して・・・) [#lbafc18e] **とりあえず同じ分割数で荷重を100Nずつ増やした時の引張試験 [#peba59db] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/tetsu.png -長さ:1m -幅:10cm -高さ:10cm -測定結果 ,荷重(N),伸び(m),ひずみ(m),応力(N/m^2),ヤング率(GPa) ,100,4.797E-8,4.797E-8,10000,208.4 ,200,9.593E-8,9.593E-8,20000,208.4 ,300,1.439E-7,1.439E-7,30000,208.4 ,400,1.919E-7,1.919E-7,40000,208.4 ,500,2.398E-7,2.398E-7,50000,208.5 ,600,2.878E-7,2.878E-7,60000,208.4 ,700,3.358E-7,3.358E-7,70000,208.4 ,800,3.837E-7,3.837E-7,80000,208.5 ,900,4.317E-7,4.317E-7,90000,208.5 ,1000,4.797E-7,4.797E-7,100000,208.5 ,1100,5.276E-7,5.276E-7,110000,208.5 ,1200,5.756E-7,5.756E-7,120000,208.5 ,1300,6.236E-7,6.236E-7,130000,208.5 ,1400,6.715E-7,6.715E-7,140000,208.5 ,1500,7.195E-7,7.195E-7,150000,208.5 ,1600,7.675E-7,7.675E-7,160000,208.5 ,1700,8.154E-7,8.154E-7,170000,208.5 ,1800,8.634E-7,8.634E-7,180000,208.5 ,1900,9.114E-7,9.114E-7,190000,208.5 ,2000,9.593E-7,9.593E-7,200000,208.5 -荷重と変位の関係(縦軸:荷重 横軸:変位) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/tetsug.png -ヤング率:207.47GPa --データをとった20個の平均をとる -相対誤差:0.713% **次に荷重を変えずに分割数のみを変えたときの引張試験 [#rcc9dc7b] -荷重:100N -応力:10000N/m^2 -ヤング率(理論値):206GPa ,分割数,伸び(m),ひずみ(m),ヤング率(GPa) ,8,4.810E-8,4.810E-8,207.9 ,16,4.828E-8,4.828E-8,207.1 ,20,4.831E-8,4.831E-8,206.9 ,22,4.832E-8,4.832E-8,206.9 -無次元化したグラフ(縦軸:無次元化した値 横軸:分割数) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/tetsu2.png -相対誤差:0.437% *鉄の引張試験(弾塑性解析) [#cb99257e] **とりあえず1MN～3MNまで1MNずつ増やし、3MN～4MNまで0.1MNずつ増やしたときの弾塑性解析 [#ecf978c1] -ファイル名:danso -鉄の種類:SS400 -降伏応力:400N/mm^2 -鉄の寸法 --長さ:1m --幅:10cm --高さ:10cm -断面積:100cm^2=10000mm^2 -降伏時の荷重:4MN -荷重は面荷重(FACE LOADを使う) ,荷重(MN),伸び(m),ひずみ(m),ヤング率(GPa),相対誤差(%) ,1.0,4.836E-4,4.836E-4,206.7,0.340 ,2.0,9.671E-4,9.671E-4,206.8,0.388 ,3.0,1.451E-3,1.451E-3,206.7,0.340 ,3.5,1.692E-3,1.692E-3,206.8,0.388 ,3.6,1.741E-3,1.741E-3,206.7,0.340 ,3.7,1.789E-3,1.789E-3,206.8,0.388 ,3.8,1.838E-3,1.838E-3,206.7,0.340 ,3.9,1.897E-3,1.897E-3,205.5,0.243 ,4.0,3.658E-3,3.658E-3 ※ここから線形の関係が成立しない,109.3,46.94 -応力とひずみの関係(縦軸:応力 横軸:ひずみ) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/danso1.png →これだと、bilinearモデルに想定するのか・・・ **次に3.9MN～4.0MNまで0.01MNずつ増やしたときの弾塑性解析 [#pc817f17] -ファイル名:danso -鉄の種類:SS400 -降伏応力:400N/mm^2 -鉄の寸法等 --長さ:1m --幅:10cm --高さ:10cm -断面積:100cm^2=10000mm^2 -降伏時の荷重:4MN -荷重は面荷重(FACE LOADを使う) ,荷重(MN),伸び(m),ひずみ(m),ヤング率(GPa),相対誤差(%) ,1.0,4.836E-4,4.836E-4,206.7,0.340 ,2.0,9.671E-4,9.671E-4,206.8,0.388 ,3.0,1.451E-3,1.451E-3,206.7,0.340 ,3.5,1.692E-3,1.692E-3,206.8,0.388 ,3.6,1.741E-3,1.741E-3,206.7,0.340 ,3.7,1.789E-3,1.789E-3,206.8,0.388 ,3.8,1.838E-3,1.838E-3,206.7,0.340 ,3.9,1.897E-3,1.897E-3,205.5,0.243 ,3.91,1.906E-3,1.906E-3,205.1,0.436 ,3.92,1.900E-3,1.900E-3,206.3,0.146 ,3.93,1.906E-3,1.906E-3,206.1,0.0485 ,3.94,1.914E-3,1.914E-3,205.6,0.194 ,3.95,1.923E-3,1.923E-3,205.4,0.291 ,3.96,1.915E-3,1.915E-3,206.8,0.388 ,3.97,1.949E-3,1.949E-3,203.6,1.16 ,3.98,1.954E-3,1.954E-3,203.7,1.12 ,3.99,1.989E-3,1.989E-3,200.6,2.62 ,3.99999,3.623E-3,3.623E-3,110.4 ,4.0,3.658E-3,3.658E-3 ※ここから線形の関係が成立しない,109.3,46.94 -応力とひずみの関係(縦軸:応力 横軸:ひずみ) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/danso3.png →降伏点近くはほとんど変化しない **4.0MN以上の荷重をかけた時 [#c5a8a0ee] -ファイル名:danso -鉄の種類:SS400 -降伏応力:400N/mm^2 -鉄の寸法等 --長さ:1m --幅:10cm --高さ:10cm -断面積:100cm^2=10000mm^2 -降伏時の荷重:4MN -荷重は面荷重(FACE LOADを使う) ,荷重(MN),伸び(m),ひずみ(m) ,4.0,3.658E-3,3.658E-3 ,4.1,7.122E-1,7.122E-1 ,4.2,1.423,1.423 ,4.3,2.138,2.138 ,4.4,2.853,2.853 ,4.5,3.569,3.569 -荷重と変位の関係(縦軸:荷重 横軸:変位) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/kajuheni.png -応力とひずみの関係(縦軸:応力 横軸:ひずみ) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/danso4.png →降伏後は比例関係になっている *引張降伏応力とせん断降伏応力 [#qc37093a] -引張降伏応力とせん断降伏応力の関係は$\sigma_{Y}=\sqrt{3}\tau_{Y}$となる。 -摩擦試験で作ったモデルを使って考える --ファイル名:jis2 --節点番号 ---上:127 ---下:1104 **鋼と鋼の場合 [#w4e9ddfe] -滑らせた材料:鋼材 -台:鋼材 ,荷重(N),変位(m),せん断応力$\tau_{xz}$,ひずみ ,0.1N,5.016E-13,6.076E-6,2.960E-11 ,0.2N,1.003E-12,1.214E-5,5.893E-11 ,0.3N,1.505E-12,1.821E-5,8.840E-11 ,0.4N,1.999E-12,2.428E-5,1.179E-10 ,0.5N,2.508E-12,3.035E-5,1.473E-10 ,0.6N,3.279E-12,3.453E-5,1.676E-10 ,0.7N,4.872E-12,4.175E-5,2.026E-10 ,0.8N,6.762E-12,4.961E-5,2.408E-10 ,0.9N,8.932E-12,5.804E-5,2.817E-10 ,1.0N,1.240E-10,6.198E-5,3.009E-10 ,1.1N,4.074E-10,1.083E-4,5.257E-10 ,1.2N,6.143E-10,1.302E-4,6.320E-10 ,1.3N,5.000E-9,8.791E-5,4E-9 ,1.31N,5.000E-9,8.742E-5,4E-9 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/jis4.png **木と鋼の場合 [#zc34142e] -滑らせた材料:木材 -台:鋼材 ,荷重(N),変位(m),せん断応力$\tau_{xz}$,ひずみ ,0.1N,7.738E-13,6.093E-6,1.105E-6 ,0.2N,1.468E-12,1.219E-5,2.031E-6 ,0.3N,2.202E-12,1.828E-5,3.046E-6 ,0.4N,2.937E-12,2.437E-5,4.061E-6 ,0.5N,3.671E-12,3.047E-5,5.078E-6 ,0.6N,4.405E-12,3.656E-5,6.093E-6 ,0.7N,5.139E-12,4.265E-5,7.108E-6 ,0.8N,5.873E-12,4.875E-5,8.125E-6 ,0.9N,6.607E-12,5.484E-5,9.14E-6 ,1.0N,2.182E-10,6.192E-5,1.032E-5 ,1.1N,0.0851,5.716E-5,9.526E-6 ,応力(MPa),変位(m),鋼材と木材のズレ(m),せん断応力σzx(MPa),Von Mises応力,相当塑性ひずみ, ,0.5,1.577E-08,0,0.256,0.954043,0 ,0.6,1.893E-08,0,0.307,0.9958988,0 ,0.7,2.208E-08,0,0.358,1.04355,0 ,0.8,2.523E-08,0,0.409,1.09624,0 ,0.9,2.839E-08,0,0.461,1.15337,0 ,1.0,3.074E-06,3.044E-06,0.501,1.19659,0 ,1.1,1.573E-05,1.570E-05,0.519,1.22852,0 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/jis5.png **鋼と木の場合 [#f723a5ef] -滑らせた材料:鋼材 -台:木材 ,荷重(N),変位(m),せん断応力$\tau_{xz}$,ひずみ ,0.1N,1.205E-10,7.944E-6,3.8E-8 ,0.2N,2.409E-10,1.580E-5,7.6E-8 ,0.3N,3.614E-10,2.383E-5,1.15E-7 ,0.4N,4.819E-10,3.178E-5,1.54E-7 ,0.5N,6.023E-10,3.972E-5,1.92E-7 ,0.6N,7.228E-10,4.767E-5,2.31E-7 ,0.7N,8.506E-10,5.434E-5,2.63E-7 ,0.8N,1.0E-9,6.260E-5,3.03E-7 ,0.9N,1.0E-9,7.051E-5,3.42E-7 ,1.0N,9.0E-9,8.418E-5,4.08E-7 ,1.01N,1.5E-8,9.596E-5,4.65E-7 ,1.02N,1.5E-8,8.578E-5,4.16E-7 ,1.03N,1.8E-8,8.872E-5,4.3E-7 ,1.04N,2.3E-8,9.108E-5,4.42E-7 ,1.05N,2.8E-8,9.332E-5,4.53E-7 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/mo.png *間に摩擦係数の代わりになる材料を挟んだ場合の変位 [#te033b5f] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/model1.png -ファイル名:hasamu3 -上に載せている材料:木材 -鋼材と木材の間に挟んだ材料:木材 **分割数128の場合 [#j10b4d22] ,荷重(N),節点番号30762での変位(m) ,0.1,3.9814E-13, ,0.2,7.96281E-13 ,0.3,1.19442E-12 ,0.4,1.59256E-12 ,0.5,1.99221E-12 ,0.6,2.39171E-12 ,0.7,2.79377E-12 ,0.8,3.2E-12 →0.9Nあたりから解析ができなくなった。 -荷重と変位の関係(縦軸:荷重 横軸:変位) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/128.png →比例・・・ **分割数16の場合 [#x8cbe2e3] 分割数128では途中から解析が不可能になったので分割数を少なくしてみた。 ,荷重,節点番号1184での変位(m) ,0,0 ,0.1,3.93612E-13 ,0.2,7.87224E-13 ,0.3,1.18084E-12 ,0.4,1.57472E-12 ,0.5,1.9715E-12 ,0.6,2.37079E-12 ,0.7,2.77595E-12 ,0.8,3.21053E-12 ,0.9,3.72077E-12 ,1.0,3.85327E-12 ,1.1,3.86749E-12 ,1.2,3.93149E-12 ,1.3,3.924E-12 →1.2Nが限界 -荷重と変位の関係(縦軸：荷重 横軸：変位) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/as.png **分割数8の場合 [#c0033042] 分割数16でも1.2Nで限界だった。さらに分割数を少なくしてみる。 ,荷重(N),節点番号624での変位(m) ,0.1,3.541E-13 ,0.2,7.081E-13 ,0.3,1.062E-12 ,0.4,1.416E-12 ,0.5,1.770E-12 ,0.6,2.152E-12 ,0.7,2.532E-12 ,0.8,2.933E-12 ,0.9,3.132E-12 ,1.0,3.603E-12 ,1.1,3.999E-12 ,1.2,4.378E-12 ,1.3,4.769E-12 ,1.4,5.153E-12 ,1.5,5.533E-12 ,1.6,5.953E-12 ,1.7,6.336E-12 ,1.8,6.720E-12 ,1.9,7.101E-12 ,2.0,7.482E-12 ,3.0,1.127E-11 ,4.0,1.504E-11 ,5.0,1.881E-11 →動いた -荷重と変位の関係(縦軸:荷重 横軸:変んだ場合の変位) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/8.png →やっぱり比例関係・・・ **分割数8の場合 [#ad007736] 鋼材を挟んでみると・・・ -降伏応力:1.00528E-4 ,荷重(N),節点番号624での変位(m) ,0.1,3.273E-13 ,0.2,6.546E-13 ,0.3,9.807E-13 ,0.4,1.294E-12 ,0.5,1.595E-12 ,0.6,1.976E-12 -荷重と変位の関係(縦軸:荷重 横軸:変位) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/8te.png *接触解析の場合 [#pa9fe17f] **分割数128の場合 [#j8e0c0c7] -解析モデル http://www.str.ce.akita-u.ac.jp/~gotouhan/oyama/masatu128.png --上からの荷重:2N -台(鋼材)の寸法等 --長さ:60cm --幅:24cm --高さ:1cm --ヤング率:206GPa --ポアソン比:0.3 -試験体(木材)の寸法 --長さ:20cm --幅:8cm --高さ:10cm --${E_{x}}$=${6GPa}$,${E_{y}}$=${\frac{E_{x}}{25}}$,${E_{z}}$=${\frac{E_{x}}{25}}$ --${n_{xy}}$=${0.4}$,${n_{yz}}$=${0.016}$,${n_{zx}}$=${0.016}$ --${G_{xy}}$=${400}$,${G_{yz}}$=${400}$,${G_{zx}}$=${400}$ ,荷重(N),変位(m),鋼材と木材のズレ(m),せん断応力${\sigma_{zx}}$(MPa) ,0.5,1.971E-12,0,3.199E-05 ,0.6,2.366E-12,0,3.839E-05 ,0.7,2.760E-12,0,4.479E-05 ,0.8,3.154E-12,0,5.118E-05 ,0.9,3.549E-12,0,5.759E-05 ,1.0,2.484E-10,2.446E-10,6.328E-05 ,1.1,1.717E-09,1.713E-09,6.571E-05 -せん断応力と変位の関係(縦軸：せん断応力${\sigma_{zx}}$ 横軸：変位(m)) http://www.str.ce.akita-u.ac.jp/~gotouhan/oyama/jis128.png **分割数16の場合 [#hbb0d618] ,荷重(N),変位,せん断応力$\sigma_{zx}$(MPa) ,0.1,0,6.357E-6 ,0.2,0,1.271E-5 ,0.3,0,1.907E-5 ,0.4,0,2.543E-5 ,0.5,0,3.181E-5 ,0.6,0,3.818E-5 ,0.7,0,4.445E-5 ,0.8,0,5.092E-5 ,0.9,0,5.730E-5 ,0.91,0,5.793E-5 ,0.92,0,5.856E-5 ,0.93,0,5.920E-5 ,0.94,0,5.984E-5 ,0.95,0,6.048E-5 ,0.96,0,6.115E-5 ,0.97,9.495152E-11,6.171E-5 ,0.98,3.3654032E-10,6.166E-5 ,0.99,1.41E-7,4.805E-5 ,1.0,2.39E-7,4.450E-5 ,1.1,0.00357797,0.000125 -荷重と変位の関係(縦軸：荷重 横軸：変位) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/aa.png -せん断応力と変位の関係(縦軸：せん断応力${\sigma_{zx}}$ 横軸：変位(m)) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/aaa.png *摩擦試験モデルと挟んだモデルについて [#u3f03ed1] -のっけている材料を基準に考える -摩擦試験モデル→masatsusikentenp -挟んだモデル→hasamutemp -tempはMeshを切っていない -分割数 --x方向:任意 --y方向:6 --z方向:8 --挟んだモデルの挟んでいる物体のz方向の分割数は2 *弾塑性解析(POWER LOWにした時) [#qae98445] -COEFFICIENT B=2の時 ,応力(MPa),変位(m),ひずみ,ヤング率(GPa) ,100,4.836E-4,4.836E-4,206.8 ,200,9.671E-4,9.671E-4,206.8 ,235,1.151E-4,1.151E-4,204.1 ,250,9.251E-1,9.251E-1, ,300,4.502,4.502, ,350,8.079,8.079, ,400,11.66,11.66, -COEFFICIENT B=5の時 ,応力(MPa),変位(m),ひずみ,ヤング率(GPa) ,100,4.836E-4,4.836E-4,206.8 ,200,9.671E-4,9.671E-4,206.8 ,235,1.138E-3,1.138E-3,206.5 ,250,7.106E-1,7.106E-1, ,300,4.228,4.228, ,350,7.864,7.864, ,400,11.44,11.44, -COEFFICIENT B=10の時 ,応力(MPa),変位(m),ひずみ,ヤング率(GPa) ,100,4.836E-4,4.836E-4,206.8 ,200,9.671E-4,9.671E-4,206.8 ,235,1.137E-3,1.137E-3,206.6 ,250,7.106E-1,7.106E-1, ,300,4.228,4.228, ,350,7.864,7.864, ,400,11.44,11.44, -COEFFICIENT B=30の時 ,応力(MPa),変位(m),ひずみ,ヤング率(GPa) ,100,4.836E-4,4.836E-4,206.8 ,200,9.671E-4,9.671E-4,206.8 ,235,1.137E-3,1.137E-3,206.6 ,250,7.106E-1,7.106E-1, ,300,2.499,2.449, ,350,7.864,7.864, ,400,11.44,11.44, -COEFFICIENT B=30の時 ,応力(MPa),変位(m),ひずみ,ヤング率(GPa) ,100,4.836E-4,4.836E-4,206.8 ,200,9.671E-4,9.671E-4,206.8 ,235,1.157E-3,1.157E-3,206.6 ,250,7.106E-1,7.106E-1, ,300,2.499,2.449, ,350,7.864,7.864, ,400,11.44,11.44, -荷重と変位の関係(縦軸:荷重 横軸:変位) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/aaaa.png *TABLE機能を用いた弾塑性解析 [#t67b5520] -MATERIAL PROPATYの中にあるTABLE機能を用いて、応力と塑性ひずみを入力すれば、完全弾塑性体になる...? -傾きは自分で与える -与えた傾き、応力から塑性ひずみを逆算し、入力 -弾性域の傾きは206000 **傾き2とした場合 [#wacf35d1] -与えた塑性ひずみ ,応力(MPa),塑性ひずみ ,235,0 ,235.1,5.05705E-2 ,235.2,1.005705E-1 ,235.3,1.505705E-1 ,235.4,2.005705E-1 ,235.5,2.505705E-1 -解析結果 ,応力(MPa),変位(m),ひずみ ,100,4.836E-4,4.836E-4 ,200,9.671E-4,9.671E-4 ,235,2.083E-3,2.083E-3 ,235.1,6.229E-3,6.229E-3 ,235.2,1.258E-2,1.258E-2 ,235.3,1.856E-2,1.856E-2 ,235.4,2.433E-2,2.433E-2 ,235.5,3.015E-2,3.015E-2 -荷重と変位の関係(縦軸:荷重 横軸:ひずみ) 傾き:17.81451527 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/235-4.png **傾き20とした場合 [#rad990f4] -与えた塑性ひずみ ,応力(MPa),塑性ひずみ ,235,0 ,235.1,5.05705E-3 ,235.2,1.005705E-2 ,235.3,1.505705E-2 ,235.4,2.005705E-2 ,235.5,2.505705E-2 -解析結果 ,応力(MPa),変位(m),ひずみ ,100,4.836E-4,4.836E-4 ,200,9.671E-4,9.671E-4 ,235,1.316E-3,1.775E-3 ,235.1,3.829E-3,3.829E-3 ,235.2,6.453E-3,6.453E-3 ,235.3,8.930E-3,8.930E-3 ,235.4,1.143E-2,1.143E-2 ,235.5,1.394E-2,1.394E-2 -荷重と変位の関係(縦軸:荷重 横軸:変位) 傾き;39.57574798 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/235-3.png **傾き200とした場合 [#t4821ec6] -与えた塑性ひずみ ,応力(MPa),塑性ひずみ ,235,0 ,235.1,5.05705E-4 ,235.2,1.005705E-3 ,235.3,1.505705E-3 ,235.4,2.005705E-3 ,235.5,2.505705E-3 -解析結果 ,応力(MPa),変位(m),ひずみ ,100,4.836E-4,4.836E-4 ,200,9.671E-4,9.671E-4 ,235,1.316E-3,1.316E-3 ,235.1,1.738E-3,1.738E-3 ,235.2,2.119E-3,2.119E-3 ,235.3,2.508E-3,2.508E-3 ,235.4,2.902E-3,2.902E-3 ,235.5,3.298E-3,3.298E-3 -応力と変位の関係(縦軸:応力 横軸:変位) 傾き:252.2704339 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/235.png **傾き2000とした場合 [#zcef5465] -与えた塑性ひずみ ,応力(MPa),塑性ひずみ ,235,0 ,235.1,5.05705E-5 ,235.2,1.005705E-4 ,235.3,1.505705E-4 ,235.4,2.005705E-4 ,235.5,2.505705E-4 -解析結果 ,応力(MPa),変位(m),ひずみ ,100,4.836E-4,4.836E-4 ,200,9.671E-4,9.871E-4 ,235,1.183E-3,1.183E-3 ,235.1,1.225E-3,1.225E-3 ,235.2,1.270E-3,1.270E-3 ,235.3,1.315E-3,1.315E-3 ,235.4,1.361E-3,1.361E-3 ,235.5,1.406E-3,1.406E-3 -応力とひずみの関係(縦軸:応力 横軸:ひずみ) 傾き:2242.152466 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/235-1.png **傾き20000とした場合 [#v6717f0a] -与えた塑性ひずみ ,応力(MPa),塑性ひずみ ,235,0 ,235.1,5.05705E-6 ,235.2,1.005705E-5 ,235.3,1.505705E-5 ,235.4,2.005705E-5 ,235.5,2.505705E-5 -解析結果 ,応力(MPa),変位(m),ひずみ ,100,4.836E-4,4.836E-4 ,200,9.671E-4,9.671E-4 ,235,1.151E-3,1.151E-3 ,235.1,1.150E-3,1.150E-3 ,235.2,1.155E-3,1.155E-3 ,235.3,1.160E-3,1.160E-3 ,235.4,1.170E-3,1.165E-3 ,235.5,1.170E-3,1.170E-3 -応力とひずみの関係(縦軸:応力 横軸:ひずみ) 傾き:26315.78947 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/235-2.png **まとめ [#t8f30010] -傾きを決め、そこから逆算し塑性ひずみを求め、与えたことにより、各々の塑性ひずみは変わる。 -塑性ひずみが変わったことにより、傾きは変わる。 -応力とひずみの関係(縦軸:応力 横軸:ひずみ) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/235ga.png -上のグラフの拡大版(縦軸:応力 横軸:ひずみ) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/235ga1.png *弾塑性材料を挟んだ場合の水平変位 [#c6af7b88] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/masamugg.png -摩擦試験の条件で台と材料の間に挟んだ場合の水平変位 -摩擦係数に相当する弾塑性材料を挟む -変位を見る場所は何点かとって接触解析の解析結果と一番近い箇所を取ることとする *弾塑性材料(挟んだ材料) [#gbc7b530] -ポアソン比:0.3 -235MPaで降伏 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/hasamugg1.png -塑性して動いた！ -進む面に荷重をかけ過ぎると、傾いたりするみたい... -荷重は変えずに、左右に配分すると傾かなくなる。 --1:4とか1:5みたいにかければいい感じ。 --できるだけ大きい荷重は後ろ(押す方向)にかけてあげれば良い。 **ヤング率206MPaの場合 [#tdac7216] -Plasticityに降伏応力を入れた場合(TABLE機能は使用していない) ,応力(MPa),先端の点(1)での変位(m),真ん中の点(2204)での変位(m),後ろの点(6236)での変位(m) ,235,3.11384E-4,3.00343E-4,2.89997E-4 ,240,3.18009E-4,3.06733E-4,2.96139E-4 ,250,3.31259E-4,3.19513E-4,3.08478E-4 ,260,3.4451E-4,3.32294E-4,3.20817E-4 ,270,5.19406E-3,5.23197E-3,5.27268E-3 -上と同じ条件でポアソン比0.4とした場合 ,応力(MPa),先端の点(1)での変位(m),真ん中の点(2204)での変位(m),後ろの点(6236)での変位(m) ,235,3.34033E-4,3.22853E-4,3.12506E-4 ,240,3.4114E-4,3.29722E-4,3.19155E-4 ,250,3.55354E-4,3.4346E-4,3.32453E-4 ,260,3.70568E-4,3.57199E-4,3.45751E-4 -Plasticityに降伏応力を入れない場合(弾性?としての解析) ,応力(MPa),先端の点(1)での変位(m),真ん中の点(2204)での変位(m),後ろの点(6236)での変位(m) ,235,3.11384E-4,3.00343E-4,2.8997E-4 ,240,3.18009E-4,3.06733E-4,2.96139E-4 ,250,3.331259E-4,3.19513E-4,3.08478E-4 ,260,3.4451E-4,3.32294E-4,3.20817E-4 ,270,3.5776E-4,3.45075E-4,3.33156E-4 **ヤング率20MPaにした時 [#nd9367c4] -Plasticityに降伏応力を入れた場合(TABLE機能は使用していない) ,応力(MPa),先端の点(1)での変位(m),真ん中の点(2204)での変位(m),後ろの点(6236)での変位(m) ,235,3.07785E-3,3.06637E-3,3.05518E-3 ,240,3.14333E-3,3.13161E-3,3.12019E-3 ,250,3.27431E-3,3.26209E-3,3.25019E-3 ,260,3.40528E-3,3.39258E-3,3.3802E-3 ,270,3.53625E-3,3.52306E-3,3.51021E-3 ,280,3.66722E-3,3.65354E-3,3.64022E-3 ,290,3.79819E-3,3.78403E-3,3.777022E-3 ,300,3.92917E-3,3.91451E-3,3.90023E-3 *挟まないモデル [#jce5f07d] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/hasamanai.png -接合部分の点は変形しない(動かない) -荷重を加える場所によって変形は変わる -後はヤング率をいじって理想の変形を探す... *横方向に一様に荷重をかけた場合 その1 [#gaa52195] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/sendan1.png -降伏応力:100MPa -長さ:200mm=0.2m -幅:80mm=0.08m -高さ:1mm=0.001m -相当応力$\bar{\sigma}$を求める式 $\bar{\sigma}=\sqrt{\frac{3}{2}}\sqrt{S_{xx}^2+S_{yy}^2+S_{zz}^2+2(S_{xy}^2+S_{yz}^2+S_{zx}^2)}$ $S_{xx}=\sigma_{xx}-\sigma_{m}$ $S_{yy}=\sigma_{yy}-\sigma_{m}$ $S_{zz}=\sigma_{zz}-\sigma_{m}$ $S_{xy}=\sigma_{xy}$ $S_{yz}=\sigma_{yz}$ $S_{zx}=\sigma_{zx}$ $\sigma_{m}=\frac{1}{3}(\sigma_{xx}+\sigma_{yy}+\sigma_{zz})$ **全面に荷重をかけた場合 [#o86cce0c] ***真ん中の点(0.1m,0.04m,0.001m) [#cd1556fa] ,荷重(N),$\sigma_{Mises}$,相当応力,$\sigma_{xx}$,$\sigma_{yy}$,$\sigma_{zz}$,$\tau_{xy}$,$\tau_{yz}$,$\tau_{zx}$,変位(m) ,1000,0.0860258,0.086016067,0,0,0,0,1.34224E-18,0.04966140,6.29E-10 ,5000,0.430129,0.430080339,0,5.18696E-15,8.10463E-17,0,0,0.248307,3.14584E-9 ,10000,0.860258,0.860160679,0,1.037E-14,1.620E-16,0,3.13985E-18,0.496614,6.29167E-9 ,15000,1.29039,1.29024102,0,0,1.29674E-15,0,2.148E-17,0.744921,9.4375E-9 ,20000,1.72052,1.720321359,0,0,0,0,1.514E-17,0.993228,1.25833E-8 ,25000,2.15064,2.150393039,8.299E-14,0,6.4837E-16,0,4.29518E-17,1.24153,1.573E-8 ,26000,2.23667,2.236424003,0,0,0,0,2.269E-17,1.2912,1.63583E-8 ,29000,2.49475,2.494464932,0,0,0,0,2.269E-17,1.44018,1.82458E-8 ,30000,2.58077,2.580478575,0,0,0,0,8.347E-17,1.48984,1.8875E-8 -荷重と変位の関係(縦軸:変位 横軸:荷重) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/so.png ***端の点(0.2m,0m,0.001m) [#f5d42cd7] -応力が集中している ,荷重(N),$\sigma_{Mises}$,相当応力$\bar{\sigma}$,$\sigma_{xx}$,$\sigma_{yy}$,$\sigma_{zz}$,$\tau_{xy}$,$\tau_{yz}$,$\tau_{zx}$,変位(m) ,1000,0.701445,0.7134096016,0.544447,0.085576,-0.079059,-0.0129,0.00276,0.255,5.232-9 ,10000,7.01445,7.12810292,5.44447,0.855756,-0.79059,-0.129121,0.0275558,2.54534,5.232E-8 ,15000,10.5217,10.692151098,8.1667,1.28363,-1.18589,-0.193681,0.0413337,3.818,7.84864E-8 ,20000,12.2935,12.540298909,9.96516,2.07722,-1.06716,-0.200825,0.0464743,4.48041,1.08583E-7 ,25000,11.1995,11.42868098,8.53547,3.6565,-0.334745, -0.034878,0.0211037,4.8786,1.75359E-7 ,30000,10,10.218137513,4.51615,3.26422,-1.53686,-0.118612,-0.00761243,4.95781,4.65926E-7 **端に荷重をかけない場合 [#l228d234] ***応力が大きい点(8082) [#v6ef6038] ,荷重(N),$\sigma_{Mises}$,相当応力$\bar{\sigma}$,$\sigma_{xx}$,$\sigma_{yy}$,$\sigma_{zz}$,$\tau_{xy}$,$\tau_{yz}$,$\tau_{zx}$,変位(m) ,1000,0.329444,0.43778551,-0.350768,-0.0535491,0.028944,0.00173891,-0.000481659,0.0225921,-6.05511E-10 ,5000,1.64722,2.188877546,-1.75384,-0.267745,0.14472,0.00169455,-0.0024083,0.11296,-3.02755E-9 ,10000,3.29444,4.277855102,-3.50768,-0.535491,0.28944,0.0173891,-0.0048165,0.225921,-6.05511E-9 ,15000,4.94166,6.566782584,-5.26152,-0.803236,0.434161,0.0260837,-0.00722489,0.338881,-9.08266E-9 ,20000,6.58888,8.755709782,-7.01536,-1.07098,0.578881,0.0347782,-0.0096331,0.451841,-1.21102E-8 ,25000,8.23609,10.94463824,-8.7692,-1.33873,0.723601,0.0434728,-0.0120415,0.564802,-1.51378E-8 ,30000,9.88331,13.13351663,-10.523,-1.60647,0.868321,0.0521673,-0.0144498,0.677762,-1.81653E-8 ,35000,12.5305,15.32143528,-12.276,-1.87422,1.01304,0.0608619,-0.0168581,0.790722,-2.11929E-8 ***真ん中の点 [#d74617a8] ,荷重(N),$\sigma_{Mises}$,相当応力$\bar{\sigma}$,$\sigma_{xx}$,$\sigma_{yy}$,$\sigma_{zz}$,$\tau_{xy}$,$\tau_{yz}$,$\tau_{zx}$,変位(m) ,1000,0.139547,0.139524532,0,0,0,0,0,0.0805603,1.0134E-9 ,5000,0.697734,0.697673529,0,0,0,0,0,0.402802,5.06701E-9 ,10000,1.39547,1.395345327,0,0,0,0,0,0.805603,1.0134E-8 ,15000,2.0932,2.093010196,0,0,0,0,0,1.2084,1.5201E-8 ,20000,2.79094,2.790697582,0,0,0,0,0,1.61121,2.0268E-8 ,25000,3.48867,3.488367647,0,0,0,0,0,2.01401,2.5335E-8 ,30000,4.1864,4.186037712,0,0,0,0,0,2.41681,3.0402E-8 ,35000,4.88414,4.883707778,0,0,0,0,0,2.81961,3.54691E-8 ,50000,6.97734,6.976735294,0,0,0,0,0,4.02802,5.06701E-8 ,60000,8.37973,8.378986307,0,0,0,0,0,4.83761,6.08377E-8 *横方向に一様に荷重をかけた場合 その2 [#n6df3660] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/sendan.png -せん断によって変形を起こすかをcheck -$\sigma=\sqrt{3}\tau$の関係からせん断応力を出す -せん断応力に断面積をかけて荷重とする -Point Loadで解析する ,応力(MPa),荷重(MN),変位(m) ,100,5.77350269,1.323E-4 ,150,8.660254038,1.984E-4 ,160,9.237604307,2.210E-4 ,170,9.814954576,2.494E-4 ,180,10.39230485,2.943E-4 ,190,10.9695511,3.678E-4 ,200,11.54700538,9.527E-4 -応力と変位の関係(縦軸:応力 横軸:変位) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/sita1.png *側面に硬い材料を挟めた場合 [#k0787238] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/soku.png -弾塑性材料 --ヤング率:206GPa --降伏応力:235MPa -側面の硬い材料 --ヤング率:2060GPa -ファイル名:soku ,応力(MPa),変位(m) ,100,5.782E-5 ,200,1.157E-4 ,235,1.317E-4 ,300,2.303E-4 ,500,2.891E-4 ,1000,5.782E-4 -1000MPaを超え始めると、底面は固定されているから動かない。 -上面だけが動く(下図は2000MPaの時の挙動) http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/2000.png *等方性でヤング率を木材のヤング率にした場合の引張試験 [#h4f53868] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/moy.png -ヤング率:6GPa -降伏応力:10MPa ,応力(MPa),変位(m),ヤング率(GPa) ,1,1.661E-4,6.02 ,5,8.307E-4,6.02 ,9,1.496E-3,6.02 ,10,4.711E-3,2.12 ,10.01,3.367E-1,0.0297 ,10.02,7.141E-1,0.0140 ,10.03,1.088,0.00922 *間に載せている材料と同じヤング率の材料を挟んだ場合 [#pf588fd1] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/moy2.png -ヤング率:6GPa -挟んだ材料の降伏応力:10MPa ,応力(MPa),変位(m) ,1,2.968E-5 ,5,1.495E-4 ,6,1.804E-4 ,7,2.116E-4 ,8,2.452E-4 ,8.1,2.487E-4 ,8.2,2.523E-4 ,8.3,2.559E-4 ,8.4,2.595E-4 ,8.5,2.630E-4 ,8.6,2.665E-4 *接触解析 [#zc94cc55] http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/se.png **木と鋼材の接触解析 [#e4dde26e] -上の材料:木材 -下の材料:鋼材 -上の番号:2204 -下の番号:14142 -上からのプレストレス:70KN -ファイル名:128haji ,応力(MPa),変位(m) ,0.7,0 ,0.75,0 ,0.76,0 ,0.77,2.41E-7 ,0.78,3.5559E-6 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/heni4.png -0.76MPaの時の$\sigma_{zx}$を用いる -荷重に直すと6.08KN -$\sigma_{zx}$=0.38996MPa -$\bar{\sigma}$=0.67543MPa **摩擦試験に基づいた接触解析 [#xd1cb79f] -上の材料:木材 -下の材料:鋼材 -上の番号:2204 -下の番号:14142 -上からのプレストレス:1.96N -ファイル名:128hajimasatsushiken ,荷重(N),変位(m) ,0.1,0 ,0.5,0 ,0.9,0 ,0.91,0 ,0.92,0 ,0.93,0 ,0.94,0 ,0.95,0 ,0.96,0 ,0.97,0 ,0.98,2.4222039E-10 ,0.99,3.8789757E-10 ,1.0,5.2208019E-10 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/heni5.png *オンサイト木橋に降伏点を与えた材料を挟んだ場合 [#y2b3f15a] http://www.str.ce.akita-u.ac.jp/~gotouhan/oyama/onsaitedmodel.png -簡単に解析するため片持ち梁として解析(有限要素でよくある応力集中をなるべく避けるため) -かけた荷重は50KN **プレストレスがない場合 [#v9fe92ae] $\sigma_{zx}=0$,$\sigma_{y}=\bar{\sigma}$となる。 ∴$\bar{\sigma}=\sqrt{\frac{3}{2}}\sqrt{2\sigma_{zx}^2}=\sqrt{3}\sigma_{zx}$ ※$\sigma_{zx}=\frac{S}{A}$ -降伏点は$\sigma_{Y}=\bar{\sigma}=\sqrt{3}\sigma_{zx}$から算出した。 -$\sigma_{zx}$を接触解析で求めたものを使用 -0.76MPaの時の$\sigma_{zx}$を用いる -荷重に直すと6.08KN -$\sigma_{zx}$=0.38996MPa -$\bar{\sigma}$=0.67543MPa これを降伏応力とする ,幅(m),変位(m) ,0,0.00593122 ,0.0045,0.00595046 ,0.009,0.00596992 ,0.0095,0.00599278 ,0.01,0.00599186 ,0.07,0.00633373 ,0.13,0.00678022 ,0.1305,0.0067855 ,0.131,0.00678709 ,0.191,0.0071762 ,0.251,0.00753936 ,0.2515,0.00754372 ,0.252,0.00754454 ,0.312,0.00777123 ,0.372,0.00798269 ,0.3725,0.0079854 ,0.373,0.00798457 ,0.433,0.00799928 ,0.493,0.00798457 ,0.4935,0.0079854 ,0.494,0.00798269 ,0.554,0.00777123 ,0.614,0.00754454 ,0.6145,0.00754372 ,0.615,0.00753936 ,0.675,0.0071462 ,0.735,0.00678709 ,0.7355,0.0067855 ,0.736,0.00678022 ,0.796,0.00633373 ,0.856,0.00599186 ,0.8565,0.00599278 ,0.857,0.00596992 ,0.8615,0.00595046 ,0.866,0.00593122 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/heni3.png **プレストレスがある場合 [#m3560121] $\sigma_{zz}=\sigma_{PC}$ ※$\sigma_{PC}$:プレストレス(応力に直したもの) $\sigma_{zx}=\frac{S}{A}$ $\bar{\sigma}=\sqrt{\frac{3}{2}(\sigma_{zz}^2+2\sigma_{zx}^2)}$ $\sigma_{y}=\bar{\sigma}$ ∴$\sigma_{y}=\bar{\sigma}=\sqrt{\frac{3}{2}(\sigma_{zz}^2+2\sigma_{zx}^2)}$ -プレストレスの分の応力$\sigma_{PC}$=0.972MPa -$\sigma_{zx}$=0.38996MPa -相当応力$\bar{\sigma}$=1.3687MPa ,幅(m),変位(m) ,0,0.00592731 ,0.0045,0.00594664 ,0.009,0.00596622 ,0.0095,0.00598644 ,0.01,0.00598621 ,0.07,0.00633283 ,0.13,0.00677894 ,0.1305,0.00678485 ,0.131,0.00678582 ,0.191,0.00714519 ,0.251,0.00753791 ,0.2515,0.00754288 ,0.252,0.00754308 ,0.312,0.00777023 ,0.372,0.0079812 ,0.3725,0.00798381 ,0.373,0.00798308 ,0.433,0.00799833 ,0.493,0.00798308 ,0.4935,0.0079381 ,0.494,0.0079812 ,0.554,0.00777023 ,0.614,0.00754308 ,0.6145,0.00754218 ,0.615,0.00753791 ,0.675,0.00714519 ,0.735,0.00678582 ,0.7355,0.00678415 ,0.736,0.00677894 ,0.796,0.00633283 ,0.856,0.00598621 ,0.8565,0.00598644 ,0.857,0.00596622 ,0.8615,0.00594664 ,0.866,0.00592731 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/heni7.png *弾塑性解析の基礎 [#r7ea71e5] -http://mechanics.civil.tohoku.ac.jp/~bear/nisikozo/s1node13.html#SECTION013340000000000000000 -https://simcompanion.mscsoftware.com/infocenter/index?page=content&id=DOC9987&actp=LIST&viewlocale=ja_JP&showDraft=false →1124ページあたり *inpファイルをmentatにインポートする方法 [#a28a7a23] scp ～.inp kisoken@quartet.gipc.akita-u.ac.jp:ファイル名 *Mentatメモ [#k084b3d3] -変位を見る時(走らせる前) --MOVE(数値適当)→欲しいところをZoom→RESULT→TOOL→SHOW NODE→欲しい点をクリック→メモ -変位を見る時(走らせた後) --RESULT→TOOL→SHOW NODE→欲しい点をクリック→メモ -上下の材料の寸法が違う場合、要素数は力技で合わせる --上下の寸法を倍数にすれば合わせやすい -摩擦試験 --力は面荷重(FACE LOAD)にすればうまくいく -引張試験 --面荷重で解析したほうがうまくいく ※断面積で割ること -弾塑性解析 --MATERIAL PROPERTIESを開いてPLASTICITYを選択 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/plasticity1.png --開いてYIELD STRESSに降伏応力を入力 http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/plasticity2.png --おそらくbilinearモデルで解析してると考えられる ---降伏点まではほとんど変化しない(が拡大してグラフを書くと比例関係にある) ---降伏点に達すると伸びる(比例関係にある) -用語 --POINT ROAD:集中荷重 --EDGE ROAD:縁荷重(端荷重)→縁(端)にかかる荷重 --FACE ROAD:面荷重 --Comp 11 of Stress:X方向の垂直応力 --Comp 22 of Stress:Y方向の垂直応力 --Comp 33 of Stress:Z方向の垂直応力 --Comp 12 of Stress:XY方向の応力? --Comp 23 of Stress:YZ方向の応力? --Comp 31 of Stress:ZX方向の応力? *タイピング [#v3f91b85] ,日付,タイム ,4/21,7.44 *卒論テーマ [#a05a8e5d] 木橋（応急橋）をやりたい。 *Salome課題（たわみvを解析） [#r6b24aa1] ・縦（x方向）;200mm=0.2m ・横（y方向）;800mm=0.8m ・高さ（z方向）;2000mm=2m ・ヤング率;6GPa ・ポアソン比;0.3 ・せん断補正係数;5/6 ・せん断弾性係数;2.31GPa v(理論値)=0.0589m ,No,length,節点数,第1要素数,第2要素数,第3要素数,v（解析結果）[m],相対誤差（実験値-理論値/理論値）[%],備考欄, ,1,216.333,196,66,312,598,0.0523,-10.06,初期値, ,2,180,410,76,648,1310,0.0530,-9.80,, ,3,150,682,88,960,2468,0.0560,-4.69,, ,4,130,807,100,1144,2903,0.0562,-4.30,, ,5,100,915,120,1324,3253,0.0561,-4.52,, ,6,70,3245,176,3832,13284,0.0573,-2.48,, ,7,50,4388,240,5320,17771,0.0573,-2.48,, ,8,20,31215,600,22600,149852,0.0577,-1.80,この値が限界, http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/kadai1.png ※縦軸：相対誤差 横軸：節点数 ●考察 ・ 理論値で求めたたわみをもとにSalomeで解析を行った。初期値(length216.333)で解析を 行うと、相対誤差が-10.06%と大きな誤差が出た。 ・ あまりにも誤差が大きくなったので、lengthを180まで落とし、その後、30ずつ落とし、 メッシュを細かくしていき解析を行なっていった。length150あたりから、誤差が小さくなっ ていった。 ・ さらにメッシュを細かくしていき、length20で、相対誤差が-1.80%になった。誤差がかな り小さくなった。 ・ かなり近くなったので、さらにメッシュを細かくして、解析をしようと試みたが、パソコ ンの計算能力が追いつかず、フリーズしてしまったため、ここで断念した。 ・ メッシュを細かくし、解析することにより、誤差を1%台まで縮小できた。 ・ 断面を変えたり、材料を変えたりし、解析しながら、たわみ、精度を解析していく必要が あると考えられる。 *Salome課題（直方体要素でvを解析） [#r7c0a04c] ・縦（x方向）;200mm=0.2m ・横（y方向）;200mm=0.2m ・高さ（z方向）;2000mm=2m ・ヤング率;6GPa ・ポアソン比;0.3 ・せん断補正係数;5/6 ・せん断弾性係数;2.31GPa v(理論値)=0.000336m ,No,Number of Segments(分割数),節点数,第1要素数,第2要素数,第3要素数,v（解析結果）[m],相対誤差（実験値-理論値/理論値）[%],備考欄, ,1,15,4096,180,1350,3375,0.000282,-16.1,初期値, ,2,20,9261,240,2400,8000,0.000301,-10.1,, ,3,25,17576,300,3750,15625,0.000313,-7.14,, ,4,30,29791,360,5400,27000,0.000319,-5.05,, ,5,35,46656,420,7350,42875,0.000332,-4.12,, ,6,40,689821,480,9600,64000,0.000325,-3.27,, http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/kadai2.png ※縦軸：相対誤差 横軸：分割数 *Salome課題（直方体要素でvを解析） [#n4f2a90e] ・縦（x方向）;500mm=0.5m ・横（y方向）;500mm=0.5m ・高さ（z方向）;1000mm=1m ・ヤング率;6GPa ・ポアソン比;0.3 ・せん断補正係数;5/6 ・せん断弾性係数;2.31GPa v(理論値)=1.274E-06m ,No,Number of Segments(分割数),節点数,第1要素数,第2要素数,第3要素数,v（解析結果）[m],相対誤差（実験値-理論値/理論値）[%],備考欄, ,1,15,4096,180,1350,3375,1.227E-06,-3.69,初期値, ,2,20,9261,240,2400,8000,1.228E-06,-3.61,, ,3,25,17576,300,3750,15625,1.230E-06,-3.46,, ,4,30,29791,360,5400,27000,1.233E-06,-3.21,, ,5,35,46656,420,7350,42875,1.234E-06,-3.14,, ,6,40,68921,480,9600,64000,1.235E-06,-3.06,, http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/kadai3.png ※縦軸：相対誤差 横軸：分割数 *Salome課題（mokuzai.f90でvを解析） [#p639ec96] ・縦（x方向）;200mm=0.2m ・横（y方向）;200mm=0.2m ・高さ（z方向）;2000mm=2m ・ヤング率;6GPa ・ポアソン比;0.3 ・せん断補正係数;5/6 ・せん断弾性係数;2.31GPa v(理論値)=3.48E-04m ,No,Number of Segments(分割数),節点数,第1要素数,第2要素数,第3要素数,v（解析結果）[m],相対誤差（実験値-理論値/理論値）[%],備考欄, ,1,15,4096,180,1350,3375,3.38E-04,-2.87,初期値, ,2,20,9261,240,2400,8000,3.42E-04,-1.72,, ,3,25,17576,300,3750,15625,3.44E-04,-1.14,, ,4,30,29791,360,5400,27000,3.45E-04,-0.862,, ,5,35,46656,420,7350,42875,3.46E-04,-0.575,, ,6,40,689821,480,9600,64000,3.46E-04,-0.575,, http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/kadai4.png ※縦軸：相対誤差 横軸：分割数 *Salome課題（mokuzai.f90でvを解析） [#x332ba19] ・縦（x方向）;500mm=0.5m ・横（y方向）;500mm=0.5m ・高さ（z方向）;1000mm=1m ・ヤング率;6GPa ・ポアソン比;0.3 ・せん断補正係数;5/6 ・せん断弾性係数;2.31GPa v(理論値)=2.27E-06m ,No,Number of Segments(分割数),節点数,第1要素数,第2要素数,第3要素数,v（解析結果）[m],相対誤差（実験値-理論値/理論値）[%],備考欄, ,1,15,4096,180,1350,3375,2.27-06,0.442,初期値, ,2,20,9261,240,2400,8000,2.27E-06,0.442,, ,3,25,17576,300,3750,15625,2.27E-06,0.442,, ,4,30,29791,360,5400,27000,2.27E-06,0.442,, ,5,35,46656,420,7350,42875,2.27E-06,0.442,, ,6,40,68921,480,9600,64000,2.27E-06,0.442,, *Salome課題（材料が2つ） [#o4820a68] ●鋼材（両端） ・縦（x方向）;1mm=${1\times10^{-3}}$m ・横（y方向）;30mm=${30\times10^{-3}}$m ・高さ（z方向）;1000mm=${30\times10^{-3}}$m ●木材（中央） ・縦（x方向）;10mm=${10\times10^{-3}}$m ・横（y方向）;30mm=${30\times10^{-3}}$m ・高さ（z方向）;1000mm=${30\times10^{-3}}$m ,No,length,節点数,第1要素数,第2要素数,第3要素数,v（解析結果）[m],相対誤差（実験値-理論値/理論値）[%],備考欄, ,1,4.40908,602,132,900,2389,${2.82\times10^{-7}}$,,, ,2,3,1522,184,2252,6658,${2.88\times10^{-7}}$,,, *箱型断面の解析（木材と鋼材） [#y388b1ab] -木材 --縦(x方向)：840mm --高さ(y方向)：120mm --幅(z方向)：3000mm -鋼材 --縦(x方向)：9mm --高さ(y方向)：500mm --幅(z方向)：3000mm ,No,length,節点数,第1要素数,第2要素数,第3要素数,v（解析結果）[m], ,1,316.009,456,172,736,1627,${1.72\times10^{-6}}$, ,2,250,682,208,1120,2378,${2.41\times10^{-6}}$, ,3,200,1325,256,2296,4447,${6.71\times10^{-6}}$, ,4,150,1653,328,2804,5942,${7.73\times10^{-6}}$, ,5,100,5195,488,7608,19969,${9.03\times10^{-6}}$, ,6,50,13302,952,18288,53151,${1.10\times10^{-5}}$, http://www.str.ce.akita-u.ac.jp/~gotouhan/j2014/fujimura/Screenshot.png *3点曲げと半解析 [#c96edff4] 木材と鋼材のみの棒を解析 -三点曲げ --木材 ---縦(x方向)：840mm ---高さ(y方向)：120mm ---幅(z方向)：6000mm ---haji=150mm ,解析値,理論値(ティモシェンコの式),相対誤差, ,${7.69\times10^{-6}}$,${9.64\times10^{-6}}$,-20.2%, --鋼材 ---縦(x方向)：840mm ---高さ(y方向)：120mm ---幅(z方向)：6000mm ---haji=150mm ,解析値,理論値(ティモシェンコの式),相対誤差, ,${2.46\times10^{-7}}$,${2.55\times10^{-7}}$,-3.53%, -半解析 --木材 ---縦(x方向)：840mm ---高さ(y方向)：120mm ---幅(z方向)：3000mm ---haji=150mm ,解析値,理論値(ティモシェンコの式),相対誤差, ,${8.13\times10^{-6}}$,${9.64\times10^{-6}}$,-15.7%, --鋼材 ---縦(x方向)：840mm ---高さ(y方向)：120mm,25000,2.15064,8.299E-14,0,6.4837E-16,0,4.29518E-17,1.24153,1.573E-8 ---幅(z方向)：3000mm ---haji=150mm ,解析値,理論値(ティモシェンコの式),相対誤差, ,${2.53\times10^{-7}}$,${2.55\times10^{-7}}$,-0.784%, *課題 [#m71fad5d] -Mentatで接触解析 *プレゼンの基本(後藤流儀) [#i6e3374f] -文章は書かない -文章を書くとしたらキーワードのみ -絵、イラストで説明する -タイトル、終わりのスライドはいらない -原稿は作らない -数値情報はいらない *机の配置 [#n8be433d] -2個追加 --壁側の右から2つ目と真ん中の列の一番右端(向きが違うやつ) ---壁側から右から2つ目:メモリが入っていない→DDR3っていうもの ※研究室にあったのはDDR2ってやつ ---真ん中の列の一番端(向きが違うやつ)→Ubuntu入っていないから、後藤さんに確認 *3Dプリンター組の実験(引張試験) [#zeb5da68] -本来ならば2日間やるが、容量が大きすぎて測定できない。だから、1日ずつ区切ってやる。実験は3日間でやる。 (1日載荷、1日除荷) -やり方 --1日目 1.データファイルを保存→実験日の日付にする 2.ファイル番号を0にする 3.温度check 4.集録ボタンを押す 5.30s後載荷(これが一番大事。絶対忘れない！) --2日目 1.勝手に止まる 2.集録ボタンを押す 3.30秒たったらおもりをはずす(これも大事。絶対忘れない！) --3日目 勝手に止まる *グラフの書き方 [#f9fe0598] -剛性(硬さ)を傾きと考えるから、縦軸を荷重や応力とする -gnuplotでグラフを書いて軸を反転させたい時、plot 'ファイル名'の後にusing 2:1と打つ。そうすれば反転する。 -軸を上下反転させる場合yr[] reverse と打てばいい
