柴田真杜の卒論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
柴田真杜
**1月15日実験と解析 [#s0238563]
ケーブルの断面、ヤング率が一様で解析を行ったところ、ケー...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
分割数を100,40,20,10と変更しても変わらない。~
実験と同じくケーブルをアラミドとばねに分けて解析を行った。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ばねのヤング率を求め、アラミドのヤング率を計算で求めた。~
ケーブル張力の変化量が小さすぎ~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
実験と解析で振動数が一致しないため、ケーブルを取り除いた...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
結果は5Hzであり理論値も5Hzであったことからケーブル要素が...
実験はiphoneのカメラで撮影をしたが撮影した動画のfpsが平均...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**ケーブルをたるませてから引っ張る [#ycb3dafc]
モデルと結果~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
死荷重と同じ大きさ(24.5N$\times$2)より大きい力を上向きに...
3本モデルでケーブルが破断したときのケーブル張力と、２本モ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**4本モデル、6本モデル、8本モデルの固定部応力の比較 [#a20...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**ケーブル断面を上下で変えた [#h2ede71c]
初期張力はそのままで、ケーブル断面を上下で変更することで...
cable1,2の断面積：11.64156428mm2~
cable3,4の断面積：4.656625711mm2~
初期応力を一致させることで挙動が正常になった。~
*ケーブルのヤング率をすべて100MPa、張力はcable,2を50N、ca...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*ケーブルのヤング率を全て150MPa、張力はcable1,2を50N、cab...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**HHTでも挙動がおかしくなった [#j92e2b3e]
ケーブルのヤング率をすべて100MPa、張力はケーブル1,2を50N...
分割数100 ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
HHT法でも張力がケーブルによって異なると挙動がおかしくなる。
cable4の張力~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**主塔の最大変位と固定部の最大応力の関係を2本モデルと4本...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**主塔先端の最大変位と主塔固定部の最大応力の関係をcable要...
ケーブルの初期張力を10N,20N,30N,40N,50Nにしたときの主塔先...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
時間積分アルゴリズムによる違いはあまりない。pou_d_em要素...
マルチファイバー要素とcable要素(HHT-oui)で幾何学非線形の...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブルの初期張力が上がるに連れグラフの傾きが変わったこ...
**cable要素のnewmark-βとHHT、マルチファイバー要素でcable...
変位と固定部の応力を比較した時、初期張力のとき50Nで大きな...
そこで初期張力50Nのときの破断時の各ケーブルの張力の変化を...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**11月26日 [#ofb1b2ec]
newmark-β法とHHT法のOUI,NONの3つの解析方法で、最大先端変...
予想として、newmark-β法では非線形問題が不安定になりやすい...
**陽解法と陰解法 [#f728f4ad]
[[陽解法と陰解法の違いについて:http://jikosoftcom.blog25....
陽解法は衝撃など動的な現象を得意とし、陰解法は静的な現象...
ケーブルの破断問題では陽解法の方がいい?
陽解法ではtimestep<=2/構造物の最大固有振動数とする必要が...
**11月26日 [#q3e3d99c]
HHT法で解析を行った時、減衰が小さいと破断後5波目くらいの...
残ったケーブルの張力がどのように変化しているのか調べた。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
NON、OUIどちらもpoint1で張力が破断直後ではなく5波目くらい...
**HHT法、MODI_EQUI-OUIのαの値について~ [#b2672f52]
MODE_EQUIーOUIのときのαを変えた解析結果を比較すると違いは...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
αの値が影響するのはMODI_EQUIーNONのときであり、αの絶対値...
**11月25日 [#m77337ea]
HHT法にすることで減衰が大きくかかった。そこで主塔に今まで...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
多少の違いはあったが主塔の減衰をなくしても減衰はかかって...
**HHT法のMODI_EQUI-OUIとNONの違いについて [#u039d716]
salomeのマニュアルにはOUIは完全なHHT法、NONは修正平均加速...
この違いを明らかにするため、αの値を変えずOUIとNONを比較し...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
OUIにすることでNONよりも減衰が小さくなっている。~
OUIの破断後の挙動が急に変化している部分がある。~
**newmark-β法とHHT法について [#u8446c9f]
newmark-β法:[[newmark-β法の説明:https://ocw.u-tokyo.ac.jp...
物理量として変位、変位速度、加速度を考慮した 1ステップ（...
HHT法:動的応答では、低次の大域的な応答が重要であることが...
時間ステップ i, i + 1 の寄与を人工的なパラメータ α で調整...
salomeでは(-0.3<=α<=0)?確かではない、、~
salomeではHHT法の選択時にMODI_EQUI-OUI,NONを選択するとこ...
この違いはあまりよくわかっていないが、OUIは完全なHHT法、N...
**11月21日 [#z637e90b]
時間積分アルゴリズムの選択は、数値解法の精度や安定性の根...
これまではnewmark-β法を用いた解析を行ってきたが、非線形解...
HHT法による解析をすることで分割数を増やしても安定した結果...
解析モデル~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
,,①,②,③,④
,newmark-β,○,☓,○,☓
,HHT法,○,○,○,○
cable3破断のように安定しやすい（高周波の振動が少ない)モデ...
cable1破断のような安定しにくい(高周波の振動が多い)モデル...
newmark-β法とHHT法どちらも正常な挙動を得られた③において先...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
どちらも正常な結果が得られている③ではnewmark-β法、HHT法に...
HHT法で正常な結果が得られた④の先端変位~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
減衰が大きくかかりすぐに収束している。この原因として、HHT...
この結果は主塔に減衰をかけたものである。今後は主塔に減衰...
破断直後の変位よりそれより少し遅れた変位が最大になってい...
④の収束変位と静的解析による静的変位を比較した。~
収束変位:24.24mm~
静的変位:22.44mm~
収束変位と静的変位が一致していない。~
この原因が幾何学非線形を考慮したため衝撃の影響が現れたの...
**11月13日 [#l53b00d7]
下モデルで破断ケーブル1(18)と破断ケーブル3(26)において主...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**cable要素とマルチファイバー要素の比較 [#s9699d04]
２本モデルの実験をもとに、cable要素とマルチファイバー要素...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
２本モデルにおいては、cable要素とマルチファイバー要素に大...
**11月11日 [#q77810b9]
ケーブル要素にしたときの分割数が先端変位、張力に及ぼす影...
結果として分割数は先端変位、張力に影響を及ぼさないことが...
破断モデル~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
分割数を1、2、10、50、100にした時の先端変位とケーブル2の...
*cable1破断の場合 [#oeade346]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*cable3破断の場合 [#xc73ed66]
先端変位
...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**10月28日 [#r2b03967]
ケーブル主塔positiveの動的解析による収束変位と静的解析に...
,モデル,解析,変位
,2本モデル,静的解析,18.1808
,２本モデル,動的解析,18.3353
,,,
,4本モデル(47),静的解析,44.3
,4本モデル(47),動的解析,44.3
,,,
,4本モデル(38),静的変位,18.1894
,4本モデル(38),動的変位,収束しない
4本モデル(38)の動的解析 減衰定数0.3と大きくしたが収束して...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
47モデルでは一致しているがその他は一致していないため現状...
**10月25日 [#y91753ee]
主塔にも幾何学非線形を考慮する必要があるのか。~
まずは2本モデルで実験値、ケーブルのみpositive、ケーブル＋...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
収束する値が実験値と近いのはケーブルのみpositiveのときで...
静的解析を行い釣り合ったときの変位（静的変位）は、ケーブ...
したがって収束する値は18.335に近いはずだがケーブル主塔pos...
静的解析で釣り合い時の横方向変位、縦方向変位をみた~
,,横方向変位,縦方向変位
,ケーブルのみpositive,13.4464mm,0.0082mm
,ケーブル主塔positive,18.1894mm,0.24mm
主塔にも幾何学非線形を考慮することで変位は変化したが、縦...
**10月22日 [#y645c0a6]
2本モデルでcable要素を使用し、解析を行った。ケーブルの分...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
2本モデルにおいてケーブルの分割数は影響がなかった。~
次に4本モデルにおいてcable要素を使用して解析を行ったとこ...
cable1と2の材料特性、3と4の材料特性がそれぞれ同じためこれ...
4本モデル(破断ケーブル3）だと正常に回るようになったが、４...
**cable要素 [#gc41eef3]
ケーブルに圧縮力が生じないようにする方法を調べたところ、...
cable要素では圧縮時のヤング率を設定することができ、今回は...
cable要素の設定を進める上で幾何学非線形をいれることになる...
**幾何学非線形 [#d7b616f5]
幾何学非線形を考慮したい理由：幾何学非線形を考慮すること...
また、変形にともなった剛性の変化を考慮する必要があると考...
**10月22日 [#i72b4d18]
幾何学非線形positiveとnegativeの先端変位を静的解析と計算...
,,静的解析,計算
,positive,7.8mm,7.6mm
,negative,8.67mm,8.42mm
幾何学非線形を入れたほうが先端変位が小さくなることを確認...
**幾何学非線形の入れ方は正しいのか [#xdff88ba]
2本モデルで幾何学非線形を入れたとき(positive)と入れなかっ...
先端変位とケーブル張力~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
幾何学非線形をいれた場合、破断後のケーブル張力が収束せず...
幾何学非線形の入れ方を間違っている可能性がある。~
**10月11日 幾何学非線形を考慮した解析をしてみる。 [#p9ec...
その理由としては、これまでの解析ではケーブルの初期張力が...
幾何学非線形を考慮することで、初期張力が大きくなる<と、引...
まず、幾何学的非線形を考慮していないもの(negative)と考慮...
[[条件:https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/s...
先端変位~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断時のケーブル張力の変化量~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
動き
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
幾何学非線形を考慮したためケーブルの引張の剛性が高くなり...
またケーブルが曲がりにくくなり張力の変化量も小さくなった?~
結果に影響が出たことから、解析方法の見直しが必要。~
これまで２本モデルでは実験と解析一致していたが、幾何学非...
**10月9日 [#y87cb38a]
4本モデルで初期張力がケーブル破断時の張力に与える影響をわ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブル破断時の張力の変化量に違いがでるのは、ケーブル1,2...
ケーブル1,2と3,4で初期張力が異なるものでも張力の変化量が...
**10月7日 [#gf270028]
ケーブル破断時の張力への影響を初期張力、ヤング率、破断す...
[[条件:https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/s...
主塔の材料特性~
寸法:5mm $\times$ 10mm $\times$ 900mm~
密度:7.668$N/mm^{2}$ $\times 10^{-9}$~
ケーブルの材料特性~
直径:3.85mm
密度:3.13$N/mm^{2}$ $\times 10^{-10}$~
主塔には全てに同じ減衰を与える（Salomeの材料特性の部分のα...
ケーブル番号~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断時のケーブル張力を読み取るpoint~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
まずは初期張力の条件を変え、それ以外の条件を同じにしたも...
比較するものの初期張力を1にする。（9では初期張力20Nなので...
・２本モデル
''9と13、10と14、12と16''を比較
,,初期張力,ヤング率,サグ
,9,20N,50MPa,無
,13,50N,50MPa,無
,,,,
,10,20N,50MPa,有
,14,50N,50MPa,有
,,,,
,12,20N,100MPa,有
,16,50N,100Mpa,有
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
この結果より２本モデルにおいては、初期張力は破断後のケー...
・４本モデル~
''17と21、18と22、20と24、26と30、28と32、51と55、59と63...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
1〜4のケーブルで初期張力がそろっているものは(17など）初期...
しかし1〜4のケーブルで初期張力がそろっていないものは（34...
**9月19日 6本モデル [#r893c93b]
6本モデル ケーブルのヤング率はすべて50N/mm2に設定~
張力 上:28.70N 中:20.09N 下:11.48N~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断１~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断２~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
4本モデルの破断１と同じ挙動~
破断３~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
４本モデルの破断２と同じ挙動~
**9月26日 張力の入れ方は同じでヤング率を変えてみる~ [#i59...
張力 上:28.71N 下:8.61Nは同じまま、ヤング率を上ケーブル80...
張力の変化量~
,ue_hidari,ue_migi,sita_hidari,sita_migi,(上の張力÷下の張...
,4.608NN,-4.608N,-,-2.201,2.094
ケーブルのヤング率が変わることで破断時のケーブル張力の変...
**9月18日 ヤング率を同じにする [#we683b9d]
*破断1 [#zbaea865]
ヤング率を全てのケーブルで50N/mm2と同じにしたときの塔の動...
張力 上:28.71N 下:8.61N~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
張力 上:14.35N 下:20.10N~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ヤング率を同じにしても動きに違いはなかった。~
張力の変化量~
,,ue_hidari,ue_migi,sita_hidari,sita_migi,(上の張力÷下の...
,張力 上:28.71N 下:8.61Nのとき,3.71N,-3.70N,-,-3.50N,1.06
,張力 上:14.35N 下:20.10Nのとき,8.65N,-8.65N,-,-8.15N,1.06
,張力 上:14.3526N 下:14.3526Nのとき,5.826N,6.18N,-,6.17N...
ケーブルのヤング率を揃えることで張力の変化量比に影響が出...
ヤング率を100N/mm2に揃えたときの張力の変化量比は、1.03と...
ヤング率を揃えると張力の変化量比は1.0に近づくと考えられる...
*破断2 [#ve53e7c8]
ヤング率を50N/mm2に揃えた時、張力の変化量比は1.06になった...
ヤング率を100N/mm2に揃えた時、張力の変化量比は1.11になっ...
このことから、破断2の場合はヤング率を揃えると1.11に近づく...
**9月17日 張力の大きさを上下逆にする[#d441503d]
9月11日の実験モデルでケーブル張力の大小を上下逆にしてみた~
張力 上:14.4081N 下:27.0556N~
破断１~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
張力の大小を逆にしても違いはなかった。~
破断２~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
張力の大小を逆にしても動きに違いはなかった。~
破断１と２の破断後のケーブルの張力~
破断１~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断２~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
張力の変化量~
,,ue_hidari,ue_migi,sita_hidari,sita_migi,(上の張力÷下の...
,破断１,14.7217N,-14.7136N,-,-6.6275N,2.218
,破断２,-,-12.8582N,5.4554N,-5.4553N,2.357
張力の大きさを上下で逆にしても、上下での張力の変化量の比...
**9月11日の実験 [#d67e2b24]
実験モデル~
張力 ケーブル上:24.0161N ケーブル下:15.7814N
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
今までは先端変位を測定してDAFを求めてきた。しかし本数を増...
今後は破断時の塔の動きを追うことで、ケーブル本数を増やし...
モデル破断１の塔の動きは次のようになった。(ケーブル要素は...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブルあり（scale1とscale30)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
初期位置に戻ってきたときの塔の形~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
モデル破断2の塔の動きは次のようになった。(ケーブル要素は...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブルあり(scale1とscale10)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
初期位置に戻ってきたときの塔の形~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断１と２の破断後のケーブルの張力~
破断１~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断２~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断時の張力の変化量~
,,ue_hidari,ue_migi,sita_hidari,sita_migi,(上の張力÷下の...
,破断１,8.5869N,-8.5872N,-,-3.8658N,2.221
,破断２,-,-21.43226N,9.0947N,-9.09228N,2.358
ケーブル１本モデル（scale foctor 100）~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブル２本モデル(scale factor 400)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブル１本、２本、４本のときで動きが異なることから、本...
**8月1日の実験と解析の結果 [#k9f0f58b]
*ケーブル１本の時 [#d61893f5]
解析は健全時の変位がほぼ一致するように温度荷重を与えた。~
実験と解析の結果~
,,静的,動的,健全時,DAF
,実験,0,30.01,57.99,1.932
,解析,0,30.00,59.20,1.973
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
実験と解析のグラフ~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブル破断時の周期が解析の方が速い。~
*ケーブル2本の時 [#id80de6e]
解析は静的変位がほぼ一致するように温度荷重を与えた。~
初めは温度荷重によりケーブル張力を求め、片方のケーブルを...
次に、千代岡さんから教えていただいた２段階解析をした。~
2段階解析では温度荷重はできない？エラーがでてしまった。
そのため、ケーブル２本とも「あたかもあったかも」の方法を...
この方法で減衰はかかった。（この方法を①とする。）~
実験と解析の結果~
,,静的,動的,健全時,DAF
,実験,18.34,34.54,0,1.884
,解析,18.40,35.43,0,1.926
実験と解析のグラフ~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブル破断時の周期が実験の方が速い。~
DYNA_VIBRAを使用し、片方を温度荷重、もう片方を同じ張力で...
つまり、温度荷重と点載荷の両方を使用するときはDYNA_VIBRA...
DYNA_NON_LINEを使用して片方温度荷重、もう片方を点載荷にし...
DYNA_VIBRAを使用しているときはcalc_modeとcalc_amor_modal...
そこで材料設定の部分で減衰に関わるα、βを設定し減衰をかけ...
（この理由はケーブルのみにα、βを設定していたから）~
このときの実験と解析のグラフは次のようになった。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
このグラフでは周期のズレは小さくなっている。このことから...
方法①と②のモデル~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
方法②を解析の途中でtime stepを変える解析してみた。(0s〜19...
またα、βは梁のみに設定して解析を行った。~
このときの実験と解析のグラフは次のようになった。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
,,静的,動的,健全時,DAF
,解析,18.40,35.79,0,1.945
**片持ち梁の動的解析について [#e65d6f36]
[[FrontlSTRチュートリアルP85〜:https://manual.frontistr.c...
これによると片持ち梁のDAFの理論値は1.9753となる。~
断面10mm*50mm,長さ500mmの片持ち梁モデルで動的解析を行って...
まずはステップ荷重を梁の先端に載荷する解析を行った。しか...
Exception utilisateur levee mais pas interceptee. ...
! Les bases sont fermees. ...
! Type de l'exception : NonConvergenceError ...
! ...
! Arrêt pour cause d'absence de convergence avec le...
! dans l'algorithme non-linéaire de Newton. ...
! La base globale est sauvegardée. Elle contient le...
! l'arrêt. ...
! ...
! Conseils : ...
! - Augmentez ITER_GLOB_MAXI. ...
! - Réactualisez plus souvent la matrice tangente. ...
! - Raffinez votre discrétisation temporelle. ...
! - Essayez d'activer la gestion des événements (dé...
! exemple) dans la commande DEFI_LIST_INST.
次に「あたかもあったかも」を用いてtimestep0.001で解析を行...
理論値と考えていた1.9753にはならなかった。~
**塔の長さ、境界条件とDAFの関係 [#w1cc0437]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
,長さ(mm),片持のDAF,単純のDAF
,100,1.258,1.0962
,200,1.7225,1.36
,300,1.9082,1.6528
,400,1.9609,1.8328
,500,1.9791,1.9148
,600,1.9878,1.9527
,700,1.9884,1.9701
,800,1.9915,1.9789
,900,1.9944,1.9858
,1000,1.9947,1.9889
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
・境界条件：片持ばり~
塔の長さを変えながら先端に水平方向の荷重を与えた。~
塔の長さを長くすることで先端が境界条件を考えない自由な動...
・境界条件：単純梁~
DRYのみを自由にした単純梁モデルを作り、梁の長さを変えなが...
片持ち梁と同様に長くすることで中央が境界条件を考えない自...
・両方を比べてみて
片持ち梁の方が同じ長さ時のDAFが大きく、また長さの伸びに対...
しかし、固定部からの距離で見る（片持ちL=100mmと単純L=200m...
この理由を考えながら、境界条件がどのようにDAFに影響するか...
**ケーブルの角度とDAFの関係 [#c71d2581]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
塔の長さ、荷重は同じにして、ケーブルの角度のみを変化させD...
,,case1,case2,case3,case4
,DAF,1.88059,1.88059,1.88059,1.88059
この結果からケーブルの角度はDAFに影響しないことが分かった...
**7月2日の実験 [#me6f7b2c]
・２本の場合(減衰なし)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断していないケーブルがあるため、初期位置より破断したケ...
**７月２日の実験 [#vef2bead]
7月2日の実験ではタワーにケーブルを１本取り付けたものと２...
・１本の場合~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
１本のときは上のような状態からケーブルを切った。~
実験と解析の先端変位のグラフは次のようになり、近い値とな...
しかし、解析の方が周期が早く見える。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
time stepを0.01から0.005にすることで周期の問題は解決した...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
DAFは1.897となり、理論値に近づいた。~
time stepを0.005から0.001にすることでDAFが1.919となりさら...
この原因はまだわかっていない。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
実験と解析結果~
,,動的変位,静的変位,DAF,DAF理論値,
,実験,31.76mm,17.70mm,1.794,1.963,
,解析step0.01,32.84mm,17.66mm,1.859,1.963,
,解析steo0.005,33.50mm,17.66mm,1.897,1.963,
,解析step0.001,33.90mm,17.66mm,1.919,1.963,
**7月10日 [#z401a45c]
７月２日に行った実験の再現を解析で行っている。~
ケーブル部分にプレストレスを与えるため、温度荷重のかけか...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
温度荷重によってプレストレスが入り、塔がたわむことを確認...
今後はたわんだ状態からケーブルを破断させたときの挙動を確...
**減衰について [#f0235647]
これまでの減衰定数の計算の考え方は間違っていた。~
https://wwwra.meijo-u.ac.jp/labs/ra007/murata/pdf/textboo...
このページを参考にした。~
$δ=\displaystyle \frac{Y_1}{Y_2}$で対数減衰率δを求めるが...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
この減衰定数を計算すると0.00477と与えたかった0.005に近い...
・減衰定数を0.05に設定し、解析を行った。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
このグラフからDAFを求めると1.850で理論値のDAF=1.855に近い...
・減衰定数0.1に設定した場合~
理論値が1.730なのに対し、解析から求めたDAFは1.725となった...
**7月3日 [#s396f2bf]
構造関係メモから$DAF=1-\cos\omega t$と表すことが出来るこ...
減衰についてだが、現在自分の与えたい減衰定数を反映するが...
減衰定数0.005で与えたいところを、変位のグラフから計算する...
**3本のばねモデルでの実験結果 [#jf53c381]
,,実験値,解析値
,3本静的,32mm,32.48mm
,破断時(動的),68mm,65.15mm
,2本静的,47mm,48.72mm
--動的増幅率 $DAF={\displaystyle \frac{S_{dyn}-S_{0}}{S_{...
---$S_{dyn}$: 動的解析の最大応答値
---$S_{stat}$: 静的解析の応答値
---$S_{0}$: 初期値
健全時を原点にとり、動的変位と静的変位の比をとったものと...
健全時を原点にとったものを$DAF_{健}$、自然長を原点にとっ...
構造関係メモから自然長を原点にとったDAFは$DAF={1-\display...
tは破断開始から最大変位にいたるまでの時間を代入する。~
解析、実験の結果は次のようになった。~
,,,$DAF_{健}$,$DAF_{自}$,$DAF_{理}$,$DAF_{自理}$,
,,解析,2.01,1.337,2.00,1.333,
,,実験,2.40,1.447,2.00,,
**6月26日 [#td51e133]
固有振動数を解析で求め、レイリー減衰で用いるα、βを計算し...
減衰定数0.005~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
減衰定数0.005よりもかかっている気がする、、
解析方法を見直す必要あり
**6月24日 [#zc7a53dc]
実験した3本モデルを減衰を考慮して解析してみた。~
グラフを見ると減衰しすぎな気がする。なにか設定をミスって...
減衰0.5%のとき~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
減衰1%のとき~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
減衰5%のとき~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**6月19日 [#k44de0d4]
断面、縦22mm、横70mmで高いところから紐+バネで吊るした解析...
これを踏まえ断面、縦70mm、横22mmの縦長にし、低いところか...
今回の実験の結果からばね定数等を求め解析を行い、実験と解...
,,実験値,解析値
,3本静的,32mm,32.48mm
,破断時,68mm,65.15mm
,2本静的,47mm,48.72mm
動的変位と2本の静的変位の比から求めるDAFは次のようになっ...
,,実験値,解析値
,DAF,1.447,1.337
今回のバネの本数が破断により3本から2本になったため、DAFの...
DAFを衝撃係数問題から求めると理論値は２になり、実際に解析...
**6月11日 [#e33f392d]
5月31日に実験をIHIや他大学の教授に見せ意見をいただいた。~
初めの実験モデルでは複数のパラメータにより、調べたい対称...
そこで簡単なモデルから徐々に進めていくことにした。~
初めは3本のケーブルで鋼材を吊るしたモデルで、真ん中のケー...
千代岡さんに"あたかもあったかも"を教えていただき、salome...
3本のケーブルで真ん中の1本を破断させたときの時間と変位の...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
この結果に加え、ケーブルが2本のときの静的変位を調べること...
（DAFの導き方は千代岡さんの卒論を参考にした。）~
次にケーブルの剛性を変化させることでDAFの違いを求めようと...
そこでケーブルの本数を増やし、ケーブルの破断する数を変化...
ケーブルの本数を9本にして真ん中の1本を破断させたときの時...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
DAFは1.177となった。~
今後は3本破断、5本破断、7本破断の解析を進めていく。~
**5月29日 [#uab9dea7]
衝撃を与える時間を変えてみた
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
マルチファイバーに衝撃を与えた時の先端変位
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**シェル要素の片持ちばり [#mcfd08a0]
[[http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin...
[[http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin...
このページを参考にして解析した
AFFE_CARA_ELEMのわからない用語を調べる必要あり
EXCENTREMENT_FO： 法線方向で、格子と平均曲面の間の距離を...
[[AFFE_CARE_ELMEの説明:https://code-aster.org/doc/v14/fr/...
**ビーム要素の片持ちばり [#d39d49cd]
点を固定できるのか:
lineをグループ化で固定、載荷の点を作ったがAFFE_CHAR_MECA...
GROUP_MAを選択していたため点が出てこなかった、GROUP_NOを...
固定面をどうやってつくるのか:
長方形フェースを固定したい面に付ける方法を考えたがうまく...
**Midas [#ge1b6317]
https://in.midasit.co.jp/wp-content/uploads/2020/07/const...
https://in.midasit.co.jp/wp-content/uploads/2020/07/const...
*参考資料 [#u9b8ef7d]
[[code-asterの説明:https://code-aster.org/doc/v15/en/inde...
[[code-aster質問ページ:https://forum.code-aster.org/publi...
[[構造要素の概要:https://www.opencae.or.jp/wp-content/upl...
[[減衰について:https://www.cybernet.co.jp/ansys/learning/...
[[面内と面外の説明:https://bonperson-civil.com/in-out-pla...
[[cable軸力の計算?:https://code-aster.org/doc/v15/en/man_...
*12月8日 [#m89d8b57]
サンドイッチ梁
,メッシュの長さ,要素数,変位,相対誤差,計算者
,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,佐々木
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*11月29日 [#u96a8789]
異方性一次 [#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,佐々木
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
等方性二次 [#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,佐々木
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*11月24日 [#h4d56322]
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.7,1455234,0.422484,0.01388,安藤
,0.8,142973,0.422570,0.01409,安藤
,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,佐々木
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*11月17日 [#gb489580]
,メッシュ長さ,要素数,変位,相対誤差,計算者
,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,0.458,服部
,5,494,3.8580375,0.422,服部
,6,581,2.50682,62.416,梶原
,7,133,1.41225,78.827,梶原
,8,78,1.2887175,80.68,工藤
,9,72,1.2879925,80.69,工藤
,10,60,1.14344,82.85,佐々木
,11,65,1.23124,81.154,佐々木
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*11月10日 [#z1f8fbe7]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*10月27日 [#dcab77d9]
コマンドを学んだ。
pwd:今いる場所
ls:今いる場所のフォルダ名ファイル名等を見る
mkdir:フォルダの作成
cd:ファイル間の移動（ディレクトリーの移動)
gedeit ファイル名 .tex & :テキストファイルを開く
vi:ファイルを開く
Esc : q:開いたファイルを閉じる
cp 元のファイル名 新しいファイル名:元のファイルをコピー
cat:ファイルの中を確認
rm:ファイルの削除、cd .. :1個前に戻る
rmdir:フォルダを削除
終了行:
#contents
柴田真杜
**1月15日実験と解析 [#s0238563]
ケーブルの断面、ヤング率が一様で解析を行ったところ、ケー...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
分割数を100,40,20,10と変更しても変わらない。~
実験と同じくケーブルをアラミドとばねに分けて解析を行った。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ばねのヤング率を求め、アラミドのヤング率を計算で求めた。~
ケーブル張力の変化量が小さすぎ~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
実験と解析で振動数が一致しないため、ケーブルを取り除いた...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
結果は5Hzであり理論値も5Hzであったことからケーブル要素が...
実験はiphoneのカメラで撮影をしたが撮影した動画のfpsが平均...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**ケーブルをたるませてから引っ張る [#ycb3dafc]
モデルと結果~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
死荷重と同じ大きさ(24.5N$\times$2)より大きい力を上向きに...
3本モデルでケーブルが破断したときのケーブル張力と、２本モ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**4本モデル、6本モデル、8本モデルの固定部応力の比較 [#a20...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**ケーブル断面を上下で変えた [#h2ede71c]
初期張力はそのままで、ケーブル断面を上下で変更することで...
cable1,2の断面積：11.64156428mm2~
cable3,4の断面積：4.656625711mm2~
初期応力を一致させることで挙動が正常になった。~
*ケーブルのヤング率をすべて100MPa、張力はcable,2を50N、ca...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*ケーブルのヤング率を全て150MPa、張力はcable1,2を50N、cab...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**HHTでも挙動がおかしくなった [#j92e2b3e]
ケーブルのヤング率をすべて100MPa、張力はケーブル1,2を50N...
分割数100 ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
HHT法でも張力がケーブルによって異なると挙動がおかしくなる。
cable4の張力~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**主塔の最大変位と固定部の最大応力の関係を2本モデルと4本...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**主塔先端の最大変位と主塔固定部の最大応力の関係をcable要...
ケーブルの初期張力を10N,20N,30N,40N,50Nにしたときの主塔先...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
時間積分アルゴリズムによる違いはあまりない。pou_d_em要素...
マルチファイバー要素とcable要素(HHT-oui)で幾何学非線形の...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブルの初期張力が上がるに連れグラフの傾きが変わったこ...
**cable要素のnewmark-βとHHT、マルチファイバー要素でcable...
変位と固定部の応力を比較した時、初期張力のとき50Nで大きな...
そこで初期張力50Nのときの破断時の各ケーブルの張力の変化を...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**11月26日 [#ofb1b2ec]
newmark-β法とHHT法のOUI,NONの3つの解析方法で、最大先端変...
予想として、newmark-β法では非線形問題が不安定になりやすい...
**陽解法と陰解法 [#f728f4ad]
[[陽解法と陰解法の違いについて:http://jikosoftcom.blog25....
陽解法は衝撃など動的な現象を得意とし、陰解法は静的な現象...
ケーブルの破断問題では陽解法の方がいい?
陽解法ではtimestep<=2/構造物の最大固有振動数とする必要が...
**11月26日 [#q3e3d99c]
HHT法で解析を行った時、減衰が小さいと破断後5波目くらいの...
残ったケーブルの張力がどのように変化しているのか調べた。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
NON、OUIどちらもpoint1で張力が破断直後ではなく5波目くらい...
**HHT法、MODI_EQUI-OUIのαの値について~ [#b2672f52]
MODE_EQUIーOUIのときのαを変えた解析結果を比較すると違いは...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
αの値が影響するのはMODI_EQUIーNONのときであり、αの絶対値...
**11月25日 [#m77337ea]
HHT法にすることで減衰が大きくかかった。そこで主塔に今まで...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
多少の違いはあったが主塔の減衰をなくしても減衰はかかって...
**HHT法のMODI_EQUI-OUIとNONの違いについて [#u039d716]
salomeのマニュアルにはOUIは完全なHHT法、NONは修正平均加速...
この違いを明らかにするため、αの値を変えずOUIとNONを比較し...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
OUIにすることでNONよりも減衰が小さくなっている。~
OUIの破断後の挙動が急に変化している部分がある。~
**newmark-β法とHHT法について [#u8446c9f]
newmark-β法:[[newmark-β法の説明:https://ocw.u-tokyo.ac.jp...
物理量として変位、変位速度、加速度を考慮した 1ステップ（...
HHT法:動的応答では、低次の大域的な応答が重要であることが...
時間ステップ i, i + 1 の寄与を人工的なパラメータ α で調整...
salomeでは(-0.3<=α<=0)?確かではない、、~
salomeではHHT法の選択時にMODI_EQUI-OUI,NONを選択するとこ...
この違いはあまりよくわかっていないが、OUIは完全なHHT法、N...
**11月21日 [#z637e90b]
時間積分アルゴリズムの選択は、数値解法の精度や安定性の根...
これまではnewmark-β法を用いた解析を行ってきたが、非線形解...
HHT法による解析をすることで分割数を増やしても安定した結果...
解析モデル~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
,,①,②,③,④
,newmark-β,○,☓,○,☓
,HHT法,○,○,○,○
cable3破断のように安定しやすい（高周波の振動が少ない)モデ...
cable1破断のような安定しにくい(高周波の振動が多い)モデル...
newmark-β法とHHT法どちらも正常な挙動を得られた③において先...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
どちらも正常な結果が得られている③ではnewmark-β法、HHT法に...
HHT法で正常な結果が得られた④の先端変位~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
減衰が大きくかかりすぐに収束している。この原因として、HHT...
この結果は主塔に減衰をかけたものである。今後は主塔に減衰...
破断直後の変位よりそれより少し遅れた変位が最大になってい...
④の収束変位と静的解析による静的変位を比較した。~
収束変位:24.24mm~
静的変位:22.44mm~
収束変位と静的変位が一致していない。~
この原因が幾何学非線形を考慮したため衝撃の影響が現れたの...
**11月13日 [#l53b00d7]
下モデルで破断ケーブル1(18)と破断ケーブル3(26)において主...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**cable要素とマルチファイバー要素の比較 [#s9699d04]
２本モデルの実験をもとに、cable要素とマルチファイバー要素...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
２本モデルにおいては、cable要素とマルチファイバー要素に大...
**11月11日 [#q77810b9]
ケーブル要素にしたときの分割数が先端変位、張力に及ぼす影...
結果として分割数は先端変位、張力に影響を及ぼさないことが...
破断モデル~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
分割数を1、2、10、50、100にした時の先端変位とケーブル2の...
*cable1破断の場合 [#oeade346]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*cable3破断の場合 [#xc73ed66]
先端変位
...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**10月28日 [#r2b03967]
ケーブル主塔positiveの動的解析による収束変位と静的解析に...
,モデル,解析,変位
,2本モデル,静的解析,18.1808
,２本モデル,動的解析,18.3353
,,,
,4本モデル(47),静的解析,44.3
,4本モデル(47),動的解析,44.3
,,,
,4本モデル(38),静的変位,18.1894
,4本モデル(38),動的変位,収束しない
4本モデル(38)の動的解析 減衰定数0.3と大きくしたが収束して...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
47モデルでは一致しているがその他は一致していないため現状...
**10月25日 [#y91753ee]
主塔にも幾何学非線形を考慮する必要があるのか。~
まずは2本モデルで実験値、ケーブルのみpositive、ケーブル＋...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
収束する値が実験値と近いのはケーブルのみpositiveのときで...
静的解析を行い釣り合ったときの変位（静的変位）は、ケーブ...
したがって収束する値は18.335に近いはずだがケーブル主塔pos...
静的解析で釣り合い時の横方向変位、縦方向変位をみた~
,,横方向変位,縦方向変位
,ケーブルのみpositive,13.4464mm,0.0082mm
,ケーブル主塔positive,18.1894mm,0.24mm
主塔にも幾何学非線形を考慮することで変位は変化したが、縦...
**10月22日 [#y645c0a6]
2本モデルでcable要素を使用し、解析を行った。ケーブルの分...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
2本モデルにおいてケーブルの分割数は影響がなかった。~
次に4本モデルにおいてcable要素を使用して解析を行ったとこ...
cable1と2の材料特性、3と4の材料特性がそれぞれ同じためこれ...
4本モデル(破断ケーブル3）だと正常に回るようになったが、４...
**cable要素 [#gc41eef3]
ケーブルに圧縮力が生じないようにする方法を調べたところ、...
cable要素では圧縮時のヤング率を設定することができ、今回は...
cable要素の設定を進める上で幾何学非線形をいれることになる...
**幾何学非線形 [#d7b616f5]
幾何学非線形を考慮したい理由：幾何学非線形を考慮すること...
また、変形にともなった剛性の変化を考慮する必要があると考...
**10月22日 [#i72b4d18]
幾何学非線形positiveとnegativeの先端変位を静的解析と計算...
,,静的解析,計算
,positive,7.8mm,7.6mm
,negative,8.67mm,8.42mm
幾何学非線形を入れたほうが先端変位が小さくなることを確認...
**幾何学非線形の入れ方は正しいのか [#xdff88ba]
2本モデルで幾何学非線形を入れたとき(positive)と入れなかっ...
先端変位とケーブル張力~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
幾何学非線形をいれた場合、破断後のケーブル張力が収束せず...
幾何学非線形の入れ方を間違っている可能性がある。~
**10月11日 幾何学非線形を考慮した解析をしてみる。 [#p9ec...
その理由としては、これまでの解析ではケーブルの初期張力が...
幾何学非線形を考慮することで、初期張力が大きくなる<と、引...
まず、幾何学的非線形を考慮していないもの(negative)と考慮...
[[条件:https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/s...
先端変位~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断時のケーブル張力の変化量~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
動き
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
幾何学非線形を考慮したためケーブルの引張の剛性が高くなり...
またケーブルが曲がりにくくなり張力の変化量も小さくなった?~
結果に影響が出たことから、解析方法の見直しが必要。~
これまで２本モデルでは実験と解析一致していたが、幾何学非...
**10月9日 [#y87cb38a]
4本モデルで初期張力がケーブル破断時の張力に与える影響をわ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブル破断時の張力の変化量に違いがでるのは、ケーブル1,2...
ケーブル1,2と3,4で初期張力が異なるものでも張力の変化量が...
**10月7日 [#gf270028]
ケーブル破断時の張力への影響を初期張力、ヤング率、破断す...
[[条件:https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/s...
主塔の材料特性~
寸法:5mm $\times$ 10mm $\times$ 900mm~
密度:7.668$N/mm^{2}$ $\times 10^{-9}$~
ケーブルの材料特性~
直径:3.85mm
密度:3.13$N/mm^{2}$ $\times 10^{-10}$~
主塔には全てに同じ減衰を与える（Salomeの材料特性の部分のα...
ケーブル番号~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断時のケーブル張力を読み取るpoint~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
まずは初期張力の条件を変え、それ以外の条件を同じにしたも...
比較するものの初期張力を1にする。（9では初期張力20Nなので...
・２本モデル
''9と13、10と14、12と16''を比較
,,初期張力,ヤング率,サグ
,9,20N,50MPa,無
,13,50N,50MPa,無
,,,,
,10,20N,50MPa,有
,14,50N,50MPa,有
,,,,
,12,20N,100MPa,有
,16,50N,100Mpa,有
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
この結果より２本モデルにおいては、初期張力は破断後のケー...
・４本モデル~
''17と21、18と22、20と24、26と30、28と32、51と55、59と63...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
1〜4のケーブルで初期張力がそろっているものは(17など）初期...
しかし1〜4のケーブルで初期張力がそろっていないものは（34...
**9月19日 6本モデル [#r893c93b]
6本モデル ケーブルのヤング率はすべて50N/mm2に設定~
張力 上:28.70N 中:20.09N 下:11.48N~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断１~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断２~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
4本モデルの破断１と同じ挙動~
破断３~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
４本モデルの破断２と同じ挙動~
**9月26日 張力の入れ方は同じでヤング率を変えてみる~ [#i59...
張力 上:28.71N 下:8.61Nは同じまま、ヤング率を上ケーブル80...
張力の変化量~
,ue_hidari,ue_migi,sita_hidari,sita_migi,(上の張力÷下の張...
,4.608NN,-4.608N,-,-2.201,2.094
ケーブルのヤング率が変わることで破断時のケーブル張力の変...
**9月18日 ヤング率を同じにする [#we683b9d]
*破断1 [#zbaea865]
ヤング率を全てのケーブルで50N/mm2と同じにしたときの塔の動...
張力 上:28.71N 下:8.61N~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
張力 上:14.35N 下:20.10N~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ヤング率を同じにしても動きに違いはなかった。~
張力の変化量~
,,ue_hidari,ue_migi,sita_hidari,sita_migi,(上の張力÷下の...
,張力 上:28.71N 下:8.61Nのとき,3.71N,-3.70N,-,-3.50N,1.06
,張力 上:14.35N 下:20.10Nのとき,8.65N,-8.65N,-,-8.15N,1.06
,張力 上:14.3526N 下:14.3526Nのとき,5.826N,6.18N,-,6.17N...
ケーブルのヤング率を揃えることで張力の変化量比に影響が出...
ヤング率を100N/mm2に揃えたときの張力の変化量比は、1.03と...
ヤング率を揃えると張力の変化量比は1.0に近づくと考えられる...
*破断2 [#ve53e7c8]
ヤング率を50N/mm2に揃えた時、張力の変化量比は1.06になった...
ヤング率を100N/mm2に揃えた時、張力の変化量比は1.11になっ...
このことから、破断2の場合はヤング率を揃えると1.11に近づく...
**9月17日 張力の大きさを上下逆にする[#d441503d]
9月11日の実験モデルでケーブル張力の大小を上下逆にしてみた~
張力 上:14.4081N 下:27.0556N~
破断１~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
張力の大小を逆にしても違いはなかった。~
破断２~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
張力の大小を逆にしても動きに違いはなかった。~
破断１と２の破断後のケーブルの張力~
破断１~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断２~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
張力の変化量~
,,ue_hidari,ue_migi,sita_hidari,sita_migi,(上の張力÷下の...
,破断１,14.7217N,-14.7136N,-,-6.6275N,2.218
,破断２,-,-12.8582N,5.4554N,-5.4553N,2.357
張力の大きさを上下で逆にしても、上下での張力の変化量の比...
**9月11日の実験 [#d67e2b24]
実験モデル~
張力 ケーブル上:24.0161N ケーブル下:15.7814N
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
今までは先端変位を測定してDAFを求めてきた。しかし本数を増...
今後は破断時の塔の動きを追うことで、ケーブル本数を増やし...
モデル破断１の塔の動きは次のようになった。(ケーブル要素は...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブルあり（scale1とscale30)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
初期位置に戻ってきたときの塔の形~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
モデル破断2の塔の動きは次のようになった。(ケーブル要素は...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブルあり(scale1とscale10)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
初期位置に戻ってきたときの塔の形~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断１と２の破断後のケーブルの張力~
破断１~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断２~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断時の張力の変化量~
,,ue_hidari,ue_migi,sita_hidari,sita_migi,(上の張力÷下の...
,破断１,8.5869N,-8.5872N,-,-3.8658N,2.221
,破断２,-,-21.43226N,9.0947N,-9.09228N,2.358
ケーブル１本モデル（scale foctor 100）~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブル２本モデル(scale factor 400)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブル１本、２本、４本のときで動きが異なることから、本...
**8月1日の実験と解析の結果 [#k9f0f58b]
*ケーブル１本の時 [#d61893f5]
解析は健全時の変位がほぼ一致するように温度荷重を与えた。~
実験と解析の結果~
,,静的,動的,健全時,DAF
,実験,0,30.01,57.99,1.932
,解析,0,30.00,59.20,1.973
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
実験と解析のグラフ~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブル破断時の周期が解析の方が速い。~
*ケーブル2本の時 [#id80de6e]
解析は静的変位がほぼ一致するように温度荷重を与えた。~
初めは温度荷重によりケーブル張力を求め、片方のケーブルを...
次に、千代岡さんから教えていただいた２段階解析をした。~
2段階解析では温度荷重はできない？エラーがでてしまった。
そのため、ケーブル２本とも「あたかもあったかも」の方法を...
この方法で減衰はかかった。（この方法を①とする。）~
実験と解析の結果~
,,静的,動的,健全時,DAF
,実験,18.34,34.54,0,1.884
,解析,18.40,35.43,0,1.926
実験と解析のグラフ~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
ケーブル破断時の周期が実験の方が速い。~
DYNA_VIBRAを使用し、片方を温度荷重、もう片方を同じ張力で...
つまり、温度荷重と点載荷の両方を使用するときはDYNA_VIBRA...
DYNA_NON_LINEを使用して片方温度荷重、もう片方を点載荷にし...
DYNA_VIBRAを使用しているときはcalc_modeとcalc_amor_modal...
そこで材料設定の部分で減衰に関わるα、βを設定し減衰をかけ...
（この理由はケーブルのみにα、βを設定していたから）~
このときの実験と解析のグラフは次のようになった。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
このグラフでは周期のズレは小さくなっている。このことから...
方法①と②のモデル~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
方法②を解析の途中でtime stepを変える解析してみた。(0s〜19...
またα、βは梁のみに設定して解析を行った。~
このときの実験と解析のグラフは次のようになった。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
,,静的,動的,健全時,DAF
,解析,18.40,35.79,0,1.945
**片持ち梁の動的解析について [#e65d6f36]
[[FrontlSTRチュートリアルP85〜:https://manual.frontistr.c...
これによると片持ち梁のDAFの理論値は1.9753となる。~
断面10mm*50mm,長さ500mmの片持ち梁モデルで動的解析を行って...
まずはステップ荷重を梁の先端に載荷する解析を行った。しか...
Exception utilisateur levee mais pas interceptee. ...
! Les bases sont fermees. ...
! Type de l'exception : NonConvergenceError ...
! ...
! Arrêt pour cause d'absence de convergence avec le...
! dans l'algorithme non-linéaire de Newton. ...
! La base globale est sauvegardée. Elle contient le...
! l'arrêt. ...
! ...
! Conseils : ...
! - Augmentez ITER_GLOB_MAXI. ...
! - Réactualisez plus souvent la matrice tangente. ...
! - Raffinez votre discrétisation temporelle. ...
! - Essayez d'activer la gestion des événements (dé...
! exemple) dans la commande DEFI_LIST_INST.
次に「あたかもあったかも」を用いてtimestep0.001で解析を行...
理論値と考えていた1.9753にはならなかった。~
**塔の長さ、境界条件とDAFの関係 [#w1cc0437]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
,長さ(mm),片持のDAF,単純のDAF
,100,1.258,1.0962
,200,1.7225,1.36
,300,1.9082,1.6528
,400,1.9609,1.8328
,500,1.9791,1.9148
,600,1.9878,1.9527
,700,1.9884,1.9701
,800,1.9915,1.9789
,900,1.9944,1.9858
,1000,1.9947,1.9889
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
・境界条件：片持ばり~
塔の長さを変えながら先端に水平方向の荷重を与えた。~
塔の長さを長くすることで先端が境界条件を考えない自由な動...
・境界条件：単純梁~
DRYのみを自由にした単純梁モデルを作り、梁の長さを変えなが...
片持ち梁と同様に長くすることで中央が境界条件を考えない自...
・両方を比べてみて
片持ち梁の方が同じ長さ時のDAFが大きく、また長さの伸びに対...
しかし、固定部からの距離で見る（片持ちL=100mmと単純L=200m...
この理由を考えながら、境界条件がどのようにDAFに影響するか...
**ケーブルの角度とDAFの関係 [#c71d2581]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
塔の長さ、荷重は同じにして、ケーブルの角度のみを変化させD...
,,case1,case2,case3,case4
,DAF,1.88059,1.88059,1.88059,1.88059
この結果からケーブルの角度はDAFに影響しないことが分かった...
**7月2日の実験 [#me6f7b2c]
・２本の場合(減衰なし)~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
破断していないケーブルがあるため、初期位置より破断したケ...
**７月２日の実験 [#vef2bead]
7月2日の実験ではタワーにケーブルを１本取り付けたものと２...
・１本の場合~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
１本のときは上のような状態からケーブルを切った。~
実験と解析の先端変位のグラフは次のようになり、近い値とな...
しかし、解析の方が周期が早く見える。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
time stepを0.01から0.005にすることで周期の問題は解決した...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
DAFは1.897となり、理論値に近づいた。~
time stepを0.005から0.001にすることでDAFが1.919となりさら...
この原因はまだわかっていない。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
実験と解析結果~
,,動的変位,静的変位,DAF,DAF理論値,
,実験,31.76mm,17.70mm,1.794,1.963,
,解析step0.01,32.84mm,17.66mm,1.859,1.963,
,解析steo0.005,33.50mm,17.66mm,1.897,1.963,
,解析step0.001,33.90mm,17.66mm,1.919,1.963,
**7月10日 [#z401a45c]
７月２日に行った実験の再現を解析で行っている。~
ケーブル部分にプレストレスを与えるため、温度荷重のかけか...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
温度荷重によってプレストレスが入り、塔がたわむことを確認...
今後はたわんだ状態からケーブルを破断させたときの挙動を確...
**減衰について [#f0235647]
これまでの減衰定数の計算の考え方は間違っていた。~
https://wwwra.meijo-u.ac.jp/labs/ra007/murata/pdf/textboo...
このページを参考にした。~
$δ=\displaystyle \frac{Y_1}{Y_2}$で対数減衰率δを求めるが...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
この減衰定数を計算すると0.00477と与えたかった0.005に近い...
・減衰定数を0.05に設定し、解析を行った。~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
このグラフからDAFを求めると1.850で理論値のDAF=1.855に近い...
・減衰定数0.1に設定した場合~
理論値が1.730なのに対し、解析から求めたDAFは1.725となった...
**7月3日 [#s396f2bf]
構造関係メモから$DAF=1-\cos\omega t$と表すことが出来るこ...
減衰についてだが、現在自分の与えたい減衰定数を反映するが...
減衰定数0.005で与えたいところを、変位のグラフから計算する...
**3本のばねモデルでの実験結果 [#jf53c381]
,,実験値,解析値
,3本静的,32mm,32.48mm
,破断時(動的),68mm,65.15mm
,2本静的,47mm,48.72mm
--動的増幅率 $DAF={\displaystyle \frac{S_{dyn}-S_{0}}{S_{...
---$S_{dyn}$: 動的解析の最大応答値
---$S_{stat}$: 静的解析の応答値
---$S_{0}$: 初期値
健全時を原点にとり、動的変位と静的変位の比をとったものと...
健全時を原点にとったものを$DAF_{健}$、自然長を原点にとっ...
構造関係メモから自然長を原点にとったDAFは$DAF={1-\display...
tは破断開始から最大変位にいたるまでの時間を代入する。~
解析、実験の結果は次のようになった。~
,,,$DAF_{健}$,$DAF_{自}$,$DAF_{理}$,$DAF_{自理}$,
,,解析,2.01,1.337,2.00,1.333,
,,実験,2.40,1.447,2.00,,
**6月26日 [#td51e133]
固有振動数を解析で求め、レイリー減衰で用いるα、βを計算し...
減衰定数0.005~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
減衰定数0.005よりもかかっている気がする、、
解析方法を見直す必要あり
**6月24日 [#zc7a53dc]
実験した3本モデルを減衰を考慮して解析してみた。~
グラフを見ると減衰しすぎな気がする。なにか設定をミスって...
減衰0.5%のとき~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
減衰1%のとき~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
減衰5%のとき~
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**6月19日 [#k44de0d4]
断面、縦22mm、横70mmで高いところから紐+バネで吊るした解析...
これを踏まえ断面、縦70mm、横22mmの縦長にし、低いところか...
今回の実験の結果からばね定数等を求め解析を行い、実験と解...
,,実験値,解析値
,3本静的,32mm,32.48mm
,破断時,68mm,65.15mm
,2本静的,47mm,48.72mm
動的変位と2本の静的変位の比から求めるDAFは次のようになっ...
,,実験値,解析値
,DAF,1.447,1.337
今回のバネの本数が破断により3本から2本になったため、DAFの...
DAFを衝撃係数問題から求めると理論値は２になり、実際に解析...
**6月11日 [#e33f392d]
5月31日に実験をIHIや他大学の教授に見せ意見をいただいた。~
初めの実験モデルでは複数のパラメータにより、調べたい対称...
そこで簡単なモデルから徐々に進めていくことにした。~
初めは3本のケーブルで鋼材を吊るしたモデルで、真ん中のケー...
千代岡さんに"あたかもあったかも"を教えていただき、salome...
3本のケーブルで真ん中の1本を破断させたときの時間と変位の...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
この結果に加え、ケーブルが2本のときの静的変位を調べること...
（DAFの導き方は千代岡さんの卒論を参考にした。）~
次にケーブルの剛性を変化させることでDAFの違いを求めようと...
そこでケーブルの本数を増やし、ケーブルの破断する数を変化...
ケーブルの本数を9本にして真ん中の1本を破断させたときの時...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
DAFは1.177となった。~
今後は3本破断、5本破断、7本破断の解析を進めていく。~
**5月29日 [#uab9dea7]
衝撃を与える時間を変えてみた
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
マルチファイバーに衝撃を与えた時の先端変位
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
**シェル要素の片持ちばり [#mcfd08a0]
[[http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin...
[[http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?plugin...
このページを参考にして解析した
AFFE_CARA_ELEMのわからない用語を調べる必要あり
EXCENTREMENT_FO： 法線方向で、格子と平均曲面の間の距離を...
[[AFFE_CARE_ELMEの説明:https://code-aster.org/doc/v14/fr/...
**ビーム要素の片持ちばり [#d39d49cd]
点を固定できるのか:
lineをグループ化で固定、載荷の点を作ったがAFFE_CHAR_MECA...
GROUP_MAを選択していたため点が出てこなかった、GROUP_NOを...
固定面をどうやってつくるのか:
長方形フェースを固定したい面に付ける方法を考えたがうまく...
**Midas [#ge1b6317]
https://in.midasit.co.jp/wp-content/uploads/2020/07/const...
https://in.midasit.co.jp/wp-content/uploads/2020/07/const...
*参考資料 [#u9b8ef7d]
[[code-asterの説明:https://code-aster.org/doc/v15/en/inde...
[[code-aster質問ページ:https://forum.code-aster.org/publi...
[[構造要素の概要:https://www.opencae.or.jp/wp-content/upl...
[[減衰について:https://www.cybernet.co.jp/ansys/learning/...
[[面内と面外の説明:https://bonperson-civil.com/in-out-pla...
[[cable軸力の計算?:https://code-aster.org/doc/v15/en/man_...
*12月8日 [#m89d8b57]
サンドイッチ梁
,メッシュの長さ,要素数,変位,相対誤差,計算者
,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,佐々木
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*11月29日 [#u96a8789]
異方性一次 [#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,佐々木
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
等方性二次 [#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,佐々木
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*11月24日 [#h4d56322]
,メッシュ長さ,要素数,変位,相対誤差,計算者
,0.7,1455234,0.422484,0.01388,安藤
,0.8,142973,0.422570,0.01409,安藤
,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,佐々木
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*11月17日 [#gb489580]
,メッシュ長さ,要素数,変位,相対誤差,計算者
,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,0.458,服部
,5,494,3.8580375,0.422,服部
,6,581,2.50682,62.416,梶原
,7,133,1.41225,78.827,梶原
,8,78,1.2887175,80.68,工藤
,9,72,1.2879925,80.69,工藤
,10,60,1.14344,82.85,佐々木
,11,65,1.23124,81.154,佐々木
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*11月10日 [#z1f8fbe7]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/shibata2...
*10月27日 [#dcab77d9]
コマンドを学んだ。
pwd:今いる場所
ls:今いる場所のフォルダ名ファイル名等を見る
mkdir:フォルダの作成
cd:ファイル間の移動（ディレクトリーの移動)
gedeit ファイル名 .tex & :テキストファイルを開く
vi:ファイルを開く
Esc : q:開いたファイルを閉じる
cp 元のファイル名 新しいファイル名:元のファイルをコピー
cat:ファイルの中を確認
rm:ファイルの削除、cd .. :1個前に戻る
rmdir:フォルダを削除
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
