服部の卒論日誌
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
|
ログイン
]
開始行:
#contents
*卒論テーマ [#id2b2754]
CLT床版用の防護柵の開発（安部さんサポート）旭川の実験モデ...
2つに割れている鞘管が1つずつ順番に降伏することで、どれぐ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
↑防護柵 図面
・小川さんの論文 CLT床版への取付・交換を考慮した鋼製防護...
https://doi.org/10.11532/structcivil.70A.885
・小川さん 卒論：CLT床版に取り付けられた鋼製防護柵の性能...
・小川さん 修論：CLT 床版への取付・交換が容易な鋼製防護柵...
・安部さん：支柱交換を考慮したCLT床版用鋼製防護柵の挙動（...
*4年次 進捗[#o6991381]
**2月5日 [#k16c1e69]
! <BoucleGeometrieError> <MECANONLINE9_9> ...
! ...
! Arrêt par échec dans la boucle de point fixe sur ...
! La base globale est sauvegardée. Elle contient le...
<ループ・ジオメトリー・エラー> <MECANONLINE9_9> ！
! !
! ジオメトリの固定点ループの失敗で停止 ！
! グローバルデータベースが保存される。停止前のステ...
**1月26日 [#scf4a32b]
角柱のモデルを正12角形から正24角形に変更して解析を回して...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
前回の正12角柱と比較したものが
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
スリットサイズが正12角柱が30mmで正24角柱が24mmになってい...
この正24角柱のモデルに支柱をかぶせて拘束荷重を与えたモデ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
ただこれはエラーが出てしまった。
! <MatriceSinguliereError> <MECANONLINE9_4> ...
! ...
! Arrêt pour cause de matrice non inversible. ...
! La base globale est sauvegardée. Elle contient le...
! ...
! Conseils : ...
! - Vérifiez vos conditions aux limites. ...
! - Vérifiez votre modèle, la cohérence des unit...
! - Si vous faites du contact, il ne faut pas qu...
! ...
! - Parfois, en parallèle, le critère de détecti...
! possible de faire passer le calcul complet e...
! en le débranchant (valeur du mot-clé NPREC=-...
<MatrixSinguliereError> <MECANONLINE9_4> ！
! !
! 非可逆行列のため停止 ！
! グローバルデータベースが保存される。停止前のステ...
! !
! アドバイス: ！
! - 境界条件をチェックしてください！
! - モデルと単位の整合性をチェックする！
! - 接触を使用する場合、接触だけで構造を保持して...
- 並行して、MUMPSの特異点検出基準が悲観的すぎることがあ...
この基準を緩和する（NPRECキーワードの値を1または2増加...
切断する（NPRECキーワードの値を1）か、より少ないプロセ...
接触のみで支柱を固定しようとしたことが原因と思われる。
圧縮側の接触部分に半円の仮想材料を挟んでみる。
**1月24日 [#z01a8ccd]
以前の中間発表で用いたモデルが鞘管と下の鋼板がしっかりと...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
そこでまず次のような正12多角柱の鞘管モデルを作成して解析...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
続いてこのモデルに穴を開けずに鞘管を鋼板に一体化させたモ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
この結果がうまく行っていた為、次のような穴を開けたモデル...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
この結果からメッシュの隙間はないほうがいいと言えると思っ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
一旦中間発表時点で出ていた一つの問題はジオメトリでの作成...
円柱の結果と角柱の結果を比較してみたら
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
青が角柱 オレンジが円柱
圧縮側の鞘管の降伏のタイミングが違うのは円柱と角柱のスリ...
スリットサイズを揃えた鞘管のみの角柱のモデルでエラーが出...
この時出て解決できなかったエラー
! <EXCEPTION> <ADAPTATION_11> ...
! ...
! La valeur du pas de temps retenu < 8.91009488412...
PAS_MINIで設定できる最小値は1e-12までなのでどうすれはいい...
円柱のモデルで降伏のタイミングが早くなってしまっているよ...
**12月5日 [#f5c340cc]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
縦軸を縦軸モーメント、横軸を変位として値をそろえて比較し...
変位を揃えるために、解析結果の方の変位を支柱と鞘管の高さ...
このようなズレが生じた原因としてはまず、私の行った解析結...
次に実験では鞘管に支柱を固定するためのボルトが刺さってい...
そのため実験に近い形のモデルを作成する必要がやはりあると...
**12月4日 [#a4bf5258]
まず、接触面を平行に直す前の実際の実験と同じモデルで同じ...
<S> 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.
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
ピンクの点がデータを取った箇所
荷重863N/mm （200mm[h]×862.8 N/mm ×130mm[線荷重]=22432800...
(過去の実験値 701mm[h]×320 N/mm ×100mm[線荷重]=22432000 ...
縦軸モーメントを基準に比較をしようとしたが、よくよく考え...
**11月26日 [#y9f99f09]
<EXCEPTION> <FACTOR_11> ...
! ...
! Problème : la matrice est singulière ou presque sing...
! Lors de la factorisation de la matrice, on a renco...
! (pivot nul ou presque nul) à la ligne 406 qui corr...
問題：行列が特異またはほとんど特異である： ！
! 行列を因数分解すると、問題が発生した！
! (ピボットがゼロまたはほとんどゼロ)が406行目にあり...
! ...
! Risques et conseils : ...
! * Si la ligne correspond a un degré de liberté ph...
! de corps rigide mal bloqué. ...
! Vérifiez les conditions aux limites. ...
! Si vous faites du contact, il ne faut pas que l...
! Vérifiez également les caractéristiques matéria...
リスクとアドバイス： ！
! * その線が物理的自由度に対応している場合、それは...
の動きである可能性が高い！
! 境界条件をチェックする！
! 接触している場合、その構造は接触によってのみ「...
! また、材料の特性（ヤング率など）もチェックして...
! ...
! * Si la ligne correspond a un degré de liberté de...
! limite redondante. ...
! En particulier, il se peut que la relation liné...
! Peut-être devriez vous exclure certains noeuds ...
! (mots clés SANS_NOEUD et SANS_GROUP_NO). ...
線がラグランジュの自由度に対応する場合、これはおそらく...
冗長限界！
! 特に、その冗長な線形関係は、接触条件から来てい...
! 接触条件からいくつかの節点を除外すべきかもしれ...
! (キーワードWITHOUT_NOEUDとWITHOUT_GROUP_NO)。
! ...
! * Si le solveur utilisé est LDLT ou MULT_FRONT, v...
! car celui-ci est le seul à pouvoir factoriser l...
使用するソルバーが LDLT または MULT_FRONT の場合は、MU...
これは正定値でない行列を因数分解できる唯一のソルバーで...
! ...
! * Parfois, en parallèle, le critère de détection ...
! possible de faire passer le calcul complet en r...
! en le débranchant (valeur du mot-clé NPREC=-1) ...
並行して、MUMPSの特異点検出基準が悲観的すぎることがあ...
この基準を緩和する（NPRECキーワードの値を1または2大き...
切断（NPRECキーワードの値を1）するか、より少ないプロセ...
! ...
! * Il se peut aussi que ce phénomène soit tout à f...
! très près d'un noeud. ...
! Si le nombre de décimales perdues n'est pas tro...
! vous pouvez relancer le calcul en augmentant le...
! mot-clé NPREC du mot clé facteur SOLVEUR. ...
! Sinon, contactez l'équipe de développement.
また、亀裂が節に非常に近いところを通過する場合、X-FEM...
節に非常に近い！
! 失われた小数点以下の桁数がそれほど大きくなけれ...
小数点以下の桁数を増やして計算を再開することができます！
キーワード係数SOLVEURのキーワードNPREC！
! そうでなければ、開発チームに連絡してください。
2種類の鋼材を合わせた隙間のないモデル（岐阜高専さんの使っ...
まず、岐阜高専さんのやり方でやっていた時、接触面に共通の...
一旦今後の方針として後藤先生と相談した結果、摩擦なしの設...
**11月22日 [#a6e3681a]
<NonConvergenceError> <MECANONLINE9_7> ...
! ...
! Arrêt pour cause d'absence de convergence avec le...
! La base globale est sauvegardée. Elle contient le...
! ...
! 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é...
<NonConvergenceError> <MECANONLINE9_7> ！
! !
! ニュートンの非線形アルゴリズムにおいて、必要な反...
! グローバルデータベースが保存される。停止前のステ...
! !
! アドバイス: ！
! - ITER_GLOB_MAXI を増やす！
! - タンジェント行列の更新頻度を上げる！
! - 時間離散化を洗練させる ！
! - DEFI_LIST_INSTコマンドでイベント管理（例えば時...
conotinueを使って摩擦を考慮した解析を回していっているが、...
**11月20日 [#fb8c2e15]
前回の鞘管のモデルをdiscreatを使って解析を回していたがう...
これまでに摩擦を考慮した簡単なモデルをやったことがなかっ...
ただうまくいかず、安部さんによると現在の設定だけでは足り...
<Erreur> Échec dans le traitement du contact discret
<Action> On essaie de découper le pas de temps.
On utilise la découpe manuelle.
Découpe uniforme à partir de l'instant < 6.0000...
(soit un incrément constant de < 1.94...
Le pas de temps minimum < 1.000000000000e-20> (...
<Action><Échec> Échec dans la tentative de découper le p...
!<TraitementContactError> <MECANONLINE9_5> ...
! ...
! Arrêt pour cause d'échec lors du traitement du co...
! La base globale est sauvegardée. Elle contient le...
**11月18日 [#p51bfa2a]
摩擦を考慮するためにdiscreteを使った設定を作成中
! <A> <JEVEUX1_64> ...
! ...
! Le volume des données temporaires (objets de la bas...
! est plus de 100.00 fois supérieur au volume de donn...
! ...
! Risques et conseils : ...
! Ce déséquilibre n'a pas de conséquence sur les résu...
! certaines structures de données temporaires ont été...
! été relues. C'est le cas lorsque vous utilisez le s...
! données sont volontairement déchargées pour maximis...
! ...
! ...
! Ceci est une alarme. Si vous ne comprenez pas le sen...
! alarme, vous pouvez obtenir des résultats inattendus...
! ディスクに書き込まれた一時データ（揮発性データベー...
! !
! リスクとアドバイス： ！
! この不均衡は計算結果には何の影響もなく、単に特定の...
! !
! !
! これはアラームである。この！
アラームの意味を理解していないと、予期せぬ結果を招くこ...
このエラーと接触不良に関するエラーが消えない。
**11月13日 [#wcbb0360]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
contenueを使った設定では上記のような結果を得ることができ...
ただこの設定で摩擦を考慮できているのかがよくわからないた...
中間発表のものはdiscreteでやっていたが摩擦を考慮しない設...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
このグラフは上のparavisの結果を使って右側の鞘管の載荷線部...
**11月5日 [#w7e46ee9]
これまで使ってきた設定をノート等に調べてまとめている。
接触面が複数ある場合はzoneのアイテム数をその分増やす必要...
discreteからcontinueにかえて、それに伴う設定も変えてみて...
! <MatriceSinguliereError> <MECANONLINE9_4> ...
! ...
! Arrêt pour cause de matrice non inversible. ...
! La base globale est sauvegardée. Elle contient le...
! ...
! Conseils : ...
! - Vérifiez vos conditions aux limites. ...
! - Vérifiez votre modèle, la cohérence des unit...
! - Si vous faites du contact, il ne faut pas qu...
! ...
! - Parfois, en parallèle, le critère de détecti...
! possible de faire passer le calcul complet e...
! en le débranchant (valeur du mot-clé NPREC=-...
!非可逆行列による停止 ！
! グローバルデータベースが保存される。停止前のステ...
! !
! アドバイス: ！
! - 境界条件をチェックしてください！
! - モデルと単位の整合性をチェックする！
! - 接触を使用する場合、接触だけで構造を保持して...
! !
! - 並行して、MUMPSの特異点検出基準が悲観的すぎ...
！それでも、しばしばこの基準を緩和する（NPRECキーワード...
！切断する（NPRECキーワードの値を1にする）か、より少な...
荷重のかけすぎ、または特定の点に荷重が集中しすぎてしまっ...
一旦鋼材二枚のモデルにこの設定で試してみる
→小さいモデルだと解析で計算してくれず、大きくしたもので試...
弾塑性を消して解析を回してみたりもしたが結果は変わらなか...
岡田さんのcontinueを使用してうまく行ったhdfファイルを借り...
→DEFI_CONTACTのCONTACT_INITをNONに変更したら解析が回った...
warp by vector が何故か使えないため
https://note.com/mmer547/n/n1a0ce48712e4
これを参考に試してみる
**10月23日 [#x4acd072]
以前のモデルで接触面の設定を変えたものはうまく行ったが、...
! <TraitementContactError> <MECANONLINE9_5> ...
! ...
! Arrêt pour cause d'échec lors du traitement du co...
! La base globale est sauvegardée. Elle contient le...
**10月22日 [#ef116cfa]
前回のモデルの解析は要素数が47000程度のものだと解析に4日...
このエラーの原因がよくわからなかったため、今まで参考にし...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
→DEFI_CONTACTの設定より、
触する面と接触される面の設定でitemの数を増やす。
**10月15日 [#h730cbd5]
円柱状のモデルを作成。サイズやスリットの角度などはCtype4....
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
解析時間は26時間ほど。加えた荷重は線荷重で引張側の鞘管の...
解析結果は「接触不良によるエラー」が出てしまい解析が回ら...
詳しい原因はわからないが、おそらく今までは線が面に接触す...
そのため一旦次のモデルは実際のものとは少し違うが、スリッ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
↑のようなモデルを作成してみたため現在解析を回してみている。
**10月10日 [#a257e150]
大きさを変えても問題がないかの検証のため、200×130×9mmのモ...
線荷重を100N/mm加えた結果、解析自体は回り、時間は26時間程...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
今回は少し荷重が足りず、2つめの鋼材が降伏しきるところまで...
今度時間があるタイミングでもう一度荷重を増やして回してみ...
とりあえず次の段階として実際に用いられているモデルと同じ2...
**10月2日 [#zd02de0c]
以前までのモデルサイズの1/10での接触解析はparavisや荷重−...
まずジオメトリでの作成の際はモデル本体作成後、グループ作...
メッシュの作成時いつもどうりにメッシュサイズを設定後、サ...
その後メッシュを切る。今回は全体を1つとしてメッシュを切っ...
aster-studyでは以下を使用した。
,コマンド名 ,使用用途 ,カテゴリーの名前
,LIRE_MAILLAGE ,メッシュの選択 ,mesh
,AFFE_MODELE ,モデル形状の指定 ,model Definition
,DEFI_MATERIAU ,材料諸量の指定 ,material
,AFFE_MATERIAU ,各モデルの材料諸量の指定 ,material
,DEFI_FONCTION ,時間の指定 ,Functions and Lists
,DEFI_LIST_INST , 時間ステップ幅の自動調整機能の設定 ...
,DEFI_LIST_REEL ,時間ステップの指定 ,Functions and Lists
,AFFE_CHAR_MECA ,境界条件の指定 ,BC and Load
,DEFI_CONTACT ,摩擦なしの接触の指定 ,BC and Load
,STAT_NON_LINE ,非線形解析の設定 ,Analysis
,CALC_CHAMP ,計算の指定 ,Post processing
,IMPR_RESU ,出力の指定 ,output
接触解析メモのページにかかれていたものとほとんど使ってい...
DEFI_LIST_INSTの時間ステップ幅の自動調整機能の設定をしな...
STAT_NON_LINEでNEWTONの設定でMATRICEがTANGENTEだとエラー...
解析時間短縮のためSTAT_NON_LINEのCONVERGENCEのRESI_GLOB_M...
DEFI_LIST_INSTのADAPTATION よりEVENEMENTはSEUILで設定。EC...
DEFI_CONTACTのFORMULATIONはDISCRETEで設定した。
今回は線荷重で荷重を加えて載荷を行った。
**9月9日 [#s9ce2c6d]
おそらく前回一定に降伏してしまった原因は境界条件の設定を...
**9月7日 [#pfb5c5e0]
現在接触解析をうまくできないか試していっているが、岐阜高...
ただ、鋼材としての挙動があきらかにおかしく、一定にずっと...
バージョンの違いも有り、設定を変えていく必要があったり、...
**7月31日 現状の情報整理 [#nfbe50a3]
まず200×130×9の鋼板１枚のモデルを荷重5N/mm^2をかけ解析を...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
縦軸 荷重 横軸 鉛直方向変位
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
応力-ひずみ曲線のグラフと降伏度周辺のデータ （左 ひずみ...
続いて鋼材2枚で仮想材料を挟んだモデルは鋼材の厚さが9mmで...
載荷面の6点をプロットして荷重と鉛直方向の変位をグラフに落...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
このモデルの応力-ひずみ曲線を圧縮側の中央一点を上の鋼材と...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
左 上の鋼材の応力-ひずみ曲線 右 下の鋼材の応力-ひ...
とりあえず2つの鋼材は降伏していることが見受けられる。
ただ、荷重と変位の実際の実験結果を再現できそうではなかっ...
上記図のような載荷面側の片側のみ、計14個の7mmのバネ要素を...
バネ定数を32N/mmのモデルに荷重を10N/mm^2かけてみた結果次...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
上の鋼材に変位は見られたものの、下の鋼材変位がほとんど見...
ばね定数を80N/mmにして荷重を8N/mm^2かけてみた結果
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
ばね定数を56N/mmにして荷重を8N/mm^2かけてみた結果
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
応力歪曲線も確認してみたかったが、EPEQがバネ要素を用いて...
ここまで接触解析を用いずに接触を容易にモデル化するために...
そのために接触解析のための簡単なモデルから試してみる。
**24/7/16 [#u951dcaf]
全体にバネを設定していたのを載荷面側の一列のみに変更。
今まで降伏点を理論値のもので考えてしまっていたため、もう...
→（こうすることで仮想材料を入れていたときのものに変化がな...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
変位に対する剛性（ばね定数？）を32N/mmにし、荷重を10Nに設...
**24/7/10 [#e9d0b372]
バネの変位に対する剛性を30にしたが下の鋼材に接する程の荷...
バネの変位に対する剛性を15にすると鋼材は貫通してしまった。
バネの変位に対する剛性を18で荷重を15加えると
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
下の鋼材が少しだけ圧縮されているもののあまり変位が見られ...
**24/7/9 [#ydd08ccf]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
10×10間隔でバネ要素を断面全体に取り込んだモデルの解析を行...
バネ要素のz軸の剛性を10以上にしないと鋼材が貫通してしまい...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
↑荷重15 バネのｚ軸方向の変位に対する剛性20N/mm
**24/7/8 [#q3c2baf2]
前回2枚の鋼材の間に仮想材料の代わりにバネ要素を使ったモデ...
予想はしていたが、バネ一つだけでは上の鋼材が降伏したあと...
この解析を回す時バネ要素が鋼材と接している部分がメッシュ...
**24/7/5 [#cb4cf7c1]
今まで荷重が10までしか回っていなかったが仮想材料がヤング...
どうしても部材が同時に降伏してしまっているようで今まで以...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
**24/7/4 [#a3ca8113]
仮想材料のヤング率を0.01から0.1にしてみたが同様にすり抜け...
ヤング率を0.5、0.7にし、載荷を10にして解析を回したが、す...
***今後の方針 [#r1137367]
直方体の仮想材料を入れるだけではうまくいかなかったので次...
**24/7/3 [#v149df48]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
縦軸：荷重 横軸：変位
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
ポアソン比が0に近い場合、仮に縦に荷重を加えても横に膨らん...
仮想材料のポアソン比を0.001にし、解析を回したら、想定通り...
仮想材料をもっと柔らかいものにするために、仮想材料のヤン...
***現在の目標の再確認と問題点 [#j003f185]
現在は間の空いた2枚の部材の片側を降伏させて、荷重の様子を...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
**24/7/1,2 [#kb54db15]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
縦軸：荷重 横軸：変位
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
設定荷重が3.526Nだとあまり降伏が見られなかったため、荷重...
タイムステップの数も10分割から20分割に変更。
結果的に降伏は確認することができたものの１回分の降伏しか...
**24/6/28 [#qd4de2b5]
たわみの理論値計算がミスっていた
訂正 33.801→3.3801
仮想材料の大きさが大きく変わるため結果が変わるはず。
モデル図 鋼材 200×130×9 (上部、下部の部材）
仮想材料 200×130×3.3801 （中間の部材）
ヤング率 ポアソン比
鋼材 235 0.3
仮想材料 1 0.3
設定荷重 3.526N（面載荷 上部鋼材の断面のみ）
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
**24/6/27 [#cde629c2]
○前回のモデルの変更点
・線荷重から載荷側の上の鋼材の断面による面荷重に
・設定荷重を増やした
・理論値をしっかり仮想材料も考慮したものへ
設定荷重を増やして解析を回したが、設定荷重を理論値の3倍以...
Arrêt suite à l'échec de l'intégration de la loi de comp...
! La base globale est sauvegardée. Elle contient le...
! ...
! Conseils : ...
! - Vérifiez vos paramètres, la cohérence des un...
! - Essayez d'augmenter ITER_INTE_MAXI, ou de su...
このエラーによりITER_INTE_MAXIを大きくすればいいというこ...
設定荷重を2倍にしたときは解析が回った。しかし下の鋼材の外...
**24/6/~25 [#nf21dd51]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
スケールファクター 20倍
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
スケールファクター 25倍
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
スケールファクター 30倍
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
まず簡単なモデルで2段階の降伏が見られるグラフを作成できる...
・鋼材 200（高さ）×130（幅）×9（厚さ）
仮想材料 200×130×34.232（理論値）
ヤング率 ポアソン比
鋼材 20600 0.3
仮想材料 1 0.3
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
しかし今回出たグラフは直線になってしまった。（縦軸 荷重...
一見Paravisのスケールファクターが20倍のものをみたら良さげ...
荷重方法が今回は線荷重にして、荷重部分のグラフを作成した...
あと仮想材料をParavisで表示させない方法も調べておきたい。...
**24/6/~18 [#dd9728dd]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
左のグラフが過去の実験結果を示したグラフで右が解析結果に...
仮想材料を用いても降伏が一回しか見られないという小川さん...
**24/6/~11 [#hfd016ca]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
前回までのモデルで相当応力−相当ひずみ曲線を作成した。上の...
前回の応力−ひずみ曲線はマイナスの値を取っていたり、与えて...
Paravisで相当ひずみを確認するためにAsterStudyで設定を付け...
相当応力−相当ひずみ曲線を作成した際、スライスした箇所は固...
卒論に関しては、以前までは「Atype」と「Btype」と書かれて...
当然であったが15分程度では解析が回らず、設定時間が足りな...
**24/5/30 [#d4d56b6d]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
上記の写真が今回使用した片持ちばりの載荷されている様子の...
応力がマイナスになっているのは圧縮、プラスになっているの...
**24/5/22~28 [#qc76d78d]
先週作っていたモデルをecro_lineを用いた設定の仕方でも作成...
弾塑性解析で作ったモデルをParavisで設定する方法等を練習。...
小川さんのモデルをParavisで見る時、どのように応用できるか...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
これは載荷面の右上端の点のものだが、上記のようなグラフに...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
**24/5/16～21 [#wfdc6df1]
・やったこと
小川さんの卒論のスライドや概要の内容を読んだ。そして卒論...
解析は回ったが、弾塑性の内容についてよくわかっていなかっ...
salome-meca 弾塑性解析（2021)
https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?cmd=re...
**24/5/15(水） [#i0e06d43]
C種について
https://www.mlit.go.jp/notice/noticedata/sgml/070/7900013...
*春課題について [#jf096965]
**24/5/1（水） [#o0b3c7a2]
メモ
・10*5*700 単純1次
鉛直1次 理論値 37.15451
,メッシュ, 要素数, 解析値, 相対誤差
,0.7,604706, 37.8469, 1.864
,1,150957, 39.3168, 5.820
,2,21270, 47.7953, 28.639
,4,2609, 75.567, 103.386
鉛直2次 理論値 120.3942
,メッシュ, 要素数, 解析値, 相対誤差
,0.7, 604706, 122.609, 1.840
,1, 150957, 127.365, 5.790
,2, 21270, 154.866, 28.632
,4, 2609, 244.558, 103.131
鉛直3次 理論値 251.1547
,メッシュ, 要素数, 解析値, 相対誤差
,0.7, 604706, 255.685, -1.804
,1, 150957, 265.629, -5.763
,2, 21270, 322.8, 28.526
,4, 2609, 509.496, 102.861
水平1次 理論値 74.30902
,メッシュ, 要素数, 解析値, 相対誤差
,0.7, 604706, 74.1075, 0.271
,1, 150957, 74.9296, 0.835
,2, 21270, 77.9669, 4.923
,4, 2609, 87.7007, 18.022
水平2次 理論値 240.7885
,メッシュ, 要素数, 解析値, 相対誤差
,0.7, 604706, 239.837, 0.395
,1, 150957, 242.486, 0.705
,2, 21270, 252.306, 4.783
,4, 2609, 283.626, 17.791
水平3次 理論値 502.3093
,メッシュ, 要素数, 解析値, 相対誤差
,0.7, 604706, 499.386, 0.582
,1, 150957, 504.892, 0.514
,2, 21270, 525.213, 4.560
,4, 2609, 590.132, 17.484
**24/4/16(火) [#ieee64e9]
共振周波数解析について
https://www.muratasoftware.com/fwp/wp-content/uploads/201...
振動設計について
https://www.fml.t.u-tokyo.ac.jp/lecture/handout/DE/yurenu...
**24/4/15(月） [#ta19d4ae]
メモ
・Texで図を挿入する際に一段に2つ以上横に並べて図を貼る方法
\begin{figure}[h]
\begin{tabular}{cc}
%---- 最初の図 ---------------------------
\begin{minipage}{0.45\hsize}
\begin{center}
\includegraphics[width=サイズ]{ファイル1}
\caption{キャプション1}
\label{ラベル1}
\end{center}
\end{minipage} &
%---- 2番目の図 --------------------------
\begin{minipage}{0.45\hsize}
\centering
\includegraphics[width=サイズ]{ファイル2}
\caption{キャプション2}
\label{ラベル2}
\end{center}
\end{minipage}
%---- 図はここまで ----------------------
\end{tabular}
\end{figure}
図がかぶってしまうときは[width]のサイズを変えたりする。
今回は2つを横に並べただけで試していないがおそらく\begin{t...
参考文献
http://www.yamamo10.jp/~yamamoto/comp/latex/make_doc/inse...
・texで表を作成した際に全体のサイズを変更する方法
→\scalebox で拡大・縮小する
\begin{table}[htbp]
\begin{center}
\caption{表題}
\label{heading}
\scalebox{0.5}[0.9]{ %ココ
\begin{tabular}{|c|c|c|c|} \hline
& & & \\\hline
& & & \\\hline
\end{tabular}
}
\end{center}
\end{table}
通常のグラフを作成して\label{heading}の後に\scalebox{0.5}...
\scaleboxの後ろの｛｝で横倍率、[]で縦倍率が変更可能。ただ...
参考文献
https://blog.goo.ne.jp/shafact/e/8219383eb299d5b6bdfd1544...
・texでの矩形選択とコピペ方法
Escを押した後、Ctrl＋vで矩形選択モードになり矢印キーで選...
選択したい箇所が決まったら、ｙでコピー。
その後、ｐで貼り付けるが、この時下の列に注意。貼り付ける...
表を作成する際、メモを開いて
%\begin{figure}[h] ←これはいらない
%1 & \\
%2 & \\
%3 & \\
%4 & \\
%5 & \\
%6 & \\
%7 & \\
%8 & \\
%9 & \\
%1 & \\
的なものを用意してtexに作成する表のしたに貼り付けておくと...
他にも表を多用する際は一度メモに使用するコマンドを貼り付...
*創造工房実習 (3年次）[#o6991381]
**23/12/14(金) [#ed1d2b72]
texの基本的な書き方https://www.str.ce.akita-u.ac.jp/~goto...
**23/12/08(金) [#o5aea6ef]
,メッシュの長さ ,要素数 ,先端変位[mm] ,相対誤差 ,計算者
,0.7 ,155192 ,0.08378905246 ,15.365 ,安藤
,0.8,138808,0.08380386491,15.350,安藤
,0.9,82587,0.083707073981,15.45 ,兼田
,1.1,38671 ,0.084201207602,14.95 ,兼田
,1.2,31929 ,0.083688 ,15.466. ,柴田
,1.3,28621 ,0.083669 ,15.4857 ,柴田
,1.4,28854 ,0.08368 ,15.47 ,佐藤
,1.5,20015 ,0.084052 ,15.10 ,佐藤
,1.6,19448 ,0.0835402938 ,15.62 ,皆川
,1.7,13801 ,0.0834355098 ,15.72 ,皆川
,1.8,12528 ,0.083733 ,15.42 ,永山
,1.9 ,11769 ,0.083924 ,15.23 ,永山
,2 ,10699,0.084076876559,15.074 ,辻
,3,3579,0.08414561753,15.004 ,辻
,4,1628 ,0.082794,16.37 ,服部
,5,1016 ,0.083033,18.89,服部
,6,839 ,-0.082882 ,16.26 ,梶原
,7,554 ,-0.080871 ,18.28 ,梶原
,8,285,0.079995,-19.20 ,工藤
,9,261 ,0.078980,-20.22 ,工藤
,10,232 ,0.081911,17.26 ,佐々木
,11,208 ,0.075676,23.56,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
**23/11/29(金） [#m5e37185]
異方性１次について
自分のデータ
,メッシュ長さ,要素数,変位,相対誤差,計算者
,4,1500,0.430156,12.52,服部
,5,432,0.282968,42.45,服部
全員のデータ
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.7,144563,0.505252,2.76,安藤
,0.8,141517,0.504692,2.64,安藤
,0.9,91648,0.502595,2.216,兼田
,1.1,27160,0.489914,0.363,兼田
,1.2,24675,0.487088,0.791,柴田
,1.3,23446,0.4868010,0.995,柴田
,1.4,17738,0.485999,1.16,佐藤
,1.5,15438,0.485180,1.33,佐藤
,1.6,15900,0.483286,1.71,皆川
,1.7,12142,0.477952,2.80,皆川
,1.8,11604,0.482085,1.9554,永山
,1.9,10391,0.470887,4.2329,永山
,2,10291,0.480910,2.19,辻
,3,2328,0.431937,12.15,辻
,4,1500,0.430156,12.52,服部
,5,432,0.282968,42.45,服部
,6,356,0.3441556,30.00,梶原
,7,196,0.213934,56.49,梶原
,8,104,0.229874,53.25,工藤
,9,81,0.232308,52.75,工藤
,10,78,0.203271,58.65,佐々木
,11,63,0.222316,54.78,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
等方性２次について
自分のデータ
,メッシュ長さ,要素数,変位,相対誤差,計算者
,4,1500,0.429254,3.01,服部
,5,432,0.428170,2.75,服部
全員のデータ
,メッシュ長さ ,要素数 ,先端変位 ,相対誤差 ,計算者
,0.7,144563,0.430124,3.22,安藤
,0.8,141517,0.430132,3.22,安藤
,0.9,91648,0.430020,3.197,兼田
,1.1,27160,0.429828,3.151,兼田
,1.2,24675,0.429836,3.15,柴田
,1.3,23446,0.42974,3.13 ,柴田
,1.4,17738,0.429797,1.3 ,佐藤
,1.5,15438,0.429958,3.14,佐藤
,1.6,15900,0.429755,3.18,皆川
,1.7,12142,0.429676,3.11,皆川
,1.8,11604,0.429829,3.1507,永山
,1.9,10391,0.429684,3.1159,永山
,2,10291,0.429620,3.10,辻
,3,2328,0.429169,2.99,辻
,4,1500,0.429254,3.01,服部
,5,432,0.428170,2.75,服部
,6 ,356 ,0.428452 ,2.82 ,梶原
,7 ,196 ,0.42591 ,2.21,梶原
,8 ,104 ,0.426074 ,2.25 ,工藤
,9 ,81 ,0.425552 ,2.12 ,工藤
,10 ,78 ,0.488382 ,17.20, 佐々木
,11 ,63 ,0.423972 ,9.0534 ,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
**23/11/24(金） [#dbef5dfe]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.7,1455234,0.422484,1.388,安藤
,0.8,142973,0.422570,1.409,安藤
,0.9,91648,-0.420437,0.897,兼田
,1.1,27160,-0.405618,2.659,兼田
,1.2,24675,0.404349,2.96,柴田
,1.3,23446,0.404185,3.00,柴田
,1.4,17738,0.398604,4.34,佐藤
,1.5,15438,0.396593,4.83,佐藤
,1.6,16122,0.398212,4.44,皆川
,1.7,12026,0.393411,5.59,皆川
,1.8,11604,0.393668,5.53,永山
,1.9,10391,0.390695,6.24,永山
,2,10921,-0.395103,5.18,辻
,3,2328,-0.324762,22.06,辻
,4,1500,-0.155013,62.80,服部
,5,432,-0.065278,84.33,服部
,6,357,0.213062,48.87,梶原
,7,196,0.1019,75.55,梶原
,8,104,0.1158624,72.20,工藤
,9,81,0.1255118,69.88,工藤
,10,78,-0.07733,81.44,佐々木
,11,63,-0.1999,52.03,佐々木
**23/11/17(金) [#vcf762ab]
すいません、僕のデータの相対誤差が100倍し忘れていたので直...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
,メッシュ長さ,要素数,先端変位(4隅の平均値)[mm],相対誤差,...
,0.7,198464,6.54281,-1.91,安藤
,0.8,113812,6.5104,-2.39,安藤
,0.9,40280,6.3631525,-4.60,兼田
,1.1,30055,6.3363525,-5.00,兼田
,1.2,264667,6.3043375,5.48,柴田
,1.3,25180,6.304355,5.48,柴田
,1.4,32212,6.31612,5.31,佐藤
,1.5,17753,6.1209,8.23,佐藤
,1.6,14296,6.2044625,-6.98,皆川
,1.7,13596,6.2156625,-6.81,皆川
,1.8,2866,5.737755,-13.98,永山
,1.9,6001,5.7263625,-14.15,永山
,2,5617,5.6458525,-15.355,辻
,3,2309,5.4728755,-17.948,辻
,4,617,3.6160575,45.12,服部
,5,494,3.8580375,42.16,服部
,6,581,2.50682,-62.416,梶原
,7,133,1.41225,-78.82,梶原
,8,78,1.2887175,-80.68,工藤
,9,72,1.2879925,-80.69,工藤
,10,60,1.14344,-82.85,佐々木
,11,65,1.23124,-81.154,佐々木
**23/11/10(金) [#c94eac4c]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
**23/10/27(金) コマンドと操作確認 [#t9ee68ee]
pwd : 自分の位置
ls : すべてのファイルを確認する
mkdir : 新規ファイルの作成
ls (作成できているかの確認）
cd ファイル名 ：複数のファイルから特定のファイルに入る
pwd (現在位置の確認）
ls -a : 隠しファイル
gedit ファイル名.tex & : ファイルを作りそのファイルに入...
vi ファイル名.tex : 編集（作業スペースからファイル内容を）
Esc → : → wq : 内容保存
cp コピーしたいファイル名.tex 新しいファイル名.tex : 新...
rm ファイル名.tex : ファイルを消す
cd .. : 1つ前のフォルダに戻る（他ファイルの選択が可能に）
rmdir ファイル名 ：mkdirで作ったものの削除
マウスの両クリックor中央クリック ： 選択した文章コピーa...
*salome-meca 説明 [#z8702cc2]
http://dugi6514.odns.fr/doc/v15/en/index.php?man=commande
**salome-meca 設定に関する情報等 [#e412098d]
https://www.slideshare.net/slideshow/salomemeca-59878666/...
https://qiita.com/Jun_Tatsuno/items/3e6558c1adc4dd9499fe
PARM_THETA によって，θ の値を変えられる. 時間と空間の再分...
RESI_INTE_RELA, ITER_INTE_MAXIについて
「これらは線形研究のパラメータである。最大反復回数は、...
と、線形研究の収束のために到達すべき精度reslin が与えられ...
を使用しないことを推奨する。
CORD法では、精度や反復回数をあまり高く指定する必要はあり...
実践によると、線形研究の反復は 2、3 回で十分です。したが...
デフォルトの精度で3回の反復を求めることで満足できる。ユー...
を超える線形調査の反復を行うことはできません。
一方、MIXED法では、損傷のある問題では、数十回の反復が有効...
が有効です。」 Deepl翻訳より
サンドイッチ梁について
https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome...
AsterStudyでよく使うコマンド
,コマンド名 ,使用用途 ,カテゴリーの名前
,LIRE_MAILLAGE ,メッシュの選択 ,mesh
,AFFE_MODELE ,モデル形状の指定 ,model Definition
,DEFI_MATERIAU ,材料諸量の指定 ,material
,AFFE_MATERIAU ,各モデルの材料諸量の指定 ,material
,DEFI_FONCTION ,時間の指定 ,Functions an...
,DEFI_LIST_REEL ,時間ステップの指定 ,Functions a...
,AFFE_CHAR_MECA ,境界条件の指定 ,BC and ...
,DEFI_CONTACT ,摩擦なしの接触の指定 ,BC and Load
,STAT_NON_LINE ,非線形解析の設定 ,Analysis
,CALC_CHAMP ,計算の指定 ,Post processing
,IMPR_RESU ,出力の指定 ,output
**salome-mecaエラーについて [#sd314b3c]
Arrêt suite à l'échec de l'intégration de la loi de com...
! La base globale est sauvegardée. Elle contient le...
! ...
! Conseils : ...
! - Vérifiez vos paramètres, la cohérence des un...
! - Essayez d'augmenter ITER_INTE_MAXI, ou de su...
ITER_INTE_MAXIの値を大きくしてみる。（ただ大きくすればい...
ITER_INTE_PASを考慮するでも解決するらしい（未検討 調べた...
! <NonConvergenceError> <MECANONLINE9_7> ...
! ...
! Arrêt pour cause d'absence de convergence avec le...
! La base globale est sauvegardée. Elle contient le...
! ...
! 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é...
! ニュートンの非線形アルゴリズムにおいて、必要な反...
! グローバルデータベースが保存される。停止前のステ...
! !
! アドバイス: ！
! - ITER_GLOB_MAXI を増やす！
! - タンジェント行列の更新頻度を上げる！
! - 時間離散化を洗練させる ！
! - DEFI_LIST_INSTコマンドでイベント管理（例えば時...
! <EXCEPTION> <FACTOR_11> ...
! ...
! Problème : la matrice est singulière ou presque sing...
! Lors de la factorisation de la matrice, on a renco...
! (pivot nul ou presque nul) à la ligne 12414 qui co...
! ...
! Risques et conseils : ...
! * Si la ligne correspond a un degré de liberté ph...
! de corps rigide mal bloqué. ...
! Vérifiez les conditions aux limites. ...
! Si vous faites du contact, il ne faut pas que l...
! Vérifiez également les caractéristiques matéria...
! ...
! * Si la ligne correspond a un degré de liberté de...
! limite redondante. ...
! En particulier, il se peut que la relation liné...
! Peut-être devriez vous exclure certains noeuds ...
! (mots clés SANS_NOEUD et SANS_GROUP_NO). ...
! ...
! * Si le solveur utilisé est LDLT ou MULT_FRONT, v...
! car celui-ci est le seul à pouvoir factoriser l...
! ...
! * Parfois, en parallèle, le critère de détection ...
! possible de faire passer le calcul complet en r...
! en le débranchant (valeur du mot-clé NPREC=-1) ...
! ...
! * Il se peut aussi que ce phénomène soit tout à f...
! très près d'un noeud. ...
! Si le nombre de décimales perdues n'est pas tro...
! vous pouvez relancer le calcul en augmentant le...
! mot-clé NPREC du mot clé facteur SOLVEUR. ...
! Sinon, contactez l'équipe de développement. ...
**メモ [#u32aa2f5]
・サロメの過去の人のデータをコピーする際はRun_Caseも一緒...
・最大主応力と最小主応力について
物体に力が作用するときに方向によりせん断力がゼロになる面...
ParavisではPRIN＿１が最小主応力 PRIN＿３が最大主応力
最大主応力は引張応力、最小主応力は圧縮応力となる。
尚、最大主応力の結果表示でプラスの値を示す場合は、引張応力
マイナスを示す場合は、圧縮応力
最小主応力の結果表示でプラスの値を示す場合は、引張応力
マイナスを示す場合は、圧縮応力
↑実際にParavisで確認しながらやると理解しやすい
・最小主ひずみ、最大主ひずみについて
最大主ひずみは材料にかかる最大の引張りひずみの評価、最小...
・appariementについて
ジオメトリ上の固定点ループのパラメータ（デフォルト・パラ...
§ 3.4.1 .
ALGO_RESO_ GEOM = 'NEWTON'を選択すると、一般化ニュートン...
幾何学的非線形性は、各ニュートン反復で考慮され、マッチン...
(コンタクト・ニュートン(CONTACT NEWTON GENE CRIT. GEOM.)...
Les paramètres de la boucle de point fixe sur la géométri...
§ 3.4.1 .
Quand on choisit ALGO_RESO_ GEOM ='NEWTON', on bascule ve...
La non-linéarité géométrique est alors prise en compte à ...
géométrique (colonne CONTACT NEWTON GENE CRIT. GEOM.).
・ペナルティ要素（法）について
https://www.cybernet.co.jp/ansys/learning/glossary/penalt...
・岐阜高専さんのページ
http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-M...
・接触と摩擦について
file:///home/kouzou/%E3%83%80%E3%82%A6%E3%83%B3%E3%83%AD%...
・CLA
https://www.sciencedirect.com/science/article/pii/S004578...
終了行:
#contents
*卒論テーマ [#id2b2754]
CLT床版用の防護柵の開発（安部さんサポート）旭川の実験モデ...
2つに割れている鞘管が1つずつ順番に降伏することで、どれぐ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
↑防護柵 図面
・小川さんの論文 CLT床版への取付・交換を考慮した鋼製防護...
https://doi.org/10.11532/structcivil.70A.885
・小川さん 卒論：CLT床版に取り付けられた鋼製防護柵の性能...
・小川さん 修論：CLT 床版への取付・交換が容易な鋼製防護柵...
・安部さん：支柱交換を考慮したCLT床版用鋼製防護柵の挙動（...
*4年次 進捗[#o6991381]
**2月5日 [#k16c1e69]
! <BoucleGeometrieError> <MECANONLINE9_9> ...
! ...
! Arrêt par échec dans la boucle de point fixe sur ...
! La base globale est sauvegardée. Elle contient le...
<ループ・ジオメトリー・エラー> <MECANONLINE9_9> ！
! !
! ジオメトリの固定点ループの失敗で停止 ！
! グローバルデータベースが保存される。停止前のステ...
**1月26日 [#scf4a32b]
角柱のモデルを正12角形から正24角形に変更して解析を回して...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
前回の正12角柱と比較したものが
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
スリットサイズが正12角柱が30mmで正24角柱が24mmになってい...
この正24角柱のモデルに支柱をかぶせて拘束荷重を与えたモデ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
ただこれはエラーが出てしまった。
! <MatriceSinguliereError> <MECANONLINE9_4> ...
! ...
! Arrêt pour cause de matrice non inversible. ...
! La base globale est sauvegardée. Elle contient le...
! ...
! Conseils : ...
! - Vérifiez vos conditions aux limites. ...
! - Vérifiez votre modèle, la cohérence des unit...
! - Si vous faites du contact, il ne faut pas qu...
! ...
! - Parfois, en parallèle, le critère de détecti...
! possible de faire passer le calcul complet e...
! en le débranchant (valeur du mot-clé NPREC=-...
<MatrixSinguliereError> <MECANONLINE9_4> ！
! !
! 非可逆行列のため停止 ！
! グローバルデータベースが保存される。停止前のステ...
! !
! アドバイス: ！
! - 境界条件をチェックしてください！
! - モデルと単位の整合性をチェックする！
! - 接触を使用する場合、接触だけで構造を保持して...
- 並行して、MUMPSの特異点検出基準が悲観的すぎることがあ...
この基準を緩和する（NPRECキーワードの値を1または2増加...
切断する（NPRECキーワードの値を1）か、より少ないプロセ...
接触のみで支柱を固定しようとしたことが原因と思われる。
圧縮側の接触部分に半円の仮想材料を挟んでみる。
**1月24日 [#z01a8ccd]
以前の中間発表で用いたモデルが鞘管と下の鋼板がしっかりと...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
そこでまず次のような正12多角柱の鞘管モデルを作成して解析...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
続いてこのモデルに穴を開けずに鞘管を鋼板に一体化させたモ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
この結果がうまく行っていた為、次のような穴を開けたモデル...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
この結果からメッシュの隙間はないほうがいいと言えると思っ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
一旦中間発表時点で出ていた一つの問題はジオメトリでの作成...
円柱の結果と角柱の結果を比較してみたら
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
青が角柱 オレンジが円柱
圧縮側の鞘管の降伏のタイミングが違うのは円柱と角柱のスリ...
スリットサイズを揃えた鞘管のみの角柱のモデルでエラーが出...
この時出て解決できなかったエラー
! <EXCEPTION> <ADAPTATION_11> ...
! ...
! La valeur du pas de temps retenu < 8.91009488412...
PAS_MINIで設定できる最小値は1e-12までなのでどうすれはいい...
円柱のモデルで降伏のタイミングが早くなってしまっているよ...
**12月5日 [#f5c340cc]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
縦軸を縦軸モーメント、横軸を変位として値をそろえて比較し...
変位を揃えるために、解析結果の方の変位を支柱と鞘管の高さ...
このようなズレが生じた原因としてはまず、私の行った解析結...
次に実験では鞘管に支柱を固定するためのボルトが刺さってい...
そのため実験に近い形のモデルを作成する必要がやはりあると...
**12月4日 [#a4bf5258]
まず、接触面を平行に直す前の実際の実験と同じモデルで同じ...
<S> 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.
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
ピンクの点がデータを取った箇所
荷重863N/mm （200mm[h]×862.8 N/mm ×130mm[線荷重]=22432800...
(過去の実験値 701mm[h]×320 N/mm ×100mm[線荷重]=22432000 ...
縦軸モーメントを基準に比較をしようとしたが、よくよく考え...
**11月26日 [#y9f99f09]
<EXCEPTION> <FACTOR_11> ...
! ...
! Problème : la matrice est singulière ou presque sing...
! Lors de la factorisation de la matrice, on a renco...
! (pivot nul ou presque nul) à la ligne 406 qui corr...
問題：行列が特異またはほとんど特異である： ！
! 行列を因数分解すると、問題が発生した！
! (ピボットがゼロまたはほとんどゼロ)が406行目にあり...
! ...
! Risques et conseils : ...
! * Si la ligne correspond a un degré de liberté ph...
! de corps rigide mal bloqué. ...
! Vérifiez les conditions aux limites. ...
! Si vous faites du contact, il ne faut pas que l...
! Vérifiez également les caractéristiques matéria...
リスクとアドバイス： ！
! * その線が物理的自由度に対応している場合、それは...
の動きである可能性が高い！
! 境界条件をチェックする！
! 接触している場合、その構造は接触によってのみ「...
! また、材料の特性（ヤング率など）もチェックして...
! ...
! * Si la ligne correspond a un degré de liberté de...
! limite redondante. ...
! En particulier, il se peut que la relation liné...
! Peut-être devriez vous exclure certains noeuds ...
! (mots clés SANS_NOEUD et SANS_GROUP_NO). ...
線がラグランジュの自由度に対応する場合、これはおそらく...
冗長限界！
! 特に、その冗長な線形関係は、接触条件から来てい...
! 接触条件からいくつかの節点を除外すべきかもしれ...
! (キーワードWITHOUT_NOEUDとWITHOUT_GROUP_NO)。
! ...
! * Si le solveur utilisé est LDLT ou MULT_FRONT, v...
! car celui-ci est le seul à pouvoir factoriser l...
使用するソルバーが LDLT または MULT_FRONT の場合は、MU...
これは正定値でない行列を因数分解できる唯一のソルバーで...
! ...
! * Parfois, en parallèle, le critère de détection ...
! possible de faire passer le calcul complet en r...
! en le débranchant (valeur du mot-clé NPREC=-1) ...
並行して、MUMPSの特異点検出基準が悲観的すぎることがあ...
この基準を緩和する（NPRECキーワードの値を1または2大き...
切断（NPRECキーワードの値を1）するか、より少ないプロセ...
! ...
! * Il se peut aussi que ce phénomène soit tout à f...
! très près d'un noeud. ...
! Si le nombre de décimales perdues n'est pas tro...
! vous pouvez relancer le calcul en augmentant le...
! mot-clé NPREC du mot clé facteur SOLVEUR. ...
! Sinon, contactez l'équipe de développement.
また、亀裂が節に非常に近いところを通過する場合、X-FEM...
節に非常に近い！
! 失われた小数点以下の桁数がそれほど大きくなけれ...
小数点以下の桁数を増やして計算を再開することができます！
キーワード係数SOLVEURのキーワードNPREC！
! そうでなければ、開発チームに連絡してください。
2種類の鋼材を合わせた隙間のないモデル（岐阜高専さんの使っ...
まず、岐阜高専さんのやり方でやっていた時、接触面に共通の...
一旦今後の方針として後藤先生と相談した結果、摩擦なしの設...
**11月22日 [#a6e3681a]
<NonConvergenceError> <MECANONLINE9_7> ...
! ...
! Arrêt pour cause d'absence de convergence avec le...
! La base globale est sauvegardée. Elle contient le...
! ...
! 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é...
<NonConvergenceError> <MECANONLINE9_7> ！
! !
! ニュートンの非線形アルゴリズムにおいて、必要な反...
! グローバルデータベースが保存される。停止前のステ...
! !
! アドバイス: ！
! - ITER_GLOB_MAXI を増やす！
! - タンジェント行列の更新頻度を上げる！
! - 時間離散化を洗練させる ！
! - DEFI_LIST_INSTコマンドでイベント管理（例えば時...
conotinueを使って摩擦を考慮した解析を回していっているが、...
**11月20日 [#fb8c2e15]
前回の鞘管のモデルをdiscreatを使って解析を回していたがう...
これまでに摩擦を考慮した簡単なモデルをやったことがなかっ...
ただうまくいかず、安部さんによると現在の設定だけでは足り...
<Erreur> Échec dans le traitement du contact discret
<Action> On essaie de découper le pas de temps.
On utilise la découpe manuelle.
Découpe uniforme à partir de l'instant < 6.0000...
(soit un incrément constant de < 1.94...
Le pas de temps minimum < 1.000000000000e-20> (...
<Action><Échec> Échec dans la tentative de découper le p...
!<TraitementContactError> <MECANONLINE9_5> ...
! ...
! Arrêt pour cause d'échec lors du traitement du co...
! La base globale est sauvegardée. Elle contient le...
**11月18日 [#p51bfa2a]
摩擦を考慮するためにdiscreteを使った設定を作成中
! <A> <JEVEUX1_64> ...
! ...
! Le volume des données temporaires (objets de la bas...
! est plus de 100.00 fois supérieur au volume de donn...
! ...
! Risques et conseils : ...
! Ce déséquilibre n'a pas de conséquence sur les résu...
! certaines structures de données temporaires ont été...
! été relues. C'est le cas lorsque vous utilisez le s...
! données sont volontairement déchargées pour maximis...
! ...
! ...
! Ceci est une alarme. Si vous ne comprenez pas le sen...
! alarme, vous pouvez obtenir des résultats inattendus...
! ディスクに書き込まれた一時データ（揮発性データベー...
! !
! リスクとアドバイス： ！
! この不均衡は計算結果には何の影響もなく、単に特定の...
! !
! !
! これはアラームである。この！
アラームの意味を理解していないと、予期せぬ結果を招くこ...
このエラーと接触不良に関するエラーが消えない。
**11月13日 [#wcbb0360]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
contenueを使った設定では上記のような結果を得ることができ...
ただこの設定で摩擦を考慮できているのかがよくわからないた...
中間発表のものはdiscreteでやっていたが摩擦を考慮しない設...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
このグラフは上のparavisの結果を使って右側の鞘管の載荷線部...
**11月5日 [#w7e46ee9]
これまで使ってきた設定をノート等に調べてまとめている。
接触面が複数ある場合はzoneのアイテム数をその分増やす必要...
discreteからcontinueにかえて、それに伴う設定も変えてみて...
! <MatriceSinguliereError> <MECANONLINE9_4> ...
! ...
! Arrêt pour cause de matrice non inversible. ...
! La base globale est sauvegardée. Elle contient le...
! ...
! Conseils : ...
! - Vérifiez vos conditions aux limites. ...
! - Vérifiez votre modèle, la cohérence des unit...
! - Si vous faites du contact, il ne faut pas qu...
! ...
! - Parfois, en parallèle, le critère de détecti...
! possible de faire passer le calcul complet e...
! en le débranchant (valeur du mot-clé NPREC=-...
!非可逆行列による停止 ！
! グローバルデータベースが保存される。停止前のステ...
! !
! アドバイス: ！
! - 境界条件をチェックしてください！
! - モデルと単位の整合性をチェックする！
! - 接触を使用する場合、接触だけで構造を保持して...
! !
! - 並行して、MUMPSの特異点検出基準が悲観的すぎ...
！それでも、しばしばこの基準を緩和する（NPRECキーワード...
！切断する（NPRECキーワードの値を1にする）か、より少な...
荷重のかけすぎ、または特定の点に荷重が集中しすぎてしまっ...
一旦鋼材二枚のモデルにこの設定で試してみる
→小さいモデルだと解析で計算してくれず、大きくしたもので試...
弾塑性を消して解析を回してみたりもしたが結果は変わらなか...
岡田さんのcontinueを使用してうまく行ったhdfファイルを借り...
→DEFI_CONTACTのCONTACT_INITをNONに変更したら解析が回った...
warp by vector が何故か使えないため
https://note.com/mmer547/n/n1a0ce48712e4
これを参考に試してみる
**10月23日 [#x4acd072]
以前のモデルで接触面の設定を変えたものはうまく行ったが、...
! <TraitementContactError> <MECANONLINE9_5> ...
! ...
! Arrêt pour cause d'échec lors du traitement du co...
! La base globale est sauvegardée. Elle contient le...
**10月22日 [#ef116cfa]
前回のモデルの解析は要素数が47000程度のものだと解析に4日...
このエラーの原因がよくわからなかったため、今まで参考にし...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
→DEFI_CONTACTの設定より、
触する面と接触される面の設定でitemの数を増やす。
**10月15日 [#h730cbd5]
円柱状のモデルを作成。サイズやスリットの角度などはCtype4....
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
解析時間は26時間ほど。加えた荷重は線荷重で引張側の鞘管の...
解析結果は「接触不良によるエラー」が出てしまい解析が回ら...
詳しい原因はわからないが、おそらく今までは線が面に接触す...
そのため一旦次のモデルは実際のものとは少し違うが、スリッ...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
↑のようなモデルを作成してみたため現在解析を回してみている。
**10月10日 [#a257e150]
大きさを変えても問題がないかの検証のため、200×130×9mmのモ...
線荷重を100N/mm加えた結果、解析自体は回り、時間は26時間程...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
今回は少し荷重が足りず、2つめの鋼材が降伏しきるところまで...
今度時間があるタイミングでもう一度荷重を増やして回してみ...
とりあえず次の段階として実際に用いられているモデルと同じ2...
**10月2日 [#zd02de0c]
以前までのモデルサイズの1/10での接触解析はparavisや荷重−...
まずジオメトリでの作成の際はモデル本体作成後、グループ作...
メッシュの作成時いつもどうりにメッシュサイズを設定後、サ...
その後メッシュを切る。今回は全体を1つとしてメッシュを切っ...
aster-studyでは以下を使用した。
,コマンド名 ,使用用途 ,カテゴリーの名前
,LIRE_MAILLAGE ,メッシュの選択 ,mesh
,AFFE_MODELE ,モデル形状の指定 ,model Definition
,DEFI_MATERIAU ,材料諸量の指定 ,material
,AFFE_MATERIAU ,各モデルの材料諸量の指定 ,material
,DEFI_FONCTION ,時間の指定 ,Functions and Lists
,DEFI_LIST_INST , 時間ステップ幅の自動調整機能の設定 ...
,DEFI_LIST_REEL ,時間ステップの指定 ,Functions and Lists
,AFFE_CHAR_MECA ,境界条件の指定 ,BC and Load
,DEFI_CONTACT ,摩擦なしの接触の指定 ,BC and Load
,STAT_NON_LINE ,非線形解析の設定 ,Analysis
,CALC_CHAMP ,計算の指定 ,Post processing
,IMPR_RESU ,出力の指定 ,output
接触解析メモのページにかかれていたものとほとんど使ってい...
DEFI_LIST_INSTの時間ステップ幅の自動調整機能の設定をしな...
STAT_NON_LINEでNEWTONの設定でMATRICEがTANGENTEだとエラー...
解析時間短縮のためSTAT_NON_LINEのCONVERGENCEのRESI_GLOB_M...
DEFI_LIST_INSTのADAPTATION よりEVENEMENTはSEUILで設定。EC...
DEFI_CONTACTのFORMULATIONはDISCRETEで設定した。
今回は線荷重で荷重を加えて載荷を行った。
**9月9日 [#s9ce2c6d]
おそらく前回一定に降伏してしまった原因は境界条件の設定を...
**9月7日 [#pfb5c5e0]
現在接触解析をうまくできないか試していっているが、岐阜高...
ただ、鋼材としての挙動があきらかにおかしく、一定にずっと...
バージョンの違いも有り、設定を変えていく必要があったり、...
**7月31日 現状の情報整理 [#nfbe50a3]
まず200×130×9の鋼板１枚のモデルを荷重5N/mm^2をかけ解析を...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
縦軸 荷重 横軸 鉛直方向変位
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
応力-ひずみ曲線のグラフと降伏度周辺のデータ （左 ひずみ...
続いて鋼材2枚で仮想材料を挟んだモデルは鋼材の厚さが9mmで...
載荷面の6点をプロットして荷重と鉛直方向の変位をグラフに落...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
このモデルの応力-ひずみ曲線を圧縮側の中央一点を上の鋼材と...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
左 上の鋼材の応力-ひずみ曲線 右 下の鋼材の応力-ひ...
とりあえず2つの鋼材は降伏していることが見受けられる。
ただ、荷重と変位の実際の実験結果を再現できそうではなかっ...
上記図のような載荷面側の片側のみ、計14個の7mmのバネ要素を...
バネ定数を32N/mmのモデルに荷重を10N/mm^2かけてみた結果次...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
上の鋼材に変位は見られたものの、下の鋼材変位がほとんど見...
ばね定数を80N/mmにして荷重を8N/mm^2かけてみた結果
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
ばね定数を56N/mmにして荷重を8N/mm^2かけてみた結果
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
応力歪曲線も確認してみたかったが、EPEQがバネ要素を用いて...
ここまで接触解析を用いずに接触を容易にモデル化するために...
そのために接触解析のための簡単なモデルから試してみる。
**24/7/16 [#u951dcaf]
全体にバネを設定していたのを載荷面側の一列のみに変更。
今まで降伏点を理論値のもので考えてしまっていたため、もう...
→（こうすることで仮想材料を入れていたときのものに変化がな...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
変位に対する剛性（ばね定数？）を32N/mmにし、荷重を10Nに設...
**24/7/10 [#e9d0b372]
バネの変位に対する剛性を30にしたが下の鋼材に接する程の荷...
バネの変位に対する剛性を15にすると鋼材は貫通してしまった。
バネの変位に対する剛性を18で荷重を15加えると
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
下の鋼材が少しだけ圧縮されているもののあまり変位が見られ...
**24/7/9 [#ydd08ccf]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
10×10間隔でバネ要素を断面全体に取り込んだモデルの解析を行...
バネ要素のz軸の剛性を10以上にしないと鋼材が貫通してしまい...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
↑荷重15 バネのｚ軸方向の変位に対する剛性20N/mm
**24/7/8 [#q3c2baf2]
前回2枚の鋼材の間に仮想材料の代わりにバネ要素を使ったモデ...
予想はしていたが、バネ一つだけでは上の鋼材が降伏したあと...
この解析を回す時バネ要素が鋼材と接している部分がメッシュ...
**24/7/5 [#cb4cf7c1]
今まで荷重が10までしか回っていなかったが仮想材料がヤング...
どうしても部材が同時に降伏してしまっているようで今まで以...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
**24/7/4 [#a3ca8113]
仮想材料のヤング率を0.01から0.1にしてみたが同様にすり抜け...
ヤング率を0.5、0.7にし、載荷を10にして解析を回したが、す...
***今後の方針 [#r1137367]
直方体の仮想材料を入れるだけではうまくいかなかったので次...
**24/7/3 [#v149df48]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
縦軸：荷重 横軸：変位
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
ポアソン比が0に近い場合、仮に縦に荷重を加えても横に膨らん...
仮想材料のポアソン比を0.001にし、解析を回したら、想定通り...
仮想材料をもっと柔らかいものにするために、仮想材料のヤン...
***現在の目標の再確認と問題点 [#j003f185]
現在は間の空いた2枚の部材の片側を降伏させて、荷重の様子を...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
**24/7/1,2 [#kb54db15]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
縦軸：荷重 横軸：変位
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
設定荷重が3.526Nだとあまり降伏が見られなかったため、荷重...
タイムステップの数も10分割から20分割に変更。
結果的に降伏は確認することができたものの１回分の降伏しか...
**24/6/28 [#qd4de2b5]
たわみの理論値計算がミスっていた
訂正 33.801→3.3801
仮想材料の大きさが大きく変わるため結果が変わるはず。
モデル図 鋼材 200×130×9 (上部、下部の部材）
仮想材料 200×130×3.3801 （中間の部材）
ヤング率 ポアソン比
鋼材 235 0.3
仮想材料 1 0.3
設定荷重 3.526N（面載荷 上部鋼材の断面のみ）
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
**24/6/27 [#cde629c2]
○前回のモデルの変更点
・線荷重から載荷側の上の鋼材の断面による面荷重に
・設定荷重を増やした
・理論値をしっかり仮想材料も考慮したものへ
設定荷重を増やして解析を回したが、設定荷重を理論値の3倍以...
Arrêt suite à l'échec de l'intégration de la loi de comp...
! La base globale est sauvegardée. Elle contient le...
! ...
! Conseils : ...
! - Vérifiez vos paramètres, la cohérence des un...
! - Essayez d'augmenter ITER_INTE_MAXI, ou de su...
このエラーによりITER_INTE_MAXIを大きくすればいいというこ...
設定荷重を2倍にしたときは解析が回った。しかし下の鋼材の外...
**24/6/~25 [#nf21dd51]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
スケールファクター 20倍
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
スケールファクター 25倍
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
スケールファクター 30倍
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
まず簡単なモデルで2段階の降伏が見られるグラフを作成できる...
・鋼材 200（高さ）×130（幅）×9（厚さ）
仮想材料 200×130×34.232（理論値）
ヤング率 ポアソン比
鋼材 20600 0.3
仮想材料 1 0.3
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
しかし今回出たグラフは直線になってしまった。（縦軸 荷重...
一見Paravisのスケールファクターが20倍のものをみたら良さげ...
荷重方法が今回は線荷重にして、荷重部分のグラフを作成した...
あと仮想材料をParavisで表示させない方法も調べておきたい。...
**24/6/~18 [#dd9728dd]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
左のグラフが過去の実験結果を示したグラフで右が解析結果に...
仮想材料を用いても降伏が一回しか見られないという小川さん...
**24/6/~11 [#hfd016ca]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
前回までのモデルで相当応力−相当ひずみ曲線を作成した。上の...
前回の応力−ひずみ曲線はマイナスの値を取っていたり、与えて...
Paravisで相当ひずみを確認するためにAsterStudyで設定を付け...
相当応力−相当ひずみ曲線を作成した際、スライスした箇所は固...
卒論に関しては、以前までは「Atype」と「Btype」と書かれて...
当然であったが15分程度では解析が回らず、設定時間が足りな...
**24/5/30 [#d4d56b6d]
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
上記の写真が今回使用した片持ちばりの載荷されている様子の...
応力がマイナスになっているのは圧縮、プラスになっているの...
**24/5/22~28 [#qc76d78d]
先週作っていたモデルをecro_lineを用いた設定の仕方でも作成...
弾塑性解析で作ったモデルをParavisで設定する方法等を練習。...
小川さんのモデルをParavisで見る時、どのように応用できるか...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
これは載荷面の右上端の点のものだが、上記のようなグラフに...
https://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori2...
**24/5/16～21 [#wfdc6df1]
・やったこと
小川さんの卒論のスライドや概要の内容を読んだ。そして卒論...
解析は回ったが、弾塑性の内容についてよくわかっていなかっ...
salome-meca 弾塑性解析（2021)
https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?cmd=re...
**24/5/15(水） [#i0e06d43]
C種について
https://www.mlit.go.jp/notice/noticedata/sgml/070/7900013...
*春課題について [#jf096965]
**24/5/1（水） [#o0b3c7a2]
メモ
・10*5*700 単純1次
鉛直1次 理論値 37.15451
,メッシュ, 要素数, 解析値, 相対誤差
,0.7,604706, 37.8469, 1.864
,1,150957, 39.3168, 5.820
,2,21270, 47.7953, 28.639
,4,2609, 75.567, 103.386
鉛直2次 理論値 120.3942
,メッシュ, 要素数, 解析値, 相対誤差
,0.7, 604706, 122.609, 1.840
,1, 150957, 127.365, 5.790
,2, 21270, 154.866, 28.632
,4, 2609, 244.558, 103.131
鉛直3次 理論値 251.1547
,メッシュ, 要素数, 解析値, 相対誤差
,0.7, 604706, 255.685, -1.804
,1, 150957, 265.629, -5.763
,2, 21270, 322.8, 28.526
,4, 2609, 509.496, 102.861
水平1次 理論値 74.30902
,メッシュ, 要素数, 解析値, 相対誤差
,0.7, 604706, 74.1075, 0.271
,1, 150957, 74.9296, 0.835
,2, 21270, 77.9669, 4.923
,4, 2609, 87.7007, 18.022
水平2次 理論値 240.7885
,メッシュ, 要素数, 解析値, 相対誤差
,0.7, 604706, 239.837, 0.395
,1, 150957, 242.486, 0.705
,2, 21270, 252.306, 4.783
,4, 2609, 283.626, 17.791
水平3次 理論値 502.3093
,メッシュ, 要素数, 解析値, 相対誤差
,0.7, 604706, 499.386, 0.582
,1, 150957, 504.892, 0.514
,2, 21270, 525.213, 4.560
,4, 2609, 590.132, 17.484
**24/4/16(火) [#ieee64e9]
共振周波数解析について
https://www.muratasoftware.com/fwp/wp-content/uploads/201...
振動設計について
https://www.fml.t.u-tokyo.ac.jp/lecture/handout/DE/yurenu...
**24/4/15(月） [#ta19d4ae]
メモ
・Texで図を挿入する際に一段に2つ以上横に並べて図を貼る方法
\begin{figure}[h]
\begin{tabular}{cc}
%---- 最初の図 ---------------------------
\begin{minipage}{0.45\hsize}
\begin{center}
\includegraphics[width=サイズ]{ファイル1}
\caption{キャプション1}
\label{ラベル1}
\end{center}
\end{minipage} &
%---- 2番目の図 --------------------------
\begin{minipage}{0.45\hsize}
\centering
\includegraphics[width=サイズ]{ファイル2}
\caption{キャプション2}
\label{ラベル2}
\end{center}
\end{minipage}
%---- 図はここまで ----------------------
\end{tabular}
\end{figure}
図がかぶってしまうときは[width]のサイズを変えたりする。
今回は2つを横に並べただけで試していないがおそらく\begin{t...
参考文献
http://www.yamamo10.jp/~yamamoto/comp/latex/make_doc/inse...
・texで表を作成した際に全体のサイズを変更する方法
→\scalebox で拡大・縮小する
\begin{table}[htbp]
\begin{center}
\caption{表題}
\label{heading}
\scalebox{0.5}[0.9]{ %ココ
\begin{tabular}{|c|c|c|c|} \hline
& & & \\\hline
& & & \\\hline
\end{tabular}
}
\end{center}
\end{table}
通常のグラフを作成して\label{heading}の後に\scalebox{0.5}...
\scaleboxの後ろの｛｝で横倍率、[]で縦倍率が変更可能。ただ...
参考文献
https://blog.goo.ne.jp/shafact/e/8219383eb299d5b6bdfd1544...
・texでの矩形選択とコピペ方法
Escを押した後、Ctrl＋vで矩形選択モードになり矢印キーで選...
選択したい箇所が決まったら、ｙでコピー。
その後、ｐで貼り付けるが、この時下の列に注意。貼り付ける...
表を作成する際、メモを開いて
%\begin{figure}[h] ←これはいらない
%1 & \\
%2 & \\
%3 & \\
%4 & \\
%5 & \\
%6 & \\
%7 & \\
%8 & \\
%9 & \\
%1 & \\
的なものを用意してtexに作成する表のしたに貼り付けておくと...
他にも表を多用する際は一度メモに使用するコマンドを貼り付...
*創造工房実習 (3年次）[#o6991381]
**23/12/14(金) [#ed1d2b72]
texの基本的な書き方https://www.str.ce.akita-u.ac.jp/~goto...
**23/12/08(金) [#o5aea6ef]
,メッシュの長さ ,要素数 ,先端変位[mm] ,相対誤差 ,計算者
,0.7 ,155192 ,0.08378905246 ,15.365 ,安藤
,0.8,138808,0.08380386491,15.350,安藤
,0.9,82587,0.083707073981,15.45 ,兼田
,1.1,38671 ,0.084201207602,14.95 ,兼田
,1.2,31929 ,0.083688 ,15.466. ,柴田
,1.3,28621 ,0.083669 ,15.4857 ,柴田
,1.4,28854 ,0.08368 ,15.47 ,佐藤
,1.5,20015 ,0.084052 ,15.10 ,佐藤
,1.6,19448 ,0.0835402938 ,15.62 ,皆川
,1.7,13801 ,0.0834355098 ,15.72 ,皆川
,1.8,12528 ,0.083733 ,15.42 ,永山
,1.9 ,11769 ,0.083924 ,15.23 ,永山
,2 ,10699,0.084076876559,15.074 ,辻
,3,3579,0.08414561753,15.004 ,辻
,4,1628 ,0.082794,16.37 ,服部
,5,1016 ,0.083033,18.89,服部
,6,839 ,-0.082882 ,16.26 ,梶原
,7,554 ,-0.080871 ,18.28 ,梶原
,8,285,0.079995,-19.20 ,工藤
,9,261 ,0.078980,-20.22 ,工藤
,10,232 ,0.081911,17.26 ,佐々木
,11,208 ,0.075676,23.56,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
**23/11/29(金） [#m5e37185]
異方性１次について
自分のデータ
,メッシュ長さ,要素数,変位,相対誤差,計算者
,4,1500,0.430156,12.52,服部
,5,432,0.282968,42.45,服部
全員のデータ
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.7,144563,0.505252,2.76,安藤
,0.8,141517,0.504692,2.64,安藤
,0.9,91648,0.502595,2.216,兼田
,1.1,27160,0.489914,0.363,兼田
,1.2,24675,0.487088,0.791,柴田
,1.3,23446,0.4868010,0.995,柴田
,1.4,17738,0.485999,1.16,佐藤
,1.5,15438,0.485180,1.33,佐藤
,1.6,15900,0.483286,1.71,皆川
,1.7,12142,0.477952,2.80,皆川
,1.8,11604,0.482085,1.9554,永山
,1.9,10391,0.470887,4.2329,永山
,2,10291,0.480910,2.19,辻
,3,2328,0.431937,12.15,辻
,4,1500,0.430156,12.52,服部
,5,432,0.282968,42.45,服部
,6,356,0.3441556,30.00,梶原
,7,196,0.213934,56.49,梶原
,8,104,0.229874,53.25,工藤
,9,81,0.232308,52.75,工藤
,10,78,0.203271,58.65,佐々木
,11,63,0.222316,54.78,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
等方性２次について
自分のデータ
,メッシュ長さ,要素数,変位,相対誤差,計算者
,4,1500,0.429254,3.01,服部
,5,432,0.428170,2.75,服部
全員のデータ
,メッシュ長さ ,要素数 ,先端変位 ,相対誤差 ,計算者
,0.7,144563,0.430124,3.22,安藤
,0.8,141517,0.430132,3.22,安藤
,0.9,91648,0.430020,3.197,兼田
,1.1,27160,0.429828,3.151,兼田
,1.2,24675,0.429836,3.15,柴田
,1.3,23446,0.42974,3.13 ,柴田
,1.4,17738,0.429797,1.3 ,佐藤
,1.5,15438,0.429958,3.14,佐藤
,1.6,15900,0.429755,3.18,皆川
,1.7,12142,0.429676,3.11,皆川
,1.8,11604,0.429829,3.1507,永山
,1.9,10391,0.429684,3.1159,永山
,2,10291,0.429620,3.10,辻
,3,2328,0.429169,2.99,辻
,4,1500,0.429254,3.01,服部
,5,432,0.428170,2.75,服部
,6 ,356 ,0.428452 ,2.82 ,梶原
,7 ,196 ,0.42591 ,2.21,梶原
,8 ,104 ,0.426074 ,2.25 ,工藤
,9 ,81 ,0.425552 ,2.12 ,工藤
,10 ,78 ,0.488382 ,17.20, 佐々木
,11 ,63 ,0.423972 ,9.0534 ,佐々木
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
**23/11/24(金） [#dbef5dfe]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
,メッシュ長さ,要素数,先端変位,相対誤差,計算者
,0.7,1455234,0.422484,1.388,安藤
,0.8,142973,0.422570,1.409,安藤
,0.9,91648,-0.420437,0.897,兼田
,1.1,27160,-0.405618,2.659,兼田
,1.2,24675,0.404349,2.96,柴田
,1.3,23446,0.404185,3.00,柴田
,1.4,17738,0.398604,4.34,佐藤
,1.5,15438,0.396593,4.83,佐藤
,1.6,16122,0.398212,4.44,皆川
,1.7,12026,0.393411,5.59,皆川
,1.8,11604,0.393668,5.53,永山
,1.9,10391,0.390695,6.24,永山
,2,10921,-0.395103,5.18,辻
,3,2328,-0.324762,22.06,辻
,4,1500,-0.155013,62.80,服部
,5,432,-0.065278,84.33,服部
,6,357,0.213062,48.87,梶原
,7,196,0.1019,75.55,梶原
,8,104,0.1158624,72.20,工藤
,9,81,0.1255118,69.88,工藤
,10,78,-0.07733,81.44,佐々木
,11,63,-0.1999,52.03,佐々木
**23/11/17(金) [#vcf762ab]
すいません、僕のデータの相対誤差が100倍し忘れていたので直...
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
,メッシュ長さ,要素数,先端変位(4隅の平均値)[mm],相対誤差,...
,0.7,198464,6.54281,-1.91,安藤
,0.8,113812,6.5104,-2.39,安藤
,0.9,40280,6.3631525,-4.60,兼田
,1.1,30055,6.3363525,-5.00,兼田
,1.2,264667,6.3043375,5.48,柴田
,1.3,25180,6.304355,5.48,柴田
,1.4,32212,6.31612,5.31,佐藤
,1.5,17753,6.1209,8.23,佐藤
,1.6,14296,6.2044625,-6.98,皆川
,1.7,13596,6.2156625,-6.81,皆川
,1.8,2866,5.737755,-13.98,永山
,1.9,6001,5.7263625,-14.15,永山
,2,5617,5.6458525,-15.355,辻
,3,2309,5.4728755,-17.948,辻
,4,617,3.6160575,45.12,服部
,5,494,3.8580375,42.16,服部
,6,581,2.50682,-62.416,梶原
,7,133,1.41225,-78.82,梶原
,8,78,1.2887175,-80.68,工藤
,9,72,1.2879925,-80.69,工藤
,10,60,1.14344,-82.85,佐々木
,11,65,1.23124,-81.154,佐々木
**23/11/10(金) [#c94eac4c]
http://www.str.ce.akita-u.ac.jp/~gotouhan/j2023/hattori23...
**23/10/27(金) コマンドと操作確認 [#t9ee68ee]
pwd : 自分の位置
ls : すべてのファイルを確認する
mkdir : 新規ファイルの作成
ls (作成できているかの確認）
cd ファイル名 ：複数のファイルから特定のファイルに入る
pwd (現在位置の確認）
ls -a : 隠しファイル
gedit ファイル名.tex & : ファイルを作りそのファイルに入...
vi ファイル名.tex : 編集（作業スペースからファイル内容を）
Esc → : → wq : 内容保存
cp コピーしたいファイル名.tex 新しいファイル名.tex : 新...
rm ファイル名.tex : ファイルを消す
cd .. : 1つ前のフォルダに戻る（他ファイルの選択が可能に）
rmdir ファイル名 ：mkdirで作ったものの削除
マウスの両クリックor中央クリック ： 選択した文章コピーa...
*salome-meca 説明 [#z8702cc2]
http://dugi6514.odns.fr/doc/v15/en/index.php?man=commande
**salome-meca 設定に関する情報等 [#e412098d]
https://www.slideshare.net/slideshow/salomemeca-59878666/...
https://qiita.com/Jun_Tatsuno/items/3e6558c1adc4dd9499fe
PARM_THETA によって，θ の値を変えられる. 時間と空間の再分...
RESI_INTE_RELA, ITER_INTE_MAXIについて
「これらは線形研究のパラメータである。最大反復回数は、...
と、線形研究の収束のために到達すべき精度reslin が与えられ...
を使用しないことを推奨する。
CORD法では、精度や反復回数をあまり高く指定する必要はあり...
実践によると、線形研究の反復は 2、3 回で十分です。したが...
デフォルトの精度で3回の反復を求めることで満足できる。ユー...
を超える線形調査の反復を行うことはできません。
一方、MIXED法では、損傷のある問題では、数十回の反復が有効...
が有効です。」 Deepl翻訳より
サンドイッチ梁について
https://www.str.ce.akita-u.ac.jp/cgi-bin/pukiwiki/?Salome...
AsterStudyでよく使うコマンド
,コマンド名 ,使用用途 ,カテゴリーの名前
,LIRE_MAILLAGE ,メッシュの選択 ,mesh
,AFFE_MODELE ,モデル形状の指定 ,model Definition
,DEFI_MATERIAU ,材料諸量の指定 ,material
,AFFE_MATERIAU ,各モデルの材料諸量の指定 ,material
,DEFI_FONCTION ,時間の指定 ,Functions an...
,DEFI_LIST_REEL ,時間ステップの指定 ,Functions a...
,AFFE_CHAR_MECA ,境界条件の指定 ,BC and ...
,DEFI_CONTACT ,摩擦なしの接触の指定 ,BC and Load
,STAT_NON_LINE ,非線形解析の設定 ,Analysis
,CALC_CHAMP ,計算の指定 ,Post processing
,IMPR_RESU ,出力の指定 ,output
**salome-mecaエラーについて [#sd314b3c]
Arrêt suite à l'échec de l'intégration de la loi de com...
! La base globale est sauvegardée. Elle contient le...
! ...
! Conseils : ...
! - Vérifiez vos paramètres, la cohérence des un...
! - Essayez d'augmenter ITER_INTE_MAXI, ou de su...
ITER_INTE_MAXIの値を大きくしてみる。（ただ大きくすればい...
ITER_INTE_PASを考慮するでも解決するらしい（未検討 調べた...
! <NonConvergenceError> <MECANONLINE9_7> ...
! ...
! Arrêt pour cause d'absence de convergence avec le...
! La base globale est sauvegardée. Elle contient le...
! ...
! 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é...
! ニュートンの非線形アルゴリズムにおいて、必要な反...
! グローバルデータベースが保存される。停止前のステ...
! !
! アドバイス: ！
! - ITER_GLOB_MAXI を増やす！
! - タンジェント行列の更新頻度を上げる！
! - 時間離散化を洗練させる ！
! - DEFI_LIST_INSTコマンドでイベント管理（例えば時...
! <EXCEPTION> <FACTOR_11> ...
! ...
! Problème : la matrice est singulière ou presque sing...
! Lors de la factorisation de la matrice, on a renco...
! (pivot nul ou presque nul) à la ligne 12414 qui co...
! ...
! Risques et conseils : ...
! * Si la ligne correspond a un degré de liberté ph...
! de corps rigide mal bloqué. ...
! Vérifiez les conditions aux limites. ...
! Si vous faites du contact, il ne faut pas que l...
! Vérifiez également les caractéristiques matéria...
! ...
! * Si la ligne correspond a un degré de liberté de...
! limite redondante. ...
! En particulier, il se peut que la relation liné...
! Peut-être devriez vous exclure certains noeuds ...
! (mots clés SANS_NOEUD et SANS_GROUP_NO). ...
! ...
! * Si le solveur utilisé est LDLT ou MULT_FRONT, v...
! car celui-ci est le seul à pouvoir factoriser l...
! ...
! * Parfois, en parallèle, le critère de détection ...
! possible de faire passer le calcul complet en r...
! en le débranchant (valeur du mot-clé NPREC=-1) ...
! ...
! * Il se peut aussi que ce phénomène soit tout à f...
! très près d'un noeud. ...
! Si le nombre de décimales perdues n'est pas tro...
! vous pouvez relancer le calcul en augmentant le...
! mot-clé NPREC du mot clé facteur SOLVEUR. ...
! Sinon, contactez l'équipe de développement. ...
**メモ [#u32aa2f5]
・サロメの過去の人のデータをコピーする際はRun_Caseも一緒...
・最大主応力と最小主応力について
物体に力が作用するときに方向によりせん断力がゼロになる面...
ParavisではPRIN＿１が最小主応力 PRIN＿３が最大主応力
最大主応力は引張応力、最小主応力は圧縮応力となる。
尚、最大主応力の結果表示でプラスの値を示す場合は、引張応力
マイナスを示す場合は、圧縮応力
最小主応力の結果表示でプラスの値を示す場合は、引張応力
マイナスを示す場合は、圧縮応力
↑実際にParavisで確認しながらやると理解しやすい
・最小主ひずみ、最大主ひずみについて
最大主ひずみは材料にかかる最大の引張りひずみの評価、最小...
・appariementについて
ジオメトリ上の固定点ループのパラメータ（デフォルト・パラ...
§ 3.4.1 .
ALGO_RESO_ GEOM = 'NEWTON'を選択すると、一般化ニュートン...
幾何学的非線形性は、各ニュートン反復で考慮され、マッチン...
(コンタクト・ニュートン(CONTACT NEWTON GENE CRIT. GEOM.)...
Les paramètres de la boucle de point fixe sur la géométri...
§ 3.4.1 .
Quand on choisit ALGO_RESO_ GEOM ='NEWTON', on bascule ve...
La non-linéarité géométrique est alors prise en compte à ...
géométrique (colonne CONTACT NEWTON GENE CRIT. GEOM.).
・ペナルティ要素（法）について
https://www.cybernet.co.jp/ansys/learning/glossary/penalt...
・岐阜高専さんのページ
http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-M...
・接触と摩擦について
file:///home/kouzou/%E3%83%80%E3%82%A6%E3%83%B3%E3%83%AD%...
・CLA
https://www.sciencedirect.com/science/article/pii/S004578...
ページ名:
![[PukiWiki]](image/pukiwiki.png "[PukiWiki]")
