皆川の卒論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
[[FrontPage]]
#contents
*卒論 [#x5b6f34a]
**テーマ [#xf3b401f]
ケーブルの腐食やクラックを考慮した固有振動解析（皆川）
昨年度の山本の研究の続き．ケーブルの腐食方法やimpact...
昨年できなかった130mくらいの長いケーブルを要素分割し...
**参考資料 [#m1131b7d]
・山本さんの卒論
ケーブル腐食をくり抜くことで再現~
腐食の長さ・深さ・位置を変えて解析
・角田さんの修論
健全モデルと腐食モデルを使って、ケーブル破断時の連鎖崩壊...
C1~C14をそれぞれ破断させた時に崩壊するかの解析~
ケーブル健全時・腐食時の動的増幅率(DAF)と荷重再分配率(LRR...
*ケーブルの固有振動数の理論値 [#pdf8bc9e]
**長さ71m 直径0.315m [#ndfe1ce2]
・面固定~
1次モード 2.92969Hz~
3次モード 8.59375Hz~
5次モード 13.9648Hz~
7次モード 18.8477Hz~
・線固定~
1次モード 2.83203Hz~
3次モード 8.49609Hz~
5次モード 13.7695Hz~
7次モード 18.5547Hz~
*2/20〜 [#h56ec792]
**2/20 1次モードのsin波 Scale1.0 [#s868fbc5]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**2/28 ケーブルのみモデル 1次モード 500Nのsin波(1.66016Hz...
Scale1.0~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
中央の鉛直方向変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
*11/27〜 [#o4e79b21]
**11/27 引張実験 [#x9106e6a]
腐食レベル1 1-2 最大引張荷重64.25kN 応力1669.5MPa~
腐食レベル2 2-2 最大引張荷重64.00kN 応力1663.0MPa~
腐食レベル3 3-2 最大引張荷重63.55kN 応力1651.3MPa~
**IISさんからのアドバイス [#e874b20f]
①重力加速度を入れてみる~
②加速度と張力の関係を見てみる~
③ノッチ位置の変更~
**中間発表での質問 [#d9f25442]
①なぜ引張解析の結果で半径分の深さのノッチを入れると引張荷...
②どれくらいノッチが入った時に見つけられればよいのか(深く...
③新しい腐食発見方法はどんなものがありそうか~
**重力加速度をいれた振動解析 [#r95bd45a]
***健全モデル 解析時間50s [#a92626fb]
これまでと同じように解析時間を30s、インパクトを15sで与え...
解析時間を50s、インパクトを30sで与える時間設定にした。~
解析結果は下のようになった。卓越振動数は1次モードのみであ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**加速度の比較 100MN 重力加速度なし [#hfcbf938]
健全モデル、深さ157.5mmモデル、深さ45.00mmモデルで比較し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**30MN 重力加速度なし [#vf7662cc]
健全モデルとクラックモデルの加速度波形を比較した。重力加...
***加速度の比較 [#hd648f33]
クラックの深さでは加速度は変わらなかったが、健全とは異な...
・ノッチ深さでの比較~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
・各測定点での比較~
健全~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
ノッチあり~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
***PSDの比較 [#cdf20b98]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**ケーブルのみ 振動解析 重力加速度あり str30MN 1/22〜...
***時間設定 [#p93958f1]
0〜5秒:段階的に張力負荷~
5秒〜:重力加速度を入れる~
ケーブルを安定させる~
10〜14.99秒:ケーブルを安定させる~
15秒:インパクト~
15〜30秒:速度と加速度を測定する~
***解析結果 振動数 [#e10c2075]
健全 1次:1.66016Hz 3次:4.98047Hz 5次:8.49609Hz~
Crack15.75 1次:1.66016Hz 3次:4.98047Hz 5次:8.49609Hz~
Crack157.5 1次:1.66016Hz 3次:4.98047Hz 5次:8.49609Hz~
やはりノッチの有無・深さで振動数は変化なかった。~
***解析結果 加速度 [#o9d00940]
測定点3点ごとの応答加速度~
だいたい同じような結果が出た。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
***解析結果 インパクト後の振動の様子 [#v22ad56e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**カバーありモデル 振動解析 重力加速度あり str30MN 1/2...
***材料と時間設定 [#h993758e]
材料~
ケーブル:E=195000,nu=0.3,ρ=7.8☓10^(-9)~
仮想材料:E=50,nu=0.3,ρ=1.0☓10^(-11)~
カバー:E=1000,nu=0.3,ρ=9.0☓^(-10)~
時間設定~
0〜5秒:段階的に張力負荷~
5秒〜:重力加速度を入れる~
ケーブルを安定させる~
10〜14.99秒:ケーブルを安定させる~
15秒:インパクト~
15〜30秒:速度と加速度を測定する~
***卓越振動数 [#v04da4a2]
健全 １次:1.56Hz ３次:4.88Hz ５次:8.40Hz~
Crack15.75 １次:1.56Hz ３次:4.88Hz ５次:8.40Hz~
Crack157.5 １次:1.56Hz ３次:4.88Hz ５次:8.40Hz~
***応答加速度 [#ab11b7c2]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
***インパクト後の振動の様子 [#z266255b]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**参考 腐食モデル W69000 D15750 振動解析 重力加速度あ...
str30MN~
ノッチの深さでは卓越振動数などが変化していなかったので、...
***卓越振動数 [#m0fd5707]
1次:2.05078Hz ３次:6.15234Hz ５次:10.2539Hz~
健全モデルとの相対誤差 １次:23.5% ３次:23.5% ５次20.7%~
***応答加速度 [#dc2c71e3]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**1/31 応答加速度が変化した原因の検討 [#gadcee86]
***張力が変化したことで応答加速度が変化した? [#ndafb0df]
ノッチが深く入ることで応力集中しやすくなることは引張解析...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
ノッチが深いと変位が大きくなっていたが、これだけだと張力...
*9/20〜 [#w87deece]
**幾何学非線形を考慮した振動解析 [#v7534fc8]
振動解析を幾何学非線形を考慮して解析した。~
理論値との誤差が非常に小さくなった。(理論値 1次:2.837Hz ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
卓越振動数である2.83203Hzのsin波を与えた時の振動~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
卓越振動数である8.49609Hzのsin波を与えた時の振動~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
卓越振動数である13.7695Hzのsin波を与えた時の振動~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**クラック 幾何学非線形性 [#la89c235]
***長さ71m 半径157.5mm [#j80656c1]
①クラック深さ157.5mm~
1次 2.83203Hz 3次 8.49609Hz 5次 13.7695Hz~
②クラック深さ146.25mm~
1次 2.83203Hz 3次 8.49609Hz 5次 13.7695Hz~
③健全~
1次 2.83203Hz 3次 8.49609Hz 5次 13.7695Hz~
クラックでは固有振動数に変化はなかった。~
昨年の山本さんの腐食モデルでは振動数が変化していたことを...
**腐食レベル2 [#p9c7275f]
腐食レベル2のケーブルは下のようにサビが入っていた。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
そこでこのケーブルのメッキ部分を削り、どの程度腐食が進展...
***引張解析 [#t096d3e7]
・300mm 7mmモデル
まずはサビがある部分のみを削り、サビの深さと幅をモデルで...
引張応力はσ=1164MPaとなり、健全時と同様の結果が得られ、サ...
・71000mm 315mmモデル
次に300mmモデルの腐食深さを45倍(315mm/7mm)、幅を236.7倍(7...
引張応力はσ=1146MPaとなり、腐食による影響は見られなかった...
**6角柱モデル 2本 [#h16a7ac1]
下の図形のように面接続と線接続で同じ結果が得られるかを確...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・振動解析~
,モード,面接続 結果,線接続 結果,相対誤差%
,1次,3.125Hz,3.125Hz,0.29%
,3次,9.17969Hz,9.17969Hz,1.81%
,5次,14.7461Hz,14.7461Hz,5.36%
・引張解析~
面接続 引張応力σ=1152MPa、point71000での応力=1156MPa~
線接続 引張応力σ=1152MPa、point71000での応力=1156.1MPa~
振動解析、引張解析ともに面接続と線接続では同じ結果が得ら...
**6角柱モデル 面結合[#me8ea140]
実際のケーブルに近づけるために下のように6角柱を7本・19本...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***7本モデル [#if448901]
7本モデルでは半径60mmの円上に点を作って、底面の6角形を作...
・振動解析~
100MNの張力を与えての振動解析~
1次 3.125Hz 相対誤差1.06%~
3次 9.277Hz 相対誤差0.012%~
5次14.844Hz 相対誤差3.99%~
7本モデルの底面と同じ面積の円柱で振動解析を行ったところ、...
・引張解析~
引張応力σ=1152MPa point71000での応力=1156.8MPa~
***19本モデル [#seb4db28]
**6角柱モデル 線結合 [#s28eb61f]
・振動解析
線結合でも面結合と同じ結果になった。~
**腐食 振動解析 [#v6c834e9]
***腐食再現モデルの種類 [#dddba9f9]
幅と深さを下の図のようにする。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
まずは中央部に下にある、幅7種類と深さ6種類の計42種類の解...
・幅 ()内は長さに対する割合~
①710mm(1%) ②3550mm(5%) ③7100mm(10%) ④14200mm(20%) ⑤21300m...
・深さ ()内は半径に対する割合~
①15.75mm(10%) ②47.25mm(30%) ③78.75mm(50%) ④110.25mm(70%) ...
幅・深さ以外にも腐食の位置も固有振動数に影響してくると思...
また、ボリュームが同じであれば、えぐり方(腐食の入り方)で...
***結果 [#e6d401bb]
下のグラフは、横軸に健全時のボリュームに対する腐食モデル...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
次のグラフは、上のグラフの横軸を直径に対する腐食深さの割...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
次は腐食の位置が各振動モードの固有振動数にどのように影響...
①中央にW28400の腐食~
②両端から1000mmの位置にW14200ずつの腐食(合計でW28400の腐...
③左端から1000mmの位置にW28400の腐食~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**被覆ありモデル 振動解析 面固定[#rbede1f0]
***被覆の半径250mm 仮想材料の半径200mm [#p7eaf383]
ケーブルの半径はこれまで通り157.5mm、仮想材料の半径を200m...
材料諸元は~
ケーブル:E=195000MPa,nu=0.3,ρ=7.8☓10^(-9)ton/mm^3~
仮想材料:E=50MPa,nu=0.3,ρ=1☓10^(-11)ton/mm^3~
被覆(高密度ポリエチレン HDPE):E=1000MPa,nu=0.3,ρ=9☓^(-10)...
固有振動数の結果 ()はケーブルのみモデルとの相対誤差~
1次モード 2.73438Hz (6.67%)~
3次モード 8.20313Hz (4.55%)~
5次モード 13.3789Hz (4.20%)~
7次モード 18.0664Hz (4.15%)~
→ρが大きくなったことでケーブルのみの時よりも固有振動数が...
***被覆の半径190mm 仮想材料の半径170mm [#x22166de]
上の被覆ありモデルの被覆と仮想材料の半径のみを変えて解析...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
固有振動数の結果 ()は面固定のケーブルのみモデルとの相対...
1次モード 2.83203Hz (3.33%)~
3次モード 8.49609Hz (1.14%)~
5次モード 13.7695Hz (1.40%)~
7次モード 18.5547Hz (1.55%)~
**腐食 引張解析 [#m2b0cdda]
振動解析の時の腐食モデルをすべて同じ張力を与えた時のひず...
張力は500MPaを与えた。~
結果はえぐった部分の断面積(腐食深さ)でひずみが決まるよう...
*振動解析 [#n270fa13]
**2段階解析 [#occ4ae6c]
***1段階目に張力 2段階目にimpact [#g7002f2d]
1段階目:固定条件と張力を設定し、ASSEMBLAGEの境界条件に固...
2段階目:impactを設定し、ASSEMBLAGEの境界条件に固定条件とi...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***1段階目に張力 2段階目に固有振動解析(CALC_MODES) 由利...
今回の解析では由利橋の1/100モデルを使用した。~
1段階目:固定条件と張力を設定し、ASSEMBLAGEの境界条件に固...
2段階目:固有振動解析(CALC_MODES)~
y方向の振動数の結果は下の表のようになった。~
1次 14.00Hz~
2次 55.65Hz~
3次 124.82Hz~
4次 222.89Hz~
5次 346.32Hz~
***1段階目に張力 2段階目にimpact 由利橋の1/100モデル [#d...
今回の解析は由利橋の1/100モデルで行った。~
①1段階目・2段階目ともにDYNA_VIBRAを使い、2段階目はETAT_IN...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
②1段階目STAT_NONLINE、2段階目DYNA_NONLINEを使い解析をおこ...
結果は上のグラフと同じ振動数がでた。
**9/16 振動解析 PRES_REP [#af4e71ac]
指定面にPRESを与えるコマンドPRES_REPを使って解析を行った。~
SIZZは与えたい張力分が入っていた。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**9/12 振動解析 PRE_EPSI [#y21e4eb2]
PRE_EPSIについては3Dでも適用することができたが、σ=Eεで100...
PRE_EPSIでは正の値と負の値を両方を与えたが、どちらの解析...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**9/10 振動解析 PRE_SIGM [#q7a525ce]
今回の解析ではPRE_SIGMを使って解析した。(DYNA_VIBRA)~
PRE_SIGMについては、impactなしで解析したところ与えたい張...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**長さ1000mm、半径0.5mmの健全モデル [#c0a583be]
固定の仕方を線固定にした。(張力は面載荷で与えた。)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
固有振動数の理論値 f=6.37Hz 解析時間0〜30秒 細長比354...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
,時間,行うこと
,0〜5秒,張力(1N)をかける。
,5〜14.99秒,ケーブルを安定させる。
,15秒,10Nの力を与える。
,15.01〜30秒,卓越した振動を検出する。
**由利橋のケーブルモデル 張力0 [#jadf6985]
張力0のモデルにインパクトを与えて、その結果をFFT処理した...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**由利橋のケーブル 長さ71000mm、半径157.5mmの健全モデル...
固有振動数の理論値 f=2.85Hz 解析時間0〜100秒 細長比25...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
,時間,行うこと
,0〜5秒,張力(100MN)をかける。
,5〜14.99秒,ケーブルを安定させる。
,15秒,100kNの力を与える。
,15.01〜100秒,卓越した振動を検出する。
結果~
Mesh50 ボリューム169538~
中央部の振動数 ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
波形を書かせるとおかしな波形となり、モードの推定ができな...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***6/24 内部張力が、与えている張力と等しいかの解析 [#nc22...
中央部(point35500)と張力を与えた部分(point71000)の15秒〜1...
中央部(point35500)~
平均値99.6MN 最大値104.5MN 最小値98.7MN~
張力を与えた部分(point71000)~
平均値・最大値・最小値ともに95.6MN~
***インパクト1000kNに変更 [#ofefb2d1]
上の解析のインパクトを100kN→1000kNに変更して解析を行った。~
インパクトを与えた後(16.20s)の100kNの時のモードと1000kNの...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***6/21 張力とインパクトをいろいろ変えて解析 [#ldc16f87]
①張力100MN、インパクト100kN~
②張力100MN、インパクト1000kN~
③張力200MN、インパクト100kN~
④張力200MN、インパクト1000kN~
どの解析でも下のような同じ振動数が出力された。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**由利橋のケーブル(長さ71000mm、半径157.5mm)に振動を与え...
***両端に振動を与える解析 [#xaa17730]
両端に理論値である2.85Hzの振動を与える解析~
ParaVISでモードを見ると下のようになり、微細な振動は見られ...
両端を固定しており、その部分に振動を入れているからだと思...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***中心に振動を与える解析 [#s13d0a3e]
長さ71000mm、張力100MNの弦の理論値である2.85Hzを中心の点...
モードは以下のようになった。~
振動を与えた後~
16.41秒~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
16.58秒~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
17.94秒~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
中心に1.16Hzの振動を与える解析を行った。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
振動を与える直前(14.99秒)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
振動を与えた直後(15.5秒)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
中心に3.2Hzの振動を与える解析を行った。~ [#t2afbbe4]
振動を与える直前(14.99秒)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
振動を与えた後は下の2つの位相が繰り返された。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
中心に1.16Hzと3.2Hzの振動を与えたときの比較 [#j7e3c02b]
振動を与えた点(中心部)のy方向変位~
y方向の最大変位~
1.16Hz:11070.7mm 3.2Hz:4837.51mm~
y方向の最小変位~
1.16Hz:-10744.2mm 3.2Hz:-4925.62mm~
***両端から5mmの位置に振動を与える解析 [#ue6d276e]
振動を与える位置のモデル~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
両端から5mmの位置に0.122Hzの振動を与える解析を行った。~
振動モードは下のように、１次モードとなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
両端から5mmの位置に1.16Hzの振動を与える解析を行った。~[#b...
振動モードは下のように、３次モードとなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
両端から5mmの位置に3.2Hzの振動を与える解析を行った。~[#t5...
振動モードは下のように、５次モードとなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**7/30 張力の与え方を単位体積あたりの荷重に変えて解析 [#z...
張力の与え方は、今までの解析では面載荷(FORCE_FACE)で与え...
FFT処理をすると以下のような結果となった。(理論値は2.85Hz)~
張力を100MN→1000MNに変更して解析したが、下の結果と同じに...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**6/28 ビーム要素(温度荷重ではない) [#ne6248ee]
ビーム要素で由利橋健全ケーブル(長さ71000mm、半径157.5mm)...
梁の曲げ振動数の理論値と等しくなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**強制変位 [#n235dd68]
100MNの張力を与えたときに475mmz方向に変位していたので、今...
強制変位させることでも応力が入ることがわかったが、場所ご...
***インパクトを与える解析 [#db9e6d09]
これまで通りインパクトを与えて、それで得た結果をFFT処理し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***静的解析 [#t0385b50]
強制変位させた場合の応力の入り方が初期張力のように入って...
1次モード 0.13Hz~
2次モード 0.51Hz~
3次モード 1.15Hz~
4次モード 2.05Hz~
5次モード 3.20Hz~
張力が0の時とほぼ同じ結果となった。~
**両側引張 [#zd1f0a59]
これまでの解析では片側のみを引っ張って張力を与えていたが...
結果を見ると、片側引張では完全固定部分で張力が小さくなっ...
両側引張では全体的にほぼ均一な張力が入っていたが、結果は...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**インパクトの与え方をじっくり与えた [#j5bb89cb]
これまでのインパクトの与え方は15sで一気に与えていたが、今...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**解析モデルを由利橋の1/100にした(長さ710mm、直径3.15mm) ...
***ホールドあり [#b7e1b619]
振動後に面がz方向に垂直でないことで与えたい張力が与えられ...
右の表がインパクトを与えてMATLABで処理した時の結果(理論値...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***周波数応答解析 [#zdcabd3d]
CENTER:周波数応答解析 面固定
,モード数,結果,理論値,理論値との誤差,梁の曲げ振動,曲げ振...
,1次,21.9Hz,28.6Hz,23.4%,19.2Hz,14.1%
,2次,71.2Hz,57.1Hz,24.7%,62.1Hz,14.7%
,3次,147.8Hz,85.7Hz,72.5%,129.6Hz,14.0%
CENTER:周波数応答解析 線固定
,モード数,結果,理論値,理論値との誤差,梁の曲げ振動,曲げ振...
,1次,14.1Hz,28.6Hz,50.7%,19.2Hz,26.6%
,2次,56.2Hz,57.1Hz,1.6%,62.1Hz,9.5%
,3次,125.6Hz,85.7Hz,46.6%,129.6Hz,3.1%
**由利橋 長さ71000mm、半径157.5mmの中央腐食モデル(長さ71...
固有振動数の理論値 f=2.85Hz 解析時間0〜100秒~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
載荷線と腐食部分~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
*引張解析・引張実験 [#f72db340]
**引張解析 [#h726c276]
クラックの解析は専用の解析方法があるらしいので、それでお...
クラックは三角柱でくり抜いて表現した。(下の図はクラック角...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***健全モデル [#y73b55e6]
このモデルでは
***クラック角度 90° クラック深さ157.5[#je549823]
このモデルでは445MPaまで引っ張ることができた。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***9/27 長さ300mm 直径7mmモデル 角度一定 クラック幅0.4 ...
角度一定でクラックの深さを0.25mm刻みで引張解析を行った。...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
応力集中の様子が下の図(クラック深さ3.5mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
上の①〜⑤の解析のz方向変位と荷重の関係を表したものが下のグ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***10/7 長さ71000mm 直径157.5mmモデル 角度一定 クラック...
今回の引張解析は、9/27の引張解析のクラックと角度を同じに...
角度一定、深さを11.25mm刻みで解析した結果が下のグラフ・表...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**引張実験 [#u3aa1cad]
***8/28 健全の鋼棒・腐食3の鋼棒 [#g746217e]
①健全の鋼棒~
まずは8/22の実験と同様に、下のジグで掴む部分をネジ状に溝...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・ひずみゲージから求めたヤング率~
ひずみゲージNo1：E=197GPa~
ひずみゲージNo2：E=204GPa~
ひずみゲージNo3：E=203GPa~
ひずみゲージNo4：E=221GPa~
②健全の鋼棒~
次に上下のジグで掴む部分を平面にして、実験をおこなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・ひずみゲージから求めたヤング率~
ひずみゲージNo1：E=245GPa~
ひずみゲージNo2：E=198GPa~
ひずみゲージNo3：E=191GPa~
ひずみゲージNo4：E=220GPa~
③腐食レベル3~
最後に上下のつかみ部のメッキを剥がした腐食レベル3の鋼棒を...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・ひずみゲージから求めたヤング率~
ひずみゲージNo1：E=239GPa~
ひずみゲージNo2：E=226GPa~
ひずみゲージNo3：E=215GPa~
ひずみゲージNo4：E=191GPa~
***8/22 健全の鋼棒 [#od2d91a9]
今回の引張実験では、鋼棒の下のジグ部分をネジ状に溝を入れ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・ひずみゲージから求めたヤング率~
ひずみゲージNo1：E=203GPa~
ひずみゲージNo2：E=197GPa~
ひずみゲージNo3：E=197GPa~
ひずみゲージNo4：E=224GPa~
***7/16 腐食レベル2の鋼棒 [#d86d0401]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・ひずみから求めたヤング率~
ひずみゲージNo1：E=274GPa~
ひずみゲージNo2：E=233GPa~
ひずみゲージNo3：E=349GPa~
ひずみゲージNo4：E=269GPa~
・伸びから求めたヤング率~
E=305GPa~
何らかの原因で実験がうまくできていなかった。~
ジグの部分をうまく掴ませるように工夫したい。~
**論文まとめ [#ff52d2ef]
***振動法によるケーブル張力の実用算定式 [#gf30995b]
曲げ剛性やサグ、傾斜角の影響を考慮した張力算定式~
曲げ剛性の影響を表すパラメーターが小さい時(曲げ剛性が大き...
(サグや傾斜角の影響を表すパラメーター)＜3.0の時は、対称振...
***斜材点検用非破壊検査装置の開発と運用 自走式点検装置 [#...
斜ケーブルの変状調査を目的とした自走式点検装置の機能・硬...
装置の特徴~
①カメラを4台設置することで全周を同時に調査可能~
②保護管上のクラックケージを0.1mmまで判読可能~
③保護管内の鋼線のうち外周部分の破断の有無の調査可能~
④鋼線1本の破断から調査可能~
⑤地上のパソコンから遠隔操作可能~
⑥照明を装備しているため夜間の点検可能~
利点~
①点検者の近接目視と同じくらいの精度で点検可能~
②地上からの遠隔操作で安全に点検可能~
③高所作業車を使用する方法よりも安価~
④交通規制の規模も低減可能~
⑤保護管に変状があった場合、同時に斜ケーブルの変状の有無も...
**幾何学的非線形性の参考サイト [#o7f2f199]
https://www.cybernet.co.jp/ansys/learning/glossary/kikaga...
**進捗 [#f2f072c8]
***9/10〜青木班ゼミまで [#p8847e9f]
9/10 PRE_SIGMコマンドを使って解析した。静的解析だと張力が...
結果は1次0.097Hz、3次1.17Hz、5次3.22Hzだった。~
9/12 imp100kN(5sから一定)、str(5sから一定)の時間設定でDYN...
9/13 impact後の応力と変位について様々な点(y=0の点)をとっ...
9/14 DEFI_FOND_FISSでクラック解析をした。1個ずつエラーを...
9/15 PRE_EPSIで解析した。正の値と負の値を入れた解析の2つ...
9/16 圧力をかけるコマンドPRES_REPを張力面にかけて解析した...
9/17 先週の水曜日に青木さんが考えていただいたモデルで解析...
Liaisonで結合させたらうまくできた。Liaisonでは結合させた...
***8/27〜9/9 [#z8ca6aed]
8/27 振動解析において、張力を与える側がz方向に固定されて...
しかし、結果は0sから固定されており、張力→固定ができなかっ...
8/28 引張実験をおこなった。①ネジ＆メッキ剥がし②平面③腐食...
8/29 昨日の引張実験の結果の確認をした。1200MPa前後で降伏...
8/30 象潟のPCラーメン箱桁橋の現場見学に行った。~
9/2 振動解析では710mmモデルでFORCE_INTERNE→CALC_MODESをお...
引張解析ではクラックの導入の仕方を調べた。今まではクラッ...
XFEMはクラックの大きさ、形、始点を指定することで自動でク...
XFEMではクラック内に線要素があるというエラーがでた。(実際...
9/3 振動解析では強制変位+100MNで張力を与えた。結果は振動...
引き続き引張解析のクラックの導入の仕方を模索している。~
9/4 クラックの導入の仕方でDEFI_FOND_FISSを使う場合、3Dで...
INI_XFEM_ELNOに関するエラーが発生したがよくわからない。~
9/5 振動解析では両端固定＆FORCE_INTERNEの結果を確認した。...
9/6 振動解析において、張力側の面がz軸方向に垂直にならない...
結果は0.195Hz、1.37Hz、3.61Hz、6.84Hz・・・となった。~
9/9 振動解析について、時間で固定条件を変更して張力をかけ...
また、PRE_SIGMコマンドについて調べ、解析を回したところ、P...
***8/8〜青木班ゼミまで [#n0462cf2]
8/8 千代岡さんが見つけた2段階解析を少し変更して(千代岡さ...
今回の解析は1段階目に固定条件と張力を設定(loadとfunction)...
また、引張試験のジグの噛ませ方を調べたがまだ良い方法が見...
8/9 所々変更しながら2段階解析を行った。~
昨日行った2段階解析の結果をFFT処理した。結果は結果欄に書...
次に1段階目に張力をかけ、2段階目に固有振動解析(CALC_MODES...
今日解析をして、千代岡さんと話していて疑問に思っていたETA...
そのため私の解析では使えないかもしれない。(私の解析は1段...
8/10 1段階目の解析をSTAT_NONLINEに変えて解析を行った。~
また由利橋の1/100モデルでETAT_INIT内のINST_INITに0を入力...
現在のところで、ETAT_INITについてわかっていることとしては...
1段階目の解析の時間の途中から2段階目の解析を始めるときはE...
8/19 由利橋(71000mm)モデルでは容量が足りなくなってしまう...
また、引張解析について調べ、由利橋(71000mm)の健全モデルの...
8/20 由利橋の1/100モデルの解析結果を再確認した。DYNA_VIBR...
STAT_NONLINE(張力)→DYNA_NONLINE(impact)では1段階目の張力...
また、引張解析のCrackの与え方をいろいろ試してみた。三角柱...
8/21 引張解析についていろいろと試した。降伏応力とひずみの...
8/22 引張実験をおこなった。鋼棒を掴む部分を①ネジのように...
①ネジのように溝を入れる方法では、溝は機械を使ってピッチを...
③平たくする方法では0.3mm分削った。~
今日の引張実験では①②の鋼棒を使って実験をおこなった。~
ジグの部分は滑っていないようだった。(今回の実験では20kNま...
また、引張解析では90°クラックの解析を行った。引張応力は、...
また、引張解析では90°クラックの解析を行った。引張応力は、...
8/23 昨日の引張実験の結果をまとめた。途中までしか実験がで...
引張解析については結果を見るときにどの点で見ればよいかを...
引張解析のモデルは片持ち式でおこなっている。~
***8/1〜8/7 [#nd0852d8]
8/1 impactの位置を変えて解析したが、そもそも中心部以外にi...
結果は今までと変わらなかった。~
8/2 初期応力を与える良い方法がないか探した。PRE_SIGM(おそ...
PRE_EPSI(初期ひずみ)での解析を今度行ってみる。(初期ひずみ...
千代岡さんと2段階解析について確認した。まずは出力した結果...
8/3 単位体積あたりの荷重(FORCE_INTERNE)についてもう一度確...
z方向変位が入れたい張力分の1/2の変位しか入らなかった。~
8/7 PRE_EPSIについて解析を行った。結果は張力が入っていな...
PRE_EPSI(初期ひずみ)の与え方が違っていたかもしれないので...
***7/25〜7/31 [#f452332d]
7/25 及川さんからいただいたケーブル要素の例をもとに解析を...
モデルを作るところから頑張りたい。~
7/26 及川さんからいただいたケーブル要素のモデル(ケーブル4...
・ケーブル4本にするとMUMPSのエラーが発生する。(ten1〜4はX...
・プレート&ケーブル1本のモデルにするとエラーが発生する。~
・ケーブルに質点を与えた解析~
→張力が圧縮として入っていた。また、1次モードの振動数が縦...
7/28 振動した時に張力部分がz方向に垂直でないことで張力が...
Meshの設定の仕方に問題があるのかモデルが一体とならずに貫...
7/29 連成解析の方法を探した。参考資料は見つけることができ...
7/30 周波数応答解析を行ってみた。解析を少しでも早くするた...
結果は結果欄にある。~
7/31 FORCE_INTERNE(単位体積あたりの荷重)を使って張力を与...
FORCE_INTERNEを使った解析は下の3つを行ったが、全て同じ結...
・張力面をz方向フリーにして100MN分の単位体積荷重を与えた...
・張力面をz方向フリーにして1000MN分の単位体積荷重を与えた...
・両面を完全固定にして100MN分の単位体積荷重を与えた解析~
***7/18〜7/24 [#l26bd7fc]
7/18 英語の発表の準備を行った。~
7/19 オープンキャンパスで使う振り子の調整をした。~
7/21 温度荷重の解析をもう一度行ってみた。張力が入っている...
片側のみに応力を与えるのではなく、両側に応力を与える解析...
***7/11〜7/17 [#a9e78916]
7/12 引張試験の準備を行った。~
7/15 今まで行ってきた解析結果を確認して、今後の解析で張力...
・応力として線載荷で与えると張力がうまく入っていなかった...
・ビーム要素では梁の曲げ振動の値となった。(0.28Hz、1.50Hz...
・強制変位(475mm)では25〜50Hzに多く卓越していた。また、張...
・内部張力(FORCE_INTERNE)では0.51Hz、1.16Hz、3.70Hzとなっ...
・張力なしでは
7/16 引張試験を行った。使った鋼棒は腐食レベルが2-1のもの。~
破断させることはできなかったが、応力-ひずみのグラフを作成...
ジグの部分がうまくつかめていなかったので、どのようにした...
7/17 電子顕微鏡の使い方を学んだ。元素まで見ることができる...
昨日行った実験で改善すべき点であったジグの部分を修正した。~
及川さんと青木さんが、ジグが回らないようにすることとジグ...
***7/4〜7/10 [#j6070258]
7/4 ネットで見つけた資料を参考にしながら温度荷重の解析を...
7/8 前回(4日)の温度荷重の解析を修正しながら再び解析を行っ...
AnalysisにDYNA_VIBRAではなく、DYNA_NON_LINEを使うと解析が...
7/9 今まで行っていた温度荷重の解析では100℃上昇させていた...
この557.5℃は、今まで行っていた張力(100MN)を与える解析でpo...
7/10 強制変位させることで張力が入ることが分かったので、強...
***6/26〜7/3 [#i59cb769]
6/26 今週はケーブル要素、ビーム要素、温度荷重の解析を行っ...
千代岡さんから温度荷重やビーム要素について教わった。また...
及川さんいただいた、ケーブル要素について解析を回してみた...
6/27 昨日、千代岡さんから教えていただいたソリッド要素に内...
結果は上に書いているが、やはり理論値とは異なる値となった。~
ケーブル要素のエラーについても調べていきたい。~
7/1 青木さんがもらってきていただいたケーブルの引張試験の...
いただいたケーブルはCASE1(健全)とCASE5(削孔+塩水注水)~
7/2 ケーブル要素の解析がうまく回らなかった原因を考えたが...
やはり、「nddlリストにモデルがありません」というエラーが...
7/3 今後は温度荷重と初期ひずみについて進めていく。~
***6/20〜6/26 [#yffc04b9]
6/20 由利橋のケーブルモデルの両端から5mmの位置に0.122Hzを...
結果は1次モードとなった。~
6/21 中央部に張力がしっかりと入っているかを確認するため由...
また、張力を2倍の200MNに変えた解析と左右で逆側に振動を与...
今日の解析では張力を2倍にしても振動数が√2倍になっていなか...
インパクトは今まで100kNを与えていたがモードを見るとあまり...
6/23 張力0の解析を行い、理論値と比較した。~
結果は上に書いたが、張力を与えた場合よりも理論値が近くな...
6/24 モデル中央部に張力分だけの応力が入っているかを調べた...
また、今まで線固定にして解析を行っていたが、張力を与える...
明日は温度荷重について学びたいとおもう。~
***6/13〜6/19 [#w090b26f]
6/13 振動(sin波)を与える解析を及川さんに教えていただいた。
最初に1.16Hzのsin波を与えて解析を行った。与えたsin波がsin...
そこで50000sin(2π1.16t)として50kNの振動を与えた。結果は明...
6/14 50kNの振動を両端に与えて解析を回したがうまく結果が見...
そこで振動を両端ではなく、中心の点に1.16Hzの振動を与えて...
結果は上の(由利橋 健全モデルに1.16Hzの振動を与える解析)に...
6/15 昨日の解析を3.2Hzにして解析を行った。~
結果は上の(由利橋 健全モデルに3.2Hzの振動を与える解析)に...
1.16Hzの振動を与えた場合と3.2Hzの振動を与えた場合のy方向...
また、由利橋(長さ71000mm、半径157.5mm)の固有振動数の理論...
6/18 振動させる位置を両端から5mmの位置に変えて、1.16Hz,3....
1.16Hzでは３次モード、3.2Hzでは５次モードとなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***5/30~6/12 [#ceb33306]
5/30 真杜くんの実験のサポート~
5/31 実験サポート~
6/3 100mm,半径0.5mmの解析の結果をまとめようとしたが、パ...
6/4 パソコンが重く、結果が見れなかったため振動を検出する...
,時間,行うこと
,0〜5秒,張力(1N)をかける。
,5〜14.99秒,ケーブルを安定させる。
,15秒,10Nの力を与える。
,15.01〜30秒,卓越した振動を検出する。
6/5 上記の時間で解析したがパソコンがフリーズしてしまい結...
6/6 山本さんの解析では長さ50m・半径75mmで行っており、今...
そこで張力を3.2MN→1N、中央部への線載荷を100KN→10Nとして解...
結果を得ることができた。明日、青木さんにFFT処理をお願いす...
6/7 長さ1000mm、半径0.5mmの中央部に100mm腐食を与えたモデ...
長さ100mm、半径0.5mmの中央部に10mm腐食を与えたモデルでは...
細長比によるエラーかもしれない。~
長さ1000mm半径0.5mmモデルの細長比は3544.91~
長さ100mm半径0.5mmモデルの細長比は354.5~
6/8 由利橋のケーブル(長さ71000mm、半径157.5mm)の健全モデ...
月曜日に青木さんにFFT処理を依頼する。~
6/9 由利橋のケーブル(長さ71000mm、半径157.5mm)の中央腐食...
月曜日に青木さんにFFT処理を依頼する。~
6/10 FFT処理の行い方を学んだ。理論値とズレが生じたため、...
6/11 由利橋の健全モデルの解析値(3.2Hz)と理論値(2.85Hz)に...
最初の解析では、面載荷に張力(100MN)を円の断面積で割った値...
どちらの解析の結果もf=3.2Hzとなっていたため、他の考えられ...
***5/15〜5/29 [#w9e80a4f]
・山本さんの卒論の振動解析を解析してみる。~
山本さんの解析では測点を３点設定していたが、モードをより...
5/16 山本さんの卒論，角田さんの修論を読んだ．~
5/17 健全モデル、中央10m腐食、中央1m腐食のモデルを作った。~
5/20 健全モデルのAster Studyを設定し，解析を回してもらう．~
中央に線荷重を与える際は，3.2MNを断面積(半径75mmの円)で割...
5/21 昨日設定したAster Studyのコードの意味を調べた．~
5/23 健全モデルの解析を回した。~
50m 半径75mm 健全モデル(点は下の写真の上部に設置した。)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
50m 半径75mm 中央1m腐食
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**MATLAB [#bbd750c1]
1,変数→新規変数→処理したいデータを入れる~
2,FFTcable→データ名、スタート点(私の解析ではインパクト...
3,bandcable→下限値・上限値、データ名、スタート点の入力~
4,C1P3→csvファイルで波形を描く~
**Salome-meca [#ka48d877]
***Salome-mecaの使用方法 [#z8a41cac]
-[[Salomemeca使用解説法:http://opencae.gifu-nct.ac.jp/puk...
-[[code-asterの説明:https://code-aster.org/doc/v15/en/ind...
-[[code-aster質問ページ:https://forum.code-aster.org/publ...
***重力加速度の設定 [#h7a1bdbe]
材料設定で密度を設定する。~
AFFE_CHAR_MECA→PESANTEUR→GRAVITY=9800(mm/s^2)、DIRECTION=...
上の設定のDIRECTIONは、重力方向がy方向の時の設定。DIRECTI...
***引張解析 参考 [#i0ac39fd]
-[[Salome-Meca演習_弾塑性解析(2021)]]
***クラックの解析 破壊解析 [#h9b7cae4]
https://www.youtube.com/watch?v=5ivZy5rY9l8~
https://youtu.be/8fy0hmTITEQ~
https://www.opencae.or.jp/wp-content/uploads/2019/12/08-n...
https://www.kkr.mlit.go.jp/road/shintoshikenkyukai/pdf/da...
***ビーム要素 [#c879e22d]
[[Salome-Meca ビーム要素]]
***振動(sin波)の与え方 [#be3c2550]
1,csv形式でデータを作る。私の解析ではインパクトを与えた...
2,Functions and Lists→DEFI_FONCTION→VALEにcsvファイルを...
3,BC and Load→AFFE_CHAR_MECA_F→FORCE_NODALE→振動させた...
4,Analysis→DYNA_VIBRA→EXCIT→CHARGEに3を選択~
下の設定で解析をすると2重境界条件のエラーが発生した。~
1~
func1 = DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0.0, 0.0, 14.99, 0.0, 15.0, 29389.2626146241, 15....
func3 = DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0.0, 0.0, 100.0, 0.0))
load3 = AFFE_CHAR_MECA_F(
DDL_IMPO=(_F(
DX=func3,
DY=func1,
DZ=func3,
GROUP_MA=('kotei1', )
), _F(
DX=func3,
DY=func1,
GROUP_MA=('kotei2', )
)),
MODELE=model)
***SALOME エラー [#ad1f7faa]
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!...
! <S> Exception utilisateur levee mais pas interceptee...
! Les bases sont fermees. ...
! Type de l'exception : error ...
! ...
! Erreur écriture de l'enregistrement 766379 sur ...
! code retour : -4 ...
! Erreur probablement provoquée par une taille tr...
! travail. ...
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!...
及川さんからいただいた参考モデルのplateをなくして、ten1〜...
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! <S> Exception utilisateur levee mais pas interceptee...
! Les bases sont fermees. ...
! Type de l'exception : error ...
! ...
! Solveur MUMPS : ...
! Problème ou alarme dans le solveur MUMPS. ...
! Le code retour de MUMPS est : -2 ...
! ...
! Conseils : ...
! Consulter le manuel d'utilisation de MUMPS. ...
! Prévenir l'équipe de développement de Code_Aster. ...
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!...
Erreur écriture de l'enregistrement 337548 sur la base :...
! code retour : -4 ...
! Erreur probablement prov
<I>_EXIT_CODE = 256
**実験 [#yf28f66b]
***鋼棒の引張実験 [#ecabe210]
鋼棒の第一降伏点1160MPa 荷重=1160(N/mm2)☓(3.5☓3.5π)=44.6...
・ジグへの鋼棒の設置方法~
引張実験を行う際に曲げが入らないようにジグの中央に垂直に...
・ひずみ装置の操作~
最初にスタートを押し、測定開始時にもう一度スタートを押す。~
その後は一定の間隔でスタートを押していく。~
・結果の処理の仕方~
1,ひずみの表~
①ひずみゲージ1〜4の各結果に☓10^(-6)を行う。~
②各ひずみゲージの最初の値を0にするために、得られた各値を...
2,荷重と伸びの表~
①0sと実験時にスタートボタンを押した間隔ごとに荷重を抽出す...
ex)実験時に30sごとにスタートボタンを押した場合、0s、30s...
②0s時の荷重を0にするように各値を補正する。~
③補正した荷重を断面積で除して応力を求める。~
3,応力-ひずみのグラフ作成~
1で得たひずみの値と2で得た応力から、応力-ひずみのグラフを...
ひずみゲージが4つあるときは、4本の線が1つのグラフに入るよ...
4,ヤング率の計算~
・ひずみの表から計算する方法~
①2で得た最大応力の10%と40%時の応力と各ひずみゲージのひず...
②E=(σ40%-σ10%)/(ε40%-ε10%)から計算する。~
・伸びから計算する方法~
①2で得た最大応力の10%と40%時の応力と伸びを抽出する。~
②E=(f/A)☓(l/x)=(σ40%-σ10%)☓(l/(x40%-x10%)から計算する。(l...
*創造工房実習 [#e8c318e1]
***2023 12月8日 創造工房実習 [#c048ca60]
CENTER:サンドウィッチ梁
,メッシュの長さ,要素数,変位(mm),相対誤差(%),計算者
,0.7,155192,0.08378905246,15.365,安藤
,0.8,138808,0.08380386491,15.350,安藤
,0.9,82587,0.083707073981,15.45,兼田
,1.1,38671,0.084201207602,14.95,兼田
,1.2,31929,0.083688,15.466,柴田
,1.3,28621,0.083669,15.4857,柴田
,1.4,28854,0.08368,15.47,佐藤
,1.5,20015,0.084052,15.10,佐藤
,1.6,19448,0.0835402938,15.62,皆川
,1.7,13801,0.0834355098,15.72,皆川
,1.8,12528,0.083733,15.42,永山
,1.9,11769,0.083924,15.23,永山
,2,10699,0.084076876559,15.074,辻
,3,3579,0.08414561753,15.004,辻
,4,1628,0.082794,16.37,服部
,5,1016,0.083033,18.89,服部
,6,839,-0.082882,16.26,梶原
,7,554,-0.080871,18.28,梶原
,8,285,0.079995,19.20,工藤
,9,261,0.078980,20.22,工藤
,10,232,0.081911,17.26,佐々木
,11,208,0.075676,23.56,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***2023 11月29日 創造工房実習 [#r1eacb71]
異方性一次 [#p29b4460]
CENTER:異方性一次
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.7,144563,0.505252,2.76,安藤
,0.8,141517,0.504692,2.64,安藤
,0.9,91648,0.502595,2.216,兼田
,1.1,27160,0.489914,0.363,兼田
,1.2,24675,0.487088,0.791,柴田
,1.3,23446,0.4868010,0.995,柴田
,1.4,17738,0.485999,1.16,佐藤
,1.5,15438,0.485180,1.33,佐藤
,1.6,15900,0.483286,1.71,皆川
,1.7,12142,0.477952,2.80,皆川
,1.8,11604,0.482085,1.9554,永山
,1.9,10391,0.470887,4.2329,永山
,2,10291,0.480910,2.19,辻
,3,2328,0.431937,12.15,辻
,4,1500,0.430156,12.52,服部
,5,432,0.282968,42.45,服部
,6,356,0.3441556,30.00,梶原
,7,196,0.213934,56.49,梶原
,8,104,0.229874,53.25,工藤
,9,81,0.232308,52.75,工藤
,10,78,0.203271,58.65,佐々木
,11,63,0.222316,54.78,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
等方性二次 [#vb0def3f]
CENTER:等方性二次
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.7,144563,0.430124,3.22,安藤
,0.8,141517,0.430132,3.22,安藤
,0.9,91648,0.430020,3.197,兼田
,1.1,27160,0.429828,3.151,兼田
,1.2,24675,0.429836,3.15,柴田
,1.3,23446,0.42974,3.13,柴田
,1.4,17738,0.429797,1.3,佐藤
,1.5,15438,0.429958,3.14,佐藤
,1.6,15900,0.429755,3.18,皆川
,1.7,12142,0.429676,3.11,皆川
,1.8,11604,0.429829,3.1507,永山
,1.9,10391,0.429684,3.1159,永山
,2,10291,0.429620,3.10,辻
,3,2328,0.429169,2.99,辻
,4,1500,0.429254,3.01,服部
,5,432,0.428170,2.75,服部
,6,356,0.428452,2.82,梶原
,7,196,0.42591,2.21,梶原
,8,104,0.426074,2.25,工藤
,9,81,0.425552,2.12,工藤
,10,78,0.488382,17.20,佐々木
,11,63,0.423972,9.0534,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***2023 11月24日 創造工房実習 [#eab4b630]
単純梁~
,メッシュの長さ,要素数,変位[mm],相対誤差,計算者
,0.7,145234,0.422484,1.388,安藤
,0.8,142973,0.422570,1.409,安藤
,0.9,91648,0.420437,0.897,兼田
,1.1,27160,0.405618,2.659,兼田
,1.2,24675,0.404349,2.96,柴田
,1.3,23446,0.404185,3.00,柴田
,1.4,17738,0.398604,4.34,佐藤
,1.5,15438,0.396593,4.83,佐藤
,1.6,16122,0.398212,4.44,皆川
,1.7,12026,0.393411,5.59,皆川
,1.8,11604,0.393668,5.53,永山
,1.9,10391,0.390695,6.24,永山
,2,10921,0.395103,5.18,辻
,3,2328,0.324762,22.06,辻
,4,1500,0.155013,62.80,服部
,5,432,0.065278,84.33,服部
,6,357,0.213062,48.87,梶原
,7,196,0.1019,75.55,梶原
,8,104,0.1158624,72.20,工藤
,9,81,0.1255118,69.88,工藤
,10,78,0.07733,81.44,佐々木
,11,63,0.1999,52.03,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***2023 11月17日 創造工房実習 [#xa1780bc]
Salomeの使い方を学んだ。~
Geometry~
Mesh作成~
Aster Study~
Pravis~
片持梁~
,メッシュ長さ,要素数,先端変位(4隅の平均値)[mm],相対誤差,...
,0.7,198464,6.54281,1.91,安藤
,0.8,113812,6.5104,2.39,安藤
,0.9,40280,6.3631525,4.60,兼田
,1.1,30055,6.3363525,5.00,兼田
,1.2,26467,6.3043375,5.48,柴田
,1.3,25180,6.304355,5.48,柴田
,1.4,32212,6.31612,5.31,佐藤
,1.5,17753,6.1209,8.23,佐藤
,1.6,14296,6.2044625,6.98,皆川
,1.7,13596,6.2156625,6.81,皆川
,1.8,2866,5.737755,13.98,永山
,1.9,6001,5.7263625,14.15,永山
,2,5617,5.6458525,15.355,辻
,3,2309,5.4728755,17.948,辻
,4,617,3.6160575,45.12,服部
,5,494,3.8580375,42.16,服部
,6,581,2.50682,62.416,梶原
,7,133,1.41225,78.82,梶原
,8,78,1.2887175,80.68,工藤
,9,72,1.2879925,80.69,工藤
,10,60,1.14344,82.85,佐々木
,11,65,1.23124,81.154,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***2023 11月10日 創造工房実習 [#fb95bbcf]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***2023 10月27日 創造工房実習 [#fb5ee35c]
コマンドを学んだ。次回までに使いこなせるように頑張りた...
pwd 今いる場所がわかる~
ls ファイルを表示する~
mkdir ファイル名 ファイルの作成~
rmdir ファイル名 ファイルの削除~
cd ファイル名 ファイル名のところへ移動~
vi ファイル名 編集~
esc : w 保存~
esc : q 戻る~
cp コピー元 コピー先 コピーする時~
皆川恭輔
終了行:
[[FrontPage]]
#contents
*卒論 [#x5b6f34a]
**テーマ [#xf3b401f]
ケーブルの腐食やクラックを考慮した固有振動解析（皆川）
昨年度の山本の研究の続き．ケーブルの腐食方法やimpact...
昨年できなかった130mくらいの長いケーブルを要素分割し...
**参考資料 [#m1131b7d]
・山本さんの卒論
ケーブル腐食をくり抜くことで再現~
腐食の長さ・深さ・位置を変えて解析
・角田さんの修論
健全モデルと腐食モデルを使って、ケーブル破断時の連鎖崩壊...
C1~C14をそれぞれ破断させた時に崩壊するかの解析~
ケーブル健全時・腐食時の動的増幅率(DAF)と荷重再分配率(LRR...
*ケーブルの固有振動数の理論値 [#pdf8bc9e]
**長さ71m 直径0.315m [#ndfe1ce2]
・面固定~
1次モード 2.92969Hz~
3次モード 8.59375Hz~
5次モード 13.9648Hz~
7次モード 18.8477Hz~
・線固定~
1次モード 2.83203Hz~
3次モード 8.49609Hz~
5次モード 13.7695Hz~
7次モード 18.5547Hz~
*2/20〜 [#h56ec792]
**2/20 1次モードのsin波 Scale1.0 [#s868fbc5]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**2/28 ケーブルのみモデル 1次モード 500Nのsin波(1.66016Hz...
Scale1.0~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
中央の鉛直方向変位~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
*11/27〜 [#o4e79b21]
**11/27 引張実験 [#x9106e6a]
腐食レベル1 1-2 最大引張荷重64.25kN 応力1669.5MPa~
腐食レベル2 2-2 最大引張荷重64.00kN 応力1663.0MPa~
腐食レベル3 3-2 最大引張荷重63.55kN 応力1651.3MPa~
**IISさんからのアドバイス [#e874b20f]
①重力加速度を入れてみる~
②加速度と張力の関係を見てみる~
③ノッチ位置の変更~
**中間発表での質問 [#d9f25442]
①なぜ引張解析の結果で半径分の深さのノッチを入れると引張荷...
②どれくらいノッチが入った時に見つけられればよいのか(深く...
③新しい腐食発見方法はどんなものがありそうか~
**重力加速度をいれた振動解析 [#r95bd45a]
***健全モデル 解析時間50s [#a92626fb]
これまでと同じように解析時間を30s、インパクトを15sで与え...
解析時間を50s、インパクトを30sで与える時間設定にした。~
解析結果は下のようになった。卓越振動数は1次モードのみであ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**加速度の比較 100MN 重力加速度なし [#hfcbf938]
健全モデル、深さ157.5mmモデル、深さ45.00mmモデルで比較し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**30MN 重力加速度なし [#vf7662cc]
健全モデルとクラックモデルの加速度波形を比較した。重力加...
***加速度の比較 [#hd648f33]
クラックの深さでは加速度は変わらなかったが、健全とは異な...
・ノッチ深さでの比較~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
・各測定点での比較~
健全~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
ノッチあり~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
***PSDの比較 [#cdf20b98]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**ケーブルのみ 振動解析 重力加速度あり str30MN 1/22〜...
***時間設定 [#p93958f1]
0〜5秒:段階的に張力負荷~
5秒〜:重力加速度を入れる~
ケーブルを安定させる~
10〜14.99秒:ケーブルを安定させる~
15秒:インパクト~
15〜30秒:速度と加速度を測定する~
***解析結果 振動数 [#e10c2075]
健全 1次:1.66016Hz 3次:4.98047Hz 5次:8.49609Hz~
Crack15.75 1次:1.66016Hz 3次:4.98047Hz 5次:8.49609Hz~
Crack157.5 1次:1.66016Hz 3次:4.98047Hz 5次:8.49609Hz~
やはりノッチの有無・深さで振動数は変化なかった。~
***解析結果 加速度 [#o9d00940]
測定点3点ごとの応答加速度~
だいたい同じような結果が出た。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
***解析結果 インパクト後の振動の様子 [#v22ad56e]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**カバーありモデル 振動解析 重力加速度あり str30MN 1/2...
***材料と時間設定 [#h993758e]
材料~
ケーブル:E=195000,nu=0.3,ρ=7.8☓10^(-9)~
仮想材料:E=50,nu=0.3,ρ=1.0☓10^(-11)~
カバー:E=1000,nu=0.3,ρ=9.0☓^(-10)~
時間設定~
0〜5秒:段階的に張力負荷~
5秒〜:重力加速度を入れる~
ケーブルを安定させる~
10〜14.99秒:ケーブルを安定させる~
15秒:インパクト~
15〜30秒:速度と加速度を測定する~
***卓越振動数 [#v04da4a2]
健全 １次:1.56Hz ３次:4.88Hz ５次:8.40Hz~
Crack15.75 １次:1.56Hz ３次:4.88Hz ５次:8.40Hz~
Crack157.5 １次:1.56Hz ３次:4.88Hz ５次:8.40Hz~
***応答加速度 [#ab11b7c2]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
***インパクト後の振動の様子 [#z266255b]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**参考 腐食モデル W69000 D15750 振動解析 重力加速度あ...
str30MN~
ノッチの深さでは卓越振動数などが変化していなかったので、...
***卓越振動数 [#m0fd5707]
1次:2.05078Hz ３次:6.15234Hz ５次:10.2539Hz~
健全モデルとの相対誤差 １次:23.5% ３次:23.5% ５次20.7%~
***応答加速度 [#dc2c71e3]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**1/31 応答加速度が変化した原因の検討 [#gadcee86]
***張力が変化したことで応答加速度が変化した? [#ndafb0df]
ノッチが深く入ることで応力集中しやすくなることは引張解析...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
ノッチが深いと変位が大きくなっていたが、これだけだと張力...
*9/20〜 [#w87deece]
**幾何学非線形を考慮した振動解析 [#v7534fc8]
振動解析を幾何学非線形を考慮して解析した。~
理論値との誤差が非常に小さくなった。(理論値 1次:2.837Hz ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
卓越振動数である2.83203Hzのsin波を与えた時の振動~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
卓越振動数である8.49609Hzのsin波を与えた時の振動~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
卓越振動数である13.7695Hzのsin波を与えた時の振動~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**クラック 幾何学非線形性 [#la89c235]
***長さ71m 半径157.5mm [#j80656c1]
①クラック深さ157.5mm~
1次 2.83203Hz 3次 8.49609Hz 5次 13.7695Hz~
②クラック深さ146.25mm~
1次 2.83203Hz 3次 8.49609Hz 5次 13.7695Hz~
③健全~
1次 2.83203Hz 3次 8.49609Hz 5次 13.7695Hz~
クラックでは固有振動数に変化はなかった。~
昨年の山本さんの腐食モデルでは振動数が変化していたことを...
**腐食レベル2 [#p9c7275f]
腐食レベル2のケーブルは下のようにサビが入っていた。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
そこでこのケーブルのメッキ部分を削り、どの程度腐食が進展...
***引張解析 [#t096d3e7]
・300mm 7mmモデル
まずはサビがある部分のみを削り、サビの深さと幅をモデルで...
引張応力はσ=1164MPaとなり、健全時と同様の結果が得られ、サ...
・71000mm 315mmモデル
次に300mmモデルの腐食深さを45倍(315mm/7mm)、幅を236.7倍(7...
引張応力はσ=1146MPaとなり、腐食による影響は見られなかった...
**6角柱モデル 2本 [#h16a7ac1]
下の図形のように面接続と線接続で同じ結果が得られるかを確...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・振動解析~
,モード,面接続 結果,線接続 結果,相対誤差%
,1次,3.125Hz,3.125Hz,0.29%
,3次,9.17969Hz,9.17969Hz,1.81%
,5次,14.7461Hz,14.7461Hz,5.36%
・引張解析~
面接続 引張応力σ=1152MPa、point71000での応力=1156MPa~
線接続 引張応力σ=1152MPa、point71000での応力=1156.1MPa~
振動解析、引張解析ともに面接続と線接続では同じ結果が得ら...
**6角柱モデル 面結合[#me8ea140]
実際のケーブルに近づけるために下のように6角柱を7本・19本...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***7本モデル [#if448901]
7本モデルでは半径60mmの円上に点を作って、底面の6角形を作...
・振動解析~
100MNの張力を与えての振動解析~
1次 3.125Hz 相対誤差1.06%~
3次 9.277Hz 相対誤差0.012%~
5次14.844Hz 相対誤差3.99%~
7本モデルの底面と同じ面積の円柱で振動解析を行ったところ、...
・引張解析~
引張応力σ=1152MPa point71000での応力=1156.8MPa~
***19本モデル [#seb4db28]
**6角柱モデル 線結合 [#s28eb61f]
・振動解析
線結合でも面結合と同じ結果になった。~
**腐食 振動解析 [#v6c834e9]
***腐食再現モデルの種類 [#dddba9f9]
幅と深さを下の図のようにする。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
まずは中央部に下にある、幅7種類と深さ6種類の計42種類の解...
・幅 ()内は長さに対する割合~
①710mm(1%) ②3550mm(5%) ③7100mm(10%) ④14200mm(20%) ⑤21300m...
・深さ ()内は半径に対する割合~
①15.75mm(10%) ②47.25mm(30%) ③78.75mm(50%) ④110.25mm(70%) ...
幅・深さ以外にも腐食の位置も固有振動数に影響してくると思...
また、ボリュームが同じであれば、えぐり方(腐食の入り方)で...
***結果 [#e6d401bb]
下のグラフは、横軸に健全時のボリュームに対する腐食モデル...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
次のグラフは、上のグラフの横軸を直径に対する腐食深さの割...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
次は腐食の位置が各振動モードの固有振動数にどのように影響...
①中央にW28400の腐食~
②両端から1000mmの位置にW14200ずつの腐食(合計でW28400の腐...
③左端から1000mmの位置にW28400の腐食~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
**被覆ありモデル 振動解析 面固定[#rbede1f0]
***被覆の半径250mm 仮想材料の半径200mm [#p7eaf383]
ケーブルの半径はこれまで通り157.5mm、仮想材料の半径を200m...
材料諸元は~
ケーブル:E=195000MPa,nu=0.3,ρ=7.8☓10^(-9)ton/mm^3~
仮想材料:E=50MPa,nu=0.3,ρ=1☓10^(-11)ton/mm^3~
被覆(高密度ポリエチレン HDPE):E=1000MPa,nu=0.3,ρ=9☓^(-10)...
固有振動数の結果 ()はケーブルのみモデルとの相対誤差~
1次モード 2.73438Hz (6.67%)~
3次モード 8.20313Hz (4.55%)~
5次モード 13.3789Hz (4.20%)~
7次モード 18.0664Hz (4.15%)~
→ρが大きくなったことでケーブルのみの時よりも固有振動数が...
***被覆の半径190mm 仮想材料の半径170mm [#x22166de]
上の被覆ありモデルの被覆と仮想材料の半径のみを変えて解析...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2024/minakawa/...
固有振動数の結果 ()は面固定のケーブルのみモデルとの相対...
1次モード 2.83203Hz (3.33%)~
3次モード 8.49609Hz (1.14%)~
5次モード 13.7695Hz (1.40%)~
7次モード 18.5547Hz (1.55%)~
**腐食 引張解析 [#m2b0cdda]
振動解析の時の腐食モデルをすべて同じ張力を与えた時のひず...
張力は500MPaを与えた。~
結果はえぐった部分の断面積(腐食深さ)でひずみが決まるよう...
*振動解析 [#n270fa13]
**2段階解析 [#occ4ae6c]
***1段階目に張力 2段階目にimpact [#g7002f2d]
1段階目:固定条件と張力を設定し、ASSEMBLAGEの境界条件に固...
2段階目:impactを設定し、ASSEMBLAGEの境界条件に固定条件とi...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***1段階目に張力 2段階目に固有振動解析(CALC_MODES) 由利...
今回の解析では由利橋の1/100モデルを使用した。~
1段階目:固定条件と張力を設定し、ASSEMBLAGEの境界条件に固...
2段階目:固有振動解析(CALC_MODES)~
y方向の振動数の結果は下の表のようになった。~
1次 14.00Hz~
2次 55.65Hz~
3次 124.82Hz~
4次 222.89Hz~
5次 346.32Hz~
***1段階目に張力 2段階目にimpact 由利橋の1/100モデル [#d...
今回の解析は由利橋の1/100モデルで行った。~
①1段階目・2段階目ともにDYNA_VIBRAを使い、2段階目はETAT_IN...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
②1段階目STAT_NONLINE、2段階目DYNA_NONLINEを使い解析をおこ...
結果は上のグラフと同じ振動数がでた。
**9/16 振動解析 PRES_REP [#af4e71ac]
指定面にPRESを与えるコマンドPRES_REPを使って解析を行った。~
SIZZは与えたい張力分が入っていた。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**9/12 振動解析 PRE_EPSI [#y21e4eb2]
PRE_EPSIについては3Dでも適用することができたが、σ=Eεで100...
PRE_EPSIでは正の値と負の値を両方を与えたが、どちらの解析...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**9/10 振動解析 PRE_SIGM [#q7a525ce]
今回の解析ではPRE_SIGMを使って解析した。(DYNA_VIBRA)~
PRE_SIGMについては、impactなしで解析したところ与えたい張...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**長さ1000mm、半径0.5mmの健全モデル [#c0a583be]
固定の仕方を線固定にした。(張力は面載荷で与えた。)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
固有振動数の理論値 f=6.37Hz 解析時間0〜30秒 細長比354...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
,時間,行うこと
,0〜5秒,張力(1N)をかける。
,5〜14.99秒,ケーブルを安定させる。
,15秒,10Nの力を与える。
,15.01〜30秒,卓越した振動を検出する。
**由利橋のケーブルモデル 張力0 [#jadf6985]
張力0のモデルにインパクトを与えて、その結果をFFT処理した...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**由利橋のケーブル 長さ71000mm、半径157.5mmの健全モデル...
固有振動数の理論値 f=2.85Hz 解析時間0〜100秒 細長比25...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
,時間,行うこと
,0〜5秒,張力(100MN)をかける。
,5〜14.99秒,ケーブルを安定させる。
,15秒,100kNの力を与える。
,15.01〜100秒,卓越した振動を検出する。
結果~
Mesh50 ボリューム169538~
中央部の振動数 ~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
波形を書かせるとおかしな波形となり、モードの推定ができな...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***6/24 内部張力が、与えている張力と等しいかの解析 [#nc22...
中央部(point35500)と張力を与えた部分(point71000)の15秒〜1...
中央部(point35500)~
平均値99.6MN 最大値104.5MN 最小値98.7MN~
張力を与えた部分(point71000)~
平均値・最大値・最小値ともに95.6MN~
***インパクト1000kNに変更 [#ofefb2d1]
上の解析のインパクトを100kN→1000kNに変更して解析を行った。~
インパクトを与えた後(16.20s)の100kNの時のモードと1000kNの...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***6/21 張力とインパクトをいろいろ変えて解析 [#ldc16f87]
①張力100MN、インパクト100kN~
②張力100MN、インパクト1000kN~
③張力200MN、インパクト100kN~
④張力200MN、インパクト1000kN~
どの解析でも下のような同じ振動数が出力された。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**由利橋のケーブル(長さ71000mm、半径157.5mm)に振動を与え...
***両端に振動を与える解析 [#xaa17730]
両端に理論値である2.85Hzの振動を与える解析~
ParaVISでモードを見ると下のようになり、微細な振動は見られ...
両端を固定しており、その部分に振動を入れているからだと思...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***中心に振動を与える解析 [#s13d0a3e]
長さ71000mm、張力100MNの弦の理論値である2.85Hzを中心の点...
モードは以下のようになった。~
振動を与えた後~
16.41秒~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
16.58秒~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
17.94秒~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
中心に1.16Hzの振動を与える解析を行った。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
振動を与える直前(14.99秒)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
振動を与えた直後(15.5秒)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
中心に3.2Hzの振動を与える解析を行った。~ [#t2afbbe4]
振動を与える直前(14.99秒)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
振動を与えた後は下の2つの位相が繰り返された。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
中心に1.16Hzと3.2Hzの振動を与えたときの比較 [#j7e3c02b]
振動を与えた点(中心部)のy方向変位~
y方向の最大変位~
1.16Hz:11070.7mm 3.2Hz:4837.51mm~
y方向の最小変位~
1.16Hz:-10744.2mm 3.2Hz:-4925.62mm~
***両端から5mmの位置に振動を与える解析 [#ue6d276e]
振動を与える位置のモデル~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
両端から5mmの位置に0.122Hzの振動を与える解析を行った。~
振動モードは下のように、１次モードとなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
両端から5mmの位置に1.16Hzの振動を与える解析を行った。~[#b...
振動モードは下のように、３次モードとなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
両端から5mmの位置に3.2Hzの振動を与える解析を行った。~[#t5...
振動モードは下のように、５次モードとなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**7/30 張力の与え方を単位体積あたりの荷重に変えて解析 [#z...
張力の与え方は、今までの解析では面載荷(FORCE_FACE)で与え...
FFT処理をすると以下のような結果となった。(理論値は2.85Hz)~
張力を100MN→1000MNに変更して解析したが、下の結果と同じに...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**6/28 ビーム要素(温度荷重ではない) [#ne6248ee]
ビーム要素で由利橋健全ケーブル(長さ71000mm、半径157.5mm)...
梁の曲げ振動数の理論値と等しくなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**強制変位 [#n235dd68]
100MNの張力を与えたときに475mmz方向に変位していたので、今...
強制変位させることでも応力が入ることがわかったが、場所ご...
***インパクトを与える解析 [#db9e6d09]
これまで通りインパクトを与えて、それで得た結果をFFT処理し...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***静的解析 [#t0385b50]
強制変位させた場合の応力の入り方が初期張力のように入って...
1次モード 0.13Hz~
2次モード 0.51Hz~
3次モード 1.15Hz~
4次モード 2.05Hz~
5次モード 3.20Hz~
張力が0の時とほぼ同じ結果となった。~
**両側引張 [#zd1f0a59]
これまでの解析では片側のみを引っ張って張力を与えていたが...
結果を見ると、片側引張では完全固定部分で張力が小さくなっ...
両側引張では全体的にほぼ均一な張力が入っていたが、結果は...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**インパクトの与え方をじっくり与えた [#j5bb89cb]
これまでのインパクトの与え方は15sで一気に与えていたが、今...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**解析モデルを由利橋の1/100にした(長さ710mm、直径3.15mm) ...
***ホールドあり [#b7e1b619]
振動後に面がz方向に垂直でないことで与えたい張力が与えられ...
右の表がインパクトを与えてMATLABで処理した時の結果(理論値...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***周波数応答解析 [#zdcabd3d]
CENTER:周波数応答解析 面固定
,モード数,結果,理論値,理論値との誤差,梁の曲げ振動,曲げ振...
,1次,21.9Hz,28.6Hz,23.4%,19.2Hz,14.1%
,2次,71.2Hz,57.1Hz,24.7%,62.1Hz,14.7%
,3次,147.8Hz,85.7Hz,72.5%,129.6Hz,14.0%
CENTER:周波数応答解析 線固定
,モード数,結果,理論値,理論値との誤差,梁の曲げ振動,曲げ振...
,1次,14.1Hz,28.6Hz,50.7%,19.2Hz,26.6%
,2次,56.2Hz,57.1Hz,1.6%,62.1Hz,9.5%
,3次,125.6Hz,85.7Hz,46.6%,129.6Hz,3.1%
**由利橋 長さ71000mm、半径157.5mmの中央腐食モデル(長さ71...
固有振動数の理論値 f=2.85Hz 解析時間0〜100秒~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
載荷線と腐食部分~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
*引張解析・引張実験 [#f72db340]
**引張解析 [#h726c276]
クラックの解析は専用の解析方法があるらしいので、それでお...
クラックは三角柱でくり抜いて表現した。(下の図はクラック角...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***健全モデル [#y73b55e6]
このモデルでは
***クラック角度 90° クラック深さ157.5[#je549823]
このモデルでは445MPaまで引っ張ることができた。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***9/27 長さ300mm 直径7mmモデル 角度一定 クラック幅0.4 ...
角度一定でクラックの深さを0.25mm刻みで引張解析を行った。...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
応力集中の様子が下の図(クラック深さ3.5mm)~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
上の①〜⑤の解析のz方向変位と荷重の関係を表したものが下のグ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***10/7 長さ71000mm 直径157.5mmモデル 角度一定 クラック...
今回の引張解析は、9/27の引張解析のクラックと角度を同じに...
角度一定、深さを11.25mm刻みで解析した結果が下のグラフ・表...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**引張実験 [#u3aa1cad]
***8/28 健全の鋼棒・腐食3の鋼棒 [#g746217e]
①健全の鋼棒~
まずは8/22の実験と同様に、下のジグで掴む部分をネジ状に溝...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・ひずみゲージから求めたヤング率~
ひずみゲージNo1：E=197GPa~
ひずみゲージNo2：E=204GPa~
ひずみゲージNo3：E=203GPa~
ひずみゲージNo4：E=221GPa~
②健全の鋼棒~
次に上下のジグで掴む部分を平面にして、実験をおこなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・ひずみゲージから求めたヤング率~
ひずみゲージNo1：E=245GPa~
ひずみゲージNo2：E=198GPa~
ひずみゲージNo3：E=191GPa~
ひずみゲージNo4：E=220GPa~
③腐食レベル3~
最後に上下のつかみ部のメッキを剥がした腐食レベル3の鋼棒を...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・ひずみゲージから求めたヤング率~
ひずみゲージNo1：E=239GPa~
ひずみゲージNo2：E=226GPa~
ひずみゲージNo3：E=215GPa~
ひずみゲージNo4：E=191GPa~
***8/22 健全の鋼棒 [#od2d91a9]
今回の引張実験では、鋼棒の下のジグ部分をネジ状に溝を入れ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・ひずみゲージから求めたヤング率~
ひずみゲージNo1：E=203GPa~
ひずみゲージNo2：E=197GPa~
ひずみゲージNo3：E=197GPa~
ひずみゲージNo4：E=224GPa~
***7/16 腐食レベル2の鋼棒 [#d86d0401]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
・ひずみから求めたヤング率~
ひずみゲージNo1：E=274GPa~
ひずみゲージNo2：E=233GPa~
ひずみゲージNo3：E=349GPa~
ひずみゲージNo4：E=269GPa~
・伸びから求めたヤング率~
E=305GPa~
何らかの原因で実験がうまくできていなかった。~
ジグの部分をうまく掴ませるように工夫したい。~
**論文まとめ [#ff52d2ef]
***振動法によるケーブル張力の実用算定式 [#gf30995b]
曲げ剛性やサグ、傾斜角の影響を考慮した張力算定式~
曲げ剛性の影響を表すパラメーターが小さい時(曲げ剛性が大き...
(サグや傾斜角の影響を表すパラメーター)＜3.0の時は、対称振...
***斜材点検用非破壊検査装置の開発と運用 自走式点検装置 [#...
斜ケーブルの変状調査を目的とした自走式点検装置の機能・硬...
装置の特徴~
①カメラを4台設置することで全周を同時に調査可能~
②保護管上のクラックケージを0.1mmまで判読可能~
③保護管内の鋼線のうち外周部分の破断の有無の調査可能~
④鋼線1本の破断から調査可能~
⑤地上のパソコンから遠隔操作可能~
⑥照明を装備しているため夜間の点検可能~
利点~
①点検者の近接目視と同じくらいの精度で点検可能~
②地上からの遠隔操作で安全に点検可能~
③高所作業車を使用する方法よりも安価~
④交通規制の規模も低減可能~
⑤保護管に変状があった場合、同時に斜ケーブルの変状の有無も...
**幾何学的非線形性の参考サイト [#o7f2f199]
https://www.cybernet.co.jp/ansys/learning/glossary/kikaga...
**進捗 [#f2f072c8]
***9/10〜青木班ゼミまで [#p8847e9f]
9/10 PRE_SIGMコマンドを使って解析した。静的解析だと張力が...
結果は1次0.097Hz、3次1.17Hz、5次3.22Hzだった。~
9/12 imp100kN(5sから一定)、str(5sから一定)の時間設定でDYN...
9/13 impact後の応力と変位について様々な点(y=0の点)をとっ...
9/14 DEFI_FOND_FISSでクラック解析をした。1個ずつエラーを...
9/15 PRE_EPSIで解析した。正の値と負の値を入れた解析の2つ...
9/16 圧力をかけるコマンドPRES_REPを張力面にかけて解析した...
9/17 先週の水曜日に青木さんが考えていただいたモデルで解析...
Liaisonで結合させたらうまくできた。Liaisonでは結合させた...
***8/27〜9/9 [#z8ca6aed]
8/27 振動解析において、張力を与える側がz方向に固定されて...
しかし、結果は0sから固定されており、張力→固定ができなかっ...
8/28 引張実験をおこなった。①ネジ＆メッキ剥がし②平面③腐食...
8/29 昨日の引張実験の結果の確認をした。1200MPa前後で降伏...
8/30 象潟のPCラーメン箱桁橋の現場見学に行った。~
9/2 振動解析では710mmモデルでFORCE_INTERNE→CALC_MODESをお...
引張解析ではクラックの導入の仕方を調べた。今まではクラッ...
XFEMはクラックの大きさ、形、始点を指定することで自動でク...
XFEMではクラック内に線要素があるというエラーがでた。(実際...
9/3 振動解析では強制変位+100MNで張力を与えた。結果は振動...
引き続き引張解析のクラックの導入の仕方を模索している。~
9/4 クラックの導入の仕方でDEFI_FOND_FISSを使う場合、3Dで...
INI_XFEM_ELNOに関するエラーが発生したがよくわからない。~
9/5 振動解析では両端固定＆FORCE_INTERNEの結果を確認した。...
9/6 振動解析において、張力側の面がz軸方向に垂直にならない...
結果は0.195Hz、1.37Hz、3.61Hz、6.84Hz・・・となった。~
9/9 振動解析について、時間で固定条件を変更して張力をかけ...
また、PRE_SIGMコマンドについて調べ、解析を回したところ、P...
***8/8〜青木班ゼミまで [#n0462cf2]
8/8 千代岡さんが見つけた2段階解析を少し変更して(千代岡さ...
今回の解析は1段階目に固定条件と張力を設定(loadとfunction)...
また、引張試験のジグの噛ませ方を調べたがまだ良い方法が見...
8/9 所々変更しながら2段階解析を行った。~
昨日行った2段階解析の結果をFFT処理した。結果は結果欄に書...
次に1段階目に張力をかけ、2段階目に固有振動解析(CALC_MODES...
今日解析をして、千代岡さんと話していて疑問に思っていたETA...
そのため私の解析では使えないかもしれない。(私の解析は1段...
8/10 1段階目の解析をSTAT_NONLINEに変えて解析を行った。~
また由利橋の1/100モデルでETAT_INIT内のINST_INITに0を入力...
現在のところで、ETAT_INITについてわかっていることとしては...
1段階目の解析の時間の途中から2段階目の解析を始めるときはE...
8/19 由利橋(71000mm)モデルでは容量が足りなくなってしまう...
また、引張解析について調べ、由利橋(71000mm)の健全モデルの...
8/20 由利橋の1/100モデルの解析結果を再確認した。DYNA_VIBR...
STAT_NONLINE(張力)→DYNA_NONLINE(impact)では1段階目の張力...
また、引張解析のCrackの与え方をいろいろ試してみた。三角柱...
8/21 引張解析についていろいろと試した。降伏応力とひずみの...
8/22 引張実験をおこなった。鋼棒を掴む部分を①ネジのように...
①ネジのように溝を入れる方法では、溝は機械を使ってピッチを...
③平たくする方法では0.3mm分削った。~
今日の引張実験では①②の鋼棒を使って実験をおこなった。~
ジグの部分は滑っていないようだった。(今回の実験では20kNま...
また、引張解析では90°クラックの解析を行った。引張応力は、...
また、引張解析では90°クラックの解析を行った。引張応力は、...
8/23 昨日の引張実験の結果をまとめた。途中までしか実験がで...
引張解析については結果を見るときにどの点で見ればよいかを...
引張解析のモデルは片持ち式でおこなっている。~
***8/1〜8/7 [#nd0852d8]
8/1 impactの位置を変えて解析したが、そもそも中心部以外にi...
結果は今までと変わらなかった。~
8/2 初期応力を与える良い方法がないか探した。PRE_SIGM(おそ...
PRE_EPSI(初期ひずみ)での解析を今度行ってみる。(初期ひずみ...
千代岡さんと2段階解析について確認した。まずは出力した結果...
8/3 単位体積あたりの荷重(FORCE_INTERNE)についてもう一度確...
z方向変位が入れたい張力分の1/2の変位しか入らなかった。~
8/7 PRE_EPSIについて解析を行った。結果は張力が入っていな...
PRE_EPSI(初期ひずみ)の与え方が違っていたかもしれないので...
***7/25〜7/31 [#f452332d]
7/25 及川さんからいただいたケーブル要素の例をもとに解析を...
モデルを作るところから頑張りたい。~
7/26 及川さんからいただいたケーブル要素のモデル(ケーブル4...
・ケーブル4本にするとMUMPSのエラーが発生する。(ten1〜4はX...
・プレート&ケーブル1本のモデルにするとエラーが発生する。~
・ケーブルに質点を与えた解析~
→張力が圧縮として入っていた。また、1次モードの振動数が縦...
7/28 振動した時に張力部分がz方向に垂直でないことで張力が...
Meshの設定の仕方に問題があるのかモデルが一体とならずに貫...
7/29 連成解析の方法を探した。参考資料は見つけることができ...
7/30 周波数応答解析を行ってみた。解析を少しでも早くするた...
結果は結果欄にある。~
7/31 FORCE_INTERNE(単位体積あたりの荷重)を使って張力を与...
FORCE_INTERNEを使った解析は下の3つを行ったが、全て同じ結...
・張力面をz方向フリーにして100MN分の単位体積荷重を与えた...
・張力面をz方向フリーにして1000MN分の単位体積荷重を与えた...
・両面を完全固定にして100MN分の単位体積荷重を与えた解析~
***7/18〜7/24 [#l26bd7fc]
7/18 英語の発表の準備を行った。~
7/19 オープンキャンパスで使う振り子の調整をした。~
7/21 温度荷重の解析をもう一度行ってみた。張力が入っている...
片側のみに応力を与えるのではなく、両側に応力を与える解析...
***7/11〜7/17 [#a9e78916]
7/12 引張試験の準備を行った。~
7/15 今まで行ってきた解析結果を確認して、今後の解析で張力...
・応力として線載荷で与えると張力がうまく入っていなかった...
・ビーム要素では梁の曲げ振動の値となった。(0.28Hz、1.50Hz...
・強制変位(475mm)では25〜50Hzに多く卓越していた。また、張...
・内部張力(FORCE_INTERNE)では0.51Hz、1.16Hz、3.70Hzとなっ...
・張力なしでは
7/16 引張試験を行った。使った鋼棒は腐食レベルが2-1のもの。~
破断させることはできなかったが、応力-ひずみのグラフを作成...
ジグの部分がうまくつかめていなかったので、どのようにした...
7/17 電子顕微鏡の使い方を学んだ。元素まで見ることができる...
昨日行った実験で改善すべき点であったジグの部分を修正した。~
及川さんと青木さんが、ジグが回らないようにすることとジグ...
***7/4〜7/10 [#j6070258]
7/4 ネットで見つけた資料を参考にしながら温度荷重の解析を...
7/8 前回(4日)の温度荷重の解析を修正しながら再び解析を行っ...
AnalysisにDYNA_VIBRAではなく、DYNA_NON_LINEを使うと解析が...
7/9 今まで行っていた温度荷重の解析では100℃上昇させていた...
この557.5℃は、今まで行っていた張力(100MN)を与える解析でpo...
7/10 強制変位させることで張力が入ることが分かったので、強...
***6/26〜7/3 [#i59cb769]
6/26 今週はケーブル要素、ビーム要素、温度荷重の解析を行っ...
千代岡さんから温度荷重やビーム要素について教わった。また...
及川さんいただいた、ケーブル要素について解析を回してみた...
6/27 昨日、千代岡さんから教えていただいたソリッド要素に内...
結果は上に書いているが、やはり理論値とは異なる値となった。~
ケーブル要素のエラーについても調べていきたい。~
7/1 青木さんがもらってきていただいたケーブルの引張試験の...
いただいたケーブルはCASE1(健全)とCASE5(削孔+塩水注水)~
7/2 ケーブル要素の解析がうまく回らなかった原因を考えたが...
やはり、「nddlリストにモデルがありません」というエラーが...
7/3 今後は温度荷重と初期ひずみについて進めていく。~
***6/20〜6/26 [#yffc04b9]
6/20 由利橋のケーブルモデルの両端から5mmの位置に0.122Hzを...
結果は1次モードとなった。~
6/21 中央部に張力がしっかりと入っているかを確認するため由...
また、張力を2倍の200MNに変えた解析と左右で逆側に振動を与...
今日の解析では張力を2倍にしても振動数が√2倍になっていなか...
インパクトは今まで100kNを与えていたがモードを見るとあまり...
6/23 張力0の解析を行い、理論値と比較した。~
結果は上に書いたが、張力を与えた場合よりも理論値が近くな...
6/24 モデル中央部に張力分だけの応力が入っているかを調べた...
また、今まで線固定にして解析を行っていたが、張力を与える...
明日は温度荷重について学びたいとおもう。~
***6/13〜6/19 [#w090b26f]
6/13 振動(sin波)を与える解析を及川さんに教えていただいた。
最初に1.16Hzのsin波を与えて解析を行った。与えたsin波がsin...
そこで50000sin(2π1.16t)として50kNの振動を与えた。結果は明...
6/14 50kNの振動を両端に与えて解析を回したがうまく結果が見...
そこで振動を両端ではなく、中心の点に1.16Hzの振動を与えて...
結果は上の(由利橋 健全モデルに1.16Hzの振動を与える解析)に...
6/15 昨日の解析を3.2Hzにして解析を行った。~
結果は上の(由利橋 健全モデルに3.2Hzの振動を与える解析)に...
1.16Hzの振動を与えた場合と3.2Hzの振動を与えた場合のy方向...
また、由利橋(長さ71000mm、半径157.5mm)の固有振動数の理論...
6/18 振動させる位置を両端から5mmの位置に変えて、1.16Hz,3....
1.16Hzでは３次モード、3.2Hzでは５次モードとなった。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***5/30~6/12 [#ceb33306]
5/30 真杜くんの実験のサポート~
5/31 実験サポート~
6/3 100mm,半径0.5mmの解析の結果をまとめようとしたが、パ...
6/4 パソコンが重く、結果が見れなかったため振動を検出する...
,時間,行うこと
,0〜5秒,張力(1N)をかける。
,5〜14.99秒,ケーブルを安定させる。
,15秒,10Nの力を与える。
,15.01〜30秒,卓越した振動を検出する。
6/5 上記の時間で解析したがパソコンがフリーズしてしまい結...
6/6 山本さんの解析では長さ50m・半径75mmで行っており、今...
そこで張力を3.2MN→1N、中央部への線載荷を100KN→10Nとして解...
結果を得ることができた。明日、青木さんにFFT処理をお願いす...
6/7 長さ1000mm、半径0.5mmの中央部に100mm腐食を与えたモデ...
長さ100mm、半径0.5mmの中央部に10mm腐食を与えたモデルでは...
細長比によるエラーかもしれない。~
長さ1000mm半径0.5mmモデルの細長比は3544.91~
長さ100mm半径0.5mmモデルの細長比は354.5~
6/8 由利橋のケーブル(長さ71000mm、半径157.5mm)の健全モデ...
月曜日に青木さんにFFT処理を依頼する。~
6/9 由利橋のケーブル(長さ71000mm、半径157.5mm)の中央腐食...
月曜日に青木さんにFFT処理を依頼する。~
6/10 FFT処理の行い方を学んだ。理論値とズレが生じたため、...
6/11 由利橋の健全モデルの解析値(3.2Hz)と理論値(2.85Hz)に...
最初の解析では、面載荷に張力(100MN)を円の断面積で割った値...
どちらの解析の結果もf=3.2Hzとなっていたため、他の考えられ...
***5/15〜5/29 [#w9e80a4f]
・山本さんの卒論の振動解析を解析してみる。~
山本さんの解析では測点を３点設定していたが、モードをより...
5/16 山本さんの卒論，角田さんの修論を読んだ．~
5/17 健全モデル、中央10m腐食、中央1m腐食のモデルを作った。~
5/20 健全モデルのAster Studyを設定し，解析を回してもらう．~
中央に線荷重を与える際は，3.2MNを断面積(半径75mmの円)で割...
5/21 昨日設定したAster Studyのコードの意味を調べた．~
5/23 健全モデルの解析を回した。~
50m 半径75mm 健全モデル(点は下の写真の上部に設置した。)
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
50m 半径75mm 中央1m腐食
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
**MATLAB [#bbd750c1]
1,変数→新規変数→処理したいデータを入れる~
2,FFTcable→データ名、スタート点(私の解析ではインパクト...
3,bandcable→下限値・上限値、データ名、スタート点の入力~
4,C1P3→csvファイルで波形を描く~
**Salome-meca [#ka48d877]
***Salome-mecaの使用方法 [#z8a41cac]
-[[Salomemeca使用解説法:http://opencae.gifu-nct.ac.jp/puk...
-[[code-asterの説明:https://code-aster.org/doc/v15/en/ind...
-[[code-aster質問ページ:https://forum.code-aster.org/publ...
***重力加速度の設定 [#h7a1bdbe]
材料設定で密度を設定する。~
AFFE_CHAR_MECA→PESANTEUR→GRAVITY=9800(mm/s^2)、DIRECTION=...
上の設定のDIRECTIONは、重力方向がy方向の時の設定。DIRECTI...
***引張解析 参考 [#i0ac39fd]
-[[Salome-Meca演習_弾塑性解析(2021)]]
***クラックの解析 破壊解析 [#h9b7cae4]
https://www.youtube.com/watch?v=5ivZy5rY9l8~
https://youtu.be/8fy0hmTITEQ~
https://www.opencae.or.jp/wp-content/uploads/2019/12/08-n...
https://www.kkr.mlit.go.jp/road/shintoshikenkyukai/pdf/da...
***ビーム要素 [#c879e22d]
[[Salome-Meca ビーム要素]]
***振動(sin波)の与え方 [#be3c2550]
1,csv形式でデータを作る。私の解析ではインパクトを与えた...
2,Functions and Lists→DEFI_FONCTION→VALEにcsvファイルを...
3,BC and Load→AFFE_CHAR_MECA_F→FORCE_NODALE→振動させた...
4,Analysis→DYNA_VIBRA→EXCIT→CHARGEに3を選択~
下の設定で解析をすると2重境界条件のエラーが発生した。~
1~
func1 = DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0.0, 0.0, 14.99, 0.0, 15.0, 29389.2626146241, 15....
func3 = DEFI_FONCTION(
NOM_PARA='INST',
VALE=(0.0, 0.0, 100.0, 0.0))
load3 = AFFE_CHAR_MECA_F(
DDL_IMPO=(_F(
DX=func3,
DY=func1,
DZ=func3,
GROUP_MA=('kotei1', )
), _F(
DX=func3,
DY=func1,
GROUP_MA=('kotei2', )
)),
MODELE=model)
***SALOME エラー [#ad1f7faa]
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!...
! <S> Exception utilisateur levee mais pas interceptee...
! Les bases sont fermees. ...
! Type de l'exception : error ...
! ...
! Erreur écriture de l'enregistrement 766379 sur ...
! code retour : -4 ...
! Erreur probablement provoquée par une taille tr...
! travail. ...
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!...
及川さんからいただいた参考モデルのplateをなくして、ten1〜...
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! <S> Exception utilisateur levee mais pas interceptee...
! Les bases sont fermees. ...
! Type de l'exception : error ...
! ...
! Solveur MUMPS : ...
! Problème ou alarme dans le solveur MUMPS. ...
! Le code retour de MUMPS est : -2 ...
! ...
! Conseils : ...
! Consulter le manuel d'utilisation de MUMPS. ...
! Prévenir l'équipe de développement de Code_Aster. ...
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!...
Erreur écriture de l'enregistrement 337548 sur la base :...
! code retour : -4 ...
! Erreur probablement prov
<I>_EXIT_CODE = 256
**実験 [#yf28f66b]
***鋼棒の引張実験 [#ecabe210]
鋼棒の第一降伏点1160MPa 荷重=1160(N/mm2)☓(3.5☓3.5π)=44.6...
・ジグへの鋼棒の設置方法~
引張実験を行う際に曲げが入らないようにジグの中央に垂直に...
・ひずみ装置の操作~
最初にスタートを押し、測定開始時にもう一度スタートを押す。~
その後は一定の間隔でスタートを押していく。~
・結果の処理の仕方~
1,ひずみの表~
①ひずみゲージ1〜4の各結果に☓10^(-6)を行う。~
②各ひずみゲージの最初の値を0にするために、得られた各値を...
2,荷重と伸びの表~
①0sと実験時にスタートボタンを押した間隔ごとに荷重を抽出す...
ex)実験時に30sごとにスタートボタンを押した場合、0s、30s...
②0s時の荷重を0にするように各値を補正する。~
③補正した荷重を断面積で除して応力を求める。~
3,応力-ひずみのグラフ作成~
1で得たひずみの値と2で得た応力から、応力-ひずみのグラフを...
ひずみゲージが4つあるときは、4本の線が1つのグラフに入るよ...
4,ヤング率の計算~
・ひずみの表から計算する方法~
①2で得た最大応力の10%と40%時の応力と各ひずみゲージのひず...
②E=(σ40%-σ10%)/(ε40%-ε10%)から計算する。~
・伸びから計算する方法~
①2で得た最大応力の10%と40%時の応力と伸びを抽出する。~
②E=(f/A)☓(l/x)=(σ40%-σ10%)☓(l/(x40%-x10%)から計算する。(l...
*創造工房実習 [#e8c318e1]
***2023 12月8日 創造工房実習 [#c048ca60]
CENTER:サンドウィッチ梁
,メッシュの長さ,要素数,変位(mm),相対誤差(%),計算者
,0.7,155192,0.08378905246,15.365,安藤
,0.8,138808,0.08380386491,15.350,安藤
,0.9,82587,0.083707073981,15.45,兼田
,1.1,38671,0.084201207602,14.95,兼田
,1.2,31929,0.083688,15.466,柴田
,1.3,28621,0.083669,15.4857,柴田
,1.4,28854,0.08368,15.47,佐藤
,1.5,20015,0.084052,15.10,佐藤
,1.6,19448,0.0835402938,15.62,皆川
,1.7,13801,0.0834355098,15.72,皆川
,1.8,12528,0.083733,15.42,永山
,1.9,11769,0.083924,15.23,永山
,2,10699,0.084076876559,15.074,辻
,3,3579,0.08414561753,15.004,辻
,4,1628,0.082794,16.37,服部
,5,1016,0.083033,18.89,服部
,6,839,-0.082882,16.26,梶原
,7,554,-0.080871,18.28,梶原
,8,285,0.079995,19.20,工藤
,9,261,0.078980,20.22,工藤
,10,232,0.081911,17.26,佐々木
,11,208,0.075676,23.56,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***2023 11月29日 創造工房実習 [#r1eacb71]
異方性一次 [#p29b4460]
CENTER:異方性一次
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.7,144563,0.505252,2.76,安藤
,0.8,141517,0.504692,2.64,安藤
,0.9,91648,0.502595,2.216,兼田
,1.1,27160,0.489914,0.363,兼田
,1.2,24675,0.487088,0.791,柴田
,1.3,23446,0.4868010,0.995,柴田
,1.4,17738,0.485999,1.16,佐藤
,1.5,15438,0.485180,1.33,佐藤
,1.6,15900,0.483286,1.71,皆川
,1.7,12142,0.477952,2.80,皆川
,1.8,11604,0.482085,1.9554,永山
,1.9,10391,0.470887,4.2329,永山
,2,10291,0.480910,2.19,辻
,3,2328,0.431937,12.15,辻
,4,1500,0.430156,12.52,服部
,5,432,0.282968,42.45,服部
,6,356,0.3441556,30.00,梶原
,7,196,0.213934,56.49,梶原
,8,104,0.229874,53.25,工藤
,9,81,0.232308,52.75,工藤
,10,78,0.203271,58.65,佐々木
,11,63,0.222316,54.78,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
等方性二次 [#vb0def3f]
CENTER:等方性二次
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.7,144563,0.430124,3.22,安藤
,0.8,141517,0.430132,3.22,安藤
,0.9,91648,0.430020,3.197,兼田
,1.1,27160,0.429828,3.151,兼田
,1.2,24675,0.429836,3.15,柴田
,1.3,23446,0.42974,3.13,柴田
,1.4,17738,0.429797,1.3,佐藤
,1.5,15438,0.429958,3.14,佐藤
,1.6,15900,0.429755,3.18,皆川
,1.7,12142,0.429676,3.11,皆川
,1.8,11604,0.429829,3.1507,永山
,1.9,10391,0.429684,3.1159,永山
,2,10291,0.429620,3.10,辻
,3,2328,0.429169,2.99,辻
,4,1500,0.429254,3.01,服部
,5,432,0.428170,2.75,服部
,6,356,0.428452,2.82,梶原
,7,196,0.42591,2.21,梶原
,8,104,0.426074,2.25,工藤
,9,81,0.425552,2.12,工藤
,10,78,0.488382,17.20,佐々木
,11,63,0.423972,9.0534,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***2023 11月24日 創造工房実習 [#eab4b630]
単純梁~
,メッシュの長さ,要素数,変位[mm],相対誤差,計算者
,0.7,145234,0.422484,1.388,安藤
,0.8,142973,0.422570,1.409,安藤
,0.9,91648,0.420437,0.897,兼田
,1.1,27160,0.405618,2.659,兼田
,1.2,24675,0.404349,2.96,柴田
,1.3,23446,0.404185,3.00,柴田
,1.4,17738,0.398604,4.34,佐藤
,1.5,15438,0.396593,4.83,佐藤
,1.6,16122,0.398212,4.44,皆川
,1.7,12026,0.393411,5.59,皆川
,1.8,11604,0.393668,5.53,永山
,1.9,10391,0.390695,6.24,永山
,2,10921,0.395103,5.18,辻
,3,2328,0.324762,22.06,辻
,4,1500,0.155013,62.80,服部
,5,432,0.065278,84.33,服部
,6,357,0.213062,48.87,梶原
,7,196,0.1019,75.55,梶原
,8,104,0.1158624,72.20,工藤
,9,81,0.1255118,69.88,工藤
,10,78,0.07733,81.44,佐々木
,11,63,0.1999,52.03,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***2023 11月17日 創造工房実習 [#xa1780bc]
Salomeの使い方を学んだ。~
Geometry~
Mesh作成~
Aster Study~
Pravis~
片持梁~
,メッシュ長さ,要素数,先端変位(4隅の平均値)[mm],相対誤差,...
,0.7,198464,6.54281,1.91,安藤
,0.8,113812,6.5104,2.39,安藤
,0.9,40280,6.3631525,4.60,兼田
,1.1,30055,6.3363525,5.00,兼田
,1.2,26467,6.3043375,5.48,柴田
,1.3,25180,6.304355,5.48,柴田
,1.4,32212,6.31612,5.31,佐藤
,1.5,17753,6.1209,8.23,佐藤
,1.6,14296,6.2044625,6.98,皆川
,1.7,13596,6.2156625,6.81,皆川
,1.8,2866,5.737755,13.98,永山
,1.9,6001,5.7263625,14.15,永山
,2,5617,5.6458525,15.355,辻
,3,2309,5.4728755,17.948,辻
,4,617,3.6160575,45.12,服部
,5,494,3.8580375,42.16,服部
,6,581,2.50682,62.416,梶原
,7,133,1.41225,78.82,梶原
,8,78,1.2887175,80.68,工藤
,9,72,1.2879925,80.69,工藤
,10,60,1.14344,82.85,佐々木
,11,65,1.23124,81.154,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***2023 11月10日 創造工房実習 [#fb95bbcf]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/minakawa/...
***2023 10月27日 創造工房実習 [#fb5ee35c]
コマンドを学んだ。次回までに使いこなせるように頑張りた...
pwd 今いる場所がわかる~
ls ファイルを表示する~
mkdir ファイル名 ファイルの作成~
rmdir ファイル名 ファイルの削除~
cd ファイル名 ファイル名のところへ移動~
vi ファイル名 編集~
esc : w 保存~
esc : q 戻る~
cp コピー元 コピー先 コピーする時~
皆川恭輔
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
