川村の卒論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
*卒論日誌 [#lc8ddf7f]
,日付,時間帯,作業時間,内容,立合
,4/7,10:00～11:00,1h,卒論,
,5/21,14:00～16:30,2.5h,ゼミ課題（サロメ）,
,5/23,16:00～18:30,2.5h,ゼミ課題（サロメ）,
,5/26,14:00～17:00,3h,ゼミ課題（サロメ）,
,5月合計,,8h,,
,6/1,16:00～18:30,2.5h,ゼミ課題（サロメ）,
,6/2,13:00～18:00,5h,ゼミ課題（サロメ）,
,6/3,15:00～18:00,3h,ゼミ課題（サロメ）,
,6/7,15:00～18:00,3h,ゼミ課題,
,6/8,13:00～20:00,5.5h,ゼミ課題,
,6/13,13:30～17:00,3.5h,ゼミ課題,
,6/15,16:30～20:00,3.5h,ゼミ課題,
,6/18,14:30～16:00,1.5h,ゼミ課題,
,6/20,17:00～18:00,1h,ゼミ課題,
,6/21,16:00～18:00,2h,ゼミ課題,
,6/27,13:30～16:00,2.5h,ゼミ課題、外論,
,6/28,13:00～16:00,3h,ゼミ課題,
,6/29,13:00～17:00,4h,ゼミ課題,
,6月合計,,40h(48h),,
,7/4,15:00～18:00,3h,ゼミ課題,
,7/5,13:00～15:30,2.5h,ゼミ課題,
,7/10,13:00～17:00,4h,ゼミ課題,
,7/11,15:00～18:00,3h,ゼミ課題,
,7/12,13:00～19:30,6.5h,ゼミ課題,
,7/13,13:00～16:30,3.5h,ゼミ課題,
,7/19,13:00～17:00,4h,ゼミ課題,
,7/20,13:30～15:30,2h,ゼミ課題,
,7月合計,,28.5h(76.5h),,
,8/2,13:00～17:30,4.5h,中間,
,8/17,14:30-17:30,3h,中間,
,8/20,15:30-18:00,2.5h,中間,
,8/23,14:00～19:30,5.5h,中間（エラーなおらない）,
,8/28,14:00～19:30,5.5h,中間,
,8/31,15:00～18:30,3.3h,中間（エラー直った）,
,8月合計,,24.5h(101),,
,9/4,14:00～20:30,6.5h,中間,
,9/6,13:00～18:00,5h,中間と実験手伝い,
,9/7,13:00～16:00,3h,中間,
,9/12,14:30～18:30,4h,中間,
,9/13,14:00～18:00,4h,中間,
,9/14,13:00～17:30,4.5h,vi復習,
,9/15,13:00～19:00,4.5h,データ編集,
,9/16,14:30～20:30,4h,概要書く,
,9/19,15:00～21:00,6h,概要書く,
,9/20,11:00～19:00,8h,概要、スライド,
,9/21,15:00～18:30,3.3h,概要、スライド,
,9/26,14:00～15:30,1.5h,論文読む,
,9月合計,,54.5h(155.5h),,
,10/10,13:00~16:30,3.5h,接触解析練習,
,10/11,15:00~19:00,4h,接触解析練習,
,10/12,13:30~18:00,4.5ｈ,接触解析練習と論文調べる,
,10/16,13:00~18:30,5.5h,接触解析,
,10/18,13:00~18:30,5.5h,接触解析（エラー出まくる）と摩擦...
,10/19,13:00~16:00,3h,接触解析と摩擦練習,
,10/23,13:00~17:30,4.5h,接触解析,
,10/24,13:30~17:30,4h,摩擦解析練習,
,10/25,14:30~18:30,4h,摩擦解析練習,
,10/26,14:00~19:30,5.5h,摩擦解析練習,
,10/31,13:30~17:00,3.3h,バネについて検証,
,10月合計,,47.5h(203h),,
,11/1,14:30~18:30,4h,バネについて,
,11/2,13:00~17:00,4h,バネについて,
,11/7,15:00~18:00,3h,バネについて,
,11/8,13:00~17:00,4h,バネについて,
,11/9,13:00~16:30,3.5h,バネ,
,11/13,13:00~17:30,4.5h,摩擦,
,11/14,13:00~16:30,3.5h,摩擦,
,11/15,13:30~18:30,5h,摩擦,
,11/16,11:30~15:30,4h,摩擦,
,11/17,16:00~18:30,2.5h,二面摩擦,
,11/20,13:00~18:00,5h,二面摩擦,
,11/21,13:00~15:00,2h,二面摩擦（収束しない）,
,11/22,13:00~18:00,5h,二面摩擦,
,11/23,17:00~21:30,4.5h,二面摩擦,
,11/27,17:00~18:00,1h,一面摩擦(計算速度遅い),
,11/28,13:00~17:00,4h,まだ計算終わらない,
,11/29,14:30~,3h,１面摩擦と２面摩擦,
,11/30,13:00~20:30,7.5h,1面摩擦と2面摩擦,
,11月合計,,86.5h(289.5h),,
,12/5,14:00~,4h,salome,
,12/6,13:00~,4.5h,salome（収束しない＆結果安定せず）,
,12/7,13:00~,8h,salome,
,12/12,13:00~16:30,3.5h,salome(ボルト付き),
,12/13,13:00~16:00,3h,salome,
,12/14,13:00~17:00,4h,salome,
,12/16,22:30~25:00,2.5h,salome,
,12/19,14:30~16:30,2h,salome,
,12/21,14:00~,4h,salomeと概要,
,12/24,15:00~,5h,概要,
,12/25,16:30~24:00,7.5h,スライド,
,12月合計,,48h(337.5h),,
,1/6,17:00~18:00,1h,salome,
,1/9,13:00~17:00,4h,salome(ミーゼス応力とか),
,1/10,11:30~14:00,2.5h,salome,
,1/11,17:00~20:00,3h,salome,
,1/12,14:00~15:30,1.5h,salome,
,1/15,15:00~18:00,3h,salome,
,1/16,13:00~17:00,4h,salome,
,1/17,15:00~18:30,3.5h,支部概要,
,1/18,15:00~,3h,支部概要,
,1/19,12:00~,3h,支部概要,
,1/22,13:00~18:30,5.5h,支部概要,
,1/23,13:30~18:30,5h,支部概要,
,1/24,13:00~,2h,支部概要,
,1/25,14:30~,2h,スライド,
,1/29,13:30~,3h,salome,
,1/30,13:00~,4h,salomeとスライド,
,1/31,16:00~18:30,2.5h,スライド,
,1月合計,,51.5h(389h),,
,2/1,14:00~19:00,5h,スライド,
,2/2,13:00~,2.5h,スライド,
,2/5,13:00~,4.5h,スライド・salome,
,2/6,12:00~,4h,発表練習,
,2/7,13:00~,3.5h,スライド,
,2/8,14:30~19:00,4.5h,スライド,
,2/9,13:00~16:30,3.5h,スライド,
,2/13,14:00~18:30,4.5h,スライド,
,2/14,14:00~22:00,5h,スライド,
,2/15,13:00~19:00,5h,スライド,
,2月合計,,42h(431h),,
*5/8課題結果 [#y6602ec0]
10mm×10mm断面の角材（ヤング率：6GPa, ポアソン比：0.3）の...
荷重は面荷重で100N.
スパン長さは、100mmを分担。
理論値は初等梁＝0.4167mm、ティモシェンコ梁＝0.4297mm
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
<一次要素>
,length,8,4,2,1
,たわみ(mm),0.3549,0.3549,0.4177,0.4177
,相対誤差(理論値：初等梁),14.83%,14.83%,1.20%,1.20%
,相対誤差(理論値：ティモシェンコ梁),17.41%,17.41%,4.19%,4...
<二次要素>
,length,8,4,2,1
,たわみ(mm),0.4338,0.4338,0.4338,0.4338
,相対誤差(理論値：初等梁),4.10%,4.10%,4.10%,4.10%
,相対誤差(理論値：ティモシェンコ梁),0.95%,0.95%,0.95%,0.95%
*5/15課題結果 [#y139767b]
木材（異方性）についてサロメでたわみを見た。
スパンは100mm。荷重は線荷重で100N。
理論値は初等梁＝0.4167mm、ティモシェンコ梁＝0.4917mm
<一次要素>
,length,8,4,2,1
,たわみ(mm),0.4327,0.4699,0.5002,0.5361
,相対誤差(理論値：初等梁),3.84%,12.77%,20.04%,28.65%
,相対誤差(理論値：ティモシェンコ梁),12.00%,4.43%,1.735,9....
<二次要素>
,length,8,4,2,1
,たわみ(mm),0.5378,0.5485,0.5555,0.5711
,相対誤差(理論値：初等梁),29.06%,31.63%,33.31%,37.05%
,相対誤差(理論値：ティモシェンコ梁),9.38%,1.16%,1.30%,1.61%
*5/22課題結果 [#m8109b73]
サンドイッチ梁（縦に1mm幅の鋼材を入れた梁）を作ってたわ...
木材のヤング率等は前回の課題と同じ。鋼材はE=206GPa、ν=0....
また、レンクスは１、メッシュは2次要素で切った。
たわみの理論値は0.09398mmとなった。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
,,たわみ(mm),相対誤差
,木材を等方性とした時,0.1119,19.07%
,木材を異方性とした時,0.1708,81.74%
<木材を等方性としたとき>（E=6GPa、ν=0.4）
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
<木材を異方性としたとき>
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
誤差が大きくなってしまったのでせん断係数も考慮した式を使...
*5/29課題結果 [#ub253aa2]
海老さんから卒論のデータをもらい、同じようにサロメで結果...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
同じように結果を出すことができた。Eficasのいじり方がよく...
*6/5課題結果 [#y1f6b79b]
前回の課題のデータを使い、サロメと手計算それぞれで求めた...
サロメでひずみが出せなかったので、応力からひずみを出した。
,,手計算,サロメ
,応力(N/m^2),3.6393 10^6(木),5.8131 10^6
,,97.438 10^6(鋼),
,ひずみ,4.73 10^-4,7.56 10^-4(木)
,,,0.282 10-4(鋼)
*6/12課題結果 [#z5e80853]
6/22のサンドイッチ梁において、サロメと手計算それぞれで求...
,,手計算,サロメ
,ひずみ,5.639 10^-4,5.2550 10^-4
,応力(N/mm^2),116.1634(鋼)/3.3834(木),115.5825(鋼)/3.0491...
*6/19課題結果 [#j5f53223]
前回の課題でのひずみの分布の仕方のグラフ。5.2～5.8mmは鋼...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
それから、サロメで求めたひずみにヤング率をかけて、サロメ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
*6/26課題結果 [#f90cdb70]
前回と同じサンドイッチ梁の木材部分において、
&link(この行列,http://www.str.ce.akita-u.ac.jp/~gotou/kou...
,,行列によって求めた応力,サロメで求めた応力
,0mm,2.8686,2.8688
,2.5mm,3.0521,3.0661
,8mm,3.0970,3.0967
小数点以下2桁あたりまで一致した。
*7/3課題結果 [#c72fb25b]
近藤さんからもらったデータを使い、プレストレス木箱桁橋の...
,,荷重(MN/mm^2),ひずみ(10^-4)
,対傾構あり積雪あり,0.0136759,4.630721
,対傾構あり積雪なし,0.0076566,2.592559
,対傾構なし積雪あり,0.0135336,4.589838
,対傾構なし積雪なし,0.0075144,2.548455
対傾構あり積雪ありの時のひずみの様子。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
対傾構あり積雪なしの時のひずみの様子。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
対傾構なし積雪ありの時のひずみの様子。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
対傾構なし積雪なしの時のひずみの様子。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
これらのひずみと荷重の関係を表にしてみた。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
対傾構ありのほうが対傾構なしに比べてややひずみが大きかった
*中間発表 [#mc9f109b]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
のエラーがなおらなくて半月以上たってしまったので、試しに...
（固定しすぎていたら3方向固定から2方向固定にするorバネは...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
みたいになった。ネジ穴周りは
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
穴周りで色が変わるかと思ったけど、この場合は一色であった...
やっとエラーが直った。メッシュの切り方を変えたら上手く行...
-添接板のネジ穴2列
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-添接板のネジ穴3列
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
ネジ穴2列と3列のときの添接板のネジ穴周りの応力を比べた。...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
穴3列のほうが応力が大きかった。
$\sigma$=$E$$\epsilon$ より穴6つの方がひずみも大きくなる...
*10/6ゼミ [#w4b38148]
次回までに
-摩擦を考慮したモデルで解析する
-道路橋示方書にあった論文を探して読む
今週できたことなど
-まず本を読んで接触解析を試しにやってみた(下図)。本に書い...
-論文は取り寄せすることになり、早くて来週届くそうだ。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
*10/13ゼミ [#mc34d049]
次回までに
-自分がやってみたいモデルで接触解析をやる。必要があれば海...
今週できたことなど
-自分がやってみたいモデルに接触解析を行ってみたが、エフィ...
-摩擦の解析の練習も少し手を付けたけど、エフィカスのエラー...
-文献は無事に取り寄せできた
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-アラームはあるが結果は見ることができた。しかし、バネによ...
-荷重を二段階に分けるのは上手く行ったがやはりバネ定数を小...
*10/20ゼミ [#k0d3346a]
次回までに
-去年の先輩が解いた摩擦解析の問題をもらって解く
-余裕があったら文献読む
今週出来たことなど
-去年の4年の菊池さんの摩擦のモデルを解いてみた
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
菊池さんのデータを応用して柔らかいゴムではなく弱いバネを...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
*10/27ゼミ [#q6b74772]
次回までに
-論文少しでも読む！！！！
-高校の物理と対応した摩擦解析の問題がほしい。今度は壁から...
-http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin=...
今週出来たことなど
-前回のモデルにおいて、時間の分割数を100くらいにして動き...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-下図の上の物体の面の真ん中（線2本が交わっている所）にバ...
,バネ(x/y/z),変位(引っ張った方向),最大のひずみ,最大の応力
,0/0/0,1.439,15.6471,490
,10/10/10,1.200,15.6331,500
,10/0/0,1.1998,15.6334,490
,50/50/50,1.199,15.6382,490
,50/0/0,1.199,15.6394,490
,100/100/100,1.200,15.6446,490
,100/0/0,1.199,15.6471,490
今週できたことなど
-先週のモデルで間違いなどあったのでまたバネ定数などを変え...
,バネ定数(N/mm),静止摩擦係数,動き出す時の荷重(N)(手計算),...
,10,0.4,40,0.5,1,4.8077
,100,0.4,40,0.5,0.1,0.4981
,10,1,動かない,0.5,動かない,4.8077
,バネ無し,0.4,40,0.5,？,125
手計算とsalomeでの結果が大きく異なる結果となった。また、s...
*11/10ゼミ [#va29bacb]
次回までに
-摩擦係数の入れ方の確認を早急にやる。最初から滑り出さない...
11/13 菊地さんのデータの時間を1000分割してみが、やはり...
11/15 ついに摩擦を取り入れた計算が上手く行った。マスタ...
菊地さんのK2のデータのスレーブ面とマスター面を逆にすると...
,静止摩擦係数,動き出す時の荷重(N)(手計算),動き出した時の...
,0.4,40,40,149.494
,0.6,60,60,99.4478
*11/17ゼミ [#kfa8326d]
次回までに
-軸力を示方書で確認
-添接板での一面摩擦と二面摩擦
-ボルトをつけたモデルでやってみる
今週出来たことなど
-添接板に軸力が均等にかかっている場合の滑り出し方の比較の...
--一面摩擦 添接板二枚の上面に165kN×2(ボルトの軸力)をかけ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
--2面摩擦 こちらも真ん中の母材が飛び出してしまった（Salom...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
滑り出し方をグラフにしてみた。左のグラフが滑りだす瞬間ま...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-ボルトをつけたモデルは収束ができなくて結果が出なかった
*11/24ゼミ [#k4d6f0f2]
次回までに
-添接板に軸力が理想的（均等）に働いた時の滑り出し方の、一...
-ボルト付きの解析を引き続きやる
今週できたことなど
-１面摩擦を添接板１枚分で行った。
--境界条件など・・・軸力165kNを面載荷してから横へ引っ張っ...
--結果・・・引っ張る力がサロメでは66834Nの時動いた。手計...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-2面摩擦を行った。
--境界条件など・・・上下に軸力をかけて引っ張った。境界条...
--結果・・・引っ張る力が56880Nの時に滑りだした。前回やっ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-ボルト付きでもやってみた。計算は上手く行ったが全く滑り出...
*12/1ゼミ [#m0bd456b]
次回までに
-引き続き解析する
今週できたことなど
-引っ張る力を100分割にした一面摩擦の解析が上手くいった
-①一面摩擦、軸力×１のとき
--手計算・・・引っ張る力が165kN×0.4＝66000Nのとき滑り出す
--salome・・・引っ張る力が66834Nのとき滑り出す。手計算と...
-②一面摩擦、軸力×２のとき
--手計算・・・引っ張る力が165kN×2×0.4＝132000Nのとき滑り...
--salome・・・引っ張る力が136867.5Nのとき滑りだす。手計算...
-二面摩擦を行った。引っ張る力が135090Nのときに滑りだした...
-この二面摩擦のモデルにボルトをつけたモデルも解析した。な...
*12/8ゼミ [#x30c4047]
次回までに
-ボルト付き出やってみる
-概要書き始める
-論文を読みなおす
今週できたことなど
-収束判定を厳しめだと収束しない。あとボルト穴まわりに荷重...
-概要はちょっと手を付けた
*12/15ゼミ [#c30d4fe4]
次回までにと今後やること
-解ける所まででよい
-年明けに穴2つ3つの場合と全体をやってみる。圧縮側は穴を潰す
今週できたことなど
-ボルトをつけたモデルで滑り出しは見られなかったが、引っ張...
*中間発表後 [#o90f23f9]
次回までに
-母材のミーゼス応力の分布を見る
-母材が滑った後のモデル（母材とボルトが接触している）を作...
今週できたことなど
-母材が滑った後のモデルの様子
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-ミーゼス応力の分布の仕方
--滑りだす前、ボルトを締めた直後
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
--滑りだす前、母材を203MPaで引っ張ったとき
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
--滑りだした後のモデル（ボルトを締めた直後）
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
ボルトを締めただけで鋼材が降伏しそうだ。同じような条件で...
-リラクセーションを考慮した軸力をかけたのもやってみたが降...
*1/26ゼミ [#aa9a6483]
次回までに
-ボルト3列のモデルの解析
今週できたことなど
-フランジ部分の半解析(ボルト3本)をした。左から座金（最大6...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-スライドに画像並べてて違いがあまりないと感じたのでグラデ...
*支部用タイトル [#rfca77ef]
「プレストレス木箱桁橋のボルト接合部の数値モデル化」
*卒論タイトル [#ycae6e6c]
「プレストレス木箱桁橋のボルト接合部の挙動」
*ボルトメモ [#k989eb55]
-&link(座金の硬度,http://www.nssmc.com/product/constructi...
-&link(鋼の「硬度とヤング率」、「硬度と強度」の関,http://...
-座金の硬度は35$\sim$45HRC
-F10Tの硬度は27$\sim$38HRC
-座金の引張強さの推定：$引張強さ=3.2\times HRC\times 9.8$...
-座金の降伏応力の推定：$降伏応力=0.5\times 引張強さ$(MPa)...
*salomeメモ [#yfb9da81]
-RESI_GLOB_RELAの値が小さいほうが手計算とあっている気がす...
-参考にしていたサイト
--http://salome-meca.cocolog-nifty.com/blog/
--https://sites.google.com/site/codeastersalomemeca/
-11月らへんに計算した2面摩擦モデルでcommファイルの覚え書き
--(
DEBUT();
kou=DEFI_MATERIAU(ELAS=_F(E=206000.0, #ヤング率（単位：MP...
NU=0.3,),); #ポアソン比
MIAL=LIRE_MAILLAGE(FORMAT='MED',);
MIAL=MODI_MAILLAGE(reuse =MIAL,
MAILLAGE=MIAL,
ORIE_PEAU_3D=_F(GROUP_MA=('men2','men1...
newmesh=CREA_MAILLAGE(MAILLAGE=MIAL, ＃バネについて
CREA_POI1=_F(NOM_GROUP_MA='bane',
GROUP_NO=('bane1','ban...
MODE=AFFE_MODELE(MAILLAGE=newmesh, ＃「MIAL」も選択でき...
AFFE=(_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='3D',),
_F(GROUP_MA='bane', ＃CREA_MAILLA...
PHENOMENE='MECANIQUE',
MODELISATION='DIS_T',),),);
MATE=AFFE_MATERIAU(MAILLAGE=newmesh, ＃材料の種類。材料...
AFFE=_F(TOUT='OUI',
MATER=kou,),);
sktdn=AFFE_CARA_ELEM(MODELE=MODE,
DISCRET=_F(CARA='K_T_D_N',
GROUP_MA='bane',
VALE=(0.1,0.1,0.1,),),);...
contact=DEFI_CONTACT(MODELE=MODE,
FORMULATION='DISCRETE',
FROTTEMENT='COULOMB', ＃摩擦係数を...
ITER_GEOM_MAXI=100000, ＃
ZONE=(
_F(GROUP_MA_MAIT='men2', ＃マスター...
GROUP_MA_ESCL='men1', ＃スレーブ...
ALGO_CONT='PENALISATION', ＃ペナ...
E_N=206000.0, ＃ヤング率
COULOMB=0.4, ＃摩擦係数
ALGO_FROT='PENALISATION',
E_T=20600.0,), ＃ヤング率の10分...
_F(GROUP_MA_MAIT='men3',
GROUP_MA_ESCL='men4',
ALGO_CONT='PENALISATION',
E_N=206000.0,
COULOMB=0.4,
ALGO_FROT='PENALISATION',
E_T=20600.0,),),);
kotei=AFFE_CHAR_MECA(MODELE=MODE, ＃固定面。＃載荷、固定...
DDL_IMPO=_F(GROUP_MA=('kotei1','kote...
DX=0.0, ＃２方向固定に...
DZ=0.0,),);
osae=AFFE_CHAR_MECA(MODELE=MODE,
FORCE_FACE=(_F(GROUP_MA='joumen',
FY=20.24025,), ＃面載...
_F(GROUP_MA='kotei',
FY=-20.24025,),),);
tikara=AFFE_CHAR_MECA(MODELE=MODE,
FORCE_FACE=_F(GROUP_MA='saika', ＃...
FZ=200.0,),);
tikara_f=DEFI_FONCTION( ＃折れ線グラフの&link(イメ...
NOM_PARA='INST',
VALE=(0.0 ,0.0 , ＃VALE=(ｘ座...
10.0 ,0.0 ,
20.0 ,1.0 ,),); ＃10~20ｓで...
osae_f=DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0.0 ,0.0 ,
10.0 ,1.0 ,
20.0 ,1.0 ,),); ＃0~10ｓで徐...
kankaku=DEFI_LIST_REEL(DEBUT=0.0, ＃時間の分...
INTERVALLE=(_F(JUSQU_A=10.0,
NOMBRE=1,), ＃0~10...
_F(JUSQU_A=20.0,
NOMBRE=20,),),); ...
resu=STAT_NON_LINE(MODELE=MODE,
CHAM_MATER=MATE,
CARA_ELEM=sktdn,
EXCIT=(_F(CHARGE=kotei,), ＃はじめか...
_F(CHARGE=osae,
FONC_MULT=osae_f,),
_F(CHARGE=tikara,
FONC_MULT=tikara_f,),),
CONTACT=contact,
COMPORTEMENT=_F(RELATION='ELAS',),
INCREMENT=_F(LIST_INST=kankaku,),
NEWTON=_F(REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=0.01, ...
ITER_GLOB_MAXI=300,), ...
SOLVEUR=_F(SYME='OUI',),);
resu=CALC_CHAMP(reuse =resu,
MODELE=MODE,
CHAM_MATER=MATE,
RESULTAT=resu,
CONTRAINTE=('SIGM_ELNO','SIGM_NOEU',), ＃...
DEFORMATION=('EPSI_ELNO','EPSI_NOEU',), ...
CRITERES=('SIEQ_ELNO','SIEQ_NOEU',),); ...
IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=newmesh,
RESULTAT=resu,
NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','EPSI...
FIN();
#CHECKSUM:b2f4a2740b714f58c842db0021057ff8 -:FIN CHECKSUM
--)
-バネ付きモデルのときのcommファイルの覚え書き
--(
・
・
・
contact=DEFI_CONTACT(MODELE=MODE,
FORMULATION='CONTINUE',
FROTTEMENT='COULOMB',
ALGO_RESO_GEOM='POINT_FIXE',
REAC_GEOM='AUTOMATIQUE',
ITER_GEOM_MAXI=100,
RESI_GEOM=0.1,
ALGO_RESO_CONT='POINT_FIXE',
ITER_CONT_MAXI=100,
ALGO_RESO_FROT='NEWTON', #「ALGO_RE...
RESI_FROT=0.5,
ZONE=(
_F(GROUP_MA_MAIT='men1',
GROUP_MA_ESCL='boruto_m',
CONTACT_INIT='INTERPENETRE',
ALGO_CONT='PENALISATION',
COEF_PENA_CONT=206005.0,
COULOMB=0.4,
SEUIL_INIT=0.1, #これを入れる...
ALGO_FROT='PENALISATION',
COEF_PENA_FROT=20600.0,),
--)
終了行:
#contents
*卒論日誌 [#lc8ddf7f]
,日付,時間帯,作業時間,内容,立合
,4/7,10:00～11:00,1h,卒論,
,5/21,14:00～16:30,2.5h,ゼミ課題（サロメ）,
,5/23,16:00～18:30,2.5h,ゼミ課題（サロメ）,
,5/26,14:00～17:00,3h,ゼミ課題（サロメ）,
,5月合計,,8h,,
,6/1,16:00～18:30,2.5h,ゼミ課題（サロメ）,
,6/2,13:00～18:00,5h,ゼミ課題（サロメ）,
,6/3,15:00～18:00,3h,ゼミ課題（サロメ）,
,6/7,15:00～18:00,3h,ゼミ課題,
,6/8,13:00～20:00,5.5h,ゼミ課題,
,6/13,13:30～17:00,3.5h,ゼミ課題,
,6/15,16:30～20:00,3.5h,ゼミ課題,
,6/18,14:30～16:00,1.5h,ゼミ課題,
,6/20,17:00～18:00,1h,ゼミ課題,
,6/21,16:00～18:00,2h,ゼミ課題,
,6/27,13:30～16:00,2.5h,ゼミ課題、外論,
,6/28,13:00～16:00,3h,ゼミ課題,
,6/29,13:00～17:00,4h,ゼミ課題,
,6月合計,,40h(48h),,
,7/4,15:00～18:00,3h,ゼミ課題,
,7/5,13:00～15:30,2.5h,ゼミ課題,
,7/10,13:00～17:00,4h,ゼミ課題,
,7/11,15:00～18:00,3h,ゼミ課題,
,7/12,13:00～19:30,6.5h,ゼミ課題,
,7/13,13:00～16:30,3.5h,ゼミ課題,
,7/19,13:00～17:00,4h,ゼミ課題,
,7/20,13:30～15:30,2h,ゼミ課題,
,7月合計,,28.5h(76.5h),,
,8/2,13:00～17:30,4.5h,中間,
,8/17,14:30-17:30,3h,中間,
,8/20,15:30-18:00,2.5h,中間,
,8/23,14:00～19:30,5.5h,中間（エラーなおらない）,
,8/28,14:00～19:30,5.5h,中間,
,8/31,15:00～18:30,3.3h,中間（エラー直った）,
,8月合計,,24.5h(101),,
,9/4,14:00～20:30,6.5h,中間,
,9/6,13:00～18:00,5h,中間と実験手伝い,
,9/7,13:00～16:00,3h,中間,
,9/12,14:30～18:30,4h,中間,
,9/13,14:00～18:00,4h,中間,
,9/14,13:00～17:30,4.5h,vi復習,
,9/15,13:00～19:00,4.5h,データ編集,
,9/16,14:30～20:30,4h,概要書く,
,9/19,15:00～21:00,6h,概要書く,
,9/20,11:00～19:00,8h,概要、スライド,
,9/21,15:00～18:30,3.3h,概要、スライド,
,9/26,14:00～15:30,1.5h,論文読む,
,9月合計,,54.5h(155.5h),,
,10/10,13:00~16:30,3.5h,接触解析練習,
,10/11,15:00~19:00,4h,接触解析練習,
,10/12,13:30~18:00,4.5ｈ,接触解析練習と論文調べる,
,10/16,13:00~18:30,5.5h,接触解析,
,10/18,13:00~18:30,5.5h,接触解析（エラー出まくる）と摩擦...
,10/19,13:00~16:00,3h,接触解析と摩擦練習,
,10/23,13:00~17:30,4.5h,接触解析,
,10/24,13:30~17:30,4h,摩擦解析練習,
,10/25,14:30~18:30,4h,摩擦解析練習,
,10/26,14:00~19:30,5.5h,摩擦解析練習,
,10/31,13:30~17:00,3.3h,バネについて検証,
,10月合計,,47.5h(203h),,
,11/1,14:30~18:30,4h,バネについて,
,11/2,13:00~17:00,4h,バネについて,
,11/7,15:00~18:00,3h,バネについて,
,11/8,13:00~17:00,4h,バネについて,
,11/9,13:00~16:30,3.5h,バネ,
,11/13,13:00~17:30,4.5h,摩擦,
,11/14,13:00~16:30,3.5h,摩擦,
,11/15,13:30~18:30,5h,摩擦,
,11/16,11:30~15:30,4h,摩擦,
,11/17,16:00~18:30,2.5h,二面摩擦,
,11/20,13:00~18:00,5h,二面摩擦,
,11/21,13:00~15:00,2h,二面摩擦（収束しない）,
,11/22,13:00~18:00,5h,二面摩擦,
,11/23,17:00~21:30,4.5h,二面摩擦,
,11/27,17:00~18:00,1h,一面摩擦(計算速度遅い),
,11/28,13:00~17:00,4h,まだ計算終わらない,
,11/29,14:30~,3h,１面摩擦と２面摩擦,
,11/30,13:00~20:30,7.5h,1面摩擦と2面摩擦,
,11月合計,,86.5h(289.5h),,
,12/5,14:00~,4h,salome,
,12/6,13:00~,4.5h,salome（収束しない＆結果安定せず）,
,12/7,13:00~,8h,salome,
,12/12,13:00~16:30,3.5h,salome(ボルト付き),
,12/13,13:00~16:00,3h,salome,
,12/14,13:00~17:00,4h,salome,
,12/16,22:30~25:00,2.5h,salome,
,12/19,14:30~16:30,2h,salome,
,12/21,14:00~,4h,salomeと概要,
,12/24,15:00~,5h,概要,
,12/25,16:30~24:00,7.5h,スライド,
,12月合計,,48h(337.5h),,
,1/6,17:00~18:00,1h,salome,
,1/9,13:00~17:00,4h,salome(ミーゼス応力とか),
,1/10,11:30~14:00,2.5h,salome,
,1/11,17:00~20:00,3h,salome,
,1/12,14:00~15:30,1.5h,salome,
,1/15,15:00~18:00,3h,salome,
,1/16,13:00~17:00,4h,salome,
,1/17,15:00~18:30,3.5h,支部概要,
,1/18,15:00~,3h,支部概要,
,1/19,12:00~,3h,支部概要,
,1/22,13:00~18:30,5.5h,支部概要,
,1/23,13:30~18:30,5h,支部概要,
,1/24,13:00~,2h,支部概要,
,1/25,14:30~,2h,スライド,
,1/29,13:30~,3h,salome,
,1/30,13:00~,4h,salomeとスライド,
,1/31,16:00~18:30,2.5h,スライド,
,1月合計,,51.5h(389h),,
,2/1,14:00~19:00,5h,スライド,
,2/2,13:00~,2.5h,スライド,
,2/5,13:00~,4.5h,スライド・salome,
,2/6,12:00~,4h,発表練習,
,2/7,13:00~,3.5h,スライド,
,2/8,14:30~19:00,4.5h,スライド,
,2/9,13:00~16:30,3.5h,スライド,
,2/13,14:00~18:30,4.5h,スライド,
,2/14,14:00~22:00,5h,スライド,
,2/15,13:00~19:00,5h,スライド,
,2月合計,,42h(431h),,
*5/8課題結果 [#y6602ec0]
10mm×10mm断面の角材（ヤング率：6GPa, ポアソン比：0.3）の...
荷重は面荷重で100N.
スパン長さは、100mmを分担。
理論値は初等梁＝0.4167mm、ティモシェンコ梁＝0.4297mm
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
<一次要素>
,length,8,4,2,1
,たわみ(mm),0.3549,0.3549,0.4177,0.4177
,相対誤差(理論値：初等梁),14.83%,14.83%,1.20%,1.20%
,相対誤差(理論値：ティモシェンコ梁),17.41%,17.41%,4.19%,4...
<二次要素>
,length,8,4,2,1
,たわみ(mm),0.4338,0.4338,0.4338,0.4338
,相対誤差(理論値：初等梁),4.10%,4.10%,4.10%,4.10%
,相対誤差(理論値：ティモシェンコ梁),0.95%,0.95%,0.95%,0.95%
*5/15課題結果 [#y139767b]
木材（異方性）についてサロメでたわみを見た。
スパンは100mm。荷重は線荷重で100N。
理論値は初等梁＝0.4167mm、ティモシェンコ梁＝0.4917mm
<一次要素>
,length,8,4,2,1
,たわみ(mm),0.4327,0.4699,0.5002,0.5361
,相対誤差(理論値：初等梁),3.84%,12.77%,20.04%,28.65%
,相対誤差(理論値：ティモシェンコ梁),12.00%,4.43%,1.735,9....
<二次要素>
,length,8,4,2,1
,たわみ(mm),0.5378,0.5485,0.5555,0.5711
,相対誤差(理論値：初等梁),29.06%,31.63%,33.31%,37.05%
,相対誤差(理論値：ティモシェンコ梁),9.38%,1.16%,1.30%,1.61%
*5/22課題結果 [#m8109b73]
サンドイッチ梁（縦に1mm幅の鋼材を入れた梁）を作ってたわ...
木材のヤング率等は前回の課題と同じ。鋼材はE=206GPa、ν=0....
また、レンクスは１、メッシュは2次要素で切った。
たわみの理論値は0.09398mmとなった。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
,,たわみ(mm),相対誤差
,木材を等方性とした時,0.1119,19.07%
,木材を異方性とした時,0.1708,81.74%
<木材を等方性としたとき>（E=6GPa、ν=0.4）
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
<木材を異方性としたとき>
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
誤差が大きくなってしまったのでせん断係数も考慮した式を使...
*5/29課題結果 [#ub253aa2]
海老さんから卒論のデータをもらい、同じようにサロメで結果...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
同じように結果を出すことができた。Eficasのいじり方がよく...
*6/5課題結果 [#y1f6b79b]
前回の課題のデータを使い、サロメと手計算それぞれで求めた...
サロメでひずみが出せなかったので、応力からひずみを出した。
,,手計算,サロメ
,応力(N/m^2),3.6393 10^6(木),5.8131 10^6
,,97.438 10^6(鋼),
,ひずみ,4.73 10^-4,7.56 10^-4(木)
,,,0.282 10-4(鋼)
*6/12課題結果 [#z5e80853]
6/22のサンドイッチ梁において、サロメと手計算それぞれで求...
,,手計算,サロメ
,ひずみ,5.639 10^-4,5.2550 10^-4
,応力(N/mm^2),116.1634(鋼)/3.3834(木),115.5825(鋼)/3.0491...
*6/19課題結果 [#j5f53223]
前回の課題でのひずみの分布の仕方のグラフ。5.2～5.8mmは鋼...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
それから、サロメで求めたひずみにヤング率をかけて、サロメ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
*6/26課題結果 [#f90cdb70]
前回と同じサンドイッチ梁の木材部分において、
&link(この行列,http://www.str.ce.akita-u.ac.jp/~gotou/kou...
,,行列によって求めた応力,サロメで求めた応力
,0mm,2.8686,2.8688
,2.5mm,3.0521,3.0661
,8mm,3.0970,3.0967
小数点以下2桁あたりまで一致した。
*7/3課題結果 [#c72fb25b]
近藤さんからもらったデータを使い、プレストレス木箱桁橋の...
,,荷重(MN/mm^2),ひずみ(10^-4)
,対傾構あり積雪あり,0.0136759,4.630721
,対傾構あり積雪なし,0.0076566,2.592559
,対傾構なし積雪あり,0.0135336,4.589838
,対傾構なし積雪なし,0.0075144,2.548455
対傾構あり積雪ありの時のひずみの様子。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
対傾構あり積雪なしの時のひずみの様子。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
対傾構なし積雪ありの時のひずみの様子。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
対傾構なし積雪なしの時のひずみの様子。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
これらのひずみと荷重の関係を表にしてみた。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
対傾構ありのほうが対傾構なしに比べてややひずみが大きかった
*中間発表 [#mc9f109b]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
のエラーがなおらなくて半月以上たってしまったので、試しに...
（固定しすぎていたら3方向固定から2方向固定にするorバネは...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
みたいになった。ネジ穴周りは
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
穴周りで色が変わるかと思ったけど、この場合は一色であった...
やっとエラーが直った。メッシュの切り方を変えたら上手く行...
-添接板のネジ穴2列
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-添接板のネジ穴3列
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
ネジ穴2列と3列のときの添接板のネジ穴周りの応力を比べた。...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
穴3列のほうが応力が大きかった。
$\sigma$=$E$$\epsilon$ より穴6つの方がひずみも大きくなる...
*10/6ゼミ [#w4b38148]
次回までに
-摩擦を考慮したモデルで解析する
-道路橋示方書にあった論文を探して読む
今週できたことなど
-まず本を読んで接触解析を試しにやってみた(下図)。本に書い...
-論文は取り寄せすることになり、早くて来週届くそうだ。
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
*10/13ゼミ [#mc34d049]
次回までに
-自分がやってみたいモデルで接触解析をやる。必要があれば海...
今週できたことなど
-自分がやってみたいモデルに接触解析を行ってみたが、エフィ...
-摩擦の解析の練習も少し手を付けたけど、エフィカスのエラー...
-文献は無事に取り寄せできた
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-アラームはあるが結果は見ることができた。しかし、バネによ...
-荷重を二段階に分けるのは上手く行ったがやはりバネ定数を小...
*10/20ゼミ [#k0d3346a]
次回までに
-去年の先輩が解いた摩擦解析の問題をもらって解く
-余裕があったら文献読む
今週出来たことなど
-去年の4年の菊池さんの摩擦のモデルを解いてみた
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
菊池さんのデータを応用して柔らかいゴムではなく弱いバネを...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
*10/27ゼミ [#q6b74772]
次回までに
-論文少しでも読む！！！！
-高校の物理と対応した摩擦解析の問題がほしい。今度は壁から...
-http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin=...
今週出来たことなど
-前回のモデルにおいて、時間の分割数を100くらいにして動き...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-下図の上の物体の面の真ん中（線2本が交わっている所）にバ...
,バネ(x/y/z),変位(引っ張った方向),最大のひずみ,最大の応力
,0/0/0,1.439,15.6471,490
,10/10/10,1.200,15.6331,500
,10/0/0,1.1998,15.6334,490
,50/50/50,1.199,15.6382,490
,50/0/0,1.199,15.6394,490
,100/100/100,1.200,15.6446,490
,100/0/0,1.199,15.6471,490
今週できたことなど
-先週のモデルで間違いなどあったのでまたバネ定数などを変え...
,バネ定数(N/mm),静止摩擦係数,動き出す時の荷重(N)(手計算),...
,10,0.4,40,0.5,1,4.8077
,100,0.4,40,0.5,0.1,0.4981
,10,1,動かない,0.5,動かない,4.8077
,バネ無し,0.4,40,0.5,？,125
手計算とsalomeでの結果が大きく異なる結果となった。また、s...
*11/10ゼミ [#va29bacb]
次回までに
-摩擦係数の入れ方の確認を早急にやる。最初から滑り出さない...
11/13 菊地さんのデータの時間を1000分割してみが、やはり...
11/15 ついに摩擦を取り入れた計算が上手く行った。マスタ...
菊地さんのK2のデータのスレーブ面とマスター面を逆にすると...
,静止摩擦係数,動き出す時の荷重(N)(手計算),動き出した時の...
,0.4,40,40,149.494
,0.6,60,60,99.4478
*11/17ゼミ [#kfa8326d]
次回までに
-軸力を示方書で確認
-添接板での一面摩擦と二面摩擦
-ボルトをつけたモデルでやってみる
今週出来たことなど
-添接板に軸力が均等にかかっている場合の滑り出し方の比較の...
--一面摩擦 添接板二枚の上面に165kN×2(ボルトの軸力)をかけ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
--2面摩擦 こちらも真ん中の母材が飛び出してしまった（Salom...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
滑り出し方をグラフにしてみた。左のグラフが滑りだす瞬間ま...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-ボルトをつけたモデルは収束ができなくて結果が出なかった
*11/24ゼミ [#k4d6f0f2]
次回までに
-添接板に軸力が理想的（均等）に働いた時の滑り出し方の、一...
-ボルト付きの解析を引き続きやる
今週できたことなど
-１面摩擦を添接板１枚分で行った。
--境界条件など・・・軸力165kNを面載荷してから横へ引っ張っ...
--結果・・・引っ張る力がサロメでは66834Nの時動いた。手計...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-2面摩擦を行った。
--境界条件など・・・上下に軸力をかけて引っ張った。境界条...
--結果・・・引っ張る力が56880Nの時に滑りだした。前回やっ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-ボルト付きでもやってみた。計算は上手く行ったが全く滑り出...
*12/1ゼミ [#m0bd456b]
次回までに
-引き続き解析する
今週できたことなど
-引っ張る力を100分割にした一面摩擦の解析が上手くいった
-①一面摩擦、軸力×１のとき
--手計算・・・引っ張る力が165kN×0.4＝66000Nのとき滑り出す
--salome・・・引っ張る力が66834Nのとき滑り出す。手計算と...
-②一面摩擦、軸力×２のとき
--手計算・・・引っ張る力が165kN×2×0.4＝132000Nのとき滑り...
--salome・・・引っ張る力が136867.5Nのとき滑りだす。手計算...
-二面摩擦を行った。引っ張る力が135090Nのときに滑りだした...
-この二面摩擦のモデルにボルトをつけたモデルも解析した。な...
*12/8ゼミ [#x30c4047]
次回までに
-ボルト付き出やってみる
-概要書き始める
-論文を読みなおす
今週できたことなど
-収束判定を厳しめだと収束しない。あとボルト穴まわりに荷重...
-概要はちょっと手を付けた
*12/15ゼミ [#c30d4fe4]
次回までにと今後やること
-解ける所まででよい
-年明けに穴2つ3つの場合と全体をやってみる。圧縮側は穴を潰す
今週できたことなど
-ボルトをつけたモデルで滑り出しは見られなかったが、引っ張...
*中間発表後 [#o90f23f9]
次回までに
-母材のミーゼス応力の分布を見る
-母材が滑った後のモデル（母材とボルトが接触している）を作...
今週できたことなど
-母材が滑った後のモデルの様子
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-ミーゼス応力の分布の仕方
--滑りだす前、ボルトを締めた直後
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
--滑りだす前、母材を203MPaで引っ張ったとき
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
--滑りだした後のモデル（ボルトを締めた直後）
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
ボルトを締めただけで鋼材が降伏しそうだ。同じような条件で...
-リラクセーションを考慮した軸力をかけたのもやってみたが降...
*1/26ゼミ [#aa9a6483]
次回までに
-ボルト3列のモデルの解析
今週できたことなど
-フランジ部分の半解析(ボルト3本)をした。左から座金（最大6...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2017/kawamura/...
-スライドに画像並べてて違いがあまりないと感じたのでグラデ...
*支部用タイトル [#rfca77ef]
「プレストレス木箱桁橋のボルト接合部の数値モデル化」
*卒論タイトル [#ycae6e6c]
「プレストレス木箱桁橋のボルト接合部の挙動」
*ボルトメモ [#k989eb55]
-&link(座金の硬度,http://www.nssmc.com/product/constructi...
-&link(鋼の「硬度とヤング率」、「硬度と強度」の関,http://...
-座金の硬度は35$\sim$45HRC
-F10Tの硬度は27$\sim$38HRC
-座金の引張強さの推定：$引張強さ=3.2\times HRC\times 9.8$...
-座金の降伏応力の推定：$降伏応力=0.5\times 引張強さ$(MPa)...
*salomeメモ [#yfb9da81]
-RESI_GLOB_RELAの値が小さいほうが手計算とあっている気がす...
-参考にしていたサイト
--http://salome-meca.cocolog-nifty.com/blog/
--https://sites.google.com/site/codeastersalomemeca/
-11月らへんに計算した2面摩擦モデルでcommファイルの覚え書き
--(
DEBUT();
kou=DEFI_MATERIAU(ELAS=_F(E=206000.0, #ヤング率（単位：MP...
NU=0.3,),); #ポアソン比
MIAL=LIRE_MAILLAGE(FORMAT='MED',);
MIAL=MODI_MAILLAGE(reuse =MIAL,
MAILLAGE=MIAL,
ORIE_PEAU_3D=_F(GROUP_MA=('men2','men1...
newmesh=CREA_MAILLAGE(MAILLAGE=MIAL, ＃バネについて
CREA_POI1=_F(NOM_GROUP_MA='bane',
GROUP_NO=('bane1','ban...
MODE=AFFE_MODELE(MAILLAGE=newmesh, ＃「MIAL」も選択でき...
AFFE=(_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='3D',),
_F(GROUP_MA='bane', ＃CREA_MAILLA...
PHENOMENE='MECANIQUE',
MODELISATION='DIS_T',),),);
MATE=AFFE_MATERIAU(MAILLAGE=newmesh, ＃材料の種類。材料...
AFFE=_F(TOUT='OUI',
MATER=kou,),);
sktdn=AFFE_CARA_ELEM(MODELE=MODE,
DISCRET=_F(CARA='K_T_D_N',
GROUP_MA='bane',
VALE=(0.1,0.1,0.1,),),);...
contact=DEFI_CONTACT(MODELE=MODE,
FORMULATION='DISCRETE',
FROTTEMENT='COULOMB', ＃摩擦係数を...
ITER_GEOM_MAXI=100000, ＃
ZONE=(
_F(GROUP_MA_MAIT='men2', ＃マスター...
GROUP_MA_ESCL='men1', ＃スレーブ...
ALGO_CONT='PENALISATION', ＃ペナ...
E_N=206000.0, ＃ヤング率
COULOMB=0.4, ＃摩擦係数
ALGO_FROT='PENALISATION',
E_T=20600.0,), ＃ヤング率の10分...
_F(GROUP_MA_MAIT='men3',
GROUP_MA_ESCL='men4',
ALGO_CONT='PENALISATION',
E_N=206000.0,
COULOMB=0.4,
ALGO_FROT='PENALISATION',
E_T=20600.0,),),);
kotei=AFFE_CHAR_MECA(MODELE=MODE, ＃固定面。＃載荷、固定...
DDL_IMPO=_F(GROUP_MA=('kotei1','kote...
DX=0.0, ＃２方向固定に...
DZ=0.0,),);
osae=AFFE_CHAR_MECA(MODELE=MODE,
FORCE_FACE=(_F(GROUP_MA='joumen',
FY=20.24025,), ＃面載...
_F(GROUP_MA='kotei',
FY=-20.24025,),),);
tikara=AFFE_CHAR_MECA(MODELE=MODE,
FORCE_FACE=_F(GROUP_MA='saika', ＃...
FZ=200.0,),);
tikara_f=DEFI_FONCTION( ＃折れ線グラフの&link(イメ...
NOM_PARA='INST',
VALE=(0.0 ,0.0 , ＃VALE=(ｘ座...
10.0 ,0.0 ,
20.0 ,1.0 ,),); ＃10~20ｓで...
osae_f=DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0.0 ,0.0 ,
10.0 ,1.0 ,
20.0 ,1.0 ,),); ＃0~10ｓで徐...
kankaku=DEFI_LIST_REEL(DEBUT=0.0, ＃時間の分...
INTERVALLE=(_F(JUSQU_A=10.0,
NOMBRE=1,), ＃0~10...
_F(JUSQU_A=20.0,
NOMBRE=20,),),); ...
resu=STAT_NON_LINE(MODELE=MODE,
CHAM_MATER=MATE,
CARA_ELEM=sktdn,
EXCIT=(_F(CHARGE=kotei,), ＃はじめか...
_F(CHARGE=osae,
FONC_MULT=osae_f,),
_F(CHARGE=tikara,
FONC_MULT=tikara_f,),),
CONTACT=contact,
COMPORTEMENT=_F(RELATION='ELAS',),
INCREMENT=_F(LIST_INST=kankaku,),
NEWTON=_F(REAC_ITER=1,),
CONVERGENCE=_F(RESI_GLOB_RELA=0.01, ...
ITER_GLOB_MAXI=300,), ...
SOLVEUR=_F(SYME='OUI',),);
resu=CALC_CHAMP(reuse =resu,
MODELE=MODE,
CHAM_MATER=MATE,
RESULTAT=resu,
CONTRAINTE=('SIGM_ELNO','SIGM_NOEU',), ＃...
DEFORMATION=('EPSI_ELNO','EPSI_NOEU',), ...
CRITERES=('SIEQ_ELNO','SIEQ_NOEU',),); ...
IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=newmesh,
RESULTAT=resu,
NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','EPSI...
FIN();
#CHECKSUM:b2f4a2740b714f58c842db0021057ff8 -:FIN CHECKSUM
--)
-バネ付きモデルのときのcommファイルの覚え書き
--(
・
・
・
contact=DEFI_CONTACT(MODELE=MODE,
FORMULATION='CONTINUE',
FROTTEMENT='COULOMB',
ALGO_RESO_GEOM='POINT_FIXE',
REAC_GEOM='AUTOMATIQUE',
ITER_GEOM_MAXI=100,
RESI_GEOM=0.1,
ALGO_RESO_CONT='POINT_FIXE',
ITER_CONT_MAXI=100,
ALGO_RESO_FROT='NEWTON', #「ALGO_RE...
RESI_FROT=0.5,
ZONE=(
_F(GROUP_MA_MAIT='men1',
GROUP_MA_ESCL='boruto_m',
CONTACT_INIT='INTERPENETRE',
ALGO_CONT='PENALISATION',
COEF_PENA_CONT=206005.0,
COULOMB=0.4,
SEUIL_INIT=0.1, #これを入れる...
ALGO_FROT='PENALISATION',
COEF_PENA_FROT=20600.0,),
--)
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
