佐々木の卒論日誌

日誌†

3年†

4年†

創造工房†

10/26†

11/2課題†

手計算 6.667

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2018/sasaki/zu4.pdf

11/9課題†

120mm

手計算 -0.41667

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2018/sasaki/zuzu.pdf

100mm

100mmのモデルは、120mmの時と比べ手計算との誤差がとても小さくなっている。これは、120mmよりも大きく変位しているためであると考えるため、120mmのモデルよりも壊れやすいのではないかと考える。

11/16課題†

120mm

手計算 0.41667

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2018/sasaki/zu.pdf

張出しがない方が誤差が大きくなった。また、穴を開けて荷重を加えたほうが、メッシュサイズ１の時などの正確である値においては、誤差が小さくなった。

11/30課題†

手計算 -0.49167

http://www.str.ce.akita-u.ac.jp/~gotouhan/j2018/sasaki/zuzuzu.pdf

今までで一番全体的に誤差が小さくなった。メッシュサイズが、1,1.3,1.5の時が理論値との誤差が小さくなり、正確な値であると考える。

春休み課題†

モデルは作り終えたが、メッシュを切る際、Compare - Equal to が出てこなかったが、無視してそのままメッシュを切った。 また、２材料のポアソン比とヤング率の入れ方がわからないが、とりあえず計算してみたがやはりエラー。

エラーメッセージ 　　<S>ユーザー例外が発生しましたが、インターセプトされませんでした。 ！

  ！基地は閉鎖されています。 ！
  ！例外の種類：エラー！
  ！ ！
  ！上記の10個のメッシュはモデルに属していません。
  ！それでも、それらはキーワードfactor：FORCE_ARETE！で割り当てられています。

【4/15】

kou,mokuとグループを作り計算したがエラー。 fuseする前に載荷線の穴を空けた後に結合させたモデルを作ったが、こちらもエラー。

S>ユーザー例外が発生しましたが、インターセプトされませんでした。 ！

  ！基地は閉鎖されています。 ！
  ！例外の種類：エラー！
  ！ ！
  ！上記の10個のメッシュはモデルに属していません。
  ！それでも、それらはキーワードfactor：FORCE_ARETE！で割り当てられています。

【4/17】 model definitionのmodelのAt Least Oneの中のFinite elementのGroup of elementの選択を、mokuとkouだけではなく、saikaとkoteiも追加したら、計算が正常に出来たが、これでちゃんとできているのかは不明。

【4/22】 kouの部分を木材のヤング率などにして計算したところ、変位は-0.3286であり、以前計算した異方性材料は、0.5ぐらいであった。また、鋼材に直して計算したところ、変位は、-0.0957であった。

【5/15】 比べると、ポアソン比も違い、fieldmatもsaikaやkoteiを選択するかしないかが異なった。 　

【5/20】 120mmの梁をポアソン比0.016で作ってみたところ、mokuだけの時、変位は平均-0.504234となり、前回のティモシェンコ梁の変位は。-0.5060769091であるため、計算は正しいと考えられる。

また、kouの値をいれて計算したところ、平均-0.106392454となった。前回の異方性材料の時の理論値は、-0.49167であり、mokuだけのものと比べると誤差は2.6%であった。

しかし、今回鋼材の部分も含めて求めた理論値は-0.623411となり、kouもいれたものとの誤差は、82%だった。明らかに何か間違ってますね。

【5/29】 手計算が間違えていたため計算しなおしたところ、理論値は、-0.1178812742となり、kouも入れたものとの誤差は、9.7％となった。

【6/3】 せん断弾性係数を79230にして計算したところ、変位は-0.10505となり、誤差は10.9％となった。

卒論†

6/3ゼミ†

まず、モード解析に載っているYoutubeのことがsalome2018でできるのかをやってみる。 また、モード解析のページを分かったことややり方を更新する。

6/5†

Youtubeを見ながらモデルを作ってはみたが、何をやっているのか自分でもわからないので、解明していきたい。また、Paravisまで（Youtubeでいうと27分ぐらいまで）やったのだが、Time/modeの選択をmodeにするとmodel_DEPL以外は私のモデルだとうまく色がつかなかった。

あと、Asterstadyにおいて、output→Rasultの際、Youtubeはunnamedを選択していたが、私のはmodeしか選択できなかった。

6/10ゼミ†

モード解析の立体の方を作成してみる。やり方を意味は分からなくてもいいからとりあえず、前回の面要素や今回の立体要素の手順をメモ。

6/11†

立体要素のYoutubeのモデルを作ってみた。計算も動画通りに行うと上手く行き、Paravisも面要素のものと組み合わせると何をしているかは分からないが、振動しているのか変位しているような動きが見られた。次回、手順を確認しながらモード解析メモにまとめる。

