☐salomeで斜張橋モデル(2D・3D)を作成する (21/05/21)
☐ 今までと同様のモデルで弾塑性解析(梁モデル、ソリッドモデル・バイリニア)(21/05/21)
☑ ケーブルの軸応力を求められるようにし、分布荷重モデルの解析が可能であるかの確認。(21/05/14)
☑ 梁モデル(温度変化)のケーブル部分のノードを3つ→2つにする。節点をヒンジ設定にする。(21/05/13)
☑ salomeでケーブル要素を用いたモデル(簡易的)に対して解析可能であるのかを確認する。(21/04/26)
日付 | 時間帯 | 作業時間 | 内容 | 立会 |
5/2 | 11:00~14:30 | 3h | ケーブル解析 | |
5/6 | 10:00~18:00 | 6h | ケーブル解析 | |
5/7 | 10:00~18:00 | 3h | ケーブル解析 | |
5/10 | 10:00~18:00 | 6h | ゼミ ケーブル解析 | 後藤さん 青木さん 石黒さん |
5/11 | 10:00~18:00 | 6h | ケーブル解析 | |
5/12 | 10:00~18:00 | 6h | ゼミ ケーブル解析 | 後藤さん 青木さん 石黒さん |
5/13 | 10:00~18:00 | 6h | ケーブル解析 | |
5/17 | 10:00~18:00 | 6h | ゼミ ケーブル解析 | 後藤さん 青木さん 石黒さん |
5/18 | 10:00~18:00 | 6h | ケーブル解析 | |
5/19 | 10:00~18:00 | 6h | ゼミ ケーブル解析 | 後藤さん 青木さん 石黒さん |
5/20 | 11:00~18:00 | 5h | ケーブル解析 | |
5/21 | 11:00~18:00 | 5h | ケーブル解析 | |
5/24 | 10:00~18:00 | 6h | ゼミ 弾塑性解析 | 後藤さん 青木さん |
5/27 | 10:00~18:00 | 6h | 弾塑性解析 | |
5/28 | 10:00~18:00 | 6h | 弾塑性解析 | |
5/31 | 10:00~18:00 | 6h | ゼミ 弾塑性解析 | 後藤さん 青木さん |
ソリッドと梁要素(ケーブル)のモデルで、温度応力を与えた解析はまわることが確認できた。
・たわみ0.03373mm
・鉛直荷重1kN
・ケーブルの温度変化-200℃
・ソリッドの断面10cm×10cmの正方形、ケーブル半径10cmの円形
梁モデル(片持ち梁をケーブルで吊った構造)の弾塑性解析は色々試して入るがうまくいっていない
5/17(木)~
熱応力・梁要素(ケーブル)片持ち梁【1m・分布荷重】たわみ-3.4429mm
熱応力・梁要素(ケーブル)片持ち梁【10m・集中荷重】たわみ-4.15603mm
5/13(木)~
・ケーブル要素を用いたモデルで今まで解析を試みていたが、一旦
モデルを梁要素で作成し、熱応力でケーブルプレストレスを再現することにした。
スケール10倍
・集中荷重1000N
・温度−200℃
・先端のたわみ量 -4.15603mm
・ケーブルのはりのノードが3つになってしまっている
→メッシュの設定の Numbwr of segments でセグメント数に1を入力することで解決
・先端のたわみ量 0.03mm(ケーブルプレストレスなし)
5/6(木)~
・直近の課題はケーブル要素を用いた簡易的なモデルを作成して青木さんのFEM解析の結果と比較する
片持ち梁(分布) | 片持ち梁(先端集中) | 片持ち梁(分布・ケーブルあり) | 片持ち梁(先端荷重・ケーブルあり) |
○ | ○ | × | × |
・固定面に載荷面が触れていても問題なく実行できた
・しかし、ケーブルの設置点と載荷部分が重なるとエラーが生じた。
! les 3116 mailles imprimées ci-dessus n'appartiennent pas au modèle ! ! et pourtant elles ont été affectées dans le mot-clé facteur : FORCE_FACE
・対策として、
(1)載荷面をケーブルの接地点から1mmほどずらして設定し直す
(2)ケーブルはケーブル要素、四面体のボックスは3Dで分けて設定していたのを、ケーブル・ボックス含めて全体を3Dに設定且つ
ケーブルはケーブル要素に設定するに変更(ケーブルとボックスの接地点を解析に含まれるようにするため)
・次のエラーメッセージ
! Erreur utilisateur : ! ! Vous voulez contraindre le ddl DRY sur un ensemble de noeuds, ! ! Mais ce ddl n'existe sur aucun de ces noeuds. !
(ノードが含まれていない?)
↓
・メッシュの切り方を3Dだけ→3D且つ、アルゴリズム(補完)2D、1Dまで設定しメッシュを切る
・ケーブルの固定点のBC&Loadの設定を変えてみる
・梁だけのモデル(ケーブル要素含む)、ケーブル単体のモデルを作成
→いずれも同じエラーメッセージだった
・ケーブル要素は非線形でしか使えない可能性がある(線形では使えない)
・非線形モデルを作成して試してみる
・梁とBARRE要素を用いてモデルを作成してみる
→BARRE要素と梁を用いたモデルで解析を回してみたが同じエラーが出た
! Erreur utilisateur : ! ! Vous voulez contraindre le ddl DRY sur un ensemble de noeuds, ! ! Mais ce ddl n'existe sur aucun de ces noeuds. !
5/2(月)
・(100,100,1000)で集中荷重も試みた
→載荷面を上面に指定
・問題なく解析できた
・やはり(100,100,10000)という寸法で解析すること(片持ち梁)に無理があったのか?..
日付 | 時間帯 | 作業時間 | 内容 | 立会 |
4/1 | 10:00~18:00 | 6h | 3年生課題 | |
4/5~9 | 10:00~17:00 | 6h | 3年生課題 | |
4/21 | 10:00~18:00 | 6h | salomeゼミ ケーブル要素解析 | 後藤さん 青木さん 石黒さん |
4/22 | 10:00~17:00 | 6h | ケーブル要素解析 | |
4/26 | 10:00~18:00 | 6h | ゼミ ケーブル要素解析 | 後藤さん 青木さん |
4/27 | 10:00~18:00 | 6h | ケーブル要素解析 | |
4/28 | 10:00~17:00 | 5h | ケーブル要素解析 | |
4/30 | 11:00~18:00 | 4h | ケーブル要素解析 |
4/30(金)
モデルが(100,100,10000)だったので(100,100,1000)にしてみて試してみた。
↑10cm,10cmの断面に対して、10mが長すぎるのではないかと予想
・集中荷重(先端)
・10mにしたらうまくいった
4/28(水)
片持ち梁のモデルの解析ができていないので優先して取り組むことにした。
モデルの荷重のかけ方は等分布荷重だったので、載荷面をボックスの上面に指定して解析を回していた。
→これを線載荷に変更して同様に試みた。
! <S> Exception user raised but not interceptee. ! ! The bases are fermees. ! ! Type of the exception: error ! ! ! ! Solver MUMPS: ! ! The solution of the linear system is too vague: ! ! Computed error: 1.58999e-05 ! ! Acceptable error: 1e-06 (RESI_RELA) ! ! ! ! Advices: ! ! One can increase the value of the key word SOLVER/RESI_RELA. !
・26日のエラーとほとんど同じ
・同じモデルで集中荷重(先端)で試した
→うまくいかなかった
4/26(月)
salomeで、ケーブル要素を用いた簡易的なモデルを作成し解析を回してみた。
↑青木さんがやられていた解析プログラムで算出した値と比較する
! Exception user raised but not interceptee. ! ! The bases are fermees. ! ! Type of the exception: error ! ! ! ! Solver MUMPS: ! ! The solution of the linear system is too vague: ! ! Computed error: 5.59007e-05 ! ! Acceptable error: 1e-06 (RESI_RELA) ! ! ! ! Advices: ! ! One can increase the value of the key word SOLVER/RESI_RELA. !・収束判定基準がおかしい..?
! <S> Exception user raised but not interceptee. ! ! The bases are fermees. ! ! Type of the exception: error ! ! ! ! les 9 mailles imprimées ci-dessus n'appartiennent pas au modèle ! ! et pourtant elles ont été affectées dans le mot-clé facteur : FORCE_FACE !
1. 【集中荷重・1m】(hari4)
・Geometryの作成:新しいエンティティ->基本オブジェクト->点->(0,0,0),(0,0,1000)作成
・基本オブジェクト->線->Vertex1、2->Line1作成
・点->(0,-1000,0)->Vertex2、3->Line2作成
・Fuse->Line1、Line2->Partition->Fuse1、Line1、Line2
・Line1グループ作成->node(0,0,0)-kotei->node(0,0,1000)-saika
・Line2グループ作成->node(0,-1000,0)-kotei2
・Mesh1:Line1->ジオメトリのグループ作成->ジオメトリ-Line1->node-kotei-saika
・Mesh2:Line2->Wire Discretisaition->Number of Segments->セグメント数1->ジオメトリのグループ作成->ジオメトリ-Line2->node-kotei2
・Compound_Mesh->Mesh1、2
・Astestudy:メッシュ->Read a mesh->format-Med
・Material:Model Definition->Assign finite element->Finite element->Everywhere->Mechanic->POU_D_E
・Model Definition->AFFE_CARA_ELEM->POUTRE->RECTANGLE->Group-Line1->CARA-HY-HZ->Value-100-100
・Material:Define a material->Linear isotropic elasticity->ヤング率200000->ポアソン比0.3
・Material:Define a material->Linear isotropic elasticity->ヤング率200000->ポアソン比0.3->Thermal expansion 1.17e-0.5
・BC and Load:EnforceDOF->Group node-kotei->Laison-ENCASTRE->DX~DRZ-o
・EnforceDOF->Group node-kotei2->Laison-ENCASTRE->DX~DRZ(DRX以外)-o
・FORCE_NODLE->Group node-saika->FY-1000
・Post Processing->CREA_CHAMP->Model-AFEE_MODELE->TYPE_CHAM-NOEU_TEMP_R->OPERATION-Assingment->
Group element-Line2->Value- -200->NOM_CMP-TEMP
・Material:Assigm a material->Model-model(AFEE_MODELE)->Marerial assigment->Group-Line1->mater->
Marerial assigment->Group-Line2->mater0->External state-Temperature->Reference value-20->Group-Line2->Field-unnamed5->
・Analysis->Material field->Structural element->Model->Load=Load->Solver->Method-MUMPS->RESI_RELA-0.01
・Post processing CALC_CHAMP->SIGM_NOEU->MODEl-model->Material field->fieldmat->Structural element
・Output->Set output result->Format-Med-> Result-DEPL-SIGM_NOEU
1m.comm
(commファイル:一度保存してAsterstadyでエクスポートすることによって利用可能。hdfファイルではないため、Geometry等は反映されない)
2. 【集中荷重・10m】(10m)
1.集中荷重1mとほとんど同じ。異なる点は、Geometryで10mのモデルを作ることとAstestudyの
SolverでNPREC->-1、RESI_RELA->20にする(収束基準・条件を緩和する)
3. 【分布荷重・1m】(bump2)
1.集中荷重1mとほとんど同じ。異なる点は、BC and Loadで載荷のオプションを変更することのみ。
FORCE_POUTRE=_F(FY=1.0, GROUP_MA=('saika', )), MODELE=model)
❌線形モデル