★有限要素法(Finite Element Method)★
入出力部をいちから作るのはたいへんだが、計算部だけならそうでもないでしょ >>4 なんでそういう商用アプリを使うかな。 FreeMATなりOctaveを使えよ! オーソドックスにC言語でやるのは、かなりきついよね? なんで、そういう発想になるのかな? 物理が、分からないと、何やっても、ダメだって!!! 差分法!!!!! abaqus student がインストールできないんですが 誰かわかりませんか?涙 Spectral Element Method とか FEM なんかGPLでソース配布してるプロジェクトが 少なからずある >>15 普通にインストールできたけど…? 当方の環境はWindows XP + Abaqus Student Editon 6.9(書籍の付録) Abaqus Studentってのはフリーなの?有料ならサポがあるよな >>21 本に付録してるだけなのか? 広める気あるのかないのか・・・・w >>22 個人でAbacus使おうなんていうヒトはいないでしょwww ちなみにAbacusっぽいやつなら CALCULIX http://www.calculix.de/ http://www.bconverged.com/products.php (FreeのWindows版あり) へどうぞ。 >>15 1万円でAbaqusが使えると聞いて、私も購入しました。 ですが私もインストールできません。 windows7だから? ttp://www.google.co.jp/search?hl=&q=%83K%83X%83%89%83C%83e%83B%83%93%83O よほど大規模解析でもしないのならFlexPDEなるソフトはどう? 評価版(要素数や支配方程式の数の制限あり)ならフリーでwindows上で動くよ。 easy-σとかいうの大学で使った 原理はよくわからんけどとりあえず結果は出た 使う側は原理なんて考えなくていいよ。 イレギュラーな結果を見定める時に多少知ってればいい程度だ。 3次元図化したいのですが、フリーのソフトでオススメありますか? gnuplot でfEMの結果を表示でキルンでしょうか? 単純なグラフならよさそうだけど。 >>41 質問を小出しにするやつにFEMをやる資格はないと思え 図のインパクトで考えるなら やっぱParaViewだわな 100万要素ぐらいでもサクサク動いて くれるソフトありますか? そこまで行ったらソフトの良し悪しじゃないもんなぁ 粒子法でGPUバリバリって方があってるかもねw フリー、有料どちらでもいいんですが、粒子法のおすすめ ソフトありますか? C言語とかで、地道にFEMを計算している人っています? 層化の矢倉の婆が朝から異常行動繰り返して迷惑しているんですが 日に日にエスカレートしてまったくやめる気配がまったくありません FEMの場合は、どう考えてもFortranが多いでしょう >>57 みたいに何を解くのか なにを求めるのかも書かずに アホのような質問されても 誰も答えないっしょ。 今、有限要素法の勉強をちまちまやっているのだけど、難しくてさっぱり分からんw Octaveでプログラムを書いているのだけど、境界条件を入れた全体剛性行列のdetがゼロになって解けん・・・ 正六面体の一次要素一つの形状を作って片面を全固定(123自由度)、反対面に荷重をかけているのだけど 正六面体だと解けない理由ってのが有るのかな?まあプログラムの何処かが悪いのだろうけど 何から調べていけば良いのかアドバイスお願いします >>62 そんな大それたもの読めません(本気なのかー!) >>66 2次元の要素をまだ実装していないのです。申し訳有りません。 実装するにしてもまだ時間がかかりそうです 使える要素は四面体と六面体の一次要素です。 四面体も変なのですがまだ細かくつめていません (時間はかかると思いますが、2次元問題の方も作ってみます) Microsoft Office Excelでゆうげんようそほうはかのうですか? >>79 コマンド一覧は http://www.usennet.ne.jp/ ~kiyoshi/manual/CalculiX.html つうかその程度自分で見つけられない人間って・・・。 まぁまぁ Caliculixに関する情報は結構イッパイあるから、 全く初めての人にとっては、どこから読めばいいのか さっぱりわからん場合もあるので。 自分で見つけたれば、ググればいいし、 教えてあげたい人は教えればいいでしょ。 並列FEMのフリープログラムって どこかでdownloadできますか? >>79 2ch病なのかなあ コマンド一覧 を コマンドー覧 と誤読したorz Gmshから読み込んで3次元のtet4のコードを書いているんだが、paraview可視化の部分がうまくいかない。 手動で座標とconnectivityを書たテストは、うまく出ているので書き方は間違ってないと思うのだが? でいま積んでる。知ってるヤツいる? 電波テロ装置の戦争(始) 魂は幾何学、コピー出来る公安はサリンオウム信者の子供を40歳まで社会から隔離している オウム信者が地方で現在も潜伏している それは新興宗教を配下としている公安の仕事だ 発案で盗聴器を開発したら霊魂が寄って呼ぶ来た <電波憑依> スピリチャル全否定なら江原三輪氏、高橋佳子大川隆法氏は、幻聴で強制入院矛盾する日本宗教と精神科 <コードレス盗聴> 2004既に国民20%被害250〜700台数中国工作員3〜7000万円2005ソウルコピー2010ソウルイン医者アカギ絡む<盗聴証拠> 今年5月に日本の警視庁防課は被害者SDカード15分を保持した有る国民に出せ!!<創価幹部> キタオカ1962年東北生は二十代で2人の女性をレイプ殺害して入信した創価本尊はこれだけで潰せる<<<韓国工作員鸛<<<創価公明党 <テロ装置>>東芝部品)>>ヤクザ<宗教<同和<<公安<<魂複<<官憲>日本終Googl検索 947 名前:彼氏いない歴774年[sage] 投稿日:2010/07/05(月) 15:45:01 ID:7xkp916Q 女子大生風二人。電車にて。 A「最近まとめサイト見てんだけどさー」 B「なんのまとめだろう」 A「コピペのまとめサイトだよ。わかる?」 B「あーー…わかる」 A「わかってなさそうだなー。ま、いいけど。モジョの人たちが聞いた会話のコピペが多くて」 B「モジョ…?」 A「超面白いよ。一回モジョ板も見てみ」 B「の…奇妙な冒険?」 A「は?まあとにかくさ、どっかでモジョの人が私たちの会話聞いたりしてんだよ」 B「ネアンデルタールの次あたり…?」 A「何言ってんの?つーかさ、私たちの会話がコピペのネタになってるとか超素敵じゃん!皆に教えたいじゃん!あ私ここで降りるね!また明日ー!」 B「(手を振りながら)モジョとは…なんじょ」 お嬢ちゃん、後ろ後ろ。 >> 102 自己解決すた。vtkのconnectivityのノード番号って0から始めなきゃいけないのね strand7 ってソフト使わなきゃいけないんだけど、 使い方を簡単にまとめたサイトとか誰か知りませんか? delFEMってソフト使ってみようと思ったけど全然コンパイル通んないですね・・ これ推奨の環境ってMACなの? んー一応VSとかminGWとか使ってみてるんだけど エラー出まくる・・これはソース直して使うのがデフォ? 直せばコンパイル通すくらいはできると思うが、、 コンパイル通らない状態で公開されてるって認識で良いのだろうかってことが知りたい。 亀田も真っ青!前人未到の犯罪5階級制覇や!! 違法駐車王 塩ホ容疑者 器物損壊王 仕置容疑者 飲酒運転ほう助王 死沖容疑者 無法バーベキュー王 Jおき容疑者 児童虐待王 SALTWOOD容疑者 1.2.0でやったら動きました〜♪ お世話になりましたノシ 衝突解析でっきるっかな。 http://peacelandia.blogspot.com/2011/01/blog-post_1275.html 2011年1月24日 田原総一郎の講演料 去年の話になりますが、10月19日に弘前市で行われた「原子燃料サ イクル意見交換会」に関して、安藤晴美県議が情報開示請求し、その結 果、わかったことです。 企画は、資源エネルギー庁と青森県の共催で、(株)RABサービスに 委託されて行われており、委託料は129万6750円。その内訳は、 ジャーナリストの田原総一郎氏の講演料が110万円。交通費8万円。司会料3万円。管理手数料(食費含む)2万5千円。 前に話題になってる、 Abaqus Student Editon 6.9(本のおまけのやつ)ってmacでも動きますか? 対応OSがわかれば教えて下さい。 Abaqusの固有振動数抽出がまったくうまくいかん 単位もN mm t 系でやってんのに メッシュをめちゃくちゃ細かくしたら少しだけ理論値に近づくが・・・ Abaqusだろうがどれでも 基本的に数値は大差ない。 変ならモデルに妙な接合部があるとか 材料データ間違ってる。 単純な片持ち梁なんだけど、どこかおかしいのかな。また明日みてみるか スケッチで50*50*1000の梁を作って、密度7.8e-9、剛性206e3、ポアソン比0.3 を入れてるだけなんだけど・・・ メッシュ図をつけたもの http://www.42ch.net/UploaderSmall/source/1326465524.png ちなみに理論値は (1.875/1000)^2*sqrt((206000*50*50^3/12)/(7.8e-9*50*50))/(2*pi) =41.504 Hz だから十分あっていると思う。 あ、念のためにだけど、要素は二次要素(*ELEMENT, TYPE=C3D20, ELSET=Eall ) にしたほうがいいよん 固有値解析のポスト処理について質問です。 諸々の事情によりParaViewを使いたいと考えています。 変形表示はWarpByVectorでできそうなのですが、 これをさらにクネクネアニメーションさせたいです。 どんなやり方がありますか? >>125-128 レスありがとうございます! 今日(きのう)もう一度、最初からモデルを作って解析をすると、 Mode1 41.439 Mode2 41.448 Mode3 256.72 Mode4 256.78 Mode5 706.19 Mode6 706.34 ※ 断面:50*50[mm]正方形 長さ:1000[mm] 密度:7.8e-9 [t/mm3] 剛性:206e3 [N/mm] Poisson:0.3 となって、理論値に非常に近くなりました。どうやらメッシュの切り方等に問題があったようです・・・ 単位系ばかりに気を取られていて しかしAbaqus、メッシュを少し荒目にするとMode1が30Hz位になって誤差が大きくなりますね。 静的な応力解析の時は、メッシュの大きさが大きくても、そこそこの精度は出るのに。 詳しく説明をしていただき、本当にありがとうございました。 ソフトの使い方だけでなく有限要素法についても、さらに勉強しようと思います。 メッシュが一次要素になってない? 二次じゃないとおかしな値になるよ。 >>133 単なる非定常計算結果のアニメーションは できるんですが、固有値解析のように モードシェイプに対応する変位データが ひとつしかない時に、変形前と変形後の ノード座標値を(マイナス変位側まで) ParaView上で補完して、物体が振動して いるようなアニメーションがしたいのです。 FEMAPなどで固有値解析結果を表示するのと 同じことを、ParaViewでできないかなぁと。 変位データは自分で作成しているので、 連番ファイルで補完データを書き出すことも 不可能ではないのですが、数十個ある モードシェイプ情報に対して全て補完データを 出力すると連番ファイルが1000個近くなるので、 何かうまいやり方があれば助かります。 >>132 1次要素って、2次要素より「曲がらない」=「硬い」ので、周波数としては上がる傾向に あるみたいです。 4要素だと、2次要素なら43Hzですが、1次要素(C3D8)にすると77Hzになっちゃいます。 神道教育が起源の、創価学会から、はじまった公明党には、票を入れられません。なぜなら、日本は、国教国家ではなく、世俗国家 なのです。これは、ほぼ違憲であり、戦後の、反省の立場から、公明党が政権をとることを想定すると、神道国家に回帰するため、 アジアや旧連合諸国にたいする国際的なプレゼンスや国際関係を著しく悪化するものと思われます。また、別の観点から、創価大学を持っているといえども、「学会」と言う名称も、国際的には、本来の学会の意味でないため、非常に不適切な団体です。 創価学会との無関係を訴える公明党ですが、フジテレビの選挙報道時は、「創価学会のみなさん」と不自然なくらい必ず、放送します。 ブッシュドクトリンに協調した小泉自民・公明は、「テロとの戦い」を表明するも「オウム真理教」に破防法を適用できない失敗をし、野党に追いやられました。 つまり、宗教団体は、特に、世俗国家においては、政党を持つべきではありません。 小泉元総理は、公明党と公約を守って、「自民党をぶっこわす」だけでなく、原発も ぶっこわしたのです。これを小泉の奇跡と言いましょう。世襲制にも反対しましょう。 質問させていただきます。 Abaqusで、粘着挙動を選択しますとデフォルトでペナルティー法による粘着となっています。 このペナルティー法というのは、どういったものなのでしょうか。 調べてみますと、FEM解析における接触の基礎としてのペナルティー法や、計算が速い代わりに食い込みを許すということ等、おおまかなことは分かったのですが、 このペナルティー法でA部材とB部材を粘着させると、完全に一体として接着されたものとみるのでしょうか? 具体的には、A部材をStep1で伸ばしておき、そこにStep2でB部材を張り付けるということをしたいのですが、張り付けに、ペナルティー法の粘着を使用するとどうなるのかがわかりません。 ヘルプぐらい嫁。 購入先のサポートに問い合わせろ。 すいません、この中でFemap with NXNastranで接触解析をしたことがある人はいますか!? SOL601(NX Adv Nonlin)で二次元の弾性体半球同士の接触解析を強制変位を使ってしようとしているのですが、Connectionsが上手く働いていないのさ、両物体が思い切り貫入するような結果になってしまいます。収束も最後までしません。 Slide Lineを使った二次元解析ではなく、ConnectとSOL601を使って接触解析をするなら、プロパティやリージョンをどのように設定すれば良いのでしょうか!?アドバイス下さい! 板要素のDKMQ要素の理論について知っている人いますか? 離散kirchhoff板理論(DKT,DKQ)とひずみ仮定(MITC8)を 拡張された理論であるのは分かったのですが、 ここから先はさっぱり分かりません。 なにか参考になるソースはありますか? 今見ている文献は A new discrete Kirchhoff-Mindlin element based on Mindlin-Reissner plate theory and assumed shear strain fields. PartU 有限要素法の理論的な部分に詳しい著書を教えてください >>146 ありがとうございます その本が図書館にあるか調べてみます もし、アマゾンなどで簡単に購入できる本があれば教えてください FEMはFORTRAN。 例外は認めない! なんつったってAnsysだってFORTRANだ。 有限要素法はね、 個体の荷重ー変位曲線を割り出すのに使います <Sale>中国語版の書籍<IT|3D|NN|理系|上古音などなど>(格安) 280円より http◎//lang-8○com/194279/journals/1457398/ 致違いで申し訳ございません。 どこにもこれから質問するものに関連するものがなかったのでここで質問させていただきます。 今現在研究で(大学生4年です)解析をしています。 解析手法としては線形加速度法という陰解法を使用していますが、 この手法の欠点としては、ステップ間隔を短くしなければなりません。 実際研究で使う地震波は、100秒もあり、扱う対象のものは 固有円振動数だけで330個、その中の最小のものより1/60程度にしなければ なりません。そうしなければ、解は拡散してしまうからです。 そこで、0.0005秒間隔で100秒までとなると20万ステップしなければなりません。 構造解析のソフトウェアではそれにどれだけいったいかかるのか 図ったことはありませんが、最適化ですのでこれを少なくとも100回は回していかなければ なりません。線形加速度法で扱う行列は、330行列。。。 この前0.0005秒間隔で20000ステップやるのですら、 1時間かかていました。。 20万ステップですから単純計算で 一回終えるのに10時間。これが最低でも100回回さなければならないため、 1000時間かかってしまいます。 そこで、解析手法を変えようと思っているのですが、 どの方法が適切でしょうか?何卒ご協力お願い申し上げます >>158 >最適化ですので 何をどう最適化したいのか、エスパーしようにも どうしたものやら。 地盤?それとも構造物? 真面目に変位量の時間履歴を取りたいのなら 愚直に解くしかないんじゃないの? それをどうやったら他の手法で置き換えられるのか、 それはむしろ君自身が研究テーマとして取り組むべき 問題なんじゃないかと思うのだが。 いま京の研究公募してるから、それに応募してみたら? 今の手法でも結構頑張れちゃうかもよ? >>158 時間が限られてるなら最適化の項目を減らしなさい。 どうしてもやる必要があるなら数年掛かろうがやりなさい。 >>158 そのまえに0.005秒のデータってあんの? 俺は15万要素くらいで60秒解析してたけど、実際0.05秒とかでも、大して結果かわんなかったぞ そもそも地震関係やっていて そんな周波数を解析するって時点で無能君。 いや、こいつのテーマは 様々な地震動に対して最適な構造を見つけるってことじゃねえのか 最適化ってそういうことだと思ったんだが つーか書いてて思ったが1000時間って大したことないよね パソコンがやるんだし。連続して解析できるんだろ? 2台使えば500時間だ。一ヶ月もかからねえぞ そういえば、数年前にノートPCで解析してるペーパーがあった。 防衛大だったと思う。飛行機とか衝突してたから ただ、終わりのほうに、ノートPCでも3ヶ月かければ解析は可能だって書いてて、 こいつは何を伝える気なんだと思ってワロタ >>166 良かった! 心配でずっと眠れなかったよ。 OpenFOAMのStressAnalysisに詳しい方は居ませんか? 集中荷重の指定方法とか何故繰り返し計算を強いられるのかとか教えてください 分かっていない事が多すぎるのでソースを読めはご勘弁を >>169 OpenFOAMでStressAnalysisをする時点で、アレ。 calculix とか、そーいうの使うほうがいいかと。 >>170 StressAnalysisはオマケみたいなものなのですかね・・・ Calculixも平行して為してみているのですがプリ部分でまだ試行錯誤していますw FreeCADでSTLファイル作成→netgenでメッシュ作成→cgxで荷重、境界条件作成→ccxで計算→cgxで結果を確認 と言う感じにやってますが、単一部材のケースはわりとすんなり行くのですが複数部材のケースがどうしたら良いのか チラ裏でし >>171 アセンブリは素直にcgxでモデリングする方がてっとりばやいかと。 pythonが使えれば ttp://dip28p.web.fc2.com/calculix/netgen2calculix/index.html っていうのが使えるらしい。 どなたか塩谷豊という人物をご存じありませんか? どんな情報でも結構です どうかよろしくお願いいたします >>172 cgxか〜複雑なモデリングには向いてなさそうだけど使ってみるよ netgenまでは複数部材は考慮されているみたいで、サンプルで提供されている 複数部材を定義されたgeoファイルを読み込み、abaqus形式で吐き出すと 要素の集合名が分けられて出てきます CAD系で複雑なモデリングをしていくのが理想に思えるけどここらが自分の限界なのかな STEP形式ファイルは複数の部材を考慮した出力って出来るのでしょうか? あ、FreeCADで複数部材のモデリング方法が分からないからダメか・・・チラ裏でした >>174 Ass’yの難しさとして、 □■ A B 部材Aと部材Bが一つの面でくっついているとすると、 CGXみたくモデラーで作るときは共通平面ということで作っていけばいいけど CADからもっていくとAとBの共通平面は別々の平面になっているはず。 そこでAとBをくっつけなきゃならないんだが、たいてい、その別々の平面を 接触(固着)させて解析させることになる。 ここがちょっと難しいところですね。 モデリングでは境界条件が難しいよね。 平面と平面が面で完全に接合している物体なんて現実的ではないし。 割り切るしかないんだろうね…。 複数部材でできたモデル(片持ち梁)をうまく作ることが出来たからメモ代りにカキコ 1) FreeCADで梁を2つ(接触した形)作成し、2つ選択した状態でSTEP214形式でexport 2) gmshでSTEPファイルをオープンする。一旦geoファイルに保存しそれで開き直す。 3) Geometry:Elementary entities→Coherenceで隣接したジオメトリの節点を結合する。 4) Geometry:Physical groups→Add→Volumeで2つPhysical groupを作成する。 5) Mesh:3Dでメッシュ作成、メニューのFile→Save asでAbaqus INPを選択して適当な名前でメッシュファイルを作る 6) cgx -c メッシュファイル名 で開き境界条件、荷重を作る。適当にccx用のインプットファイルを作成して解析、結果確認 分かっている人には当たり前の事かもしれなかったけど、これを分かるまで大分時間がかかったよ 複雑な形状だと節点が重ならなかったり上手く行かないんだろうな・・・ >>177 >3) Geometry:Elementary entities→Coherenceで隣接したジオメトリの節点を結合する。 お、これは知らなかった。一回やってみよ。さんくすです。 四面体一次要素の形状関数を教えてください。とても複雑になる様な気がしてます・・・ 要素の体積を A として,要素内の点が与えられたとき,対応する節点の対面と その点が作る四面体の体積をA1,A2,A3,A4 として N1=A1/A, N2=A2/A, N3=A3/A, N4=A4/A となるらしいのです 節点をP1,P2,P3,P4としてmaxia風に書くと A: 1/6*determinant( transpose(matrix((P2-P1), (P3-P1), (P4-P1))) ); 要素内に与えられる点をP0とすると A1: 1/6*determinant( transpose(matrix((P2-P0), (P3-P0), (P4-P0))) ); A2: 1/6*determinant( transpose(matrix((P1-P0), (P3-P0), (P4-P0))) ); A3: 1/6*determinant( transpose(matrix((P1-P0), (P2-P0), (P4-P0))) ); A4: 1/6*determinant( transpose(matrix((P1-P0), (P2-P0), (P3-P0))) ); となって複雑に・・・何かが間違えている > 169 OpenFOAMってVOF法みたいので構造も有限体積で計算していたはず。 弱連成もできると聞いたことあるので、有限要素のソルバーもあるかも知れんが。 詳しくは開発元が書いた論文よめ http://onlinelibrary.wiley.com/doi/10.1002/nme.1409/abstract 機械部品の強度解析するのに、FEMのソフトを導入しようと思うのですが、 ANSYS、NASTRANなどを検討していますが、何かお勧めはありますか? メインは線形で、非線形も少しやろうかと思っています。 運用実績が無いなら、まわりで使っているソルバー又は聞ける人がいるソルバー なんちって 固有値計算プログラムが上手く動かん、固有ベクトルがぐちゃぐちゃだ 固有値はハウスホルダー+QR法で何とかoctaveのeig関数の結果と一部一致したけど (ゼロに持っていく値が振動して収束しない時があるんだよね…) 三重対角行列の固有ベクトルを求めて(怪しい実装)それをハウスホルダーの副産物を使って逆変換(とても怪しい実装) 見事に撃沈、誰かエスパー求む >固有ベクトルがぐちゃぐちゃだ メッシュがグチャグチャなのが原因と予想だな >>184 片持ち梁、片方を3成分固定、片方を集中荷重をあたえた時にそれらしい結果が出たデータを使ってます 要素は8節点6面体の一次要素です。メッシュがぐちゃぐちゃとは要素座標とかそこらの考え方がぐちゃぐちゃしているのか・・・ 剛性マトリクスを作る時についでで質量マトリクスを作っているからそこらが怪しいと思っているんだけど 積分点の位置で形状関数を使って質量を計算してはいけないのかな? 仕事とか学生とかではないから焦ってはいないんだけど、資料が無くて分からない事が多いんですよね 暇な時に少し知恵を貸してくれると幸いです >>182 まずCalculixみたいなフリーのやつ入れて、勉強したほうがいいと思う。 Salome-Mecaが英語だったら、こっちがオススメなんだけど、フランス語なんで日本人には辛い。 3DCADは何かつかってますか? SolidWorksだったら、SolidWorksアドオンのCAE(COSMOS?)を載せたほうがいいと思うんだけどさ。 >>186 レス有難う calculixは一応入れてソースを見てるんだけど・・・コンパイルが悪いのか固有振動数計算が解けないんですよね (*ERROR in arpack: info=-3 # of converged eigenvalues=70 とか出ます) 3DCADはFreeCADを使ってますが、形状がそんな複雑でないのとメッシャーが四面体要素中心なのでテストデータは手で作ってます 使える要素で静解析の結果が一番良かったのが六面体一次要素なのでそれを中心に進めています 使っているプラットフォームがdebian(squeeze 32bit)なのと貧乏なので商用の物は使えないのですが SolidWorksですかちょっと調べてみます >>182 済まない、レスを混同してしまいました>>187 は無かったことに >>187 構造解析なら、メッシュ数おとしてでも、2次要素を使ったほうが絶対いいです。 絶対っていうのは言いすぎかw だけど、自分の経験からいうと、メッシュ数がテキトーでも二次要素だとかなり精度よく解けます。 >>187 二次要素とか気にしなくて良いANSYS DesignSpaceにしとけば? 非線形は必要になったらライセンスアップグレードすれば良い。 固有振動数レベルだとANSYS、NASTRAN、ABAQUSどれも結果変わらないし。 >>189 六面体の2次要素はまだ実装していないのです >>190 ,191 すみません、レスを混同して変なレスをしています テトラは2次でないと全然駄目だけど 6面体は一次で十分です 固有値グチャグチャは気になるけど 固定が駄目なのかなー >>193 固定か〜固定している所の剛体行列と質量行列をどのように扱ったら良いか分からないから固定条件無しでやっています まず剛体モードが綺麗に直行で出てくれないし何か変だ〜 今は固有ベクトルの計算は逆反復法らしきものでやっています (上手く収束しないんだよねと言うかほとんど収束しない、振動してたり反復の上限一万回に引っかかってたりするよ) >>194 手作りだからね・・・まずは一次要素からかなと思って作ってます (ちょこちょこ弄っていて大分汚くなってきたよw) >>193 いやいやいやいや、 6面体だろうが、なんだろうが、2次要素と1次要素じゃ精度が全然違う。 普通の構造解析で積極的に1次要素を使う理由は、全く無い。 ↓ http://www.42ch.net/UploaderSmall/source/1344995786.png >>198 わからない。 学問板なので、理論的に説明してくれ。 >>197 その差は、特に長手方向に粗いのが原因っぽい メッシュ細かくしてきゃ、どっちも値大きくなって、双方一緒っぽい値になる筈 テトラ一次は細かくして行っても、2次と全然違う値になったりするけど >>200 まぁそうなんだけどねw つまり、二次要素は、メッシュ切りがテキトーでもそれなりの精度が出る。 一次はメッシュの細かさが、モロ、精度に出るので、メッシュ密度とか形状に非常に神経を使わなきゃならんので 実務的には向いていないという考えです。 >>201 テトラの2次と1次は 3割4割違ってくるので 比べるとマシだけど ちょっと差があるのは気になるね ウーム >>197 あ、書き忘れていた、これの理論値は -0.75mmです。 つまり 理論値 -0.75 1次6面体(8節点)要素 -0.72 2次6面体((20節点)要素 -0.63 >>203 う、逆だった 理論値 -0.75 2次6面体((20節点)要素 -0.72 1次6面体(8節点)要素 -0.63 >>205 その「倍」っていうのはどこから来る数値なのか? たとえばこういう単純ハリの場合、一次要素で解析するには、 縦横高さで何等分、アスペクト比をいくつにすればいい、 っていう目安ってありますか? 「倍」というのが辺の長さを指しているのであれば、 3次元だから2の3乗で8倍のメッシュを必要とするという理解で良い? メッシュ数が「倍」であれば、辺の長さを2の立方根で割れば良い? 皆遊んでいる様なので質問を・・・ 固有値、固有ベクトルの計算をする時の固定条件の入れ方を教えて頂けないでしょうか? ググっているのですがそれらしい文献を見つけられずに困っています。 詳しく書けと言われそうなので少し書いておきますが長いので面倒なら読み飛ばしてください 固有値、固有ベクトルの計算プログラムを作っているのですが、今までは固有ベクトルがぐちゃぐちゃだったのを 逆反復法で計算することでそれっぽい(と言ってもまだ怪しい)ものが出るようになったのですがoctaveで プログラムしているのでとても遅いです。そこで比較的規模が大きくても良さそうなサブスペース法で組むことに したのですがどうも小行列を作るとそれがコレスキー分解出来ないのです境界条件を設定していなく、剛体モード を含んでいるのが原因ではないかなと思っています。その他にヤコビ法で小行列の固有値、固有ベクトルを求めると サブスペース法の繰り返し計算が収束しない様で困っています。エラーが出ない行列を手作業で与えてやると ビルトインのeigの結果と比較して符号が違う固有ベクトルがでていました。ハウスホルダー+QR法+逆反復法で 出した結果はeigと比べてパターンは違いますが符号が違う固有ベクトルが出ていますがサブスペース法は収束します ここらもエスパーレスしてくれると有難いです テストデータは次の感じです ・方持ち梁 ・六面体一次要素(要素数1 と5これが主 と 625これは時間が足りなくて計算できない) ・拘束条件無し(なんかこれがいけないっぽい) >>208 その通りだと思います。 ただ、ちょっと複雑な形状をきると、アスペクト比1を目指してメッシュを切ると 要素数が非常に多くなることがあるので、1次要素を使ったほうが計算時間が かかってしまう事はあるようです。 >>210 スライドできる状態=摩擦ゼロ状態いうのは、厳密に起こり得ないので、 条件に対し解を求めるいう点では合ってるけど スライドできるいう条件設定が、果たして現実的か? そこはかなり難しい 例えば土木や建築の分野だと 地盤をモデル化したり、剛性体を挟んだり 色々しないと、単品+拘束やすべり条件設定では、合わない >>210 変形の計算(静計算って言うのかな?)はその通りだと思うのだけど、固有値、固有モード計算と同じなんですか? 変形計算の時って [K]{U}={F}で固有値計算の時は[K]{U}=λ[M]{U}ですよね。これをどう節点変位ゼロと 絡めたら良いのか分からないです。頭が悪くて申し訳ない >>213 固有振動の場合は固定点がなくても解が出るはず。 収束計算させるんじゃなくて、ガウス吐き出し法とかで解けば? >>214 境界条件がある固有値、固有ベクトル計算も知りたいので宜しくお願いします。 >215 固有値=フリーでの計算ばかり 固定など条件設定を伴ったものは見た事がない気もします できないのかも知れません >>216 出来ないってこたーないだろ。 AnsysだろうがNastranだろうが、出来ないCAEを見たことねぇぞ そもそも何の研究をしようとしてるのか わけわからん。 こんな固まった分野をいまさら自作してどうすんの? アホ? もしJAXA系や茨城日立でやってんなら 税金の無駄遣いはやめてくれ。 研究じゃないです、ただのサンデープログラマ 解析ソフトって実際何やっているのかブラックボックスで使っているじゃない どんな処理をやっているのか気になるじゃない Calculixのソース見てもさっぱり理解できなかったから小さい所でちまちま作っているわけ えらいね。そういう姿勢が大切なのかも。 オレだったらフリーか個人向けの安いのを買っていじるくらいだな。 宮脇清宮脇茂宮脇勇宮脇博宮脇進宮脇実宮脇弘宮脇正宮脇勝宮脇隆 宮脇豊宮脇昇宮脇誠宮脇明宮脇稔宮脇修宮脇勉宮脇三郎宮脇武宮脇和夫 宮脇一郎宮脇幸雄宮脇孝宮脇一男宮脇功宮脇一夫宮脇勲宮脇守宮脇健一 宮脇保宮脇秀雄宮脇正一宮脇幸男宮脇浩宮脇宏宮脇登宮脇正雄宮脇義雄 宮脇正夫宮脇秀夫宮脇敏夫宮脇光男宮脇昭宮脇徹宮脇文雄宮脇忠宮脇章 宮脇正男宮脇和男宮脇栄宮脇幸一宮脇一雄宮脇満宮脇寛宮脇良一宮脇晃 宮脇栄一宮脇勝美宮脇光雄宮脇薫宮脇剛宮脇正明宮脇正義宮脇健 宮脇幸夫宮脇英雄宮脇繁宮脇信夫宮脇一宮脇敏雄宮脇宮脇章宮脇正男 宮脇和男宮脇栄宮脇幸一宮脇一雄宮脇満宮脇寛宮脇良一宮脇晃宮脇栄一 宮脇勝美宮脇康充宮脇薫宮脇剛宮脇正明宮脇正義宮脇健宮脇幸夫宮脇英雄 宮脇繁宮脇信夫宮脇一宮脇敏雄宮脇仁宮脇健二宮脇和彦宮脇忠雄 宮脇利夫宮脇哲夫宮脇正人宮脇正美宮脇悟宮脇学宮脇武夫宮脇久雄 宮脇正治宮脇洋一宮脇文夫宮脇清一宮脇孝宮脇武雄宮脇誠一宮脇信一 宮脇猛宮脇貞夫宮脇久宮脇文男宮脇英夫宮脇忠男宮脇洋宮脇昭夫 宮脇和雄宮脇哲也宮脇彰宮脇邦夫宮脇俊夫宮脇貢宮脇健治宮脇清一 宮脇孝宮脇武雄宮脇誠一宮脇信一宮脇猛宮脇貞夫宮脇久宮脇文男 宮脇英夫宮脇忠男宮脇洋宮脇昭夫宮脇和雄宮脇哲也宮脇彰宮脇邦夫 宮脇俊夫宮脇貢宮脇健治 やっと書き込みが出来るようになったらしい・・・ 何日か書き込みできる状態ならまた書く〜 北海道札幌兵庫県三木芦屋市 北海道稚内兵庫県立花西宮市 兵庫県明石市千葉県房総兵庫県神戸市 兵庫県小野市 ★マインドコントロールの手法★ ・沢山の人が偏った意見を一貫して支持する 偏った意見でも、集団の中でその意見が信じられていれば、自分の考え方は間違っているのか、等と思わせる手法 ・不利な質問をさせなくしたり、不利な質問には答えない 誰にも質問や反論をさせないことにより、誰もが皆、疑いなど無いんだと信じ込ませる手法 ↑マスコミや、在日カルトのネット工作員がやっていること TVなどが、偏った思想や考え方に染まった人間をよく使ったり、左翼を装った人間にキチガイなフリをさせるのは、視聴者に、自分と違う考え方をする人間が世の中には大勢いるんだなと思わせる効果がある。 .. 大宮市湖南市長野市甲賀市福井市札幌市 強盗公然わいせつ児童ポルノ法放火 皆遊んでいる様なので質問を・・・ 固有値、固有ベクトルの計算をする時の固定条件の入れ方を教えて頂けないでしょうか? ググっているのですがそれらしい文献を見つけられずに困っています。 詳しく書けと言われそうなので少し書いておきますが長いので面倒なら読み飛ばしてください 固有値、固有ベクトルの計算プログラムを作っているのですが、今までは固有ベクトルがぐちゃぐちゃだったのを 逆反復法で計算することでそれっぽい(と言ってもまだ怪しい)ものが出るようになったのですがoctaveで プログラムしているのでとても遅いです。そこで比較的規模が大きくても良さそうなサブスペース法で組むことに したのですがどうも小行列を作るとそれがコレスキー分解出来ないのです境界条件を設定していなく、剛体モード を含んでいるのが原因ではないかなと思っています。その他にヤコビ法で小行列の固有値、固有ベクトルを求めると サブスペース法の繰り返し計算が収束しない様で困っています。エラーが出ない行列を手作業で与えてやると ビルトインのeigの結果と比較して符号が違う固有ベクトルがでていました。ハウスホルダー+QR法+逆反復法で 出した結果はeigと比べてパターンは違いますが符号が違う固有ベクトルが出ていますがサブスペース法は収束します ここらもエスパーレスしてくれると児童買春・児童ポルノ法違反有難いです >>231 剛体モードが出るっていうのは、つまり、逆行列が存在しない っていう意味ですかい? ((1,1) (2,2)) みたいな。 KinKi Kids・堂本剛が約4年ぶりにドラマに主演することが3日、わかった。 7月スタートのTBS・ドラマNEO枠『天魔さんがゆく』(毎週月曜 深0: 28) で、『33分探偵』シリーズ(フジテレビ系)で組んだ福田雄一監督と全く新しいホラーコメディーに挑む。 2009年3月〜4月放送の『帰ってこさせられた33分探偵』(全4話、フジテ レビ系)を最後にドラマから遠ざかっていた剛。今回、「親友の福田監督か らお声がけいただき、また面白い物が作れる、どんなものが飛び出してくる のかとワクワクしています。福田雄一と堂本剛のワールドをぜひ楽しみにし てください」と張り切っている。 同ドラマは、『コドモ警察』のヒットの記憶も新しい福田監督の脚本・ 演出による完全オリジナルドラマ。剛演じる主人公・天魔は、先祖代々 幽霊とコンタクトを取る能力を持っているが、その力に反して大の怖がりで、 「幽霊退治引き受けます」と掲げた会社の2代目という設定。米映画 『ゴースト・バスターズ』(1984年)のようだが、退治する意気地がない 天魔は、幽霊の言い分を聞き、この世に残している悔いを晴らしてあげたり、 優しく諭してあげたりしながらさまざまな幽霊退治の依頼を解決していく。 天魔の会社の女性社員に川口春奈、先代からの番頭でもある皆川猿時、 天魔を追いかけ回す近所のお姉さんに芹那、天魔と夢の中で対話する亡き 父親に佐藤二朗容疑者、天魔たちの周りを嗅ぎまわる刑事に森崎博之などが出演。 依頼人役や幽霊役のゲストにも注目だ。 patranでのsolidの四面体メッシュ作成について質問させてください. できるだけ均質な(1要素の体積や辺の長さなどが近い)メッシュを作成したいのですが どのようにすればよいでしょうか? 現在はTetmshでの自動生成を用いていますが,酷すぎて困っています. よろしくお願いいたします. お客が見える飛行機から近い場所でJALの社員数名がふざけて殴り合っていました 本人たちがじゃれあっているつもりだったのでしょうが飛行機を待っていたお客さんは JALの服を着た人たちが乱闘している喧嘩していると勘違いした それをyoutubeにのってしまってJALは大恥 JALグランドサービスはそんな会社です 俺がJALグランドサービスで働いて10ヶ月ぐらいバイトで経った頃 俺たちの課長が代わって千代っていう馬鹿課長になったんさ 前の課長はまだ部下 したの人間の言い分も聞いて器のでかさがあったけど その後きた千代って言うのがとんでもなくアホ 俺が周りの係長とかに苛められているという事実も今まで周りから不愉快な言いがかりや 暴言などを耐えていたという状態も知らず俺を悪い奴と決め付け絶対俺をやめさせると いい 毎日のように粘着し ついには俺を会社から追い出すという暴挙にでた 「お前みたいな使えない奴を長くいさせるよりかは訴えられて賠償金を払ったほうが安くすむ」という暴言も吐かれた それから「俺は今まで自分が使えないと思った奴を辞めさせてきた」と俺に豪語したこともあった 中国系列の飛行機を掃除していたとき 班長 その他数名が仕事中喋ったりふざけて遊んでいました 私はラバー(トイレ)を夢中で掃除をしていていました そのときトントン トンと私の肩を叩く人がいたので振り向いてみるとその飛行機のクルーの一人が 目を吊り上げて私の前に立っていました そして私に「あの人たち 口を動かしてばかりで仕事しない 注意してきてくれ」と私にいいました 私は急いで係長にそのことを報告しにいきました 班長はその横で口を緩ませてヘラヘラしているだけでした 他社からあいつら喋ってばかりで仕事をしないといわれても恥じもなく 適当な仕事を他社に見られても平気な JALグランドサービス このことでその飛行機会社からは信用を失い二度とJAL系列に仕事を任せないと 決断してくれることを切に願うばかりである 車の中で機内のお菓子を食べたりしていたよ あと係長が機内においてあった古新聞を一人で別の仕事をするとき車の中に持ちこんでいた これ二つともやってはいけないこと 飛行機の中の物は外に勝手に持ち出したらいけないのにもしかして係長になったら許されるとおもってんのかな? フィロソフィとかつくってもなにも役に立たないね こんなことばかりやってるからJALグランドサービスは屑会社とか馬鹿会社とか言われるんだよなw 仕事中に酒臭かったりする連中もいるよw それから仕事中ガム噛んでいたりw 正社員の質の悪さっていったらそりゃもう 機内のものは持ち出し禁止なんだけどさ 新聞を持ち出しで車の中で読んでいるクズ係長もいるしw なんでもありw あいつらの都合の悪いことは隠蔽したがるけど 俺はあのクズ会社辞めたからカキコするわw つまり俺はJALグランドサービスの敵になったわけw 密告しまくりよw JALグランドサービスは貧乏性でテーブル拭いて汚れたタオルを次の清掃場所 さらに次の清掃場所にもっていって汚れたタオルを使えっていってきた タオルからはなんともいえない臭いが漂い色もケチャップやなんやらがこびりついていて汚かった お客さまのためにを毎日言っている会社なのになにこれ? 矛盾してるよね? 青タオルはウンコがこびりついていて汚い それをキレイなタオルといいラバーのほかギャレー 床のゲロ掃除などさまざまなところで活用している バケツにそのタオルを濡らすとき気持ちが悪い バケツもすごく汚いよね バケツ洗っているところみたことがないし JALはこんな会社です 白タオルは生乾きで蒸れた悪臭がしている それでお客様がつかうテーブルを拭いているんだから あのテーブルは菌をひろげているようなもんだ テーブル拭く白タオル悪臭が酷い ラバー(トイレの便器)やトイレの床を拭く尿やウンコつきの青タオルを洗う バケツといっしょにいれてるから臭くなるんやw それでお客様がつかうテーブルを拭いているのだからJALの衛生面は最悪 こいつらは所詮お客様の姿を金としかおもっていないからそういう掃除しかできないんや 口ではお客様のためとか言ってるけどなw JALグランドサービスはパートやアルバイト、契約社員に対し若い正社員ですら横柄かつ高圧的な言動(恫喝行為など)を取ることが常態化されていて それを正せる上司が一人もいないというのが特徴の会社 具体的に証拠がないのに書き込んでると 業務妨害で訴えられちゃうよ。 calculixに軽く触れてみたけど 解析実行して.datファイルがきちんと生成されるのに(内容は変位と応力とひずみだった) cgxのコマンドに「data set」がなくて結果が読み込めないっぽい? >>249 コマンドというかcgx作業画面の空白をクリックして出てくる機能の中にdatasetがない 結果データはdatファイルを読み込むわけじゃなくてfbdファイル内にそのまま入ってるのかなと思ってfdbをメモ帳で開いてみたけど結果は入ってない ファイルの書式があればいいんだけどマニュアルにそこまで書いてあるのかな←いまここ やっぱり導入時が一番苦労するなあ 打ったコマンド書けがボケが!と言いたかったわけです。 fdbと書いている時点で怪しいかも(大分弄っていないので記憶が曖昧だけどfdbは違うんじゃないのかな) もう見つけているかもしれないけど「有限要素法のフリーのソフトウエア CalculiX を試す」でググると良いかも うちはLinux環境だからWindows特有の問題に巻き込まれているなら分からんよ ANSYSの伝熱での接触定義ってどうやるの?どうやっても断熱しちゃう >>253 Workbenchですか? Classic(APDL)ですか? 基本の講習マニュアルに載ってるレベルだろうに。 ちゃんと契約してるベンダーに聞けよ。 >>255 いやいや、ある面に熱伝達係数を設定した上で、その面の一部に接触面がある時の取扱いについては、 マニュアルに載っていないので、Quick-Aレベルだな。 >>257 コンタクトウィザードの該当面の接触プロパティの Thermal Contact Conductance が、デフォルトだと0になって断熱しているので、 この値をテキトーに入れてみてください >>258 モデリング→作成→接触ペアのところでいいですか? あと相互に熱伝導をしている状態を再現したいのですがたとえば物体A、Bが接触してるとして コンタクトAターゲットBとコンタクトBターゲットAのように二回接触定義をすればいいんでしょうか? >>259 2回やっちゃダメです。 互いに熱伝導っていうのはアタリマエで、片方に熱伝導するなんて非現実的な設定は、元々設定出来まセン。 >>260 了解しました。そうなるとThermal Contact Conductanceのテキトーな値というのはどう設定すべきなんでしょうか? 両者の熱伝導率の平均でいいんですかね? >>261 違います。 これ以上の事はノウハウなので、サポートに聞いてください。 256 僕は働かないくせにお金ほしいクズ野郎です。 288 ≫256 大丈夫。借金まみれなのにシステムトレードで不労所得目指してるクズ野郎いるよ。 http://manta.blog.jp/ 建築で使うマトリックス法による骨組み応力解析は 有限要素法っていうのでしょうか? 今日はムラタのfemtetがユーザー会だから話をきいてくるお ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。 ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。 ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。 ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。 ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。 ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。 40万(20万)という安さを考えれば、まあまあどころじゃないなぁw 使う方としては最大手の一級品と比較してああしろこうしろと言ってしまうが。 パナのCAEの達人みたいな人がキレキレのトークで絶賛しながら滅多切りしてた。 ああいう人が上手く叱咤激励してくれるとソフトも育つ。 ansysもユーザー事例発表で東芝の人が不満点をぶちまけて早く直せと発破をかけてたな。今年のユーザー会は明日だが、さて 30代後半で漫画家になろうとしている童貞ひきこもりなバカを発見 2chねらーに荒らされている。足立区に住んでいるそうだ。 http://inumenken.blog.jp/archives/7002197.html 30代後半で漫画家になろうとしている童貞ひきこもりなバカを発見 2chねらーに荒らされている。足立区に住んでいるそうだ。 http://inumenken.blog.jp/archives/7002197.html ansysスレに書いたんだけどこっちのスレでも質問。 LS dynaでS字のパイプに異方性材料を適用したいのよ。 内圧がかかる方向には高いヤング率を設定して曲げ方向には低いヤング率にする。 普通の異方性材料のカードでやるとひとつの座標系でしかベクトルを設定できないからもとから曲がりくねってるパイプにはうまく設定できない。 分かりにくくてごめん ボディを分割して部位ごとに違う材料として割り当てりゃいいだろ。 >>277 それだと座標系を山ほど作らないといけないのでは? どのぐらい方向が違う部位があるのかわからないけど それでいいんじゃね? 通常と違うことをするのに苦労をしたがらないのに変だと思うぞ ★2ch勢いランキングサイトリスト★ ◎ +ニュース ・ 2NN ・ 2chTimes ◎ +ニュース新着 ・ 2NN新着 ・ Headline BBY ◎ +ニュース他 ・ Desktop2ch ・ 記者別一覧 ◎ 全板 ・ 全板縦断勢いランキング ・ スレッドランキング総合ランキング ◎ 実況込み ・ 2勢 ・ READ2CH ・ i-ikioi ※ 要サイト名検索 ___ _ ヽo,´-'─ 、 ♪ r, "~~~~"ヽ i. ,'ノレノレ!レ〉 ☆ 日本のカクブソウは絶対に必須です ☆ __ '!从.゚ ヮ゚ノル 総務省の『憲法改正国民投票法』のURLです。 ゝン〈(つY_i(つ http://www.soumu.go.jp/senkyo/kokumin_touhyou/index.html `,.く,§_,_,ゝ, ~i_ンイノ 2014年のアンシスカンファレンスは以下で行ないます。 http://www.youtube.com/watch?v=HJKHJEJgxgo ユーザーの方もユーザーでない方もぜひご参加お待ちしております。 DelFETをフリーのwin系Cでビルドできるようにして配布すれば裾野が広がりそうだが。 Nobuyuki Umetani http://www.youtube.com/channel/UCqObO-PlZgCcz0z54GsE1YA 尋ねたいことがあります 今Abaqusで荷重の代わりに強制変位を与えて板の弾性座屈解析を行っているんだけれども この場合Abaqusが出力する固有値ってどんな意味があるんでしょう? 調べてもなかなか答えが見付からず・・・ 水溶液中のイオンの移動などを記述するネルンスト・プランク式の有限要素解析 についての調べる場合には、どんな論文集・本等を勉強すればよいか、どなたか教えていただけませんか。 >>292 遅くなりましたがありがとうございます これを読みつつなんとか解決することができました /:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::ヽ /:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::://ヽ:::::::::::::::| l:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::// ヽ::::::::::::::l l:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::/:::「'ヽ:::::::::::// ヽ:::::::::::| |::::::::::::::::::::::::::::::::::::::::::::::::::::::::::ノl:::ノ l:::::::/ ヽ::::::::| ノ:::::::::::::::::::::::::::::::::::::::::::::::::::::/ ゙゙ ノ:::/ ,,;;;;;;,, ,,,,ヽ:::::l ):::::::::::::::::::::::::::::::::::::::::::::::/ ノ/ __,'''i: ('''__):::l )::::::::::::::::::::::::::::::::::::::::::::::::::/  ̄ ̄ン:. :「 ̄`ヾ 1:::::::::::::::::::::::「 `┤l:::::::::::::::::l  ̄ , ヽ ̄ l `l:::::::::::::::::::::ヽ :l li:::::::::::::/ ヽ /´ `l | ヽ::::::::::::::::::::::\_」 lヽ::::/ .l !:-●,__ ノ / ノ:::::::::::::::::::::::::::ノ | l `゙゙ i ,,;;;;;;;;;;;;;;;;;;;;, /ヽ ,/ ヽ::::::::::::::::::::::( l l::::::::.. /.:''/´ ̄_ソ / `ヽ ヽ:::::::::::::::ヽ | l:::::::::::... /::// ̄ ̄_ソ / \ ヴッ!! ヽ:::::::\| l::::::::::::::::... / :::.ゝ` ̄ ̄/ / ヽ ヽ:::l l:::::::::::::::::::..  ̄ ̄;;'' / ヽ l l;;;;;;:::::::::::::::.....;;;;............;;;;;;''ノ l l l '''''''''''''''''''''''''''''''''''''' ̄l | | http://www.youtube.com/watch?v=z2qK2lhk9O0 ☆☆☆☆☆ ☆ 自民党、グッジョブですわ。 ☆ http://www.soumu.go.jp/senkyo/kokumin_touhyou/index.html ☆ 日本国民の皆様方、2016年7月の『第24回 参議院選挙』で、改憲の参議院議員が 3分の2以上を超えると日本国憲法の改正です。皆様方、必ず投票に自ら足を運んでください。 そして、私たちの日本国憲法を絶対に改正しましょう。☆ ρ(∂u/∂t+ u∂u/∂x+ v∂u/∂y) = ?∂p/∂x + μ∇2u^2 のxをt=0から順に計算していくにはどうすればいいの? ttp://fluid.nuae.nagoya-u.ac.jp/archives/textbook/past/nensei/Nensei2001Chapt2Print.pdf まずは有限要素法の解法の基礎を構造解析の差分法あたりからがんばればいいよ 有限要素法で定式化から解いてる本って無いですかね. 一からプログラム組んでる勢なんだけど死にたい. うおおおおお!トランザム起動!! All together now〜 All together now〜 みんなで人間革命(Yeah〜!) http://goo.gl/5YOYqX http://goo.gl/7Sastj http://goo.gl/1ondZM 後藤達俊 バックドロップ43連発 http://goo.gl/HkEPgn マ イ ン ド コ ン ト ロ ー ル の手法 ・沢山の人が、偏った意見を一貫して支持する 偏った意見でも、集団の中でその意見が信じられていれば、自分の考え方は間違っているのか、等と思わせる手法 ・不利な質問をさせなくしたり、不利な質問には答えない、スルーする 誰にも質問や反論をさせないことにより、誰もが皆、疑いなど無いんだと信じ込ませる手法 偏った思想や考え方に染まっていたり、常識が通じない人間は、頭が悪いフリをしているカルト工作員の可能性が高い 靖 国 参 拝、皇 族、国 旗 国 歌、神 社 神 道を嫌う カ ル ト 10人に一人は カ ル ト か 外 国 人 「ガ ス ラ イ テ ィ ン グ」 で 検 索 を ! アケゲ速報 募金詐欺 スパムブログ 在日 無断転載 中傷 児童買春 民主党支持 salome-meca が結構色々できるんだが、ユーザーはおらんかね やはり、理系ならC言語で有限要素法を作らないと。 有限要素法って、いわゆる変分法のことだよ。 いろいろな解析対象を細かくメッシュで区切って、力学的に解析する手法。 一度、Rubyで作ってみたいなぁ。 このスレに書き込んだ人には、有限要素法のことがよくわかる。 長岡技術科学大学の学部でも、有限要素法は教えているよ! salome-mecaっつーかcode_asterって最近アツいのか? プレミア見れない ブンデス見れない CLEL見れない 音楽聞けねえちきしょう 同和のクズ共死ねクソ共がざまあみろ気違い共 ほれ気違い共もっともっとドア閉めろ通れ それしか能のない能無し共がざまあみろ地獄に落ちろ 悔しいか、ざまあみろくたばれクソ同和 お前らの恐ろしさをもっと見せてみろ。そんなんじゃなんともねえぞ 袋とじ見たぞ 悔しいか、ざまあみろくたばれクソ同和 嫌がらせがエスカレートするのが楽しみでしょうがない 今それだけが楽しみだ 俺の生き甲斐藁 それだけ怒ってるってことだもんな藁 >>324 まぁ、それでもいいんだけどねぇ… 下手糞が組んだら何ができるものやら… (どっかのFortran90/95による有限要素法うんたら(丸善)を見ながら) >>303 昔の矢川宮崎本とか平井本とかのほうが優しいかもしれん。 平井本は入手不可能だが。 出来れば分野か、解く方程式書いたほうがわかりやすいよ。 昔2Dの弾性体を解くプログラムを作ったな こういうのを仕事にしたかったのかもしれない ばねを伸ばそうとすると、縮めようとする力がはたらきます。 縮めようとすると伸ばそうとする力がはたらきます。 このような元に戻ろうとする力を弾性力といいます。 作用⋅反作用の法則にのっとった反作用の力です。 ばねの弾性力には大きな特徴があり、それは、力の大きさが目に見えるということです。 ばねの伸び、あるいは縮みとして目に見える形で現れます。物理では珍しいことです。 伸び縮みの長さは弾性力の大きさにきれいに比例します。 (といいますか、高校物理ではそのようなばねについて考えます)。 ばね秤ばかりは、ばねの伸びの長さによって物の重さを表示する装置です。 ばねの伸び縮みの長さを x [m] 、弾性力の大きさを F [N] とすると、 F = k x という関係があります。フックの法則といいます。 k はばね定数と呼ばれる比例定数です*。ばねの伸びにくさを表します。グラフの傾きです。k が大きいと、伸ばすのに、大きい F が必要になります。 何をいまさらと言われるかもしれないが、まずは何よりも、フックの法則だ。 すべては、ここから始まる。 断面力や応力は、目に見えない内力だ しかし、バネの弾性力は目に見える だから、物理の基礎を勉強するにはバネがいい ばねを直列につなげると、伸ばしやすくなるでしょうか、伸ばしにくくなるでしょうか。 同じばねを3つ用意して、ばねが1つだけの場合と、ばねを2つ直列につなげた場合とで比較してみます。 どちらも20cm伸ばすとすると、フックの法則から、1つの場合の弾性力は F1=0.2×k [N] 。2つの場合は各10cmだけ伸びることになるから1つ当たりの弾性力は F2=0.1×k [N] 。 そうしますと、2つつなげた場合の弾性力は2倍の 2×F2=0.2×k [N] でしょうか? 直列接続のばねを伸ばしたときには各部分にまったく同じ力がはたらいています。 途中が F2[N] ならどこもかしこも F2[N] です。 ばねを伸ばして静止した状態というのは力がつり合った状態です。 ばねの各微小部分同士が同じ力で引っ張り合ってるので静止しているのです。 ミクロな視点でいえば、ばねを構成する原子たちがお互いを F2[N] で引っ張り合ってつり合って静止しているのです。 端っこでも途中でもどの部分においても各微小部分同士は同じ力で引っ張り合ってつり合って静止しています。 というわけで2つつなげた場合の弾性力は 2×F2[N] ではなくて F2=0.1×k [N] です。 ばねが1つのときの F1=0.2×k [N] の半分です。 2つつなげた方が伸ばす力がいらない、すなわち伸ばしやすくなるのです。 このことは、ばねというものは、同じ材質、同じ形状なら長い方が伸ばしやすい、つまりばね定数が小さいということを意味します。 ばねを並列につなげた場合はどうでしょうか。 重い荷物を持つ場合は2人で持てば半分の力で済みます。 2人対1人で綱引きをすれば2人の方が勝ちます。 つまり同じばねを2つ並列につなげると、1つのときに比べて倍の力が必要です。 ばねは伸びにくくなるのです。 2つで1つのばねと見なしたときのばね定数は2倍になります。 並列の場合は単純に足すだけです。 ばねは伸びにくくなります。 これはばねが太くなったことにたとえられます。 ばねは太い方が伸びにくいです。 フックの法則はバネだけに成り立つ法則ではありません。 身の回りにある物体すべてに成り立ちます。 別な言い方をすれば、「す べてのものはバネである」と言っていいのです。 考えてみるとこれはすごいことです。 非常に幅広く成り立つ普遍的な法則といえるわけです。 なぜそんなに普遍的に成り立つのかというと、物質の構成要素である原子や分子の間に働く力に「フックの法則」が成り立っているから と考えられます。 有限要素法の構造解析では、基本的にKu=fというバネの公式を用います(ここで、Kはバネ定数、uは変位、fは荷重)。 ただし、その中身はマトリクスの形式で表されていて一見複雑ですが、やっていることはバネの公式と何ら変わりがありません。 [K]{u}={f} これが、有限要素法で最終的に解く式になります。 変な括弧が付いている他はバネの公式と同一です。 [ ]はマトリクスであることを表し、{ }はベクトルであることを表しています。 [K]は剛性マトリクスと呼びます。 要素一つについての関係を表す場合には要素剛性マトリクス、解析モデル全体について重ね合わせたものを全体剛性マトリクスと呼びます。 {u}は変位ベクトルと呼び、それぞれの節点の変位になります。 これは拘束条件の入力にも用いますが、解析で求める値そのものになり、コンピュータはこれを求めるためにCPU時間のほとんどを使用します。 {f}は荷重ベクトルで荷重条件を設定した場合、各節点に対応した荷重値がここに入力されます。 この辺の詳しいところは境界条件の設定の項で詳しく説明します。 つまり、フックの法則の式を、マトリクスとベクトルを使って表記したもの、これが有限要素法の基本式 [K]{u}={f} 何を求めたいかといえば、変位ベクトル{u}を計算して求めたい。 [K]剛性マトリクス {f}は荷重ベクトル この2つが分かれば、変位ベクトル{u}も出てくる バネを引っ張るとか押すと長さが伸びたり縮みます。 伸び縮みの大きさはバネ定数に反比例し、バネに加えた荷重に比例します。 身の回りのモノはすべて、バネ、つまり弾性体と考えることができます。 変形は外力の大きさに比例し物体の剛性に反比例します。 物体の変形が外力および剛性と比例関係にある場合を線形(linear)といいます。 線形的な静的挙動を表す物体を有限要素法で表すと、[K]{u}={F}という行列方程式の問題に変換できます。 つまり、平たく言えば、剛性マトリクスは変形に対する物体の強さ 荷重ベクトルは物体にかかる荷重の大きさ そして、求めたい変位ベクトルは、変形の大きさということになる このうち、荷重ベクトルは物体にかかる荷重の大きさなのだから、これは設定次第。 境界条件の設定で決まる。 あとは、剛性マトリクスが決まればいい。 変位がわかれば、ひずみもわかる ひずみが分かれば、応力も出る だから、変位を求めることが第一目標 まずは、なんといっても、フックの法則 それが基本 中学の理科のレベルだけど、ここが一番大事 [K]{u}={f} ↑これが、フックの法則の行列表現 行列とベクトルが使われるからには、一応、高校の数学のレベルだ この式をもとにして、変位ベクトル{u}を求めるのが第一目標 節点の変位を求めたい そのために必要なのは、節点の剛性マトリクス フックの法則までは、中学生でもわかる それがなんで行列になるのかについて、中学生に説明するのは難しい それには、まず、「節点の自由度」という話をしなければならない 節点に生じる変位の方向には、水平・垂直・回転の三つがある つまり、1つの節点につき、3つの自由度がある。 でも、だからといって、節点が5つあれば、5×3で15の節点があるかっていったら、そういうわけでもない たとえば、ピン支点では、鉛直方向と水平方向に固定され、回転だけ自由になっている つまり、節点がピン支点になっている場合は、自由度が1つしかない こんな具合にして、自由度の数を数える ここでいう自由度の数、それこそが、求めるべき未知数の数ということになる でもって、自由度の数がたとえば8個だとすると、どうなるか? [K]{u}={f} のベクトル {u} と {f} のカッコ内にある要素数が、それぞれ8個になる {u}t = {u1 u2 u3 u4 u5 u6 u7 u8} {f}t = {f1 f2 f3 f4 f5 f6 f7 f8} ・・・というわけで、超初歩のお話ですた (つづく) >>358 この場合、剛性マトリックス[K]は、8行8列の行列になる。 行列を中学生にもわかるように説明するなら、 「多くの数値をまとめて書く、便利な表記法」 といったところか。 荷重ベクトル {f} = {f1 f2 f3 f4 f5 f6 f7 f8} ここで、荷重ベクトルを設定する。 2番目のところ、f2で、荷重Pがかかっているとする。 それ以外のところ、f1および、f3〜f8では、荷重がかかっていない。 そうすると、荷重ベクトルはどうなるか。 荷重ベクトル {f} = {0 P 0 0 0 0 0 0} ということになる。こんな風に、行列やベクトルの中に0が多くないと、有限要素法らしくない(笑) (転置記号は省略) コピペ 変位法とは「変位をもとめる方法」です。これはいたって明快です。 いや、そうはいっても「応力」だって出ますよね、という反論があるかもしれません。 しかし、変位が分かれば応力は自動的にもとまることになってますから、 応力は変位の後についてくるオマケみたいなものです。 変位法の実質的な計算は「変位をもとめる」ことに費やされている、 ということをまず頭に入れておいてください。 で、その「変位」とはなんのことかというと、これは「節点の変位」です。したがって、変位法で取り上げるべき最初の話題は「その節点とはいったい何なのか?」になります。 デジタル大辞泉の解説 せっ‐てん【節点】 1 構造物の骨組み部材の接合点。自由に回転する滑節点、角度の変わらない剛節点がある。 やっぱり、「骨組みの接合点」・・・これに勝る節点理解はないだろう。 機械工学ならまだしも、土木工学でやる「地盤解析」では、節点の意義が本当にわかりづらい。 単に、人為的にメッシュを切った、縦線と横線の交点という感じしかしない。 まあ、節点とはそういうものなんだけど。 でも、最初の入口段階では、骨組みの接合点で考えたい。 ノード(node) 《植物の節の意》 1 活動・組織などの中心点。集合点。 2 グラフ理論の点と線分からなる図形における、点(頂点、節点)。 路線図の駅、電気回路の素子、ワールドワイドウェブのウェブページに相当する。 3 コンピューターネットワークを構成する個々の要素または中継点。 コンピューター、ハブ、ルーター、さらにはネットワークに接続されたプリンターなどが相当する。 ひずみは、変形(形が変わる)と同じ意味です。 構造物の形が元の状態から変わらなければ、ひずみは生じません。 逆に、形が元の状態から変わると、必ずひずみが生じます。 例として、x-y面内にある正方形板を想像してみましょう。 この正方形板が、力を受けて、xy面内で移動したとします。 すると、正方形板内の各点には、すべて変位が生じます。 もうひとつ、正方形板が、力を受けて、z軸回りに回転したとしましょう。 すると、正方形板内の各点には、すべてz軸まわりの回転が生じるとともに、xy面内の変位も生じます。 以上の、変位したり回転したりした正方形板を、元の位置や向きに戻して、元の形状に重ねて見た時、ぴったりと重なれば、その正方形板には、ひずみは生じていません。要するに、”変形”していないのです。 もし、正方形板が力を受けて菱形状などに変形したとすれば、元の形状とはぴったりと重なりません。要するに、その正方形板には、ひずみ(=変形)が生じています。 要するに、ひずみとは、ある構造物の形状が、元の形から変わった、その変化量を指すのです。 力を受けて移動した後の形状が、移動前の形状とぴったり重ならなければ、ひずみが生じていることになるわけです。 ある構造物にひずみ(=変形)が生じていれば、その構造物には必ず変位も生じています。(回転は生じているとは限りません。) ある構造物に変位や回転が生じているからといって、その構造物にひずみ(=変形)が生じているとは限りません。 フランス語では”ひずみ”に対応する言葉が、日本語の”変形”に対応するdeformationとなっていて、混乱が起きにくくなっています。ドイツ語など他の欧州言語でも同様です。 ひずみは、形状が元の状態から変わった”変化量”を指すので、定義は、ある2点間の距離の変化量で表すことができます。 要するに、数式で表せば、 ひずみ=(移動後の距離−移動前の距離)/移動前の距離 などとなります。 ひずみの定義の仕方によって、分子や分母が微妙に変わってきますが、 ひずみ=長さ/長さ になることには変わりがありません。 その結果、ひずみは無次元の量になります。 無次元の量ということは、呼び名がないということにもなって、これでは不便です。 実は、角度や上記の回転も同様で、SI単位では無次元量です。しかし、これでは角度の数値を呼ぶときに不便なので、”ラジアン”などという便宜上の単位を設定しているのです。 ひずみでも、同様の工夫が行われています。 最も一般的なのが、ひずみの大きさが通常は1よりも非常に小さいことに着目して 0.001 なら、0.1% など、%で表す流儀です。 もうひとつ、主に測定関係者の間で使用されている流儀があります。 金属材料の通常の使用状態でのひずみは、10^-6のオーダーで表現すると適当な大きさで表すことができるようになることから、10^-6を”マイクロストレイン”などと呼んで、 0.001 なら、1000マイクロ(ストレイン) などと呼ぶのです。 ついでながら、もうひとつ例を。 時計の針ですが、時々刻々と回転を続けます。 要するに、時計の針は、回転と変位を生じ続けています。 言いかえれば、時計の針は、向きと位置が変わり続けます。 しかし、時計の針の形自身は変わるわけではありません。 要するに、時計の針は、時刻の変化によって変位(と回転)は生じますが、ひずみは生じません。 家電の事ならGood Price 2018にお任せください! お求めやすい価格にて販売しておりますので、是非ご覧ください! 店舗ホームページ https://goodprice2018.shop-pro.jp/ 👀 Rock54: Caution(BBR-MD5:fc5433912aa55592f73f2dda4d43bdf8) FreeFEM++ があれば、もう有限要素法のプログラムをFortranやC言語で書く必要は なくなるの? ☆★☆【神よこの者たちはもはや人間ではない悪魔であるこのような悪魔どもを一匹残らず殺してくださいお願いします】★☆★ 《超悪質!盗聴盗撮・つきまとい嫌がらせ犯罪首謀者の実名と住所/死ねっ!! 悪魔井口・千明っ!!》 【要注意!! 盗聴盗撮・つきまとい嫌がらせ犯罪工作員】 ◎井口・千明(東京都葛飾区青戸6−23−16) ※盗聴盗撮・嫌がらせつきまとい犯罪者のリーダー的存在/犯罪組織の一員で様々な犯罪行為に手を染めている 低学歴で醜いほどの学歴コンプレックスの塊/超変態で食糞愛好家である/醜悪で不気味な顔つきが特徴的である 【超悪質!盗聴盗撮・嫌がらせつきまとい犯罪者の実名と住所/井口・千明の子分たち】 @宇野壽倫(東京都葛飾区青戸6−23−21ハイツニュー青戸202) ※宇野壽倫は過去に生活保護を不正に受給していた犯罪者です/どんどん警察や役所に通報・密告してやってください A色川高志(東京都葛飾区青戸6−23−21ハイツニュー青戸103) ※色川高志は現在まさに、生活保護を不正に受給している犯罪者です/どんどん警察や役所に通報・密告してやってください 【通報先】 ◎葛飾区福祉事務所(西生活課) 〒124−8555 東京都葛飾区立石5−13−1 рO3−3695−1111 B清水(東京都葛飾区青戸6−23−19) ※低学歴脱糞老女:清水婆婆 ☆☆低学歴脱糞老女・清水婆婆は高学歴家系を一方的に憎悪している☆☆ 清水婆婆はコンプレックスの塊でとにかく底意地が悪い/醜悪な形相で嫌がらせを楽しんでいるまさに悪魔のような老婆である C高添・沼田(東京都葛飾区青戸6−26−6) ※犯罪首謀者井口・千明の子分/いつも逆らえずに言いなりになっている金魚のフン/親子孫一族そろって低能 D高橋(東京都葛飾区青戸6−23−23) E長木義明(東京都葛飾区青戸6−23−20) F若林豆腐店店主(東京都葛飾区青戸2−9−14) G肉の津南青戸店店主(東京都葛飾区青戸6−35ー2 ( https://jdsc.ai/news/ この国は変えられる AIの活用 JDSC/4418テンバガー候補) https://ja.wikipedia.org/wiki/KataGo 世界で最も強い囲碁ソフト DeepMindの論文に加え、独自の研究、強化学習の速度を向上(50倍以上)=AI上の工夫は無限 https://arxiv.org/pdf/1902.10565.pdf RPAロボパット年間24000時間の業務を削減 月の作業時間が5分の1に https://fce-pat.co.jp/case/ 3分でわかるRPAとは何か?注目される背景~今後の進化まで パーソル https://www.persol-pt.co.jp/persolrpa/rpalounge/column1/ https://www.excite.co.jp/news/article/Prtimes_2022-05-31-81631-8/ 設計・材料選定プロセスにおけるAI活用の最新動向レポートを発表。データドリブンな材料開発の裏側と、設計・材料選定プロセスにおける工程時間を55%削減するヒントについて解説 PR TIMES 0531 株式会社SUPWAT 資料DL:https://note.com/supwat/n/nb3a9ec3b318c https://spioenkopjp.blogspot.com/2022/05/blog-post_18.html Oryx Blog - ジャパン 0518 未来戦に備えよ:トルコが無人機による空戦技術の礎を築くための手法 トルコは「アクンジュ」や「クズルエルマ」といった無人戦闘機の開発に加え、いつかそれらの後継機を設計したり、先端技術を特徴とするその他の分野において働くであろう優秀な人材の確保にも入念に注意を払っています。 この国は、世界でも類を見ない規模で、子どもたちや若者の間でテクノロジー分野のあらゆるものに対する関心を高めることを通じて、その目標を達成することを試みています。 https://xtech.nikkei.com/atcl/nxt/column/18/02108/062700001/ AIは万能ではない、機械学習でできることとその限界を知る 0711 澤田千代子 日経xtech https://wired.jp/article/to-win-the-next-war-the-pentagon-needs-nerds/ AIによる戦争が現実化する時代 米軍での「高度IT人材」の不足が深刻化 WIRED 0531 (この国は変えられる AIの活用 JDSC) https://jdsc.ai/news/ ゆこゆこ JDSCのAIを導入CVRが平均2.4倍 常石造船 データサイエンティスト育成プログラム Robo Co-opで業務提携RPA×AI 弾塑性計算作ってて、修正ニュートンラプソん使ってるけど、塑性化したら一向に収束せん。詳しいアルゴリズム書いてある書籍とかない? > 彼らはよく、社会に貢献したいと口にする。 > なんでも社会悪のネトウヨを自殺に追い込むことが、社会に貢献することなんだそうで。 > イジメや嫌がらせで社会に貢献できる教師や警官になるために、あえて帰化したんであって、祖国同胞を裏切ったわけではなく、心は●●人なんだそうだ。 > > 昔は帰化すると裏切り者と呼ばれたりしたが、祖国に国籍を残したまま帰化する方法が確立された現在では、社会に貢献するためにむしろ帰化することが推奨されている。 > 拳銃所持で前科のある生粋の反日家ですら、今では普通に帰化している。 > > ●●学会などはネトウヨ認定した日本人を盗撮して、痴漢の写真だと言ってばらまいている。 > それらの写真は、集団ストーカーに使用される。 > 彼らは集団ストーカーを、[地域で子供を守る安心安全パトロール]と称している。 >>385 非線形有限要素法のための連続体力学 by Neto et al. 訳: 寺田 オレの前を歩いてた小学生の女の子が、後ろを振り向いてオレを見た瞬間に、防犯ブザーを押した read.cgi ver 07.5.5 2024/06/08 Walang Kapalit ★ | Donguri System Team 5ちゃんねる