後輩に†この日誌はあんまり親切じゃないので参考にならないかも。参考にするなら、卒論本体に解析のやり方とかちょっと書いたのでそっちを。あと菊池さんとか河原さんの日誌はとても親切なのでそっちを見てほしい。聞きたいことがあったら後藤さんに連絡先を聞いたりしてもいいと思う。作った主要なモデルのbrepファイルtate.brepと3kaku.brepはk2にあげておく。私が使ったパソコンはドアから見て一番左の奥と、右から2列目の奥から2番目です。変わるかもしれないけど。中身はごちゃごちゃしてて私もよくわからない。brepファイルをimportして使ってください。屋根材の最大主応力比率とか、菊池さんのやり方と私のやり方とで全然違うみたいで、比較するなら自分のやり方で算出しなおしたほうがいいと思う。あと荷重のかけ方もちがうみたい。異方性はsalomeでは十分にできないけど、他にできるツールを試してみてもいいと思う。salomeでエラーが出たら、1回閉じるとか、パソコン変えるとかでも何とかなる時があるので、うんざりしないで試してみると良い。縦長断面モデル三角形断面モデル改良版。斜めになって弱くなってる側面のトラスを起こして、あと下の細々してる奴を削った。近藤さんのパソコンではできた†salome2015だと出来なかった計算が、2013ではできた。今後は計算だけ前のパソコンでやります。↑2/2†asterでrunしたあと、コマンドラインは立ち上がるがエラーで止まる。下はmessファイルの中身の英訳。考えていたモデルを作成できた。実際作ってみて、スカスカしている部分に部材を追加していきたい。そのあといよいよ解析。エラーが異様に出てくる場合(ただの6面体の6面全部に出てくる時とか)は、一回閉じて開き直すと普通にできた。逆に、普通に切れてたやつが後で切れなくなることもある。パソコン替え†salome2015が入ってる奥のパソコンに移動した。体積が大きくなるとfuseやらcutができなくなるようなので、大幅に単純化できれば助かるところ。つづき†fuseに関するエラーは、結合部分の複雑さが原因っぽい。cutでいらない部分を削って解決するみたいだけど、cutでもエラー出たりする。中間発表ではとりあえず今できるものを発表する。fuseに関する謎エラー†fuseできない、もしくは勝手にcompoundになるみたいなエラー?で先に進めない。fuseできたと思ったら今度はshellになってて空洞あいてたり。新しく作ったモデルの断面。今までのエラーは、複雑だったり小さい面になったりで出てきてた感じがするので、できるだけ単純な面だけでこれを並べて作っていくつもりだが、今のうちにエラー対策をしておこうと思って、この時点でメッシュ切ってみた。compoundしたところ(下のぎざぎざ)がダメだったので、全部fuzeして解決。新しいモデルづくり†meshエラーがでないように、fuseをできるだけ使わず面を単純にしようと試みたが、やっぱり時間がめちゃくちゃかかる。中間発表12/20まで時間が全然ないので、今までどおりエラーがいっぱい出るやり方でやるけど、それも当然時間がかかる。三角形互い違いモデル(仮)の内部に直方体を通す中絶対狭い。断面の三角形は一辺8000で、通す直方体は3000×4000少なくともこの通し方だと足元めちゃくちゃ危ない軸方向に視線を向けるとこんな。2000×2000√3断面が通れる。こんな橋とか通路があったら面白そうだと思うけど、やっぱり中が狭すぎる。後で強度調べて、そのままボツにすると思う。六角形モデルと三角形たがいちがいモデル(どちらも仮)、内部の狭さが問題salomeで作ってみて、人間大の断面を持つ直方体を通して、部材にぶつかるか見る。去年の4年のデータのバックアップが河原さんの引き出しに入ってて、件のメッシュファイルもまるごと手に入った。これで自分の計算が合ってるか確かめた。3角形モデルは計算違いとかでなく単に構造的に弱いと分かった。そのデータの中にcommファイルもあったので、比較のためにも今後の設定はそれと揃えようと思う。メッシュが切れない問題については、今後も直面するだろうと思うし、解決できるに越したことはないが、とりあえず保留にしておきたい。近藤さん曰く、純粋に点と線だけからモデルを作ったなら、こういうエラーは出ないらしい。スケッチでもいくつか作る。(aptrte())っていうエラーがかなり出る。近藤さんと河原さんに教わったの試せば出来そう。応力も変位も出なかったのは面載荷なのにFORCE_ARETEって入力していたからでした。これから、剛性とか屋根材の寄与率とか出して比較とかしていよいよTeXで原稿書く。propertiesのVolumeじゃ全然あってないことが分かった。1/4モデルなのに1/1モデルより体積多かったので、少なくともどっちか違うはず。創造工房の時も体積違ってたと思う。今回から修正された滝田さんプログラム&link(unvc3d4.f90,http://www.str.ce.akita-u.ac.jp/~gotou/programoj/ccx/unvc3d4.f90)を使います。??最初なのでシンプルなモデルにしました。屋根材の最大主応力は目視でみました。菊池さんの五角形モデルと比べて高すぎるし合ってるかどうかあやしいところ。AM4:50†やっとメッシュ切れて狂喜乱舞。結局、軸方向が自動的にLになるわけじゃないということで落ち着いた。三角錐とか、六角形を並べた形の屋根付き橋のモデル作成中。先輩方に倣って、「最大主応力に対して屋根材の主応力がどの程度現れているか、比率を求める」ことで、屋根材の強度への寄与をはかる所存。斜めのモデルでやってみる†どのモデルでもx=Lになるので今度は斜めにしてみた。うまく行ってないので持ち越し。あとまた帰省するので、家ではTeXの使い方とかを調べておく。橋のモデルを考えたりもそろそろしておく。前回のL=zは計算ミスで、やり直すとx=L,y=T,z=Nになった。前々回のモデルもx=L,y=T,z=Nでつじつまが合う。今回x=20,y=30,z=10のモデルで同様にx=L,y=T,z=Nで誤差0になった。長いほうがLにならない†x=30,y=20,z=10のモデルで、x方向に圧縮したが、今度はL=zで誤差0になった。次回、引き続きポアソン比の確認。前回まで、条件が反映されてなかったり、post-proが出て来なかったりした理由が分かった。拘束条件をつけるグループ名に「x=0」とかを使っていたのが悪かったらしく、「x_0」とかに直したら万事うまくいった。L=z,T=x,N=yのつもりが、L=x,T=z,N=yになってるようだ。z方向に長いモデルなんだけども。次回、他のモデルで試してみる。反映されてなかった†まさかと思ってヤング率、ポアソン比、弾性係数全部0にしてみたけど、それまでと全く同じ結果が出た。