Paravisのやり方†

• 立体要素の動画どおりにまずはSolidcolorをmodeDEPLにかえて色を付ける
• 再生ボタンのようなやつを押すと色が変わった（はず）
• ここで動画だと開いているmedファイルを開いて読み込み出すのだが、クリックしてもダメだったため、今度は面の動画のように、左のオブジェクトインスペクター？のところにあるTime/ModeをModeにして、
• ApplyをおすとWarpByVector（電卓？のようなマークのとこにあるいかにも曲がってますよ－！ってマーク）を押せるようになったので押してみると曲がったモデルがでてきた。
• そうすると面の動画で選択できなかったsolidcolorのmodeDEPL［11］などが見られるようになるため、変位の様子？が確認できるようになった。

6/12†

構造メモにとりあえずAsterstadyのやり方をメモした。

6/17ゼミ†

• まずは固有振動数や固有振動モード、減衰定数など、ヤング率とか代入して解析解とどれくらい合うのかを計算。
• その後、創造工房の時のようにメッシュの大きさなどを変えて、どうしたら解析解に近づくか調べる。

6/17†

• Youtubeのモデルだとなんだかよく分からなかったので、&link(SALOME-Mecaの使用法解説,http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2)の10-00.pdfのモーダル解析を参考にして、
• 100mm,10mm,10mm
• ヤング率132000、ポアソン比0.343、密度8.96e-9のモデルを作って解いてみたところ、

となり、大きく異なってしまった。

• なので同じように100mm,10mm,20mmのモデルを作ったところ、

となった。

• ちなみにメッシュサイズは1で切った。ResultのNOMCHAMも入れた。
  • 使用法解説が誤差が大きくなっているのは、メッシュサイズが大きいためだと思われる。

