皆川の修論日誌

皆川の卒論日誌
FrontPage

4/21 IHIインフラシステムさん　アドバイス†

1,サグモデル
両端固定にして自重をかけることで張力導入できる。
この方法により、斜めにしても張力導入でき、固定点間距離を一定にして解析できる。

2,ケーブルを１本ではなくストランドでモデル化してはどうか？ ストランドが束なったモデルと１本の円柱モデルでの比較

4/22　神鋼鋼線さん　学んだこと†

・被覆は、直接ケーブルにポリエチレンを溶かしながら覆っていくため、ケーブルと被覆の間には空間はない。
・基本的に被覆の厚さはケーブルの種類(直径・本数)で決まっている。
・レインバイブレーションは風速が10数m/s以上で発生しづらくなる。
　→水道が破壊されてしまうから。
・張力測定に用いられる高次振動法における振動数と振動モードの対応関係は、架設時に測定しておく

サグの影響の検討†

これまでのサグは、ケーブルモデルを真っ直ぐに作成し重力を与えてサグを作っていたが、この方法だとサグによる影響が出ていない可能性がある。
そこで、モデル作成時にサグの形状に曲げて作っておき、張力を与え、その後重力を載荷することでサグを形成させるという手順で解析を行う。

ある一定のサグ比になると、サグなしモデルと比較して固有振動数が変化するという論文もあるので、まずはサグを作った状態で振動解析を行う。

サグモデルの作成†

サグモデルを作成するために、ケーブル直線モデルに張力と重力をかけてサグを作り、各点の座標を読み取る。
その座標に従ってgeometryで点を打ち、曲線状にケーブルを作成する。

直線モデルに張力&重力をかけた時のサグの大きさ†

l=113.149m、d=152mm、T=1003.69kN、ρ=7800kg/m^3、ポアソン比=0.3のときのサグの理論値は2.16m
直線モデルに張力&重力をかけた時のサグの大きさは2.25m(理論値との相対誤差は4.2%)

卒論ではSolid要素を使用して振動解析を行っていたが、今後はcable要素を使っていく予定であるためcable要素でも振動解析を行えるかを確認してみる。

直線モデルでの振動解析†

①step数0.001での振動解析

固有振動数の理論値は、1次=0.38Hz、3次=1.13Hz、5次=1.88Hz
どの卓越振動数がどのモードのものかがわからないようなグラフになっている。もしかしたら、step数の影響で雑音が入っているかも

②step数0.01での振動解析

step数0.001のときに出ていた雑音のような卓越振動数がなくなり、理論値に近い値が出力されている。
また、solid要素で解析していたときは、1次=0.391Hz、3次=1.17Hz、5次=2.05Hzであり、それにも近い値が出力されるようになった。

③直線モデルで重力かけない(サグを作らない状態で)振動解析を行った結果が下の表　step数0.01

3次と思われる卓越振動数以外は、重力ありの振動解析と同じ振動数になっている。

次はsin波をあたえて卓越振動数で共振するかを確認する。

直線モデル　sin波†

振動解析で卓越した振動数のsin波を与えて共振するか確認する。
モデルには重力をかけない状態でsin波を与えている。

0.390625Hzのsin波

1.07422Hzのsin波

1.75781Hzのsin波

この結果から、由利橋のケーブルの固有振動数は1次モード=0.390625Hz、3次モード=1.07422Hz、5次モード=1.75781Hzであることがわかった。
理論値との誤差は、1次モード=3.7%、3次モード=5.3%、5次モード=6.4%
卒論ではsolid要素(3D)で振動解析を行っていたが、cable要素(1D)でも振動解析を行えることがわかったので、次はサグモデルを作成していきたい。

由利橋　ボルト事故†

5/2〜　モデル作成†

5/2〜5/7　主塔側のプレートとボルトのモデル化完了
→エラー出ない。
5/8～5/11 桁側のプレートとボルトのモデル化完了
→エラー出ない。
5/12～　 ケーブルのモデル化作成
→接点を共有させて問題なく解析できるか確認していく。
5/22〜　　モデルのmesh
→被覆部分のmesh部分でエラーが起きてしまい色々試している。
要素数が大きくなってしまっている。

主塔側のモデル

桁側のモデル

5/29〜　モデル完成(被覆部分はmeshを切れなかったため、まずはケーブルには被覆をつけない状態)
LIAISON_MAILを使って各部材を結合したところ、エラー出ずに問題なく行うことができた。(結合の仕方はOK)
3DModelでサグを再現する設定がまだ見つけられておらずケーブルは直線にしている。

　解析†

6/1〜　重力のみ考慮
今回のモデルに重力を考慮したところ、ケーブル部分のZ方向最大変位は1.15mであった。(この点でのx方向変位は0.51mであり、xz方向の変位量は1.26m)
保護管長106.367m、張力1003.69kNのときのサグの理論値は1.95m
この点では、SIGMxx=36.8MPa、SIGMzz=6.9MPaであり、軸力は37.4MPaとなり、目標の値より小さくなっている。(1003.69kNは55.3MPa)

PRE_SIGM†

3Dでサグモデルの作成前に両端固定したまま張力導入する方法を模索
PRE_SIGMコマンドは動的解析には使えないので２段階解析で適用してみる。(1段階目に静的でPRE_SIGM、2段階目に動的解析でインパクト)
両端固定した状態の結果をFFTにかけたところ、3.91Hz、7.12Hz、11.13Hzで卓越した(sin波を与えないとモード数との関係がわからない。)
卒論のときのFORCE_FACEで与える方法のときはf1=2.93Hz,f3=8.60Hz,f5=13.96Hzであった。

全体モデルのモーダル解析†

塔の固有振動数を求めるためにモーダル解析を行った。
・面内振動
1次=1.465Hz　3次=5.786Hz　5次=15.744Hz
・面外振動
1次=0.668Hz　3次=4.137Hz　5次=11.471Hz

全体モデルの時刻歴解析†

何回解析を行っても下の一時的なエラーがたくさん出て、フリーズしてしまう。
次はメモリを下げて解析してみる。

閲覧数†

最新の100件