松田の卒論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
とりあえず理論通りの力で変形挙動を確認したい
*卒業論文 [#p2f86cec]
テーマ: 多々羅大橋の温度解析
やりたい事：ケーブルの固有振動数による張力推定→温度による...
日々の温度変化でケーブルがどう変化していくか。その後リス...
最初は↑この方向性だったが、だんだん変わってきた。最終的に...
実測平均応力変化と数値解析の平均応力変化を比較する。~
『温度変化がケーブルにどのような影響を与えるのか知る』
**冬休み前までにやりたい事 [#y8837db4]
☑実測ひずみデータがどのファイルなのか確認→すでにpython使...
・実測ひずみデータ整理して、時間-ひずみのグラフ作る~
・①サロメでサグ有モデル作る。（すでに作ったけどなんか違う...
・②軸力入れて、重力入れて解析。~
・31点の座標が、①と②で同じになるはず。近い値になればおっ...
・そのモデルに、実際の温度データを入れて、解析。~
＊最初の状態は、基準20℃ということにする。30℃にしたい時は...
・解析から、ひずみデータ得る。~
・実測ひずみデータ、解析ひずみデータ、を比較する。整合性...
・一応現状を英語でも説明できるように~
**夏休み前までやりたい事 [#a202b84f]
・DATファイルが何を示しているのか~
☑DATデータでグラフ作る~
・エクセルデータまとめて表作る~
✖・測定地点はどこ？24本あるはず。全体図にわかりやすいよう...
☑・ケーブルを一本取り出して、モデル作る。データからみてN6...
☑・★★★ケーブルに温度いれて解析する。とりあえずN6だけで。S...
✖・それが終わってから、全体のモデルつくる？そもそも今の私...
✖・全体モデルに温度入れる。
☑・温度変化によって張力がどう変化するかを知りたいの？で、...
**進捗 [#z78c72b3]
N6のケーブルの座標を資料から探し出す。解析基準点 or 図面...
けた側の解析基準点が資料に載っていない（多分）ので、図面...
図面上基準座標~
塔側(m)~
X=-447.336、
Y=-4.206、
Z=201.565~
桁側(m)~
X=-673.000、
Y=-11.839、
Z=43.499~
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/matsuda202...
N6の解析用ケーブルの長さ＝275.72m＝275720mm
ソリッド要素で解析する。被覆、仮想材料も入れてモデルを作...
青木先生と話して、ケーブル要素で解析することにした。
被覆ありでやるなら、ソリッド要素だけど、、どうするか。
結局、被覆ありでやりたいから、ソリッド要素でやることにな...
メッシュしてから結合。円柱の所からモデル作る。押し出して...
CSVとDATのデータは平均ではなくて、１分ごとや１０分ごとの...
6/17
この間ソリッド要素でやるということになったが、ソリッド要...
6/21
エクセルデータのまとめ方。~
A列をyyyy/mm/dd hh:mm:ssの表し方にする。ctrl＋１、ユーザ...
AA列には一分ごとの時刻を書く。秒のところは００にする。例...
AB２にはB列を１分ごとにした値を書く。例；AA２の式は、=AVE...
絶対に.xlsxで保存する！！！！！！！！！！！！！！！！~
62個のエクセルデータすべてやる。
6/23
ついにエクセルデータをすべて１分ごと平均して一つのファイ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/matsuda2025/62k...
6/24
「多々羅大橋計測結果の概要と整理・分析方針(案)」より、
調査では、N6だけ温湿度を計測。N6は削孔位置が多く、ケーブ...
私の研究では温度についてやるので、N6のケーブルを取り出し...
＊Nケーブルは北側（上り・走行車線）、Sケーブルは南側（下...
＊N６は生口島側から６番目のケーブル~
DATファイルのF017, F110 がN6に関係している。 ~
・DATファイルがなんのデータなのかもう一度確認。多分温度。~
・DATを、csvにして(Python 使って)、サロメに入れる。~
・DATのデータでグラフ作りたい。前にも重くて落ちたから、ど...
・「石川島播磨技報」のpdf全てよめ。
6/25
先生と話した。DATファイルにも温度のデータないかも。N1,N6...
来週のゼミまでにやりたい↓
・もらったすべてのデータから温度のデータ探す。なかったら...
・DATファイルを自分で開いてグラフにできるように。Chat Gpt...
・DATのデータをCSVにする。Chat GPT使うとできるはず。~
・「多々羅大橋計測結果の概要と整理・分析方針(案)」をもう...
・「石川島播磨技報」を読む。~
・先週の英語ゼミで選んだ英語の論文が私がやろうとしている...
・サロメでN6のケーブルモデル作って、適当な温度入れて解析...
7/1
DATデータについて、
・4チャンネルある。~
・単位はマイクロストレイン~
・おそらくひずみのデータ。ってことは、応力が分かる。σ=E×ε~
7/3
DATにそもそもなんのデータが入っているのかまだわかんない。
7/6
10mのケーブルモデル作って、サグと温度適当に入れて解析した。
(0,0,0)と(10000,0,0)でケーブル。メッシュ100。~
(0,1000,0)と(0,-1000,0)で固定線。メッシュ15。~
固定部に壁作った。これがマルチファイバー要素。たぶん。~
重力、温度25℃(プレストレス−25)にして解析。~
サグなった。~
とりあえず温度入れて適当な値で解析できた。~
ここで、疑問。~
重力→解析する→サグ出てくる→その変位？の値使ってモデルつく...
(このようなやり方でやれば良いの？)~
このやり方でやろうとしたら、重力入れてプレストレスzeroに...
温度or張力を必ず入れないといけない？だとしたら張力どうや...
柴田さんの↓参考にしてやった。~
[1]まずサグのないまっすぐなcable要素を用意する
cable要素のみの場合、解析がエラーになってしまうため、固定...
[2]cableに張力、重力を与えたときの100点の座標を読み取る。~
[3]読み取った座標をプロットし(各点をt0〜t100とする)100点...
7/9 ↑これあんまり参考にしない方がいいらしい。
サグのやり方
１本のケーブル取り出す。(何要素？)
サグを手計算で求めて、、なんだっけ？
重力かける。
点に全て軸力を入れる。
モデル完成。
温度入れる。
7/11
以下柴田さんのやり方を教えてもらって真似てみる。
サグの入れ方。
サグモデルの座標、軸力の計算方法
高さ：h、径間長：ｌ、l/2におけるサグ：f、単位長さ重量：w ...
L：ケーブル形状長、x1 y1：下端、xm ym：中央、x2 y2：上端
$ \sqrt{L^2 - h^2}= 2C \sinh\left( \frac{l}{2C} \right), ...
これらの式から$ f - \frac{h}{2} + C \cosh\left( \frac{x_m...
これを解くとC(カテナリー数)が求まる。
ここで$ x_m = C\, \mathrm{arctanh}\left(\frac{h}{L}\right...
Cを求めるここができたら下端、上端のx,y座標を求める。x座標...
軸力は $ T = w y $で各要素ごとを求める。
１．モデル作成。点たちの座標を手計算で求めて、それらを線...
２．軸力~
３．重力~
４．温度。気温をエクセルで入れる~
５．解析~
(モデル作った後)と(軸力&重力入れた後)で、点たちの座標が一...
軸力を入れて上に上がり、重力を入れて下に下がるから、一致...
まず、何メートルの長さで解析するのか決める。N6の長さで最...
N6の正確な長さ→275622mm~
点たちの座標を手計算で求める。エックセルでやる。エクセル...
ex)１点目、
塔側(m)~
(x1,y1,z1)=(-447.336,-4.206、201.565)~
桁側(m)~
(x2,y2,z2)=(-673.000,-11.839,43.499)~
h=158066mm~
l=225793mm~
L=275622mm~
石川島播磨技報に色々と有益な情報がのっていたのでもっとし...
7/16
柴田さんのモデルはだいたい１mくらいで、３０個の点。~
私がやりたいのは２７５mである。何個の点を取れば良いか？
今は、サロメでモデル作るよりも、DATファイルと温度（エクセ...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/matsuda202...
7/17
今日はあとで自分の力でPythonを使ってDATファイルをなんとか...
7/18
DATファイルについて。~
計測器番号2_N6,N1.zipには、 F001.DAT～F159.DAT と F001...
計測器番号3_S41_S42.zipには、F001.DAT～F237.DAT と F001...
・chatGPTで.DATを.csvにした。F001_multichannel_time_align...
・python開く。（リボンみたいな△のアイコン）上にある。~
・「kouzou,デスクトップ,python,dat」の中にcsvgraph_hukusu...
・「kouzou,デスクトップ,python,dat」の中にF001_multichann...
・エクスプローラーのＤＡＴの風車みたいな水色のマークにcsv...
・pythonで以下のコードを打った。~
import pandas as pd
import matplotlib.pyplot as plt
import glob
import os
# 対象CSVファイルをすべて取得
file_list = glob.glob("F*_multichannel_time_aligned.csv")
# CH名の一覧を取得（最初のファイルから）
sample_df = pd.read_csv(file_list[0])
ch_columns = [col for col in sample_df.columns if col.sta...
# CHごとにプロット（1つの画像に全体表示）
for ch in ch_columns:
plt.figure(figsize=(15, 6))
for file_path in file_list:
df = pd.read_csv(file_path)
df['Time'] = pd.to_datetime(df['Time'])
label = os.path.basename(file_path).split("_")[0]...
plt.plot(df['Time'], df[ch], label=label)
plt.xlabel("Time")
plt.ylabel("Value")
plt.title(f"Comparison of {ch} across all files")
plt.legend()
plt.grid(True)
plt.tight_layout()
# 画像保存
output_filename = f"combined_{ch}.png"
plt.savefig(output_filename)
plt.close()
print(f"{output_filename} を保存しました。")
・ターミナルで ls →F001_multichannel_time_aligned.csv F00...
・ch1だけグラフを作れた。で、PCが固まった。~
・電源切った~
・labelがいらないので、そこのコード変えた。少しは軽くなっ...
for file_path in file_list:
df = pd.read_csv(file_path)
df['Time'] = pd.to_datetime(df['Time'])
#label = os.path.basename(file_path).split("_")[0...
plt.plot(df['Time'], df[ch])#, label=label)
・ch1とch2のグラフ作れた。で、またPCが固まった。~
7/19~
・今度は、ch3とch4のグラフだけ作るように指示する。以下の...
import pandas as pd
import matplotlib.pyplot as plt
import glob
import os
# 対象CSVファイルをすべて取得
file_list = glob.glob("F*_multichannel_time_aligned.csv")
# 対象チャンネルのみ（CH_3とCH_4）
target_channels = ['CH_3', 'CH_4']
# 各チャンネルごとにプロット
for ch in target_channels:
plt.figure(figsize=(15, 6))
for file_path in file_list:
df = pd.read_csv(file_path)
df['Time'] = pd.to_datetime(df['Time'])
# サンプリングなし・labelなしで描画
plt.plot(df['Time'], df[ch])
plt.xlabel("Time")
plt.ylabel("Value")
plt.title(f"Comparison of {ch} across all files")
plt.grid(True)
plt.tight_layout()
# 画像として保存
output_filename = f"combined_{ch}.png"
plt.savefig(output_filename)
plt.close()
print(f"{output_filename} を保存しました。")
・ch3とch4だけつくって、というコードなのに、ch1から作られ...
・次は、2つ目のフォルダにある237個のDATも同じようにやる。...
•csvへの変換は終わった。~
7/20~
重いと思うから、点だけにしても固まる。
色々試したが、すべてだめだった。
研究目的：加速度データから固有振動数を抽出し、ケーブル張...
加速度データ → 固有振動数 → 張力推定 → 温度変化との整合性...
１加速度データをFFT（高速フーリエ変換）すれば、ピーク周波...
２固有振動数がわかれば、張力Tを逆算できる。~
３これを温度変化と照らし合わせて、「温度が上がると張力が...
＊FFTとは、時間ごとに変化しているデータを、**「どんな周波...
気温が上がる → ケーブルが伸びる → 張力が緩む（応力が下が...
「整合とれるか確認」とは？~
加速度データから求めた張力の変化と、温度の変化に基づいて...
つまり、、~
気温の変化に対応して、ちゃんと張力（＝固有振動数）がそれ...
を確認すること。~
7/23~
pythonでch3だけといっても、何をしてもch1がつくられる現象...
！！実行する前にコードを保存する！！
フォルダ１のF001～F159とフォルダ２のF001～F237から、それ...
7/24~
①サロメで温度れて平均応力変化＝温度変化による理論的な平均...
②温度変化による理論的な平均応力変化⇔加速度データより推定...
「温度変化による理論的な平均応力変化」いらないなら、いき...
サロメで温度れて平均応力変化⇔加速度データより推定張力（整...
7/24~
先生と話した。こんな感じで進めてく。↓
今やる事~
１．カテナリー曲線でケーブル形状を定義~
→サグを含んだ状態でのケーブル座標を手計算で求める。~
２．FEMモデルをサロメで作成~
・上記の点群を使ってケーブルジオメトリを作成~
・ケーブルに重力を作用させる~
・軸方向に「事前張力（軸力）」を入れる（手計算で求め...
３．温度荷重を与えて静的解析~
・ケース１：＋30℃~
・ケース２：－10℃~
・結果として「一番サグがある部分のｙ軸方向変位」を取...
・測定地点だけの場合と、全ての要素点に同じ温度を与え...
４．実際の温度データを使った解析~
・時系列温度データをモデルに与える。測定地点だけ~
・すべての要素点に適用して解析もやる~
・ケーブルサグの変位と応力の時間的変化を観察~
・※これは動的ではなく静的解析？多分。（時間は進めない...
５．ひずみデータとの比較~
・実測されたひずみデータとFEMの出力を比較し、~
・FEMで得られた張力または応力分布との整合性をみる。~
まとめ~
・「温度変化によってケーブルのサグが変わる」~
・それは想像よりも大きい変位かも。~
・外力（風や荷重）を受けやすい状態になる。~
・ってことは、結果として、橋全体の挙動や安全性に影響を与...
→よって、「日々の温度変化を＊＊過小評価＊＊せず、リスク要...
サグの変化→ 温度変化が構造形状を変える要因となりうる~
応力の変化→ 材料の限界に近づくリスク、疲労寿命にも影響~
実測データとの比較→ モデルの妥当性確認＋安全性評価の信頼...
日々の温度変化→ 通常は見逃されがちだが、実は構造健全性に...
7/30~
h=158.066m
l=225.793m
L=275.622m
C=3800m
x1=24689.9m
xm=24802.8m
x2=24915.7m
y1=46307.11m
ym=46385.94m
y2=46465.18m
f=0.204714m
w=929.04N/m
y2-y1=158.066
ｘ(mm) y(mm)
0 0
7526.433333 5242.494193
15052.86667 10486.80519
22579.3 15732.9332
30105.73333 20980.87842
37632.16667 26230.64106
45158.6 31482.22134
52685.03333 36735.61944
60211.46667 41990.83558
67737.9 47247.86997
75264.33333 52506.72282
82790.76667 57767.39432
90317.2 63029.88468
97843.63333 68294.19412
105370.0667 73560.32283
112896.5 78828.27103
120422.9333 84098.03891
127949.3667 89369.6267
135475.8 94643.03459
143002.2333 99918.26279
150528.6667 105195.3115
158055.1 110474.181
165581.5333 115754.8713
173107.9667 121037.3828
180634.4 126321.7157
188160.8333 131607.8701
195687.2667 136895.8463
203213.7 142185.6444
210740.1333 147477.2648
218266.5667 152770.7075
225793 158065.9728
節点力（軸力）
T1 43021157.47
T2 43026027.96
3 43030900.13
4 43035773.99
5 43040649.55
6 43045526.79
7 43050405.71
8 43055286.33
9 43060168.64
10 43065052.63
11 43069938.32
12 43074825.69
13 43079714.75
14 43084605.51
15 43089497.95
16 43094392.09
17 43099287.91
18 43104185.43
19 43109084.64
20 43113985.53
21 43118888.12
22 43123792.4
23 43128698.38
24 43133606.04
25 43138515.4
26 43143426.45
27 43148339.19
28 43153253.62
29 43158169.75
30 43163087.57
31 43168007.08
8/2~
サロメでN6モデル作った。重力と軸力入れると赤になる。
8/4~
0℃で、計算で求めたサグ状態にする。~
94.8(ケーブルkgw/m)×275.622(ケーブル長さm)×9.8(重力加速度...
256063.8629÷275622＝0.92904(N/mm)~
これで求められた0.92904をモデル全体にかける。~
重力はかけない。軸力はかける。？
モデル全体に－10℃をかけて解析する。~
モデル全体に30℃をかけて解析する。~
両方まわった!!!
とりあえずここまででプレゼンまとめる
メモ：.hdfで開く。.commで開く。cableN60804.comm
モデルはN6のケーブル。ケーブル要素でやった。両端はピン。...
10/14
中間発表後の頭の中整理。~
中間発表→加速度データによるケーブルの「張力」と、実測温度...
しかし、方法変更。~
今後→実測データの「ひずみ」と、実測温度データを入れてFEM...
★元々の方法~
・固有振動数f1をFFT解析→弦の理論式T=(2Lf1)^2μから張力推定~
・FEM解析で実測温度データいれて、ケーブル要素の平均応力を...
・両者の「変化量」や「比例関係」を比較して整合性みる~
問題点~
・弦の理論式ではサグが考慮されていない（直線弦仮定）。~
・実際のケーブルではサグによって張力分布や剛性が変わる。~
・したがって、「張力 vs 応力」の比較は前提条件が異なるた...
★修正版~
・実測ひずみデータ整理~
・FEM解析で実測温度データいれて、ケーブル要素のひずみを算...
・両者の「ひずみ（実測） vs 温度によるひずみ（解析）」を...
要するに、、、~
「実測および解析における温度—ひずみ関係の比較」をする！！！
「張力と応力の整合性を見る」から「ひずみの整合性を見る」...
「比較する物理量」が変わるが、「目的（温度応答を再現する...
むしろサグを考慮できる分、より正確で現実的な方法になる。
温度変化がケーブルの構造挙動に与える影響を、より現実的に...
実橋の“温度によるひずみ応答”を、サグを考慮したFEMモデルで...
補足：「張力」は結局「応力 × 断面積」なので、応力もひずみ...
→ より一次的で確実な指標にシフトしただけ。
ところで、ひずみデータはどのファイルにあるのか？~
csvではなかった。しかし、表３のデータ諸元ではcsvと書いて...
グラフと値的にDATっぽい。~
しかし、「多々羅大橋計測結果の概要と整理・分析方針（案）...
しかも、DATファイルから作った１〜４chのグラフと「多々羅大...
10/15
DATファイルから作った１〜４chのグラフはひずみデータだった。
10/19
今、サグ有モデルを作って、重力と軸力入れて解析したい。温...
カテナリーのサグ有ケーブルは20℃の時の形であるー１０℃とか3...
そこで、サロメの重力と軸力だけ入れて解析する時は、温度は2...
→温度＝20℃に“する”というより、“ΔT=0にする”（温度荷重を無...
10/20~
①「サグ有モデルを作って、重力と軸力だけ入れて解析して（温...
温度を入れない場合（Function and Lists→DEFI_FONCTION→VALE...
0.001℃とかにしても回らなかった。今度1℃でやってみよう。
それか、重力を等分布荷重として入れてるから、ピンの反力がg...
★そもそも「サグ有モデルを作って、重力と軸力だけ入れて解析...
・モデルが正しく初期平衡状態を再現できているか（＝初期張...
・温度などの外力を考える前に、ケーブルが“自然に釣り合った...
・温度を入れると熱ひずみも影響するため、初期形状確認が難...
②ケーブル要素でひずみを出すという件について~
サロメは20℃基準になってるから、20℃との差を考えて入力すれ...
サロメで、ケーブル要素だとひずみが出力できない。Post Proc...
σ=F/A=E・ε
軸力は出せるので、この↑式を使って、軸力からひずみを求める...
もし、サロメでひずみが出せるようだったらそっちの方が良い...
★次のゼミまでにやる事~
・「サグ有モデルを作って、重力と軸力だけ入れて解析して（...
・温度入れて解析して軸力出す。~
・そこから計算してひずみ求める。~
・軸力の入れ方、重力の入れ方を、もう一度確認して、ほんと...
10/23~
①温度の変化を0.004まで小さくしたら解析回った。0.003はダメ...
しかし、反力の値がおかしい。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/co...
10/29
どうしたら重力と軸力いれて、反力がgl/2に近い値になるのか。
10/30
そもそも反力はgl/2になれば良いのか？構造力学の式使ったけ...
反力は一旦無視して、とりあえず座標が同じになれば問題なし。~
★重力の入れ方について~
94.8(kgw/m)×275.622(ケーブルの長さｍ）×9.8(m/s^2)＝256063...
256063.8629（Ｎ）÷275622(ケーブルの長さｍｍ）＝0.92904(N/...
FORCE_NODALE FY＝ー0.92904~
しかし、解析は回るけどケーブルが上に上がったり変位が大き...
上の計算式がおかしいかもという事で、何か変えてみる。点が3...
FORCE_NODALE FY＝ー８２６０~
これでも変位が大きかったりなんかおかしい。
★軸力の入れ方について~
今まで、カテナリーの計算でつかうW単位重量が929.67042N/mで...
1 kgw = 9.80665 Nだから、~
94.8（kgw/m）=94.8×9.80665（N/m）=929.67042（N/m）
しかしケーブル図の資料に94.8kgw/mとあって、929.67042じゃ...
そしたら、当たり前だけどｗが一桁減ったから軸力も1桁減った...
で、解析回って、少し下に下がった形になった。~
考えなおしたら、94.8じゃなくて元々の929.67042N/mであって...
１０/３１にやりたいこと。
いろんなものが少しずつ値が間違っていたりして、座標も違っ...
軸力と重力の入れ方もまだ正解がよく分からない。とりあえず...
反力がgl/2になるのかというのも先生に聞いてみる。
10/31
★重力の入れ方について~
94.8(kgw/m)×9.80665(m/s^2)＝929.67042(N/m)~
929.67042(N/m)×275.6418685(m)=256256.0917(N) ２５万N~
256256.0917(N)÷３１(個)＝8266.325537(N)~
FORCE_NODALE FY＝ー8266.325537としたいところだが、
31点すべて同じ −8266 N はNGです。理由は「等分布荷重を節点...
FEMは要素と節点で計算するので、「連続的な荷重」を
節点に置き換えて（＝節点荷重）あげる必要があります。その...
・端の節点は要素が1本しかつながっていないのに、2本分の荷...
→ つまり「端の半分の分布荷重を2倍入れている」ことになる。~
・中間節点は左右の要素2本から荷重が来るのに、1本分しか入...
→ 「中央部分が軽く、端だけ重い」不自然な分布になる。~
結果：自重の重心が本来より両端寄りになり、支点反力のつり...
→ これが「端点は半分にしないとNG」の理由です。~
よって、自重の入れ方は、~
kotei1 FY=-4270.934862N~
ten2～ten30 FY=-8541.869723N~
kotei2 FY=-4270.934862N~
★軸力の入れ方~
結局、w=929.67042（N/m）で計算していくことにした。そうす...
FORCE_NODALE 点荷重 で与えると、「節点に外力（集中荷重）...
・方向のばらつきが大きくなる（サグ形状ごとに角度が違うた...
・外力として反力に現れる（内部初期張力ではない）~
・境界条件に依存して全体釣り合いが崩れる~
という理由で、解析結果が不安定になりやすいそうだ。~
そこで、N_INIT（初期軸力）をつかう。
「軸力を入れる」＝ N_INIT~
「外力を入れる」＝ FORCE_NODALE~
$ \quad H = \frac{wl^2}{8f} $
H=3360502.585N 300万N~
よって、N_INIT＝3360502.585~
ここで疑問。~
Q鉛直方向の張力はいれなくてよいのか？~
A入れなくてよい。鉛直成分は、自重とケーブルの傾きから、解...
反力FYの両端を足した値が総自重になればよい。256256.0917N
ケーブル全体の自重のつり合いから：両端の鉛直反力の合計=𝑤𝑆~
つまり、どんな形状でも 左右反力の合計 = 総自重。
（PESANTEURを使ってもFORCE_NODALEでも同じ）
11/1~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/N6...
11/18~
結局軸力と重力入れて座標が元のところに戻らなかったので、...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/N6...
↓これは一番上から、最初の座標、－30度、－10度、0.004度、1...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
11/20~
実データを入れて解析する。ひずみをみたい。
***11/25二段階解析スタート [#q450bdbc]
エクセルの温度データをサロメに入れてできる？？→できない。...
先生からのアドバイスで二段階解析をした。一段階目は０から2...
前まではSTAT_NON_LINE（静的）でやっていたが、DYNA_NON_LIN...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
二段階解析の一段階目のcomm
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/2d...
二段階解析の一段階目のcomm
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/2d...
凄い振動しているような感じでなんか嫌な予感がするが、この...
***11/26 [#j6ad4477]
昨日、5個の温度データで二段階解析できた。全部で800個（約6...
１～５６ のデータは解析回った。~
2024/10/23 14:46:04のデータを１として、2024/10/23 23:56:0...
１～１２８ 解析回らない。~
JUSQU_A=76201, PAS=127~
エラーコード
「Arrêt pour cause de matrice non inversible. ║ ║ La base...
「逆行列が存在しないため停止しました。 ║ ║ グローバルベー...
半日分が限界？
とりあえず、23日の午後データ（５６個の温度データ）での解...
温度変化とten16の座標と変位のグラフ~
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
さっきまで、56個の温度データがあるとき、PAS=55 としてい...
listr2 = DEFI_LIST_REEL(
DEBUT=1.0,
INTERVALLE=_F(JUSQU_A=33001.0,
PAS=600.0)
)
DEFI_LIST_REEL(
DEBUT=開始値,
INTERVALLE=_F(
JUSQU_A=終了値,
PAS=刻み幅
)
)
「1〜33001秒までを 600秒（10分）刻みで 56個」**「10分間隔...
ということで、PAS＝600にして、もう一度23日の午後のデータ...
***11/27 [#deae7b5e]
PAS=600で２３日の午後、解析できた。~
PAS=600で２３日の午後と２４日の午前も含めて、解析失敗。エ...
「<SolverError> <MECANONLINE9_4> ...
║ ...
║ Arrêt pour cause de matrice non inversible. ...
║ La base globale est sauvegardée. Elle contient les ...
║ ...
║ Conseils : ...
║ - Vérifiez vos conditions aux limites. ...
║ - Vérifiez votre modèle, la cohérence des unités...
║ - Si vous faites du contact, il ne faut pas que ...
║ contact. ...
║ ...
║ - Parfois, en parallèle, le critère de détection...
║ pessimiste ! Il reste néanmoins souvent ...
║ possible de faire passer le calcul complet en ...
║ la valeur du mot-clé NPREC) ou ...
║ en le débranchant (valeur du mot-clé NPREC=-1)...
║ processeurs. ...
」
PAS=600で、２４日の午前だけ、解析できた。~
PAS=600で、２４日の午後だけ、解析失敗。
<SolverError> <MECANONLINE9_4> ...
║ ...
║ Arrêt pour cause de matrice non inversible. ...
║ La base globale est sauvegardée. Elle contient les ...
║ ...
║ Conseils : ...
║ - Vérifiez vos conditions aux limites. ...
║ - Vérifiez votre modèle, la cohérence des unités...
║ - Si vous faites du contact, il ne faut pas que ...
║ contact. ...
║ ...
║ - Parfois, en parallèle, le critère de détection...
║ pessimiste ! Il reste néanmoins souvent ...
║ possible de faire passer le calcul complet en ...
║ la valeur du mot-clé NPREC) ou ...
║ en le débranchant (valeur du mot-clé NPREC=-1)...
║ processeurs.
２４日PMだけの温度データ~
１段階目comm
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/11...
２段階目comm
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/11...
２３PM~
一段階：０～25.39~
二段階：25.39～18.96~
２４AM~
一段階：０～18.96~
二段階：18.96～29.68~
２３PMと２４AMを別々に、それぞれ二段階解析した。２３PMの...
ちなみに、ひずみは結果の EFGE＿ELNO を使って、ヤング率E=1...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
***11/28 [#af6b30e9]
ここで、前に諦めた自重と軸力いれて元の座標に戻らない件に...
★★自重はGRAVITYでいれる。
★★温度変化（０→0.004）なしで、自重＋軸力いれると、元の座...
GRAVITY つかって、０→0.004で解析回ったcommファイル→
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/11...
ピンと張ったケーブル（直線）→ ten16 (112896.5, 79033)~
カテナリー計算後サグ有モデル→ ten16 (112896.5, 77271.98)...
GRAVITY（自重）＋ 各節点の FORCE_NODALE をかけてSTAT_NON...
→ 直線から見ればトータルで 1.7 + 2.6 ≒ 4.3 m サグ、みたい...
***12/2 [#gd2c39b4]
２４PMが解析できない件について。
「<Erreur> La matrice du système est singulière
<Action> On essaie de découper le pas de temps.
On utilise la découpe manuelle.
Découpe uniforme à partir de l'instant < 3.723850000000e+...
(soit un incrément constant de < 9.37...
Le nombre maximal <3> de niveaux de subdivision est attei...
<Action><Échec> Échec dans la tentative de découper le pa...
Temps CPU consommé dans ce pas de temps : 0.046 s
Nombre d'itérations de Newton : 4
Temps total intégration comportement : 0.00...
Temps total factorisation matrice : 0.00...
Temps construction second membre : 0.01...
Temps total résolution K.U=F : 0.00...
Temps assemblage matrice : 0.00...
Temps autres opérations : 0.01...
Mémoire (Mo) : 774.51 / 770.45 / 39.40 / 35.75 ...
このエラーがでたので、< 3.723850000000e+04> 秒の所を変え...
回った。２３PM、２４AM、２４PMを合わせて解析して、またさ...
今度は、２３AM～２９AMまですべて温度データを入れて、Value...
同じようにエラー。
Valueを60（1分毎）にした.
同じようにエラー
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
Valueを60（1分毎）にした.
同じようにエラー
***12/4 [#x9af0760]
・334819秒あたりでエラーが出てるので、今度はVALUEを６０か...
・しかし、これがいないような秒数の所でエラーも出たので、...
・DYNA＿NON＿LINE CONVERGENCE ITER_GLOB_MAXI＝５０だっ...
・いつものエラー~
・作戦を前に戻して、エラーが出るところのVALUEの数を小さく...
・１～37201秒 VALUE＝６００、 ～333601秒 VALUE＝６０、...
・解析回った！~
このときのcommファイル↓~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/12...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/12...
VALUEはPASのこと
PAS＝１０に所はケーブルの揺れ幅が大きすぎて変位とひずみが...
そのグラフがこれ↓
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
***12/9 [#w887599b]
12/4の結果がおかしかったので、PAS＝６０にして、エラーが出...
10月23日～27日のデータの解析結果↓
グラフ
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
このグラフの続き（10月27日～29日）を別の解析でやる。一回...
2段階解析を2回に分けてもいつものエラーが出てPASを６０から...
301000～320000秒くらいでいつもエラーでるのは、そこら辺で...
この問題は一旦無視していまある解析結果（10月27日～29日）...
***12/10~ [#j2c605e2]
ひずみは割と妥当な値がでたと思う。。
12/12~
実際のひずみとおもわれるデータとサロメから得たひずみのグ...
12/17~
先生と話して、いまあるデータで中間発表書いていく。
12/21~
スライドつくり。概要もつくり終わりたい。
12/22~
実測ひずみだと思っていたデータは平均応力変化だった。スラ...
解析ひずみ vs 実測ひずみ~
↓~
解析平均応力変化 vs 実測平均応力変化
***12/23中間発表 [#o4067663]
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
1/27~
卒論かきはじめなきゃやばい
***2/5~ [#s5634d8d]
ひずみを算出して、そこから平均応力変化を求め、実測データ...
・他の人が温度変化によるケーブルの変化の研究をやっていな...
・文献をもっと探してみて。似たのがなかったら、なぜないの...
・なぜ他の人が温度に注目していないのか、というのが大事~
・ピンピン以外をやらなかった理由~
・生のデータを分析した。これはいただいた貴重な生データ。...
・多々羅大橋のその日の天気と気温を調べておく~
・ケーブルをどうやって再現したのか、ちゃんと説明する~
・サグをどういうふうにいれたのか~
・サグ量だけでみたら温度によって変化することが分かる。こ...
・解析でサグが変化するのはみれた。温度とまあまあ似たグラ...
・橋全体を作るべきだったかも。これは今後やっていきたい~
・大量のデータ解析から始まり、研究した。~
・モデルの具体的な作り方、そして、温度入れたときの揺れる...
・他の人がサグを考慮しない理由。まだ文献調べていないので...
・ストーリーをうまくつくれ。目的、背景をもっと深く話す。...
・とにかく温度計の文献をさがせ。~
・もしかしたら、橋全体でやったら、そりゃ違うようになるよ...
・私みたいに、結果が合わない理由は沢山ある。~
・あまり脱線はしないように。~
今日、タイトル提出する
***2/16~ [#n8cd7642]
平均応力変化の平均とは、そもそもなに
2/17~
卒論の前夜というか当日の朝、平均応力変化の『平均』はいら...
2/19~
私の卒論の考察で言ったのは、拘束条件を考慮した解析の実施。~
では、今後、拘束条件をどうすればよいのか？~
ピンピンでダメなら、固定にすればいいわけではない。固定は...
まぁ、色々言わずに時間かかるけど塔桁他ケーブル作ればいい...
でも、全体モデル作ったら、そのケーブル全てにサグ考慮する...
ホンモノの多々羅を一回見てみたいような気もする
***2/20~ [#x7e6041b]
**一旦あたまの中整理 [#rf5bbb5d]
やりたいこと！！！
★★軸力と重力入れて、解析前と解析後の座標が同じになりたい...
N6について、ケーブル図より~
公称断面積 11400 mm^2 ~
単位重量PE被覆後 94.8 (kgw/m) =0.0948 (t/m) =929.67042 (...
死荷重完成時張力 516.1(t) =50580.11165(kN) ~
カテナリー計算によるケーブル形状長 276.315(m) ~
計算で求めた~
高さ ｈ=158.066m~
径間長 l=225.664m~
カテナリー Ｃ＝4417.523308~
形状長 Ｌ＝275.6418685m →これが「カテナリー計算によるケ...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
サロメ座標
x軸(mm) y軸(mm)
kotei1 0 0
ten2 7526.433333 5042.915455
ten3 15052.86667 10101.27387
ten4 22579.3 15175.08991
ten5 30105.73333 20264.37833
ten6 37632.16667 25369.15388
ten7 45158.6 30489.4314
ten8 52685.03333 35625.22573
ten9 60211.46667 40776.5518
ten10 67737.9 45943.42455
ten11 75264.33333 51125.85898
ten12 82790.76667 56323.87013
ten13 90317.2 61537.4731
ten14 97843.63333 66766.68302
ten15 105370.0667 72011.51507
ten16 112896.5 77271.98447
ten17 120422.9333 82548.1065
ten18 127949.3667 87839.89646
ten19 135475.8 93147.36972
ten20 143002.2333 98470.54169
ten21 150528.6667 103809.4278
ten22 158055.1 109164.0436
ten23 165581.5333 114534.4046
ten24 173107.9667 119920.5264
ten25 180634.4 125322.4246
ten26 188160.8333 130740.1149
ten27 195687.2667 136173.613
ten28 203213.7 141622.9347
ten29 210740.1333 147088.0959
ten30 218266.5667 152569.1124
kotei2 225793 158066
軸力[N]
1 4941134.759
2 4945823.008
3 4950525.614
4 4955242.591
5 4959973.952
6 4964719.711
7 4969479.881
8 4974254.477
9 4979043.513
10 4983847.002
11 4988664.957
12 4993497.395
13 4998344.327
14 5003205.769
15 5008081.734
16 5012972.237
17 5017877.292
18 5022796.912
19 5027731.113
20 5032679.909
21 5037643.313
22 5042621.341
23 5047614.007
24 5052621.325
25 5057643.31
26 5062679.976
27 5067731.339
28 5072797.412
29 5077878.211
30 5082973.749
31 5088084.043
salome軸力[N]
鉛直ｙ軸方向 水平ｘ軸方向
1 2750391.31 4104894.682 kotei1
2 8419.779663 -4.110606786 ten2
3 8427.776392 -4.102801731 ten3
4 8435.797621 -4.094996054 ten4
5 8443.843369 -4.087188698 ten5
6 8451.913661 -4.079379724 ten6
7 8460.008521 -4.07156985 ten7
8 8468.127973 -4.063759074 ten8
9 8476.272037 -4.055945962 ten9
10 8484.440741 -4.048133457 ten10
11 8492.634106 -4.040318667 ten11
12 8500.852157 -4.032504234 ten12
13 8509.094917 -4.024688822 ten13
14 8517.36241 -4.016871991 ten14
15 8525.65466 -4.009055587 ten15
16 8533.971693 -4.00123887 ten16
17 8542.313529 -3.993421056 ten17
18 8550.680196 -3.985603436 ten18
19 8559.071716 -3.977785623 ten19
20 8567.488115 -3.969968141 ten20
21 8575.929416 -3.962149828 ten21
22 8584.395644 -3.954332314 ten22
23 8592.886824 -3.946515071 ten23
24 8601.402979 -3.938697543 ten24
25 8609.944137 -3.93088149 ten25
26 8618.510318 -3.923064636 ten26
27 8627.101551 -3.915249141 ten27
28 8635.717859 -3.907434308 ten28
29 8644.359268 -3.899620345 ten29
30 8653.025801 -3.891806481 ten30
31 -2997901.667 -4104778.646 kotei2
*サロメ [#s104c0ef]
手動で選択する場合,コンテンツにチェックを入れた後,グルー...
Compound_Mesh_1だけを目つけてメッシュの所は目消す。そして...
Paravisにおいては、WarpByVectorを使うためには、オブジェク...
Geometry→31個の点の座標入力→31点を線で結ぶ~
Mesh作成（メッシュは点ではなくて線だけ）~
コンパウンドの作成→Compound_Mesh_1→メッシュの所に30個入れ...
Compound_Mesh_1右クリック→グループを作成→ノードはkotei1,t...
Asterstudyをやってく
*python [#w177c7bf]
サンプリング間隔→
legend→複数のデータを同じグラフに描画したときに、それぞれ...
label=“〇〇” 線につける名前（ただしこれだけでは表示され...
legend() 上でつけたlabelをグラフに表示する命令
*日記[#x22e75ea]
11/29~
11/28~
サロメ。先が見えない
11/27~
salome。半日エクセルで作ったデータ消えた。ので寝る
11/26~
salome
11/25~
二段階解析について知る。
11/21~
先生と進捗と今後どうしてくかはなし
11/1~
昼夜逆転というか朝と昼もわからなくなってきた
10/31~
今日はサロメはやってない。重力と軸力について勉強した。
10/30~
座屈の計算。サロメとか計算とか皆に助けてもらった。明日か...
10/29~
座屈のデザイン決めた。どうしたら重力と軸力いれてうまく解...
10/22~
サロメの続き
10/20~
サロメでつまって、皆に教えてもらった。先生に相談
10/17~
サロメ。クマが怖いから夜来るのはやめたい
10/15~
サロメ。いつになったらサロメになれるのか
10/14~
頭の中整理
10/7~
中間発表。夏休み全力で楽しんだのでこれからがんばる
10/4~
概要を書き始めた。
8/4~
先生に相談。１Nでもいいからかけてみる。重力と軸力なしで。...
8/2~
サロメでN6のモデル作った。軸力と重力を入れると解析回らな...
7/30~
サグケーブルのつくりかた教えてもらった。
7/27~英語の論文もう一回読む。「Temperature Effect on Tens...
Cable of Cable-Stayed Bridge」。 で、わからない単語とかな...
7/26~
気力あったらサロメやりたい。
7/25~
日本ファブテックさんの工事見学
7/24~
先生とはなした。
7/23~
英語で中間発表。みんなの発表が理解できるくらいになりたい...
DATのグラフやっと完成。石黒さん神。
7/22~
瀧上工業株式会社さんのお話聞いた。スライド作り。発表練習
7/21~スライドつくり。ねむい
7/20~
発表までにグラフが作れる気がしないから、一旦諦めてスライ...
7/19~
オーキャン。pythonやったけど欲しいグラフができあがらない...
7/18~
オーキャン準備。テント。今日はpythonの使い方を教えてもら...
7/16~
オーキャン準備。温度と湿度のグラフ作り。暑い。海いきたい
7/14~
オーキャンのポスター作り
7/13~
オーキャンの動画とポスター
7/12~
サグ手計算
7/11~
キャンパスクリーンデー。草取りした。サグのやりかた教えて...
7/10~
オーキャンの準備。ポスターとか。
7/6~
サロメ教えてもらった。感謝。一人でやってたら分からなすぎ...
7/3~
データとたたかった。何がわかんないかもわかんなくなってきた
7/1~
今の感じだと何の成果もないのでやばい。夏休みに研究室これ...
6/27~
サロメで適当にモデルつくろうとして、途中。誰かにサロメの...
6/24~
Drop boxを見てたら、石川島播磨技報のファイルがあることに...
6/23(月）~
やっとやっとエクセルまとめ終わったら、温度のデータなかっ...
6/22(日）~
エクセルやってる。嫌いになりそう。絶対pythonでやったほう...
6/21(土)~
久しぶりに何もない日で研究室きた。エクセルの件振り出しに...
6/17(火)~
明日発表のパワポ作った。
今週もゼミでいう事がなくてやばい。
6/16(月)~
IHIさんの貴重なお話をきいた。
6/15(日)~
エクセルデータをやっと全て平均できたのに、保存してその後...
6/10(火)~
エクセルのデータ平均してグラフ作るの終わらせないとそろそ...
エクセルデータについて、一つのファイルに12万行くらいある...
5/26(月）~
研究室には来た。
5/23(金)~
データやろうとしたけどやめた。
5/22(木)~
先生と進捗話した。
5/20(火)~
N6ケーブルの両端のの座標が分かったので、モデルをつかろう...
5/19(月)~
ケーブル一本取り出して解析するために、資料から座標読み取...
5/18(日)~
YEC本番。リスクアセスメントの発表。２個くらい質問されたけ...
5/12(月)〜5/16(金）~
毎日YECの発表練習。~
prof Daiaが研究室きた。YECの発表も見てもらった。~
5/11(日）~
スライドを修正して発表練習した。けっこうやばい。~
5/9(金）~
YECの原稿なんとなくできた。~
5/7(水)~
兼田さんに英語版スライドもらい、分からないところも聞いた...
たたらの方は、DATをPythonで開いてグラフにすることができた...
エクセルのデータは、今は保留。~
明日、N6ケーブルのモデルを作りたい。~
5/6(火)~
プレゼンができるようになるまで、たたらは一旦ストップで、Y...
英語の原稿作り始めた。~
ALARP領域についてもう少し詳しく調べる。~
5/5(月）~
サロメでN6のモデルを作りたい。しかし、設計図がうまく理解...
兼田さんの卒論スライドを見て、わからないところを無くす。~
YECで英語が理解できないということになりたくないので、英語...
5/3(土)~
「多々羅大橋計測結果の概要と整理・分析方針(案)」を」読ん...
調査では、N6だけ温湿度を計測。N6は削孔位置が多く、ケーブ...
私の研究では温度についてやるので、N6のケーブルを取り出し...
＊Nケーブルは北側（上り・走行車線）、Sケーブルは南側（下...
＊N６は生口島側から６番目のケーブル~
DATファイルのF017, F110 がN6に関係している。 ~
5/2(金)~
教えてもらいながら、適当にケーブル一本使ってケーブル要素...
python勉強中。~
エクセルデータの処理いついて、一旦保留。どうすればいいか...
兼田さんの卒論動画をみてどんな風に発表するか、考え中。（...
5/1(木)~
サグの実験見学した。pythonの本を読み始めた。
4/30(水)~
石黒さんと話して、結局pythonでデータをまとめて表にするこ...
たたらの細かい図を先生からもらい、正確な寸法がわかる。明...
DATファイルが文字化けする件について、ChatGPTにファイルを...
4/29(火)~
自分のパソコンだとpython入れられなかったので、研究室のや...
サロメを使って、たたらの一部分のモデルを作ってみようとし...
4/28(月)~
データの単位にkineとgalがあるのでgalに統一する？
横軸時間で縦軸を色々変えてグラフを作りたい。その為に、ま...
4/27(日)~
多々羅大橋ケーブル張力測定の調査資料を読んだ。~
HEDファイルはDATファイルの説明が書かれている？先生に聞く。~
自分の卒論とは別で兼田さんの卒論を英語で理解しようとして...
4/25(金)~
頂いた大量のデータに目を通した。DATファイルを開くと文字化...
*Risk Assessment [#eef9f3d5]
私の卒論とは関係ない。ただ自分で調べただけ。興味があると...
-兼田さんの[[卒論日誌:https://www.str.ce.akita-u.ac.jp/cg...
-ALARPは「As Low As Reasonably Practicable」の略称であり...
-[[H29年度厚生労働省委託事業 建設工業の設計段階における...
日本では、施工時の安全は主に施工業者の責任範囲であり、設...
-海外では事故が起こった時に通行止めや封鎖をしたくない。日...
Q User risk と Road Block risk を求める時に、色んなデー...
→A リスクの数字が小さくなるのは事実。数字が小さいほど安...
→Q そうなってくると、ARARP領域の下にあるリスクは、起こっ...
A
Q リスクアセスメントやるのは、結構大事そうな気がするけど...
A
Q どの国でリスクアセスメントが良く使われている？~
A
Q 日本では海外よりもリスクアセスメントが普及してないとい...
A
Q 日本ではいつから努力義務になった？いつから義務化された...
A
10/22↓
★リスクアセスメントとは、~
将来起こりうる危険（リスク）を事前に特定し、その発生の可...
★通常、3ステップで構成される。~
１リスクの特定 Risk Identification~
２リスクの分析 Risk Analysis~
３リスクの評価 Risk Evaluation~
★リスクの定量化の考え方~
リスク＝発生確率×影響の大きさ~
例：「年に1回起こるが、被害が小さい」 → 低リスク~
「100年に1回だが、崩壊する」 → 高リスク~
この評価をマトリクス化して、「どの範囲なら許容できるか（A...
★ALARP（アラープ）の原則~
「As Low As Reasonably Practicable（合理的に実行可能な限...
完全な安全を求めるのではなく、費用対効果のバランスを考え...
これはヨーロッパの安全工学・原子力・航空分野で広く用いら...
★リスクアセスメントが必要な理由は、~
・「感覚」ではなく「数値」で安全を判断するため~
・潜在的な危険を事前に発見して、災害や事故を未然に防ぐた...
・限られた予算や時間を有効に使うため。優先順位づけにより...
・説明責任（Accountability）を果たすため。社会的・法的に...
★リスクアセスメント＝危険を見つけて評価する~
リスクマネジメント＝その結果をもとに行動する~
★ここで、リスクアセスメントの問題点は、~
①定量評価の不確実性。データの不足。「100年に一度の津波」...
②評価者による主観の影響。リスクアセスメントは定量的に判断...
例：被害「重大」＝ 死亡1人？10人？経済損失10億円？~
過去データが少ない時、何年に1回発生すると想定するか？~
③リスク間の比較・優先順位づけの困難さ。異なる種類のリスク...
例：（例：落石による通行止め vs 地震による崩壊）~
④コストと安全のトレードオフ問題~
安全を高めるほどコストが上がる。安全対策を重視すると、施...
⑤社会的・心理的リスクの扱いが難しい。~
（例：新幹線の微小な地震揺れ、トンネル崩落事故後の不安）~
「感じるリスク（perceived risk）」と「実際のリスク（real ...
⑥結果の扱い方（運用）の難しさ。リスクアセスメントを実施し...
⑦時間的変化を考慮しにくい。多くのリスク評価は「ある時点の...
★
*メモ [#e50c8b22]
**YEC [#s047c458]
5/2 卒論動画見る。~
5/4 日本語版スライドで、分からないところまとめる。計算あ...
5/5 兼田さんにもらった資料全部読む。井上さんのやつも全部...
5/6 英語版概要、英語のまま理解する。英語で発表する原稿つ...
5/7 疑問に思うところを兼田さんor青木先生に聞く。~
5/8 分らないところあったら、聞く。疑問無くす。英語の原稿...
5/9-5/11 発表練習 in English ~
5/14 ゼミで発表~
5/18 本番~
**サロメ [#p92d7318]
参考になる
[[Salome-Meca ビーム要素]]
*創造工房実習 [#p2f86cec]
**2022.11.4 viとgnuplotの練習[#ja9ade9f]
授業中に作ったグラフ↓ データはてきとう
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/matsuda/b...
11/4の課題↓
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/matsuda/g...
**2022.11.11 片持ち梁 Salome-meca2020 [#ja9ade9f]
ヤング率6000
ポアソン比0.4
,メッシュの長さ,要素数(ノード),変位,相対誤差%,計算者
,0.5,59504,6.56,-1.5,千代岡
,0.6,45512,6.48774,-2.69,高井
,0.7,39075,6.54132,-2.0,関合
,0.8,13397,6.43695,-3.5,岡田
,0.9,9903,6.36315,-4.6,松田
,1.2,6256,6.3043375,-5.4,青野
,1.3,5767,6.29784,-5.6,山口
,1.4,5146,6.286015,-5.76,山本
,1.5
終了行:
#contents
とりあえず理論通りの力で変形挙動を確認したい
*卒業論文 [#p2f86cec]
テーマ: 多々羅大橋の温度解析
やりたい事：ケーブルの固有振動数による張力推定→温度による...
日々の温度変化でケーブルがどう変化していくか。その後リス...
最初は↑この方向性だったが、だんだん変わってきた。最終的に...
実測平均応力変化と数値解析の平均応力変化を比較する。~
『温度変化がケーブルにどのような影響を与えるのか知る』
**冬休み前までにやりたい事 [#y8837db4]
☑実測ひずみデータがどのファイルなのか確認→すでにpython使...
・実測ひずみデータ整理して、時間-ひずみのグラフ作る~
・①サロメでサグ有モデル作る。（すでに作ったけどなんか違う...
・②軸力入れて、重力入れて解析。~
・31点の座標が、①と②で同じになるはず。近い値になればおっ...
・そのモデルに、実際の温度データを入れて、解析。~
＊最初の状態は、基準20℃ということにする。30℃にしたい時は...
・解析から、ひずみデータ得る。~
・実測ひずみデータ、解析ひずみデータ、を比較する。整合性...
・一応現状を英語でも説明できるように~
**夏休み前までやりたい事 [#a202b84f]
・DATファイルが何を示しているのか~
☑DATデータでグラフ作る~
・エクセルデータまとめて表作る~
✖・測定地点はどこ？24本あるはず。全体図にわかりやすいよう...
☑・ケーブルを一本取り出して、モデル作る。データからみてN6...
☑・★★★ケーブルに温度いれて解析する。とりあえずN6だけで。S...
✖・それが終わってから、全体のモデルつくる？そもそも今の私...
✖・全体モデルに温度入れる。
☑・温度変化によって張力がどう変化するかを知りたいの？で、...
**進捗 [#z78c72b3]
N6のケーブルの座標を資料から探し出す。解析基準点 or 図面...
けた側の解析基準点が資料に載っていない（多分）ので、図面...
図面上基準座標~
塔側(m)~
X=-447.336、
Y=-4.206、
Z=201.565~
桁側(m)~
X=-673.000、
Y=-11.839、
Z=43.499~
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/matsuda202...
N6の解析用ケーブルの長さ＝275.72m＝275720mm
ソリッド要素で解析する。被覆、仮想材料も入れてモデルを作...
青木先生と話して、ケーブル要素で解析することにした。
被覆ありでやるなら、ソリッド要素だけど、、どうするか。
結局、被覆ありでやりたいから、ソリッド要素でやることにな...
メッシュしてから結合。円柱の所からモデル作る。押し出して...
CSVとDATのデータは平均ではなくて、１分ごとや１０分ごとの...
6/17
この間ソリッド要素でやるということになったが、ソリッド要...
6/21
エクセルデータのまとめ方。~
A列をyyyy/mm/dd hh:mm:ssの表し方にする。ctrl＋１、ユーザ...
AA列には一分ごとの時刻を書く。秒のところは００にする。例...
AB２にはB列を１分ごとにした値を書く。例；AA２の式は、=AVE...
絶対に.xlsxで保存する！！！！！！！！！！！！！！！！~
62個のエクセルデータすべてやる。
6/23
ついにエクセルデータをすべて１分ごと平均して一つのファイ...
http://www.str.ce.akita-u.ac.jp/~gotouhan/matsuda2025/62k...
6/24
「多々羅大橋計測結果の概要と整理・分析方針(案)」より、
調査では、N6だけ温湿度を計測。N6は削孔位置が多く、ケーブ...
私の研究では温度についてやるので、N6のケーブルを取り出し...
＊Nケーブルは北側（上り・走行車線）、Sケーブルは南側（下...
＊N６は生口島側から６番目のケーブル~
DATファイルのF017, F110 がN6に関係している。 ~
・DATファイルがなんのデータなのかもう一度確認。多分温度。~
・DATを、csvにして(Python 使って)、サロメに入れる。~
・DATのデータでグラフ作りたい。前にも重くて落ちたから、ど...
・「石川島播磨技報」のpdf全てよめ。
6/25
先生と話した。DATファイルにも温度のデータないかも。N1,N6...
来週のゼミまでにやりたい↓
・もらったすべてのデータから温度のデータ探す。なかったら...
・DATファイルを自分で開いてグラフにできるように。Chat Gpt...
・DATのデータをCSVにする。Chat GPT使うとできるはず。~
・「多々羅大橋計測結果の概要と整理・分析方針(案)」をもう...
・「石川島播磨技報」を読む。~
・先週の英語ゼミで選んだ英語の論文が私がやろうとしている...
・サロメでN6のケーブルモデル作って、適当な温度入れて解析...
7/1
DATデータについて、
・4チャンネルある。~
・単位はマイクロストレイン~
・おそらくひずみのデータ。ってことは、応力が分かる。σ=E×ε~
7/3
DATにそもそもなんのデータが入っているのかまだわかんない。
7/6
10mのケーブルモデル作って、サグと温度適当に入れて解析した。
(0,0,0)と(10000,0,0)でケーブル。メッシュ100。~
(0,1000,0)と(0,-1000,0)で固定線。メッシュ15。~
固定部に壁作った。これがマルチファイバー要素。たぶん。~
重力、温度25℃(プレストレス−25)にして解析。~
サグなった。~
とりあえず温度入れて適当な値で解析できた。~
ここで、疑問。~
重力→解析する→サグ出てくる→その変位？の値使ってモデルつく...
(このようなやり方でやれば良いの？)~
このやり方でやろうとしたら、重力入れてプレストレスzeroに...
温度or張力を必ず入れないといけない？だとしたら張力どうや...
柴田さんの↓参考にしてやった。~
[1]まずサグのないまっすぐなcable要素を用意する
cable要素のみの場合、解析がエラーになってしまうため、固定...
[2]cableに張力、重力を与えたときの100点の座標を読み取る。~
[3]読み取った座標をプロットし(各点をt0〜t100とする)100点...
7/9 ↑これあんまり参考にしない方がいいらしい。
サグのやり方
１本のケーブル取り出す。(何要素？)
サグを手計算で求めて、、なんだっけ？
重力かける。
点に全て軸力を入れる。
モデル完成。
温度入れる。
7/11
以下柴田さんのやり方を教えてもらって真似てみる。
サグの入れ方。
サグモデルの座標、軸力の計算方法
高さ：h、径間長：ｌ、l/2におけるサグ：f、単位長さ重量：w ...
L：ケーブル形状長、x1 y1：下端、xm ym：中央、x2 y2：上端
$ \sqrt{L^2 - h^2}= 2C \sinh\left( \frac{l}{2C} \right), ...
これらの式から$ f - \frac{h}{2} + C \cosh\left( \frac{x_m...
これを解くとC(カテナリー数)が求まる。
ここで$ x_m = C\, \mathrm{arctanh}\left(\frac{h}{L}\right...
Cを求めるここができたら下端、上端のx,y座標を求める。x座標...
軸力は $ T = w y $で各要素ごとを求める。
１．モデル作成。点たちの座標を手計算で求めて、それらを線...
２．軸力~
３．重力~
４．温度。気温をエクセルで入れる~
５．解析~
(モデル作った後)と(軸力&重力入れた後)で、点たちの座標が一...
軸力を入れて上に上がり、重力を入れて下に下がるから、一致...
まず、何メートルの長さで解析するのか決める。N6の長さで最...
N6の正確な長さ→275622mm~
点たちの座標を手計算で求める。エックセルでやる。エクセル...
ex)１点目、
塔側(m)~
(x1,y1,z1)=(-447.336,-4.206、201.565)~
桁側(m)~
(x2,y2,z2)=(-673.000,-11.839,43.499)~
h=158066mm~
l=225793mm~
L=275622mm~
石川島播磨技報に色々と有益な情報がのっていたのでもっとし...
7/16
柴田さんのモデルはだいたい１mくらいで、３０個の点。~
私がやりたいのは２７５mである。何個の点を取れば良いか？
今は、サロメでモデル作るよりも、DATファイルと温度（エクセ...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/matsuda202...
7/17
今日はあとで自分の力でPythonを使ってDATファイルをなんとか...
7/18
DATファイルについて。~
計測器番号2_N6,N1.zipには、 F001.DAT～F159.DAT と F001...
計測器番号3_S41_S42.zipには、F001.DAT～F237.DAT と F001...
・chatGPTで.DATを.csvにした。F001_multichannel_time_align...
・python開く。（リボンみたいな△のアイコン）上にある。~
・「kouzou,デスクトップ,python,dat」の中にcsvgraph_hukusu...
・「kouzou,デスクトップ,python,dat」の中にF001_multichann...
・エクスプローラーのＤＡＴの風車みたいな水色のマークにcsv...
・pythonで以下のコードを打った。~
import pandas as pd
import matplotlib.pyplot as plt
import glob
import os
# 対象CSVファイルをすべて取得
file_list = glob.glob("F*_multichannel_time_aligned.csv")
# CH名の一覧を取得（最初のファイルから）
sample_df = pd.read_csv(file_list[0])
ch_columns = [col for col in sample_df.columns if col.sta...
# CHごとにプロット（1つの画像に全体表示）
for ch in ch_columns:
plt.figure(figsize=(15, 6))
for file_path in file_list:
df = pd.read_csv(file_path)
df['Time'] = pd.to_datetime(df['Time'])
label = os.path.basename(file_path).split("_")[0]...
plt.plot(df['Time'], df[ch], label=label)
plt.xlabel("Time")
plt.ylabel("Value")
plt.title(f"Comparison of {ch} across all files")
plt.legend()
plt.grid(True)
plt.tight_layout()
# 画像保存
output_filename = f"combined_{ch}.png"
plt.savefig(output_filename)
plt.close()
print(f"{output_filename} を保存しました。")
・ターミナルで ls →F001_multichannel_time_aligned.csv F00...
・ch1だけグラフを作れた。で、PCが固まった。~
・電源切った~
・labelがいらないので、そこのコード変えた。少しは軽くなっ...
for file_path in file_list:
df = pd.read_csv(file_path)
df['Time'] = pd.to_datetime(df['Time'])
#label = os.path.basename(file_path).split("_")[0...
plt.plot(df['Time'], df[ch])#, label=label)
・ch1とch2のグラフ作れた。で、またPCが固まった。~
7/19~
・今度は、ch3とch4のグラフだけ作るように指示する。以下の...
import pandas as pd
import matplotlib.pyplot as plt
import glob
import os
# 対象CSVファイルをすべて取得
file_list = glob.glob("F*_multichannel_time_aligned.csv")
# 対象チャンネルのみ（CH_3とCH_4）
target_channels = ['CH_3', 'CH_4']
# 各チャンネルごとにプロット
for ch in target_channels:
plt.figure(figsize=(15, 6))
for file_path in file_list:
df = pd.read_csv(file_path)
df['Time'] = pd.to_datetime(df['Time'])
# サンプリングなし・labelなしで描画
plt.plot(df['Time'], df[ch])
plt.xlabel("Time")
plt.ylabel("Value")
plt.title(f"Comparison of {ch} across all files")
plt.grid(True)
plt.tight_layout()
# 画像として保存
output_filename = f"combined_{ch}.png"
plt.savefig(output_filename)
plt.close()
print(f"{output_filename} を保存しました。")
・ch3とch4だけつくって、というコードなのに、ch1から作られ...
・次は、2つ目のフォルダにある237個のDATも同じようにやる。...
•csvへの変換は終わった。~
7/20~
重いと思うから、点だけにしても固まる。
色々試したが、すべてだめだった。
研究目的：加速度データから固有振動数を抽出し、ケーブル張...
加速度データ → 固有振動数 → 張力推定 → 温度変化との整合性...
１加速度データをFFT（高速フーリエ変換）すれば、ピーク周波...
２固有振動数がわかれば、張力Tを逆算できる。~
３これを温度変化と照らし合わせて、「温度が上がると張力が...
＊FFTとは、時間ごとに変化しているデータを、**「どんな周波...
気温が上がる → ケーブルが伸びる → 張力が緩む（応力が下が...
「整合とれるか確認」とは？~
加速度データから求めた張力の変化と、温度の変化に基づいて...
つまり、、~
気温の変化に対応して、ちゃんと張力（＝固有振動数）がそれ...
を確認すること。~
7/23~
pythonでch3だけといっても、何をしてもch1がつくられる現象...
！！実行する前にコードを保存する！！
フォルダ１のF001～F159とフォルダ２のF001～F237から、それ...
7/24~
①サロメで温度れて平均応力変化＝温度変化による理論的な平均...
②温度変化による理論的な平均応力変化⇔加速度データより推定...
「温度変化による理論的な平均応力変化」いらないなら、いき...
サロメで温度れて平均応力変化⇔加速度データより推定張力（整...
7/24~
先生と話した。こんな感じで進めてく。↓
今やる事~
１．カテナリー曲線でケーブル形状を定義~
→サグを含んだ状態でのケーブル座標を手計算で求める。~
２．FEMモデルをサロメで作成~
・上記の点群を使ってケーブルジオメトリを作成~
・ケーブルに重力を作用させる~
・軸方向に「事前張力（軸力）」を入れる（手計算で求め...
３．温度荷重を与えて静的解析~
・ケース１：＋30℃~
・ケース２：－10℃~
・結果として「一番サグがある部分のｙ軸方向変位」を取...
・測定地点だけの場合と、全ての要素点に同じ温度を与え...
４．実際の温度データを使った解析~
・時系列温度データをモデルに与える。測定地点だけ~
・すべての要素点に適用して解析もやる~
・ケーブルサグの変位と応力の時間的変化を観察~
・※これは動的ではなく静的解析？多分。（時間は進めない...
５．ひずみデータとの比較~
・実測されたひずみデータとFEMの出力を比較し、~
・FEMで得られた張力または応力分布との整合性をみる。~
まとめ~
・「温度変化によってケーブルのサグが変わる」~
・それは想像よりも大きい変位かも。~
・外力（風や荷重）を受けやすい状態になる。~
・ってことは、結果として、橋全体の挙動や安全性に影響を与...
→よって、「日々の温度変化を＊＊過小評価＊＊せず、リスク要...
サグの変化→ 温度変化が構造形状を変える要因となりうる~
応力の変化→ 材料の限界に近づくリスク、疲労寿命にも影響~
実測データとの比較→ モデルの妥当性確認＋安全性評価の信頼...
日々の温度変化→ 通常は見逃されがちだが、実は構造健全性に...
7/30~
h=158.066m
l=225.793m
L=275.622m
C=3800m
x1=24689.9m
xm=24802.8m
x2=24915.7m
y1=46307.11m
ym=46385.94m
y2=46465.18m
f=0.204714m
w=929.04N/m
y2-y1=158.066
ｘ(mm) y(mm)
0 0
7526.433333 5242.494193
15052.86667 10486.80519
22579.3 15732.9332
30105.73333 20980.87842
37632.16667 26230.64106
45158.6 31482.22134
52685.03333 36735.61944
60211.46667 41990.83558
67737.9 47247.86997
75264.33333 52506.72282
82790.76667 57767.39432
90317.2 63029.88468
97843.63333 68294.19412
105370.0667 73560.32283
112896.5 78828.27103
120422.9333 84098.03891
127949.3667 89369.6267
135475.8 94643.03459
143002.2333 99918.26279
150528.6667 105195.3115
158055.1 110474.181
165581.5333 115754.8713
173107.9667 121037.3828
180634.4 126321.7157
188160.8333 131607.8701
195687.2667 136895.8463
203213.7 142185.6444
210740.1333 147477.2648
218266.5667 152770.7075
225793 158065.9728
節点力（軸力）
T1 43021157.47
T2 43026027.96
3 43030900.13
4 43035773.99
5 43040649.55
6 43045526.79
7 43050405.71
8 43055286.33
9 43060168.64
10 43065052.63
11 43069938.32
12 43074825.69
13 43079714.75
14 43084605.51
15 43089497.95
16 43094392.09
17 43099287.91
18 43104185.43
19 43109084.64
20 43113985.53
21 43118888.12
22 43123792.4
23 43128698.38
24 43133606.04
25 43138515.4
26 43143426.45
27 43148339.19
28 43153253.62
29 43158169.75
30 43163087.57
31 43168007.08
8/2~
サロメでN6モデル作った。重力と軸力入れると赤になる。
8/4~
0℃で、計算で求めたサグ状態にする。~
94.8(ケーブルkgw/m)×275.622(ケーブル長さm)×9.8(重力加速度...
256063.8629÷275622＝0.92904(N/mm)~
これで求められた0.92904をモデル全体にかける。~
重力はかけない。軸力はかける。？
モデル全体に－10℃をかけて解析する。~
モデル全体に30℃をかけて解析する。~
両方まわった!!!
とりあえずここまででプレゼンまとめる
メモ：.hdfで開く。.commで開く。cableN60804.comm
モデルはN6のケーブル。ケーブル要素でやった。両端はピン。...
10/14
中間発表後の頭の中整理。~
中間発表→加速度データによるケーブルの「張力」と、実測温度...
しかし、方法変更。~
今後→実測データの「ひずみ」と、実測温度データを入れてFEM...
★元々の方法~
・固有振動数f1をFFT解析→弦の理論式T=(2Lf1)^2μから張力推定~
・FEM解析で実測温度データいれて、ケーブル要素の平均応力を...
・両者の「変化量」や「比例関係」を比較して整合性みる~
問題点~
・弦の理論式ではサグが考慮されていない（直線弦仮定）。~
・実際のケーブルではサグによって張力分布や剛性が変わる。~
・したがって、「張力 vs 応力」の比較は前提条件が異なるた...
★修正版~
・実測ひずみデータ整理~
・FEM解析で実測温度データいれて、ケーブル要素のひずみを算...
・両者の「ひずみ（実測） vs 温度によるひずみ（解析）」を...
要するに、、、~
「実測および解析における温度—ひずみ関係の比較」をする！！！
「張力と応力の整合性を見る」から「ひずみの整合性を見る」...
「比較する物理量」が変わるが、「目的（温度応答を再現する...
むしろサグを考慮できる分、より正確で現実的な方法になる。
温度変化がケーブルの構造挙動に与える影響を、より現実的に...
実橋の“温度によるひずみ応答”を、サグを考慮したFEMモデルで...
補足：「張力」は結局「応力 × 断面積」なので、応力もひずみ...
→ より一次的で確実な指標にシフトしただけ。
ところで、ひずみデータはどのファイルにあるのか？~
csvではなかった。しかし、表３のデータ諸元ではcsvと書いて...
グラフと値的にDATっぽい。~
しかし、「多々羅大橋計測結果の概要と整理・分析方針（案）...
しかも、DATファイルから作った１〜４chのグラフと「多々羅大...
10/15
DATファイルから作った１〜４chのグラフはひずみデータだった。
10/19
今、サグ有モデルを作って、重力と軸力入れて解析したい。温...
カテナリーのサグ有ケーブルは20℃の時の形であるー１０℃とか3...
そこで、サロメの重力と軸力だけ入れて解析する時は、温度は2...
→温度＝20℃に“する”というより、“ΔT=0にする”（温度荷重を無...
10/20~
①「サグ有モデルを作って、重力と軸力だけ入れて解析して（温...
温度を入れない場合（Function and Lists→DEFI_FONCTION→VALE...
0.001℃とかにしても回らなかった。今度1℃でやってみよう。
それか、重力を等分布荷重として入れてるから、ピンの反力がg...
★そもそも「サグ有モデルを作って、重力と軸力だけ入れて解析...
・モデルが正しく初期平衡状態を再現できているか（＝初期張...
・温度などの外力を考える前に、ケーブルが“自然に釣り合った...
・温度を入れると熱ひずみも影響するため、初期形状確認が難...
②ケーブル要素でひずみを出すという件について~
サロメは20℃基準になってるから、20℃との差を考えて入力すれ...
サロメで、ケーブル要素だとひずみが出力できない。Post Proc...
σ=F/A=E・ε
軸力は出せるので、この↑式を使って、軸力からひずみを求める...
もし、サロメでひずみが出せるようだったらそっちの方が良い...
★次のゼミまでにやる事~
・「サグ有モデルを作って、重力と軸力だけ入れて解析して（...
・温度入れて解析して軸力出す。~
・そこから計算してひずみ求める。~
・軸力の入れ方、重力の入れ方を、もう一度確認して、ほんと...
10/23~
①温度の変化を0.004まで小さくしたら解析回った。0.003はダメ...
しかし、反力の値がおかしい。~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/co...
10/29
どうしたら重力と軸力いれて、反力がgl/2に近い値になるのか。
10/30
そもそも反力はgl/2になれば良いのか？構造力学の式使ったけ...
反力は一旦無視して、とりあえず座標が同じになれば問題なし。~
★重力の入れ方について~
94.8(kgw/m)×275.622(ケーブルの長さｍ）×9.8(m/s^2)＝256063...
256063.8629（Ｎ）÷275622(ケーブルの長さｍｍ）＝0.92904(N/...
FORCE_NODALE FY＝ー0.92904~
しかし、解析は回るけどケーブルが上に上がったり変位が大き...
上の計算式がおかしいかもという事で、何か変えてみる。点が3...
FORCE_NODALE FY＝ー８２６０~
これでも変位が大きかったりなんかおかしい。
★軸力の入れ方について~
今まで、カテナリーの計算でつかうW単位重量が929.67042N/mで...
1 kgw = 9.80665 Nだから、~
94.8（kgw/m）=94.8×9.80665（N/m）=929.67042（N/m）
しかしケーブル図の資料に94.8kgw/mとあって、929.67042じゃ...
そしたら、当たり前だけどｗが一桁減ったから軸力も1桁減った...
で、解析回って、少し下に下がった形になった。~
考えなおしたら、94.8じゃなくて元々の929.67042N/mであって...
１０/３１にやりたいこと。
いろんなものが少しずつ値が間違っていたりして、座標も違っ...
軸力と重力の入れ方もまだ正解がよく分からない。とりあえず...
反力がgl/2になるのかというのも先生に聞いてみる。
10/31
★重力の入れ方について~
94.8(kgw/m)×9.80665(m/s^2)＝929.67042(N/m)~
929.67042(N/m)×275.6418685(m)=256256.0917(N) ２５万N~
256256.0917(N)÷３１(個)＝8266.325537(N)~
FORCE_NODALE FY＝ー8266.325537としたいところだが、
31点すべて同じ −8266 N はNGです。理由は「等分布荷重を節点...
FEMは要素と節点で計算するので、「連続的な荷重」を
節点に置き換えて（＝節点荷重）あげる必要があります。その...
・端の節点は要素が1本しかつながっていないのに、2本分の荷...
→ つまり「端の半分の分布荷重を2倍入れている」ことになる。~
・中間節点は左右の要素2本から荷重が来るのに、1本分しか入...
→ 「中央部分が軽く、端だけ重い」不自然な分布になる。~
結果：自重の重心が本来より両端寄りになり、支点反力のつり...
→ これが「端点は半分にしないとNG」の理由です。~
よって、自重の入れ方は、~
kotei1 FY=-4270.934862N~
ten2～ten30 FY=-8541.869723N~
kotei2 FY=-4270.934862N~
★軸力の入れ方~
結局、w=929.67042（N/m）で計算していくことにした。そうす...
FORCE_NODALE 点荷重 で与えると、「節点に外力（集中荷重）...
・方向のばらつきが大きくなる（サグ形状ごとに角度が違うた...
・外力として反力に現れる（内部初期張力ではない）~
・境界条件に依存して全体釣り合いが崩れる~
という理由で、解析結果が不安定になりやすいそうだ。~
そこで、N_INIT（初期軸力）をつかう。
「軸力を入れる」＝ N_INIT~
「外力を入れる」＝ FORCE_NODALE~
$ \quad H = \frac{wl^2}{8f} $
H=3360502.585N 300万N~
よって、N_INIT＝3360502.585~
ここで疑問。~
Q鉛直方向の張力はいれなくてよいのか？~
A入れなくてよい。鉛直成分は、自重とケーブルの傾きから、解...
反力FYの両端を足した値が総自重になればよい。256256.0917N
ケーブル全体の自重のつり合いから：両端の鉛直反力の合計=𝑤𝑆~
つまり、どんな形状でも 左右反力の合計 = 総自重。
（PESANTEURを使ってもFORCE_NODALEでも同じ）
11/1~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/N6...
11/18~
結局軸力と重力入れて座標が元のところに戻らなかったので、...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/N6...
↓これは一番上から、最初の座標、－30度、－10度、0.004度、1...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
11/20~
実データを入れて解析する。ひずみをみたい。
***11/25二段階解析スタート [#q450bdbc]
エクセルの温度データをサロメに入れてできる？？→できない。...
先生からのアドバイスで二段階解析をした。一段階目は０から2...
前まではSTAT_NON_LINE（静的）でやっていたが、DYNA_NON_LIN...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
二段階解析の一段階目のcomm
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/2d...
二段階解析の一段階目のcomm
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/2d...
凄い振動しているような感じでなんか嫌な予感がするが、この...
***11/26 [#j6ad4477]
昨日、5個の温度データで二段階解析できた。全部で800個（約6...
１～５６ のデータは解析回った。~
2024/10/23 14:46:04のデータを１として、2024/10/23 23:56:0...
１～１２８ 解析回らない。~
JUSQU_A=76201, PAS=127~
エラーコード
「Arrêt pour cause de matrice non inversible. ║ ║ La base...
「逆行列が存在しないため停止しました。 ║ ║ グローバルベー...
半日分が限界？
とりあえず、23日の午後データ（５６個の温度データ）での解...
温度変化とten16の座標と変位のグラフ~
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
さっきまで、56個の温度データがあるとき、PAS=55 としてい...
listr2 = DEFI_LIST_REEL(
DEBUT=1.0,
INTERVALLE=_F(JUSQU_A=33001.0,
PAS=600.0)
)
DEFI_LIST_REEL(
DEBUT=開始値,
INTERVALLE=_F(
JUSQU_A=終了値,
PAS=刻み幅
)
)
「1〜33001秒までを 600秒（10分）刻みで 56個」**「10分間隔...
ということで、PAS＝600にして、もう一度23日の午後のデータ...
***11/27 [#deae7b5e]
PAS=600で２３日の午後、解析できた。~
PAS=600で２３日の午後と２４日の午前も含めて、解析失敗。エ...
「<SolverError> <MECANONLINE9_4> ...
║ ...
║ Arrêt pour cause de matrice non inversible. ...
║ La base globale est sauvegardée. Elle contient les ...
║ ...
║ Conseils : ...
║ - Vérifiez vos conditions aux limites. ...
║ - Vérifiez votre modèle, la cohérence des unités...
║ - Si vous faites du contact, il ne faut pas que ...
║ contact. ...
║ ...
║ - Parfois, en parallèle, le critère de détection...
║ pessimiste ! Il reste néanmoins souvent ...
║ possible de faire passer le calcul complet en ...
║ la valeur du mot-clé NPREC) ou ...
║ en le débranchant (valeur du mot-clé NPREC=-1)...
║ processeurs. ...
」
PAS=600で、２４日の午前だけ、解析できた。~
PAS=600で、２４日の午後だけ、解析失敗。
<SolverError> <MECANONLINE9_4> ...
║ ...
║ Arrêt pour cause de matrice non inversible. ...
║ La base globale est sauvegardée. Elle contient les ...
║ ...
║ Conseils : ...
║ - Vérifiez vos conditions aux limites. ...
║ - Vérifiez votre modèle, la cohérence des unités...
║ - Si vous faites du contact, il ne faut pas que ...
║ contact. ...
║ ...
║ - Parfois, en parallèle, le critère de détection...
║ pessimiste ! Il reste néanmoins souvent ...
║ possible de faire passer le calcul complet en ...
║ la valeur du mot-clé NPREC) ou ...
║ en le débranchant (valeur du mot-clé NPREC=-1)...
║ processeurs.
２４日PMだけの温度データ~
１段階目comm
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/11...
２段階目comm
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/11...
２３PM~
一段階：０～25.39~
二段階：25.39～18.96~
２４AM~
一段階：０～18.96~
二段階：18.96～29.68~
２３PMと２４AMを別々に、それぞれ二段階解析した。２３PMの...
ちなみに、ひずみは結果の EFGE＿ELNO を使って、ヤング率E=1...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
***11/28 [#af6b30e9]
ここで、前に諦めた自重と軸力いれて元の座標に戻らない件に...
★★自重はGRAVITYでいれる。
★★温度変化（０→0.004）なしで、自重＋軸力いれると、元の座...
GRAVITY つかって、０→0.004で解析回ったcommファイル→
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/11...
ピンと張ったケーブル（直線）→ ten16 (112896.5, 79033)~
カテナリー計算後サグ有モデル→ ten16 (112896.5, 77271.98)...
GRAVITY（自重）＋ 各節点の FORCE_NODALE をかけてSTAT_NON...
→ 直線から見ればトータルで 1.7 + 2.6 ≒ 4.3 m サグ、みたい...
***12/2 [#gd2c39b4]
２４PMが解析できない件について。
「<Erreur> La matrice du système est singulière
<Action> On essaie de découper le pas de temps.
On utilise la découpe manuelle.
Découpe uniforme à partir de l'instant < 3.723850000000e+...
(soit un incrément constant de < 9.37...
Le nombre maximal <3> de niveaux de subdivision est attei...
<Action><Échec> Échec dans la tentative de découper le pa...
Temps CPU consommé dans ce pas de temps : 0.046 s
Nombre d'itérations de Newton : 4
Temps total intégration comportement : 0.00...
Temps total factorisation matrice : 0.00...
Temps construction second membre : 0.01...
Temps total résolution K.U=F : 0.00...
Temps assemblage matrice : 0.00...
Temps autres opérations : 0.01...
Mémoire (Mo) : 774.51 / 770.45 / 39.40 / 35.75 ...
このエラーがでたので、< 3.723850000000e+04> 秒の所を変え...
回った。２３PM、２４AM、２４PMを合わせて解析して、またさ...
今度は、２３AM～２９AMまですべて温度データを入れて、Value...
同じようにエラー。
Valueを60（1分毎）にした.
同じようにエラー
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
Valueを60（1分毎）にした.
同じようにエラー
***12/4 [#x9af0760]
・334819秒あたりでエラーが出てるので、今度はVALUEを６０か...
・しかし、これがいないような秒数の所でエラーも出たので、...
・DYNA＿NON＿LINE CONVERGENCE ITER_GLOB_MAXI＝５０だっ...
・いつものエラー~
・作戦を前に戻して、エラーが出るところのVALUEの数を小さく...
・１～37201秒 VALUE＝６００、 ～333601秒 VALUE＝６０、...
・解析回った！~
このときのcommファイル↓~
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/12...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momoka/12...
VALUEはPASのこと
PAS＝１０に所はケーブルの揺れ幅が大きすぎて変位とひずみが...
そのグラフがこれ↓
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
***12/9 [#w887599b]
12/4の結果がおかしかったので、PAS＝６０にして、エラーが出...
10月23日～27日のデータの解析結果↓
グラフ
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
このグラフの続き（10月27日～29日）を別の解析でやる。一回...
2段階解析を2回に分けてもいつものエラーが出てPASを６０から...
301000～320000秒くらいでいつもエラーでるのは、そこら辺で...
この問題は一旦無視していまある解析結果（10月27日～29日）...
***12/10~ [#j2c605e2]
ひずみは割と妥当な値がでたと思う。。
12/12~
実際のひずみとおもわれるデータとサロメから得たひずみのグ...
12/17~
先生と話して、いまあるデータで中間発表書いていく。
12/21~
スライドつくり。概要もつくり終わりたい。
12/22~
実測ひずみだと思っていたデータは平均応力変化だった。スラ...
解析ひずみ vs 実測ひずみ~
↓~
解析平均応力変化 vs 実測平均応力変化
***12/23中間発表 [#o4067663]
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
1/27~
卒論かきはじめなきゃやばい
***2/5~ [#s5634d8d]
ひずみを算出して、そこから平均応力変化を求め、実測データ...
・他の人が温度変化によるケーブルの変化の研究をやっていな...
・文献をもっと探してみて。似たのがなかったら、なぜないの...
・なぜ他の人が温度に注目していないのか、というのが大事~
・ピンピン以外をやらなかった理由~
・生のデータを分析した。これはいただいた貴重な生データ。...
・多々羅大橋のその日の天気と気温を調べておく~
・ケーブルをどうやって再現したのか、ちゃんと説明する~
・サグをどういうふうにいれたのか~
・サグ量だけでみたら温度によって変化することが分かる。こ...
・解析でサグが変化するのはみれた。温度とまあまあ似たグラ...
・橋全体を作るべきだったかも。これは今後やっていきたい~
・大量のデータ解析から始まり、研究した。~
・モデルの具体的な作り方、そして、温度入れたときの揺れる...
・他の人がサグを考慮しない理由。まだ文献調べていないので...
・ストーリーをうまくつくれ。目的、背景をもっと深く話す。...
・とにかく温度計の文献をさがせ。~
・もしかしたら、橋全体でやったら、そりゃ違うようになるよ...
・私みたいに、結果が合わない理由は沢山ある。~
・あまり脱線はしないように。~
今日、タイトル提出する
***2/16~ [#n8cd7642]
平均応力変化の平均とは、そもそもなに
2/17~
卒論の前夜というか当日の朝、平均応力変化の『平均』はいら...
2/19~
私の卒論の考察で言ったのは、拘束条件を考慮した解析の実施。~
では、今後、拘束条件をどうすればよいのか？~
ピンピンでダメなら、固定にすればいいわけではない。固定は...
まぁ、色々言わずに時間かかるけど塔桁他ケーブル作ればいい...
でも、全体モデル作ったら、そのケーブル全てにサグ考慮する...
ホンモノの多々羅を一回見てみたいような気もする
***2/20~ [#x7e6041b]
**一旦あたまの中整理 [#rf5bbb5d]
やりたいこと！！！
★★軸力と重力入れて、解析前と解析後の座標が同じになりたい...
N6について、ケーブル図より~
公称断面積 11400 mm^2 ~
単位重量PE被覆後 94.8 (kgw/m) =0.0948 (t/m) =929.67042 (...
死荷重完成時張力 516.1(t) =50580.11165(kN) ~
カテナリー計算によるケーブル形状長 276.315(m) ~
計算で求めた~
高さ ｈ=158.066m~
径間長 l=225.664m~
カテナリー Ｃ＝4417.523308~
形状長 Ｌ＝275.6418685m →これが「カテナリー計算によるケ...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
&ref(http://www.str.ce.akita-u.ac.jp/~gotouhan/j2025/momo...
サロメ座標
x軸(mm) y軸(mm)
kotei1 0 0
ten2 7526.433333 5042.915455
ten3 15052.86667 10101.27387
ten4 22579.3 15175.08991
ten5 30105.73333 20264.37833
ten6 37632.16667 25369.15388
ten7 45158.6 30489.4314
ten8 52685.03333 35625.22573
ten9 60211.46667 40776.5518
ten10 67737.9 45943.42455
ten11 75264.33333 51125.85898
ten12 82790.76667 56323.87013
ten13 90317.2 61537.4731
ten14 97843.63333 66766.68302
ten15 105370.0667 72011.51507
ten16 112896.5 77271.98447
ten17 120422.9333 82548.1065
ten18 127949.3667 87839.89646
ten19 135475.8 93147.36972
ten20 143002.2333 98470.54169
ten21 150528.6667 103809.4278
ten22 158055.1 109164.0436
ten23 165581.5333 114534.4046
ten24 173107.9667 119920.5264
ten25 180634.4 125322.4246
ten26 188160.8333 130740.1149
ten27 195687.2667 136173.613
ten28 203213.7 141622.9347
ten29 210740.1333 147088.0959
ten30 218266.5667 152569.1124
kotei2 225793 158066
軸力[N]
1 4941134.759
2 4945823.008
3 4950525.614
4 4955242.591
5 4959973.952
6 4964719.711
7 4969479.881
8 4974254.477
9 4979043.513
10 4983847.002
11 4988664.957
12 4993497.395
13 4998344.327
14 5003205.769
15 5008081.734
16 5012972.237
17 5017877.292
18 5022796.912
19 5027731.113
20 5032679.909
21 5037643.313
22 5042621.341
23 5047614.007
24 5052621.325
25 5057643.31
26 5062679.976
27 5067731.339
28 5072797.412
29 5077878.211
30 5082973.749
31 5088084.043
salome軸力[N]
鉛直ｙ軸方向 水平ｘ軸方向
1 2750391.31 4104894.682 kotei1
2 8419.779663 -4.110606786 ten2
3 8427.776392 -4.102801731 ten3
4 8435.797621 -4.094996054 ten4
5 8443.843369 -4.087188698 ten5
6 8451.913661 -4.079379724 ten6
7 8460.008521 -4.07156985 ten7
8 8468.127973 -4.063759074 ten8
9 8476.272037 -4.055945962 ten9
10 8484.440741 -4.048133457 ten10
11 8492.634106 -4.040318667 ten11
12 8500.852157 -4.032504234 ten12
13 8509.094917 -4.024688822 ten13
14 8517.36241 -4.016871991 ten14
15 8525.65466 -4.009055587 ten15
16 8533.971693 -4.00123887 ten16
17 8542.313529 -3.993421056 ten17
18 8550.680196 -3.985603436 ten18
19 8559.071716 -3.977785623 ten19
20 8567.488115 -3.969968141 ten20
21 8575.929416 -3.962149828 ten21
22 8584.395644 -3.954332314 ten22
23 8592.886824 -3.946515071 ten23
24 8601.402979 -3.938697543 ten24
25 8609.944137 -3.93088149 ten25
26 8618.510318 -3.923064636 ten26
27 8627.101551 -3.915249141 ten27
28 8635.717859 -3.907434308 ten28
29 8644.359268 -3.899620345 ten29
30 8653.025801 -3.891806481 ten30
31 -2997901.667 -4104778.646 kotei2
*サロメ [#s104c0ef]
手動で選択する場合,コンテンツにチェックを入れた後,グルー...
Compound_Mesh_1だけを目つけてメッシュの所は目消す。そして...
Paravisにおいては、WarpByVectorを使うためには、オブジェク...
Geometry→31個の点の座標入力→31点を線で結ぶ~
Mesh作成（メッシュは点ではなくて線だけ）~
コンパウンドの作成→Compound_Mesh_1→メッシュの所に30個入れ...
Compound_Mesh_1右クリック→グループを作成→ノードはkotei1,t...
Asterstudyをやってく
*python [#w177c7bf]
サンプリング間隔→
legend→複数のデータを同じグラフに描画したときに、それぞれ...
label=“〇〇” 線につける名前（ただしこれだけでは表示され...
legend() 上でつけたlabelをグラフに表示する命令
*日記[#x22e75ea]
11/29~
11/28~
サロメ。先が見えない
11/27~
salome。半日エクセルで作ったデータ消えた。ので寝る
11/26~
salome
11/25~
二段階解析について知る。
11/21~
先生と進捗と今後どうしてくかはなし
11/1~
昼夜逆転というか朝と昼もわからなくなってきた
10/31~
今日はサロメはやってない。重力と軸力について勉強した。
10/30~
座屈の計算。サロメとか計算とか皆に助けてもらった。明日か...
10/29~
座屈のデザイン決めた。どうしたら重力と軸力いれてうまく解...
10/22~
サロメの続き
10/20~
サロメでつまって、皆に教えてもらった。先生に相談
10/17~
サロメ。クマが怖いから夜来るのはやめたい
10/15~
サロメ。いつになったらサロメになれるのか
10/14~
頭の中整理
10/7~
中間発表。夏休み全力で楽しんだのでこれからがんばる
10/4~
概要を書き始めた。
8/4~
先生に相談。１Nでもいいからかけてみる。重力と軸力なしで。...
8/2~
サロメでN6のモデル作った。軸力と重力を入れると解析回らな...
7/30~
サグケーブルのつくりかた教えてもらった。
7/27~英語の論文もう一回読む。「Temperature Effect on Tens...
Cable of Cable-Stayed Bridge」。 で、わからない単語とかな...
7/26~
気力あったらサロメやりたい。
7/25~
日本ファブテックさんの工事見学
7/24~
先生とはなした。
7/23~
英語で中間発表。みんなの発表が理解できるくらいになりたい...
DATのグラフやっと完成。石黒さん神。
7/22~
瀧上工業株式会社さんのお話聞いた。スライド作り。発表練習
7/21~スライドつくり。ねむい
7/20~
発表までにグラフが作れる気がしないから、一旦諦めてスライ...
7/19~
オーキャン。pythonやったけど欲しいグラフができあがらない...
7/18~
オーキャン準備。テント。今日はpythonの使い方を教えてもら...
7/16~
オーキャン準備。温度と湿度のグラフ作り。暑い。海いきたい
7/14~
オーキャンのポスター作り
7/13~
オーキャンの動画とポスター
7/12~
サグ手計算
7/11~
キャンパスクリーンデー。草取りした。サグのやりかた教えて...
7/10~
オーキャンの準備。ポスターとか。
7/6~
サロメ教えてもらった。感謝。一人でやってたら分からなすぎ...
7/3~
データとたたかった。何がわかんないかもわかんなくなってきた
7/1~
今の感じだと何の成果もないのでやばい。夏休みに研究室これ...
6/27~
サロメで適当にモデルつくろうとして、途中。誰かにサロメの...
6/24~
Drop boxを見てたら、石川島播磨技報のファイルがあることに...
6/23(月）~
やっとやっとエクセルまとめ終わったら、温度のデータなかっ...
6/22(日）~
エクセルやってる。嫌いになりそう。絶対pythonでやったほう...
6/21(土)~
久しぶりに何もない日で研究室きた。エクセルの件振り出しに...
6/17(火)~
明日発表のパワポ作った。
今週もゼミでいう事がなくてやばい。
6/16(月)~
IHIさんの貴重なお話をきいた。
6/15(日)~
エクセルデータをやっと全て平均できたのに、保存してその後...
6/10(火)~
エクセルのデータ平均してグラフ作るの終わらせないとそろそ...
エクセルデータについて、一つのファイルに12万行くらいある...
5/26(月）~
研究室には来た。
5/23(金)~
データやろうとしたけどやめた。
5/22(木)~
先生と進捗話した。
5/20(火)~
N6ケーブルの両端のの座標が分かったので、モデルをつかろう...
5/19(月)~
ケーブル一本取り出して解析するために、資料から座標読み取...
5/18(日)~
YEC本番。リスクアセスメントの発表。２個くらい質問されたけ...
5/12(月)〜5/16(金）~
毎日YECの発表練習。~
prof Daiaが研究室きた。YECの発表も見てもらった。~
5/11(日）~
スライドを修正して発表練習した。けっこうやばい。~
5/9(金）~
YECの原稿なんとなくできた。~
5/7(水)~
兼田さんに英語版スライドもらい、分からないところも聞いた...
たたらの方は、DATをPythonで開いてグラフにすることができた...
エクセルのデータは、今は保留。~
明日、N6ケーブルのモデルを作りたい。~
5/6(火)~
プレゼンができるようになるまで、たたらは一旦ストップで、Y...
英語の原稿作り始めた。~
ALARP領域についてもう少し詳しく調べる。~
5/5(月）~
サロメでN6のモデルを作りたい。しかし、設計図がうまく理解...
兼田さんの卒論スライドを見て、わからないところを無くす。~
YECで英語が理解できないということになりたくないので、英語...
5/3(土)~
「多々羅大橋計測結果の概要と整理・分析方針(案)」を」読ん...
調査では、N6だけ温湿度を計測。N6は削孔位置が多く、ケーブ...
私の研究では温度についてやるので、N6のケーブルを取り出し...
＊Nケーブルは北側（上り・走行車線）、Sケーブルは南側（下...
＊N６は生口島側から６番目のケーブル~
DATファイルのF017, F110 がN6に関係している。 ~
5/2(金)~
教えてもらいながら、適当にケーブル一本使ってケーブル要素...
python勉強中。~
エクセルデータの処理いついて、一旦保留。どうすればいいか...
兼田さんの卒論動画をみてどんな風に発表するか、考え中。（...
5/1(木)~
サグの実験見学した。pythonの本を読み始めた。
4/30(水)~
石黒さんと話して、結局pythonでデータをまとめて表にするこ...
たたらの細かい図を先生からもらい、正確な寸法がわかる。明...
DATファイルが文字化けする件について、ChatGPTにファイルを...
4/29(火)~
自分のパソコンだとpython入れられなかったので、研究室のや...
サロメを使って、たたらの一部分のモデルを作ってみようとし...
4/28(月)~
データの単位にkineとgalがあるのでgalに統一する？
横軸時間で縦軸を色々変えてグラフを作りたい。その為に、ま...
4/27(日)~
多々羅大橋ケーブル張力測定の調査資料を読んだ。~
HEDファイルはDATファイルの説明が書かれている？先生に聞く。~
自分の卒論とは別で兼田さんの卒論を英語で理解しようとして...
4/25(金)~
頂いた大量のデータに目を通した。DATファイルを開くと文字化...
*Risk Assessment [#eef9f3d5]
私の卒論とは関係ない。ただ自分で調べただけ。興味があると...
-兼田さんの[[卒論日誌:https://www.str.ce.akita-u.ac.jp/cg...
-ALARPは「As Low As Reasonably Practicable」の略称であり...
-[[H29年度厚生労働省委託事業 建設工業の設計段階における...
日本では、施工時の安全は主に施工業者の責任範囲であり、設...
-海外では事故が起こった時に通行止めや封鎖をしたくない。日...
Q User risk と Road Block risk を求める時に、色んなデー...
→A リスクの数字が小さくなるのは事実。数字が小さいほど安...
→Q そうなってくると、ARARP領域の下にあるリスクは、起こっ...
A
Q リスクアセスメントやるのは、結構大事そうな気がするけど...
A
Q どの国でリスクアセスメントが良く使われている？~
A
Q 日本では海外よりもリスクアセスメントが普及してないとい...
A
Q 日本ではいつから努力義務になった？いつから義務化された...
A
10/22↓
★リスクアセスメントとは、~
将来起こりうる危険（リスク）を事前に特定し、その発生の可...
★通常、3ステップで構成される。~
１リスクの特定 Risk Identification~
２リスクの分析 Risk Analysis~
３リスクの評価 Risk Evaluation~
★リスクの定量化の考え方~
リスク＝発生確率×影響の大きさ~
例：「年に1回起こるが、被害が小さい」 → 低リスク~
「100年に1回だが、崩壊する」 → 高リスク~
この評価をマトリクス化して、「どの範囲なら許容できるか（A...
★ALARP（アラープ）の原則~
「As Low As Reasonably Practicable（合理的に実行可能な限...
完全な安全を求めるのではなく、費用対効果のバランスを考え...
これはヨーロッパの安全工学・原子力・航空分野で広く用いら...
★リスクアセスメントが必要な理由は、~
・「感覚」ではなく「数値」で安全を判断するため~
・潜在的な危険を事前に発見して、災害や事故を未然に防ぐた...
・限られた予算や時間を有効に使うため。優先順位づけにより...
・説明責任（Accountability）を果たすため。社会的・法的に...
★リスクアセスメント＝危険を見つけて評価する~
リスクマネジメント＝その結果をもとに行動する~
★ここで、リスクアセスメントの問題点は、~
①定量評価の不確実性。データの不足。「100年に一度の津波」...
②評価者による主観の影響。リスクアセスメントは定量的に判断...
例：被害「重大」＝ 死亡1人？10人？経済損失10億円？~
過去データが少ない時、何年に1回発生すると想定するか？~
③リスク間の比較・優先順位づけの困難さ。異なる種類のリスク...
例：（例：落石による通行止め vs 地震による崩壊）~
④コストと安全のトレードオフ問題~
安全を高めるほどコストが上がる。安全対策を重視すると、施...
⑤社会的・心理的リスクの扱いが難しい。~
（例：新幹線の微小な地震揺れ、トンネル崩落事故後の不安）~
「感じるリスク（perceived risk）」と「実際のリスク（real ...
⑥結果の扱い方（運用）の難しさ。リスクアセスメントを実施し...
⑦時間的変化を考慮しにくい。多くのリスク評価は「ある時点の...
★
*メモ [#e50c8b22]
**YEC [#s047c458]
5/2 卒論動画見る。~
5/4 日本語版スライドで、分からないところまとめる。計算あ...
5/5 兼田さんにもらった資料全部読む。井上さんのやつも全部...
5/6 英語版概要、英語のまま理解する。英語で発表する原稿つ...
5/7 疑問に思うところを兼田さんor青木先生に聞く。~
5/8 分らないところあったら、聞く。疑問無くす。英語の原稿...
5/9-5/11 発表練習 in English ~
5/14 ゼミで発表~
5/18 本番~
**サロメ [#p92d7318]
参考になる
[[Salome-Meca ビーム要素]]
*創造工房実習 [#p2f86cec]
**2022.11.4 viとgnuplotの練習[#ja9ade9f]
授業中に作ったグラフ↓ データはてきとう
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/matsuda/b...
11/4の課題↓
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2022/matsuda/g...
**2022.11.11 片持ち梁 Salome-meca2020 [#ja9ade9f]
ヤング率6000
ポアソン比0.4
,メッシュの長さ,要素数(ノード),変位,相対誤差%,計算者
,0.5,59504,6.56,-1.5,千代岡
,0.6,45512,6.48774,-2.69,高井
,0.7,39075,6.54132,-2.0,関合
,0.8,13397,6.43695,-3.5,岡田
,0.9,9903,6.36315,-4.6,松田
,1.2,6256,6.3043375,-5.4,青野
,1.3,5767,6.29784,-5.6,山口
,1.4,5146,6.286015,-5.76,山本
,1.5
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