• 不明点
  • 使用法解説に書いてあるとおりに理論値の計算をしたのだが、
  • &link(解析解,http://mechanics.civil.tohoku.ac.jp/bear/nisikozo/s4node5.html#SECTION0410433000000000000000)は、μlが変化しているが、使用法解説は変化していない。
  • 今回は、密度(8.96e-9)\( \times \)10mm\( \times \)20mmでμl=1.792e-6を求めたため、変化する計算の仕方が分からなかった。
  • また、解析結果と理論値を使用法解説に書いてあるように対応させたが、
  • なぜこのように対応するかはわからないのでもう少し理解したい。

メモ†

• Paravisのmodes［1］－624.345などの624.345が解析結果(固有振動数)
• NMAX_FREQが固有振動数

6/24 ゼミ†

• &link(解析解,http://mechanics.civil.tohoku.ac.jp/bear/nisikozo/s4node5.html#SECTION0410433000000000000000)と&link(SALOME-Mecaの使用法解説,http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2)の
• 理論値を比べて同じか確かめる。(mの値とかμの値とか)

• Scale Factor を大きくしてApplyを押すとものすごく振動する。
• Paravisのmodes［1］－624.345などの[1]とかの1は1次とかの1。

6/24†

• &link(解析解のｍは,http://mechanics.civil.tohoku.ac.jp/bear/nisikozo/s4node5.html#SECTION0410410000000000000000)単位長さあたりの質量である。
• 密度をρとして 断面積をAとすると，m=ρAである
• そのため、μの計算は密度(ρ)\( \times \)断面積(A)で正しく、解析解と使用法解説の計算は同じであると言える。
• よって解析解のm=使用解説のμである。
• また、μ1l=1.875などの数値は、coshλcosλ+ 1 = 0
• (λ=μl)
• この式は超越方程式で代数的に解けないので数値計算(たとえば Newton-Raphson 法)を行って解くらしい。
• &link(日本機械学会宇宙工学部門,https://www.jsme.or.jp/sed/guide/dynamics5.pdf)の5－4に載っていた。
• また、3.52は1.875の2乗であった。

• μ1l=1.875,ω1=1.875^2\( \times \)√EI/ρAl^4

7/1ゼミ†

• もっと細長いモデルを作り、計算する。理論値と比べる。
• Paravisでどれが1次か2次かを見極める
• 要素数による違いを考える。
• 図面見る

7/1†

• 1次、2次は、オブジェクトインスペクターの方のVectorsをいじる。
• [0]が1次、[1]が2次、、、という感じである。
• 1次と2次は断面2次によって曲がる方向が異なる。
• modal2030でやり方を学んだので、もっと細いモデルでやる。

7/2†

• (Dx,Dy,Dz)=(1000,10,20)のモデルを作り計算してみた。
• いつも通りのAsterStudyではVERI_MODEがエラーだったため、
• VERI_MODEを選択してNoを選ぶと計算できるようになった。
• メッシュサイズ0.5,0.7,1,2,5で比べてみた。

であった。

• メッシュサイズ0.7と0.5はメモリーが足りないそうで計算できず。
• 後日また計算してみる。
• 創造工房の時のように、メッシュサイズが大きくなれば誤差も大きくなる。
• また、メッシュサイズ1,2,5の時のParavisを見比べてみると確かにちょっっっっっっとずれていた。

7/3†

• メッシュサイズ0.7は計算できた。
• メッシュサイズ0.5は計算できなかった。
• 全体的に見ると今のところメッシュサイズ0.7が一番誤差が小さい。
• 以下がエラーメッセージ
• 概要
• ステージRunCase_1_RunCase_1_Stage_1の出力メッセージ
• ステージRunCase_1_RunCase_1_Stage_1の出力メッセージ - [上/下]
• 完全なジョブ出力は、ディレクトリu '/ home / kouzou / 2019 / soturon / modal4_Files / RunCase_3 / Result-RunCase_1_RunCase_1_Stage_1 / logs'にあります。

7/8ゼミ†

• 要素数を減らす代わりに、2次要素を選んで解いてみる。(やり方は後藤資料にある)
• 曲げ以外についても解析解と合うのか。ねじれだったり。
• 減衰定数を見る方法を調べる。
• そのためにはまず、今は非減衰振動であるので、&link(SALOME-Mecaの使用法解説,http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2)の11-00.pdfとかの周波数応答の減衰有り減衰無しとかを参考にしながら、やってみる。
• 減衰振動するためには、最初から外力(sin一波とか)を与えたり、粘性係数与えたり、、などが考えられる。
• 最終的には、固有振動数や減衰定数を調べたい

7/19†

• meshを右クリックして、「convert to quadratic」をクリックし、applyすると、(中間節点を設けた)2次要素になるので粗いメッシュでも精度がよくなる。http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?Salome%2dMeca2017
• のジオメトリのグループをメッシュに反映

• ２次要素選んで解いた。メッシュサイズ1だと解けなかったので、2と1.5で解いた。
• 7/2の表に追加した。
• 2次要素で解いたほうが理論値により近かった。

7/20†

• Salome解説方の11.00減衰無しをやってみたが、
• vect=CALC_VECT_ELEM(OPTION='CHAR_MECA',
  • CHARGE=BLOCAGE,);
  • vectass=ASSE_VECTEUR(VECT_ELEM=vect,
  • NUME_DDL=NUMEDDL,);

• dynaharm=DYNA_LINE_HARM(MATR_MASS=MASSE,
  • MATR_RIGI=RIGIDITE,
  • LIST_FREQ=Lfreq,
  • EXCIT=_F(VECT_ASSE=vectass,
    • COEF_MULT=1.0,),)

• の部分が見つからなかった。

• また、DEFI_LIST_REELを入れないと値は解説法と同じになるが、
• これを入れると大きく異なってしまうため、上記の見つかっていない部分を見つける必要があると感じた。
• Lfreq=DEFI_LIST_REEL(DEBUT=100.0,
  • INTERVALLE=_F(JUSQU_A=2000.0,
    • PAS=100.0,),);

• メモ
• modal3ファイルの
• BC and Load のFOCE FACE のFzに１をいれる。
• Pre AnalysisにあるVECT_ASSEは一応てきとーに入れたが分からない。

不明点†

• CALC_VECT_ELEMが見つからない
• ASSE_VECTEURらしいものは、Pre Analysisの中にVECT_ASSEというものがあったが、中身がよく分からず、入力したいものがなかった。
• DYNA_LINE_HARMが見つからない。
  • Solverの関係らしいが、Solverというものがあるわけではなく、どこかの中にSolverがある形なので、どこのSolverか分からない。(たとえば、AnalysisのCALC_MODEのなかのSolverのような)
• IMPR_RESUがOputionのIMPR_RESU_SPというものがあったがこれも打ち込みたいものがなく別物なのか、、、
• MACRO_MATR_ASSEのCHAGE=BLOCAGEが不明。

メモ†

• MODE_ITER_SIMULTはSalome2018では、AnalysisのCALC_MODESになったらしい。
• 変動荷重を変化させる周波数の定義は、Functions and Lists　にある。

7/29ゼミ†

• 減衰定数の解き方を調べるのと並行して、
• 金峰2000年橋のモデルを作る。
• 図面を見て骨組みを作る。手すりとか接合部分は作らず、ボックス1つで良い。
• 減衰定数は、YoutubeなどでSalomemeca modal Analysisなどで調べる。
• また、前回の引き続き曲げ以外のねじれなどが解析解を調べて合うかも

8/26†

• 金峰2000年橋のモデルの円の部分は、コンパスや、パソコン上で円と合わせて円になるか確認する。
• SalomeのPravisのねじれがどうやるかもう一度確認する。

8/29†

• モデルがやっと作り始められた。
• 円環ではなく円柱を使ってアーチを作った。
• trans7(アーチ部)とFuse1(板)
• 次は支柱と横の脚？の部分を作る。

8/30†

• 昨日作成したモデルが間違っていたため、作り直し。
• 半径30000と半径28500の円柱を用いてアーチを作った。
• アーチの支間長が36900なのだが、アーチの高さが分からないので、7500にしてはみたが、支柱がほんの少しだけ長さが足りず、アーチに届かず。
• おそらく、板の厚さ600か、7500が間違っていると思う。

9/4†

• 床版の厚さは440であると思われるため、440に直したものを作成したが、当然隙間は埋まらず。
• アーチを300上にあげてみて、中心部に近い支柱は、アーチと床版に繋がったが、一番外側の支柱はまだアーチには届かず。
• アーチのところに6344と書いているのだが、これがアーチの高さなのか、、、？
• 6344にもしてみたが当然埋まらないため、どこが間違っているのか、、

9/14†

• kinpou0914
• 床版が左右に傾いているようなのでそこを修正。
• そして図面においての高さはよく分からないので支柱は床版に合わせて、高さを調節することにする。
• 支柱3は、下から6400のところにおいた。
• 次回は、支柱1の高さを調節して、その後アーチを調節する。

9/19†

• 前回に続いて、支柱1と2を作成した
• アーチを当てはめてみたが、支柱と全部繋がるようなアーチは存在せず。(床版にぶつかってしまう）
• 支柱の長さが間違っているのか、、？
• まず図面が読めない

9/26†

• &link(卒論概要,http://www.str.ce.akita-u.ac.jp/cgi-bin/gwiki/wiki.cgi?%b8%e5%c6%a3%c8%c9%a5%bc%a5%df#i26)
• モデルをとりあえず完成させた。
• 床版とアーチを基準とし、支柱の長さを調整した。
• アーチの高さは土台（500）を抜いて、6994（土台からアーチの中心部までの長さ6344＋アーチ部分の中心から外側までの長さ650）とした。
• 支柱1は1140、支柱2は887、支柱3は779長くし、
• 支柱1が6075、支柱2が3705、支柱3が2185とした。
• メッシュを切ってみたものの、ディスクの容量が足りないため？メッシュを切ることができず。また、メッシュを残した状態だと保存できなかったため、削除。
• とりあえずここまでのことを中間発表でまとめる。
• kinpou926

9/27†

• ねじれは太いほうがよく出る。
• [4]がねじれ
• Salomeの使用解説の12を中間終わったらやってみる。
• ねじれの理論値は、解析解のω3がねじれだと思われる。
• そこと比較して誤差を出す。
• 支柱の距離は、支柱1不明、支柱2は8650、支柱3は12750であったので、
• 中間発表が終わり次第モデルを直す。

9/30†

• ねじれの理論値が私の計算が間違っているのか、誤差がとても大きいのかが分からないので、後で確認する。
• とりあえず明日からはパワポを作って、Texをかく。
• 時間があればねじれの正誤性についてやる。

10/1†

• 理論値の式を直す

10/2†

• sibup2
• chuukanhappyou

10/11ゼミ†

• このまま立体モデルを作ってもメッシュを切れない可能性が高い。
• 一回支柱部分を梁要素(線)で置き換えて作成してみる
• その前に梁要素で振動解析ができるのかをやってみる。
• 梁要素ができるようなら梁の先に立体要素をくっつけてみて、解析してみる
• 梁要素で解析するなら、ヤング率やポアソン比の他に、断面2次や断面積を入れるところがあるのでは…？

10/18ゼミ†

• 梁要素について調べる

10/21†

• 1つの線ではなく、複数の線で作ってもメッシュ切れず
• FuseしてPartitionしてもメッシュ切れず。

10/24†

• hariseikou
• 両端に10,10,10のボックスをおいて、100mmの線を10分割して、
• 10本のlineをfuseして、その後両端を含めた12個の点でPartitionして、
• Partiton1と両端のBoxをCtrl押しながら同時に選択して、
• Mesh（今回はメッシュサイズ２で切ってみた）切ってメッシュを作成するとできていた
• ジオメトリのグループなどはまだ作ってない。
• もう一度同様にlineだけのモデルを作ったが、やはりメッシュは切れず。
• 片方だけにBoxをつけたモデルはメッシュを切ることができたが、
• テトラ？でしかメッシュを切ることができなかった。
• Asterでmeshを選択する際、Compound Mesh？があるようなので、
• 次はどうやってMeshをまとめられるのか調べる。

10/25ゼミ†

• lineは、1Dで、アルゴリズムがWire Discretisationでやればメッシュを切ることができる。
• メッシュサイズは、Max Sizeが選べるようになるので、それで行う。
• まずは、梁の片持ち梁のモデルを解析して、理論値との比較
• その後、今までと同様に振動解析で、理論値と比較
• 2次要素などでもやってみる

11/1ゼミ†

• 梁要素のものを解くことを頑張る
• 剛性マトリクスや伸び剛性、曲げ剛性をどこかで与えるはずなのでどこか調べる。

11/5†

• SECTIONは、Model DefinitionのAFFE_CARA_ELEMのBARREで設定できた。
• CARAとVALEがよく分からず、梁要素のページの参考にしているものと同じにしようとしても、CARAはAしか選択できず、VALEは関数のようなものだった。 －Groupの名前も選択できなかったため、とりあえずsaikaと入れてはみた。
• また、今までAFFE_MODELEのところで3Dを選択していたが、2D BARREを選択してみた。
• RHOはとりあえず7850とした。
• CHAM_MATER=materi,
• CARA_ELEM=caraelem,
• EXCIT=_F(CHARGE=charmeca,),
• INFO=2,);
• の部分はまだ見つけられない。

• とりあえず計算させてみたが、案の定計算できず。

   <S>例外ユーザーが発生しましたが、インターセプトされませんでした。！
    ！ 拠点は閉鎖されています。！
    ！ 例外のタイプ：エラー！
    ！！
    ！ GROUP_MAで要求された操作はできません！
    ！ メッシュmeshには何も含まれていないためです。！

11/8†

• test3においてのedgeやdiskなどはlineのジオメトリのグループであった。
• 今回は試しに、kotei側のlineをdisk
• saika側をedge
• その2線の間をconnectlとした。
• また、 CARA=AFFE_CARA_ELEM(MODELE=Model,

                     POUTRE=_F(SECTION='GENERALE',
                               GROUP_MA='edge',
                               CARA=('A','IY','IZ','EY','EZ','JX','RY','RZ',),
                               VALE=(4e-05,1.33333e-11,1.33333e-09,0.0,0.0,4.94237e-11,0.01,0.001,),),

• の部分は、model definition,AFFE_CARA_ELEMで選択し、定義できる事が分かった。
• Iは断面2次、Aはせん断弾性係数、Eはせん断中心についてなので今回は０、JXはねじれ定数、JYはそりねじれ定数、Rはとりあえず１。

Load=AFFE_CHAR_MECA(MODELE=Model,

                   PESANTEUR=_F(GRAVITE=9.81,
                                DIRECTION=(0.0,0.0,-1.0,),),);

• このコマンドは、BCandLoad,assign mecanical load,gravitational…のところでできた。
• test3だと0,0,-1となっているが、軸が違い－をつけると赤字になったしまったため、0,1,0とした。

• test3のコマンドをとりあえず打ち込んでみたが、計算結果はエラー

    ！ <S>例外ユーザーが発生しましたが、インターセプトではありません。 ！
     ！ ベースはフェルメです。 ！
     ！ 例外のタイプ：エラー！
     ！ ！
     ！ キーの後にキーワード「NOEUD_2」または「GROUP_NO_2」を示す必要があります！
     ！ オプションCOQ_POUのワードファクターLIAISON_ELEM！

• 私は今回3本の線で行ったため、GROUP1も2もどちらもconnectlにしてしまったため、エラーが起きたのかもしれない。

11/8ゼミ†

• 梁を解いてみる
• プログラムを使って、長方形断面のねじれ定数やそりねじれ定数を求めてみる
• まずは長方形でやってみてから、梁で。

11/8†

• group2を途中の点でやったらエラー

！ <S>例外ユーザーが発生しましたが、インターセプトではありません。 ！

  ！ベースはフェルメです。 ！
  ！例外のタイプ：エラー！
  ！ ！
  ！プログラムエラー。 ！
  ！ ！
  ！条件が満たされていない：！
  ！ .false。 ！
  ！ファイル！
  ！ /home/A21173/smeca/yamm_build/V2018.0.1_public_build/V2018_public/tools/src/！
  ！ Code_aster_stable-v136_smeca / bibfor / debug / contex_param！
  ！ F90、147行目！

memo†

• GROUP_NOは、ノード、点
• GROUP_MAは、線

11/13†

• test3ともう一度見比べたところ、修正点があったため修正したところ、以下のエラーメッセージ

！ <S>例外ユーザーが発生しましたが、インターセプトされませんでした。 ！

  ！拠点は閉鎖されています。 ！
  ！例外のタイプ：エラー！
  ！ ！
  ！問題：行列は特異またはほぼ特異です：！
  ！行列を因数分解するときに問題が発生しました！
  ！ （ゼロまたはほぼゼロのピボット）自由度に対応するライン78で！
  ！上に与えられた。 ！
  ！ ！
  ！リスクとヒント：！
  ！ *線が物理的な自由度に対応している場合、それはそうです！
  ！おそらくムーブメント！
  ！剛体がひどくブロックされています。 ！
  ！境界条件を確認してください。 ！
  ！連絡を取る場合、構造はそれを「保持」してはいけません！
  ！連絡先によって。 ！
  ！材料の特性もチェックしてください（ヤング率、...）。 ！
  ！ ！
  ！ *ラインがラグランジュの自由度に対応する場合、ラインはありません！
  ！状態の疑い！
  ！冗長制限。 ！
  ！特に、非常に豊富な線形関係が生じる可能性があります！
  ！接触条件。 ！
  ！おそらく、いくつかのノードを接触条件から除外する必要があります！
  ！ （キーワードSANS_NOEUDおよびSANS_GROUP_NO）。 ！
  ！ ！
  ！ *使用されるソルバーがLDLTまたはMULT_FRONTの場合、使用できます！
  ！ MUMPSソルバー！
  ！これはそうでない行列を因数分解できる唯一のものだからです！
  ！正の定義。 ！
  ！ ！
  ！ *並行して、MUMPSの特異点検出基準は時々！
  ！悲観的すぎる！それにもかかわらず頻繁に残っています！
  ！この基準を緩和することにより、完全な計算を渡すことができます！
  ！ （キーワードNPRECの値を1または2増やす）または！
  ！プラグを抜くこと（キーワードNPRECの値= -1）または計算を再開することにより！
  ！より少ないプロセッサで。 ！
  ！ ！
  ！ * X-FEMでは、この現象が非常に正常である可能性もあります。
  ！クラックパス！
  ！結び目に非常に近い。 ！
  ！失われた小数の数が多すぎない場合（最大10の小数）、!
  ！失われた小数の数を増やすことで計算を再開できます！
  ！許可：！
  ！キーワードファクターSOLVERのキーワードNPREC。 ！
  ！それ以外の場合は、開発チームに連絡してください。

• test3と同じコマンドファイルを作ったが、エラー。 －GROUP_NO２がよく分からない。

11/15ゼミ†

• 境界条件を入れる
• ただの片持ちで解いてみる。

11/19†

• 境界条件を入れて、ただの片持ちで解いたがエラー !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

    ! <S> Exception user raised but not interceptee.             !
    ! The bases are fermees.                                     !
    ! Type of the exception: error                               !
    !                                                            !
    !  error (S) met (S) during the checking of the assignments. !
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  ！ <S>例外ユーザーが発生しましたが、インターセプトされませんでした。！
  ！ 拠点は閉鎖されています。！
  ！ 例外のタイプ：エラー！
  ！！
  ！ 割り当ての確認中にエラーが発生しました（S）プット（S）。
  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

• &link(POUDEについて,https://sites.google.com/site/codeastersalomemeca/home/code_aster-1/bimu-youso/pou_d_e-youso)
• このリンクを見たところ、Aはせん断弾性係数ではなく、断面積と書いていたので試しにやってみたが、エラー
• koteiにDX,DY,DZだけでなく、DRX,DRY,DRZも0を入れていたので入れてみた。
• BCandloadにもconnectをedgeを選択してやってみた。（これはリンク先を見ても何をしているか分からなかったので、とりあえず）
• Analysisにもconnectを加えた
• 結果は、上記のエラーメッセージと同様にエラー
• haripoude1119

• そして今更だが、線をわけなくても1本の線でもメッシュは切れた

11/21†

• assignmentsでエラーなのかと思い、色々いじってみたが、結果同じエラーメッセージ
• haripoude1121で色々試した
• &link(SALOME-Mecaの使用法解説,http://opencae.gifu-nct.ac.jp/pukiwiki/index.php?SALOME-Meca%A4%CE%BB%C8%CD%D1%CB%A1%B2%F2%C0%E2)の16.00のシェル要素も参考にしてみる

• NOEUとELEMが違うのかも
• SIGM_NOEUなどのoutputのところが違うのかも
• BARREでもPOUDEでも色々やったが変わらず
• POUDEは節点に回転の自由度があるものらしいので、点3つにしてlineを2つに分けた。
• シェルとビームの違いを考えつつ、シェル要素の例題と比べたが、PostProcessingの部分で大幅に違うので、これがシェルとビームによる違いなのかエラーの原因なのかを考えたい
• 具体的には、CALC_CHAMP(計算_領域）で今までやっているのだが、今回のシェルの例題や他の梁要素のものはELEM（要素）を使っている
• CALA_ELEMを探したが見つからず。

• どこをいじっても同じエラーメッセージが出てくるので、コマンドの問題じゃないのかもしれない

11/22ゼミ†

• とりあえずシェル要素でとく
• 一般断面(general)ではなく、矩形断面(rectangle)でとく

11/25†

• 矩形断面を設定して、円断面のYoutubeを見て解いたところ、計算することができた
• https://www.youtube.com/watch?v=gt3Oylnj03M
• この動画と違う部分は、
• model definitionのpoutreのsectionを円形ではなく、矩形(rectangle)にする
• CARAはHYとHZを10,10で打ち込んだ
• BC&Loadで動画では、LIAISONをENCASTRE選択して固定条件を終わらせていたが、それだと計算ができなかったため、矩形は、LIAISONではなく、DX,DY,DZ,DRX,DRY,DRZに0を打ち込む
• outputのところでNOMCHAMでDEPLのみ選択

• また、普段と違う点は、
• Pravisの際に、オブジェクトインスペクターのGenerateVecotrsにチェックを入れないとWarpbyvectorは使えない
• 梁要素なので、Post Processingは不要
• また、X軸方向に1000で梁を作成

• 今回梁要素を初めて計算することができたので、立体要素と比べる

• モード解析もしてみようと思ってやってみたところ、エラー

• <EventError> <ALGELINE2_74>！

    ！ ！
    ！ VERI_MODEはそうではありません！計算されたモードの検証エラー：少なくとも1つの検証基準
    ！ respecté。 ！
    ！ ！
    ！アドバイス：！
    ！メソッドを使用したくない（リスクと危険がある！）場合は、キーを変更して計算を再開してください！
    ！キーワードVERI_MODEの下にある単語！
    ！品質基準で制約の少ない値を使用することにより、！
    ！ -または、オプションSTOP_ERREUR = 'NON'を使用します。

11/29ゼミ†

• 理論値と比較する
• 固有振動数(Natural frequency)とSalome-MecaやCoad_Aster、ビームやframeなどと合わせて検索して調べる
• シェルの振動解析と立体要素の振動解析を比較して、違いから梁の振動解析を求めてみる

11/29†

• 理論値と比べたところ、誤差はなかった
• 計算結果は、0.666なのに対し、理論値も同様に0.666だった。

12/2†

• 梁のモード解析ができた！！！！！！！！！！！！
• AnalysisのCALC_MODESで、OPTIONはPLUS_PETITE（小さいという意味だが、だいたい最初はこれを選んで解いてみようみたいなことを、モード解析のページに載っているシェル要素のYoutubeで言っていた）
• CALC FREQはNMAX_FREQ（固有振動数）をいれ、
• 梁要素では、VERI_MODEをEditし、STOP_ERREURのチェックを外し、Noにする。

• 理論値と比較したところ、手計算が間違っていそうだったので、もう一度計算する。
• また、梁の長さ1000mmでやっていたので、ねじれは確認できないと思うので、100mmのモデルを作成して、密度を変えて計算し直す。

12/4†

• 100mmの梁を作成してみたところ、誤差は以下のようになった

• また、1000mmの梁も手計算が間違っていたため、計算しなおしたところ、以下のようになった。

• となった。
• 誤差はほとんど無かった。
• メッシュの大きさは、どちらもナンバー オブ セグメント15(要素数15)

• また、メッシュの大きさ1(要素数100)でやったところ、セグメント15と同じであった。

Paravisのアニメーションのやり方†

• 表示、windows、アニメーションの概要
• WarpByVector
• フィルター、Mechanics、Normal modes animation(real)
• アニメーションの概要のコマンドのところで、Normal modes animation(real)を追加できるところがあるので（＋のとこを押す）
• そして上の方にある緑の再生ボタンを押すと再生される。

12/6ゼミ†

• 梁とソリッドの片持ちはもう比較が終わっているので(すぐできる状態なので）
• 次は、ラーメンを作って比べる
• 梁の長さは、高さ3mのラーメン。300mm✕300mmの正方形でソリッドは、梁を軸としたら150mm飛び出している状態
• 最終的には、キングポストトラスで梁とソリッドを比較する。

12/8†

• まずは、Boxの方のラーメンを解いてみたところ、1DのところでMaxsizeでメッシュを切ろうとするとどれだけ大きさを大きくしても切れなかった。
• 1Dにおいてナンバーオブセグメントでメッシュを切ると切れたので、試しにアスターに進んで計算してみたところ、今までのモーダル解析と全く同じファイルで計算できた。
• ちなみに、両足の面を固定した。
• Paravisにおいても変位の様子がよく出ていた。
• 今回は、セグメントを15で大きくメッシュを切ったので、精度はあまり良くないと思うので、もう少しセグメントの数を増やしてから、次回理論値と比較したい。

• また、梁のラーメンのモーダル解析も行ったが、今までの梁のモーダル解析のファイルだと計算できず。
• 明日改めてもう一度解いてみる。

12/9†

• 梁の方をもう一度解いてみたところ、Meshを切る際に、Maxsizeで切ると、エラーが出たが、ナンバーオブセグメントで切ると計算が成功した。
• セグメントの数は15
• この間できた梁の片持ち梁のモーダル解析のファイルで計算成功した。
• 固定は、両端の点で固定。グループオブエッジは、全てのlineを選択。
• Paravisでも立体要素と同じような変位の様子を見れた。

• 正確性は低いと思うが、セグメント15での立体要素と梁要素は以下のようになった

• 立体要素のラーメンでセグメントの数100で計算したところ、メモリーが足りなくて計算できなかった。
• メモリー5000、time 02:00:00で計算したができず。

12/15†

• 梁要素のラーメン構造のモード解析を要素数を増やして計算したところ、セグメント330（要素数990）を超えると計算できなくなった。（梁の長さは100mm）
• ソリッド要素のラーメン構造は、セグメント50で、メモリー7000、time　5:00:00では計算できず。
• 今後の予定としては、
• ねじれの理論値を出す
• ラーメン構造の正確性を高める
• ラーメンの理論値を求めて比較する。
• 概要の文量によっては、前半の振動解析についての部分を少なくする。
• 金峰2000年橋はカットした
• ラーメンのスクショや動画をのせるのもあり

12/19†

• 梁の大きさは100,20,10で、ヤング率6000、ポアソン比0.4、密度3.8e-7で統一して概要にかくことにする。
• 最初の振動解析の部分のグラフを書き直し、
• 梁の解析とも比較する
• ラーメンはまだ理論値と比較していないため、梁とソリッドが同じような振動をしていることを述べるところで終わる。

1/8†

• boxのモーダル解析を梁と同じ条件にして計算しなおした。

  メッシュサイズ

  要素数

  相対誤差

  １次曲げモード１次曲げ

  0.7	76407	0.769%
  0.9		0.922%
  0.9(２次要素)	39892	0.519%［誤差最小］
  1	26920	1.05%
  1(2次要素)		0.529%
  1.2	17966	1.223%
  1.2(2)		0.542%
  1.5	9928	1.614%
  1.5(2)		0.567%
  2	4360	2.48%
  2(2)		0.618%

  1次曲げモード2次曲げ

  0.7		-2.159%
  0.9		-2.090%
  0.9(2)		-2.267%
  1		-2.042%
  1(2)		-2.258%
  1.2		-1.962%
  1.2(2)		-2.245%
  1.5		-1.817%
  1.5(2)		-2.222
  2		-1.500%［誤差最小］
  2(2)		-2.178%

  2次曲げモード

  0.7		-3.606%
  0.9		-3.455%

  0.9(2)

  -3.857%

  1		-3.326%
  1(2)		-3.846%
  1.2		-3.153%
  1.2(2)		-3.833%
  1.5		-2.762%
  1.5(2)		-3.807%
  2		-1.898%［誤差最小］
  2(2)		-3.755%

  3次曲げ（ねじれ）

  0.7		2.5096%
  0.9		2.579%
  0.9(2)		2.389%［誤差最小］
  1		2.648%
  1(2)		2.393%
  1.2		2.725%
  1.2(2)		2.398%
  1.5		2.905%
  1.5(2)		2.410%
  2		3.346%
  2(2)		2.436%

• という結果になった。
• 曲げモード別に１番誤差が小さいメッシュサイズで梁と比較するのか、全て同じメッシュサイズで比較するのかを相談する。
• 全体的に見て、全てのモードで誤差が小さかったのは、1次要素の1とか1.2あたりだと思う

• この3000✕300✕300のラーメン構造だと、大きすぎるのか、ナンバーオブセグメントでしかメッシュを切れず、メッシュサイズではなく要素数がだいたい同じになるよう考えたところ、メッシュサイズ1に当たる要素数では、メッシュが切れなさそうである。

1/10ゼミ†

• ソリッド要素のメッシュは、切れるところまでで切る。
• メッシュの粗さも影響はあるが、ソリッド要素と骨組みで比較
• ソリッドと骨組みの結果も概要にのせる。

1/13†

• 先生から直してもらったTexファイルを確信し、ラーメン以外の部分を修正した。
• 明日明後日で、ラーメンの梁要素とソリッド要素の結果のグラフを作成し、Texに貼る。

オープンキャンパスについて†

• 折り紙は夏紀さんのパソコンに入っているkikagaku5kakuやkikagaku-5kakuが良いと思う。Paravisの動かし方は不明。
• 虹橋も夏紀さんのパソコンに入っているオープンキャンパスのもの。

最新の100件