なんでそういう商用アプリを使うかな。
FreeMATなりOctaveを使えよ！
物理が、分からないと、何やっても、ダメだって！！！
差分法！！！！！
それって、windows上で動きますか？
ttp://http://freemat.sourceforge.net/ Linux, MacOS Windows で動くらしいですね。
誰かわかりませんか？涙
とか
FEM
なんかGPLでソース配布してるプロジェクトが
少なからずある
普通にインストールできたけど…？
当方の環境はWindows XP + Abaqus Student Editon 6.9（書籍の付録）
クソガキの、巣窟！！！！！
フリーっていうか、1万円の本の付録扱い
本に付録してるだけなのか？
広める気あるのかないのか・・・・ｗ
個人でAbacus使おうなんていうヒトはいないでしょｗｗｗ
ちなみにAbacusっぽいやつなら CALCULIX
ttp://http://www.calculix.de/ ttp://http://www.bconverged.com/products.php (FreeのWindows版あり） へどうぞ。
もっと詳しく！
１万円でAbaqusが使えると聞いて、私も購入しました。
ですが私もインストールできません。
windows7だから？
ttp://www.google.co.jp/search?hl=&q=%83K%83X%83%89%83C%83e%83B%83%93%83O 評価版(要素数や支配方程式の数の制限あり)ならフリーでwindows上で動くよ。
原理はよくわからんけどとりあえず結果は出た
イレギュラーな結果を見定める時に多少知ってればいい程度だ。
東大乙
Ansysとか
gnuplotとか
単純なグラフならよさそうだけど。
質問を小出しにするやつにFEMをやる資格はないと思え
FEMのスレなんだから
やっぱParaViewだわな
品のないやつにエンジニアの資格はない。
くれるソフトありますか？
粒子法でGPUバリバリって方があってるかもねw
ソフトありますか？
日に日にエスカレートしてまったくやめる気配がまったくありません
なにを求めるのかも書かずに
アホのような質問されても
誰も答えないっしょ。
Octaveでプログラムを書いているのだけど、境界条件を入れた全体剛性行列のdetがゼロになって解けん・・・
正六面体の一次要素一つの形状を作って片面を全固定(123自由度)、反対面に荷重をかけているのだけど
正六面体だと解けない理由ってのが有るのかな?まあプログラムの何処かが悪いのだろうけど
何から調べていけば良いのかアドバイスお願いします
そんな大それたもの読めません（本気なのかー！）
他人のコードはおもしろいよ。
２次元問題だと解けるの？
2次元の要素をまだ実装していないのです。申し訳有りません。
実装するにしてもまだ時間がかかりそうです
使える要素は四面体と六面体の一次要素です。
四面体も変なのですがまだ細かくつめていません
(時間はかかると思いますが、2次元問題の方も作ってみます)
かのうです
頑張ってみるわﾉｼ
ttp://http://www.str.ce.akita-u.ac.jp/~chida/calculix.html ttp://http://freecaetester.blog62.fc2.com/ ttp://http://freecaetester.blog62.fc2.com/blog-category-1.html ここから読んでいけば、ミトコンドリアでもわかりそうだ。
ブログ作者に感謝だな。
コマンド一覧は
ttp://http://www.usennet.ne.jp/~kiyoshi/manual/CalculiX.html Caliculixに関する情報は結構イッパイあるから、
全く初めての人にとっては、どこから読めばいいのか
さっぱりわからん場合もあるので。
自分で見つけたれば、ググればいいし、
教えてあげたい人は教えればいいでしょ。
質問ある？
どこかでdownloadできますか？
さしあげるよ
PHPコードだけど
インターネッツ
2ch病なのかなあ
コマンド一覧
を
コマンドー覧 と誤読したorz
BEMってどうんな
ttp://http://www.caelinux.com/CMS/ 手動で座標とconnectivityを書たテストは、うまく出ているので書き方は間違ってないと思うのだが？
でいま積んでる。知ってるヤツいる？
魂は幾何学、コピー出来る公安はサリンオウム信者の子供を40歳まで社会から隔離している
オウム信者が地方で現在も潜伏している
それは新興宗教を配下としている公安の仕事だ
発案で盗聴器を開発したら霊魂が寄って呼ぶ来た
<電波憑依>
スピリチャル全否定なら江原三輪氏、高橋佳子大川隆法氏は、幻聴で強制入院矛盾する日本宗教と精神科
<コードレス盗聴>
2004既に国民20%被害250～700台数中国工作員3～7000万円2005ソウルコピー2010ソウルイン医者アカギ絡む<盗聴証拠>
今年５月に日本の警視庁防課は被害者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「(手を振りながら)モジョとは…なんじょ」
お嬢ちゃん、後ろ後ろ。
自己解決すた。ｖｔｋのconnectivityのノード番号って0から始めなきゃいけないのね
使い方を簡単にまとめたサイトとか誰か知りませんか？
これ推奨の環境ってMACなの？
商用使いな。
エラー出まくる・・これはソース直して使うのがデフォ？
コンパイル通らない状態で公開されてるって認識で良いのだろうかってことが知りたい。
亀田も真っ青！前人未到の犯罪５階級制覇や！！
違法駐車王 塩ホ容疑者
器物損壊王 仕置容疑者
飲酒運転ほう助王 死沖容疑者
無法バーベキュー王 Ｊおき容疑者
児童虐待王 SALTWOOD容疑者
ここできいたほうがはやいよ
アホか
お世話になりましたノシ
衝突解析でっきるっかな。
ttp://http://peacelandia.blogspot.com/2011/01/blog-post_1275.html 2011年1月24日
田原総一郎の講演料
去年の話になりますが、１０月１９日に弘前市で行われた「原子燃料サ
イクル意見交換会」に関して、安藤晴美県議が情報開示請求し、その結
果、わかったことです。
企画は、資源エネルギー庁と青森県の共催で、（株）ＲＡＢサービスに
委託されて行われており、委託料は１２９万６７５０円。その内訳は、
ジャーナリストの田原総一郎氏の講演料が１１０万円。交通費８万円。司会料３万円。管理手数料（食費含む）２万５千円。
対応OSがわかれば教えて下さい。
ttp://http://www.simulia.com/academics/student.html これとか
ttp://http://www.amazon.com/gp/aw/d/0072466855/ref=mp_bxgy_p_dp_y 単位もN mm t 系でやってんのに
メッシュをめちゃくちゃ細かくしたら少しだけ理論値に近づくが・・・
基本的に数値は大差ない。
変ならモデルに妙な接合部があるとか
材料データ間違ってる。
スケッチで50*50*1000の梁を作って、密度7.8e-9、剛性206e3、ポアソン比0.3
を入れてるだけなんだけど・・・
現状でてる結果の値も
ttp://www.42ch.net/UploaderSmall/source/1326464747.png poor man's abaqus であるところの calculixでやってみたけど、
41.57Ｈｚじゃない？
ttp://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でできそうなのですが、
これをさらにクネクネアニメーションさせたいです。
どんなやり方がありますか?
レスありがとうございます！
今日（きのう）もう一度、最初からモデルを作って解析をすると、
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位になって誤差が大きくなりますね。
静的な応力解析の時は、メッシュの大きさが大きくても、そこそこの精度は出るのに。
詳しく説明をしていただき、本当にありがとうございました。
ソフトの使い方だけでなく有限要素法についても、さらに勉強しようと思います。
二次じゃないとおかしな値になるよ。
ttp://http://d.hatena.ne.jp/kaityo/20091120/1258738842 うーん。全体を４要素という極めて少ない要素数にしても
ttp://www.42ch.net/UploaderSmall/source/1326499082.png 43Hzくらいで、まぁまぁの値になるんですが。
ttp://www.42ch.net/UploaderAnime/source/1326499532.png 単なる非定常計算結果のアニメーションは
できるんですが、固有値解析のように
モードシェイプに対応する変位データが
ひとつしかない時に、変形前と変形後の
ノード座標値を（マイナス変位側まで）
ParaView上で補完して、物体が振動して
いるようなアニメーションがしたいのです。
FEMAPなどで固有値解析結果を表示するのと
同じことを、ParaViewでできないかなぁと。
変位データは自分で作成しているので、
連番ファイルで補完データを書き出すことも
不可能ではないのですが、数十個ある
モードシェイプ情報に対して全て補完データを
出力すると連番ファイルが1000個近くなるので、
何かうまいやり方があれば助かります。
１次要素って、２次要素より「曲がらない」＝「硬い」ので、周波数としては上がる傾向に あるみたいです。
４要素だと、２次要素なら43Hzですが、１次要素（C3D8）にすると77Hzになっちゃいます。
なのです。これは、ほぼ違憲であり、戦後の、反省の立場から、公明党が政権をとることを想定すると、神道国家に回帰するため、
アジアや旧連合諸国にたいする国際的なプレゼンスや国際関係を著しく悪化するものと思われます。また、別の観点から、創価大学を持っているといえども、「学会」と言う名称も、国際的には、本来の学会の意味でないため、非常に不適切な団体です。
創価学会との無関係を訴える公明党ですが、フジテレビの選挙報道時は、「創価学会のみなさん」と不自然なくらい必ず、放送します。
ブッシュドクトリンに協調した小泉自民・公明は、「テロとの戦い」を表明するも「オウム真理教」に破防法を適用できない失敗をし、野党に追いやられました。
つまり、宗教団体は、特に、世俗国家においては、政党を持つべきではありません。
小泉元総理は、公明党と公約を守って、「自民党をぶっこわす」だけでなく、原発も
ぶっこわしたのです。これを小泉の奇跡と言いましょう。世襲制にも反対しましょう。
Abaqusで、粘着挙動を選択しますとデフォルトでペナルティー法による粘着となっています。
このペナルティー法というのは、どういったものなのでしょうか。
調べてみますと、ＦＥＭ解析における接触の基礎としてのペナルティー法や、計算が速い代わりに食い込みを許すということ等、おおまかなことは分かったのですが、
このペナルティー法でＡ部材とＢ部材を粘着させると、完全に一体として接着されたものとみるのでしょうか？
具体的には、Ａ部材をStep1で伸ばしておき、そこにStep2でB部材を張り付けるということをしたいのですが、張り付けに、ペナルティー法の粘着を使用するとどうなるのかがわかりません。
購入先のサポートに問い合わせろ。
SOL601(NX Adv Nonlin)で二次元の弾性体半球同士の接触解析を強制変位を使ってしようとしているのですが、Connectionsが上手く働いていないのさ、両物体が思い切り貫入するような結果になってしまいます。収束も最後までしません。
Slide Lineを使った二次元解析ではなく、ConnectとSOL601を使って接触解析をするなら、プロパティやリージョンをどのように設定すれば良いのでしょうか！？アドバイス下さい！
離散kirchhoff板理論(DKT,DKQ)とひずみ仮定(MITC8)を
拡張された理論であるのは分かったのですが、
ここから先はさっぱり分かりません。
なにか参考になるソースはありますか？
今見ている文献は
A new discrete Kirchhoff-Mindlin element based on Mindlin-Reissner
plate theory and assumed shear strain fields. PartⅡ
ttp://http://books.google.co.jp/books/about/%E6%9C%89%E9%99%90%E8%A6%81%E7%B4%A0%E6%B3%95%E3%81%AB%E3%82%88%E3%82%8B%E6%A7%8B%E9%80%A0%E8%A7%A3%E6%9E%90%E3%83%97.html?id=5T2HHAAACAAJ&redir_esc=y ttp://http://www.brandany.com ttp://http://www.bvbag.net ttp://http://www.lv-market.net ttp://http://www.gccbag.com ありがとうございます
その本が図書館にあるか調べてみます
もし、アマゾンなどで簡単に購入できる本があれば教えてください
なんつったってAnsysだってFORTRANだ。
ttp://http://www.youtube.com/watch?v=e7Ogui8Srr4&feature=youtu.be 個体の荷重ー変位曲線を割り出すのに使います
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時間かかってしまいます。
そこで、解析手法を変えようと思っているのですが、
どの方法が適切でしょうか？何卒ご協力お願い申し上げます
＞最適化ですので
何をどう最適化したいのか、エスパーしようにも
どうしたものやら。
地盤?それとも構造物?
真面目に変位量の時間履歴を取りたいのなら
愚直に解くしかないんじゃないの?
それをどうやったら他の手法で置き換えられるのか、
それはむしろ君自身が研究テーマとして取り組むべき
問題なんじゃないかと思うのだが。
いま京の研究公募してるから、それに応募してみたら?
今の手法でも結構頑張れちゃうかもよ?
時間が限られてるなら最適化の項目を減らしなさい。
そのまえに0.005秒のデータってあんの？
俺は15万要素くらいで60秒解析してたけど、実際0.05秒とかでも、大して結果かわんなかったぞ
そんな周波数を解析するって時点で無能君。
様々な地震動に対して最適な構造を見つけるってことじゃねえのか
最適化ってそういうことだと思ったんだが
つーか書いてて思ったが1000時間って大したことないよね
パソコンがやるんだし。連続して解析できるんだろ？
２台使えば500時間だ。一ヶ月もかからねえぞ
防衛大だったと思う。飛行機とか衝突してたから
ただ、終わりのほうに、ノートPCでも3ヶ月かければ解析は可能だって書いてて、
こいつは何を伝える気なんだと思ってワロタ
要らないよ。
チェーンで十分だ。
良かった！
心配でずっと眠れなかったよ。
集中荷重の指定方法とか何故繰り返し計算を強いられるのかとか教えてください
分かっていない事が多すぎるのでソースを読めはご勘弁を
OpenFOAMでStressAnalysisをする時点で、アレ。
calculix とか、そーいうの使うほうがいいかと。
StressAnalysisはオマケみたいなものなのですかね・・・
Calculixも平行して為してみているのですがプリ部分でまだ試行錯誤していますw
FreeCADでSTLファイル作成→netgenでメッシュ作成→cgxで荷重、境界条件作成→ccxで計算→cgxで結果を確認
と言う感じにやってますが、単一部材のケースはわりとすんなり行くのですが複数部材のケースがどうしたら良いのか
チラ裏でし
アセンブリは素直にcgxでモデリングする方がてっとりばやいかと。
pythonが使えれば
ttp://dip28p.web.fc2.com/calculix/netgen2calculix/index.html っていうのが使えるらしい。
どんな情報でも結構です
どうかよろしくお願いいたします
cgxか～複雑なモデリングには向いてなさそうだけど使ってみるよ
netgenまでは複数部材は考慮されているみたいで、サンプルで提供されている
複数部材を定義されたgeoファイルを読み込み、abaqus形式で吐き出すと
要素の集合名が分けられて出てきます
CAD系で複雑なモデリングをしていくのが理想に思えるけどここらが自分の限界なのかな
STEP形式ファイルは複数の部材を考慮した出力って出来るのでしょうか？
あ、FreeCADで複数部材のモデリング方法が分からないからダメか・・・チラ裏でした
Ａｓｓ’ｙの難しさとして、
□■
Ａ Ｂ
部材Ａと部材Ｂが一つの面でくっついているとすると、
CGXみたくモデラーで作るときは共通平面ということで作っていけばいいけど
ＣＡＤからもっていくとＡとＢの共通平面は別々の平面になっているはず。
そこでＡとＢをくっつけなきゃならないんだが、たいてい、その別々の平面を
接触（固着）させて解析させることになる。
ここがちょっと難しいところですね。
平面と平面が面で完全に接合している物体なんて現実的ではないし。
割り切るしかないんだろうね…。
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用のインプットファイルを作成して解析、結果確認
分かっている人には当たり前の事かもしれなかったけど、これを分かるまで大分時間がかかったよ
複雑な形状だと節点が重ならなかったり上手く行かないんだろうな・・・
＞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))) );
となって複雑に・・・何かが間違えている
OpenFOAMってVOF法みたいので構造も有限体積で計算していたはず。
弱連成もできると聞いたことあるので、有限要素のソルバーもあるかも知れんが。
詳しくは開発元が書いた論文よめ
ttp://http://onlinelibrary.wiley.com/doi/10.1002/nme.1409/abstract ttp://http://jp.alpha-lab.info/topic-3409-0.html ニューラルネット超入門
ANSYS、NASTRANなどを検討していますが、何かお勧めはありますか？
メインは線形で、非線形も少しやろうかと思っています。
なんちって
固有値計算プログラムが上手く動かん、固有ベクトルがぐちゃぐちゃだ
固有値はハウスホルダー＋QR法で何とかoctaveのeig関数の結果と一部一致したけど
(ゼロに持っていく値が振動して収束しない時があるんだよね…)
三重対角行列の固有ベクトルを求めて(怪しい実装)それをハウスホルダーの副産物を使って逆変換(とても怪しい実装)
見事に撃沈、誰かエスパー求む
メッシュがグチャグチャなのが原因と予想だな
片持ち梁、片方を3成分固定、片方を集中荷重をあたえた時にそれらしい結果が出たデータを使ってます
要素は8節点6面体の一次要素です。メッシュがぐちゃぐちゃとは要素座標とかそこらの考え方がぐちゃぐちゃしているのか・・・
剛性マトリクスを作る時についでで質量マトリクスを作っているからそこらが怪しいと思っているんだけど
積分点の位置で形状関数を使って質量を計算してはいけないのかな?
仕事とか学生とかではないから焦ってはいないんだけど、資料が無くて分からない事が多いんですよね
暇な時に少し知恵を貸してくれると幸いです
まずCalculixみたいなフリーのやつ入れて、勉強したほうがいいと思う。
Salome-Mecaが英語だったら、こっちがオススメなんだけど、フランス語なんで日本人には辛い。
3DCADは何かつかってますか？
SolidWorksだったら、SolidWorksアドオンのCAE（COSMOS?)を載せたほうがいいと思うんだけどさ。
レス有難う
calculixは一応入れてソースを見てるんだけど・・・コンパイルが悪いのか固有振動数計算が解けないんですよね
(*ERROR in arpack: info=-3 # of converged eigenvalues=70 とか出ます)
3DCADはFreeCADを使ってますが、形状がそんな複雑でないのとメッシャーが四面体要素中心なのでテストデータは手で作ってます
使える要素で静解析の結果が一番良かったのが六面体一次要素なのでそれを中心に進めています
使っているプラットフォームがdebian(squeeze 32bit)なのと貧乏なので商用の物は使えないのですが
SolidWorksですかちょっと調べてみます
済まない、レスを混同してしまいましたは無かったことに
構造解析なら、メッシュ数おとしてでも、２次要素を使ったほうが絶対いいです。
絶対っていうのは言いすぎかｗ
だけど、自分の経験からいうと、メッシュ数がテキトーでも二次要素だとかなり精度よく解けます。
二次要素とか気にしなくて良いANSYS DesignSpaceにしとけば？
非線形は必要になったらライセンスアップグレードすれば良い。
六面体の2次要素はまだ実装していないのです
,191
すみません、レスを混同して変なレスをしています
固有値グチャグチャは気になるけど 固定が駄目なのかなー
固定か～固定している所の剛体行列と質量行列をどのように扱ったら良いか分からないから固定条件無しでやっています
まず剛体モードが綺麗に直行で出てくれないし何か変だ～
今は固有ベクトルの計算は逆反復法らしきものでやっています
(上手く収束しないんだよねと言うかほとんど収束しない、振動してたり反復の上限一万回に引っかかってたりするよ)
手作りだからね・・・まずは一次要素からかなと思って作ってます
(ちょこちょこ弄っていて大分汚くなってきたよw)
いやいやいやいや、
６面体だろうが、なんだろうが、２次要素と１次要素じゃ精度が全然違う。
普通の構造解析で積極的に１次要素を使う理由は、全く無い。
↓
ttp://www.42ch.net/UploaderSmall/source/1344995786.png わからない。
学問板なので、理論的に説明してくれ。
その差は、特に長手方向に粗いのが原因っぽい
メッシュ細かくしてきゃ、どっちも値大きくなって、双方一緒っぽい値になる筈
テトラ一次は細かくして行っても、2次と全然違う値になったりするけど
まぁそうなんだけどねｗ
つまり、二次要素は、メッシュ切りがテキトーでもそれなりの精度が出る。
一次はメッシュの細かさが、モロ、精度に出るので、メッシュ密度とか形状に非常に神経を使わなきゃならんので
実務的には向いていないという考えです。
テトラの2次と1次は 3割4割違ってくるので 比べるとマシだけど
ちょっと差があるのは気になるね ウーム
あ、書き忘れていた、これの理論値は -0.75mmです。
つまり
理論値 -0.75
1次６面体（8節点）要素 -0.72
２次６面体（（20節点）要素 -0.63
う、逆だった
理論値 -0.75
２次６面体（（20節点）要素 -0.72
1次６面体（8節点）要素 -0.63
その「倍」っていうのはどこから来る数値なのか？
たとえばこういう単純ハリの場合、一次要素で解析するには、
縦横高さで何等分、アスペクト比をいくつにすればいい、
っていう目安ってありますか？
3次元だから2の3乗で8倍のメッシュを必要とするという理解で良い？
メッシュ数が「倍」であれば、辺の長さを2の立方根で割れば良い？
固有値、固有ベクトルの計算をする時の固定条件の入れ方を教えて頂けないでしょうか？
ググっているのですがそれらしい文献を見つけられずに困っています。
詳しく書けと言われそうなので少し書いておきますが長いので面倒なら読み飛ばしてください
固有値、固有ベクトルの計算プログラムを作っているのですが、今までは固有ベクトルがぐちゃぐちゃだったのを
逆反復法で計算することでそれっぽい(と言ってもまだ怪しい)ものが出るようになったのですがoctaveで
プログラムしているのでとても遅いです。そこで比較的規模が大きくても良さそうなサブスペース法で組むことに
したのですがどうも小行列を作るとそれがコレスキー分解出来ないのです境界条件を設定していなく、剛体モード
を含んでいるのが原因ではないかなと思っています。その他にヤコビ法で小行列の固有値、固有ベクトルを求めると
サブスペース法の繰り返し計算が収束しない様で困っています。エラーが出ない行列を手作業で与えてやると
ビルトインのeigの結果と比較して符号が違う固有ベクトルがでていました。ハウスホルダー＋QR法＋逆反復法で
出した結果はeigと比べてパターンは違いますが符号が違う固有ベクトルが出ていますがサブスペース法は収束します
ここらもエスパーレスしてくれると有難いです
テストデータは次の感じです
・方持ち梁
・六面体一次要素(要素数1 と5これが主 と 625これは時間が足りなくて計算できない)
・拘束条件無し(なんかこれがいけないっぽい)
＞・拘束条件無し(なんかこれがいけないっぽい)
その通り。
固定条件の入れ方って、単に、固定する節点の変位を0にするだけじゃないの？
ttp://http://www.jsce.or.jp/committee/struct/civil_edu_home/prog/anketo/H19_femdoc02.pdf その通りだと思います。
ただ、ちょっと複雑な形状をきると、アスペクト比１を目指してメッシュを切ると
要素数が非常に多くなることがあるので、１次要素を使ったほうが計算時間が
かかってしまう事はあるようです。
スライドできる状態＝摩擦ゼロ状態いうのは、厳密に起こり得ないので、
条件に対し解を求めるいう点では合ってるけど
スライドできるいう条件設定が、果たして現実的か？ そこはかなり難しい
例えば土木や建築の分野だと 地盤をモデル化したり、剛性体を挟んだり
色々しないと、単品＋拘束やすべり条件設定では、合わない
変形の計算(静計算って言うのかな?)はその通りだと思うのだけど、固有値、固有モード計算と同じなんですか?
変形計算の時って [K]{U}={F}で固有値計算の時は[K]{U}=λ[M]{U}ですよね。これをどう節点変位ゼロと
絡めたら良いのか分からないです。頭が悪くて申し訳ない
固有振動の場合は固定点がなくても解が出るはず。
収束計算させるんじゃなくて、ガウス吐き出し法とかで解けば？
境界条件がある固有値、固有ベクトル計算も知りたいので宜しくお願いします。
固有値＝フリーでの計算ばかり
固定など条件設定を伴ったものは見た事がない気もします できないのかも知れません
出来ないってこたーないだろ。
AnsysだろうがNastranだろうが、出来ないCAEを見たことねぇぞ
わけわからん。
アホ？
税金の無駄遣いはやめてくれ。
解析ソフトって実際何やっているのかブラックボックスで使っているじゃない
どんな処理をやっているのか気になるじゃない
Calculixのソース見てもさっぱり理解できなかったから小さい所でちまちま作っているわけ
オレだったらフリーか個人向けの安いのを買っていじるくらいだな。
バカっぽい
宮脇豊宮脇昇宮脇誠宮脇明宮脇稔宮脇修宮脇勉宮脇三郎宮脇武宮脇和夫
宮脇一郎宮脇幸雄宮脇孝宮脇一男宮脇功宮脇一夫宮脇勲宮脇守宮脇健一
宮脇保宮脇秀雄宮脇正一宮脇幸男宮脇浩宮脇宏宮脇登宮脇正雄宮脇義雄
宮脇正夫宮脇秀夫宮脇敏夫宮脇光男宮脇昭宮脇徹宮脇文雄宮脇忠宮脇章
宮脇正男宮脇和男宮脇栄宮脇幸一宮脇一雄宮脇満宮脇寛宮脇良一宮脇晃
宮脇栄一宮脇勝美宮脇光雄宮脇薫宮脇剛宮脇正明宮脇正義宮脇健
宮脇幸夫宮脇英雄宮脇繁宮脇信夫宮脇一宮脇敏雄宮脇宮脇章宮脇正男
宮脇和男宮脇栄宮脇幸一宮脇一雄宮脇満宮脇寛宮脇良一宮脇晃宮脇栄一
宮脇勝美宮脇康充宮脇薫宮脇剛宮脇正明宮脇正義宮脇健宮脇幸夫宮脇英雄
宮脇繁宮脇信夫宮脇一宮脇敏雄宮脇仁宮脇健二宮脇和彦宮脇忠雄
宮脇利夫宮脇哲夫宮脇正人宮脇正美宮脇悟宮脇学宮脇武夫宮脇久雄
宮脇正治宮脇洋一宮脇文夫宮脇清一宮脇孝宮脇武雄宮脇誠一宮脇信一
宮脇猛宮脇貞夫宮脇久宮脇文男宮脇英夫宮脇忠男宮脇洋宮脇昭夫
宮脇和雄宮脇哲也宮脇彰宮脇邦夫宮脇俊夫宮脇貢宮脇健治宮脇清一
宮脇孝宮脇武雄宮脇誠一宮脇信一宮脇猛宮脇貞夫宮脇久宮脇文男
宮脇英夫宮脇忠男宮脇洋宮脇昭夫宮脇和雄宮脇哲也宮脇彰宮脇邦夫
宮脇俊夫宮脇貢宮脇健治
NO TO ACTA!
ttp://http://www.youtube.com/watch?v=P1V2nhNIGKk 何日か書き込みできる状態ならまた書く～
北海道稚内兵庫県立花西宮市
兵庫県明石市千葉県房総兵庫県神戸市
兵庫県小野市
・沢山の人が偏った意見を一貫して支持する
偏った意見でも、集団の中でその意見が信じられていれば、自分の考え方は間違っているのか、等と思わせる手法
・不利な質問をさせなくしたり、不利な質問には答えない
誰にも質問や反論をさせないことにより、誰もが皆、疑いなど無いんだと信じ込ませる手法
↑マスコミや、在日カルトのネット工作員がやっていること
TVなどが、偏った思想や考え方に染まった人間をよく使ったり、左翼を装った人間にキチガイなフリをさせるのは、視聴者に、自分と違う考え方をする人間が世の中には大勢いるんだなと思わせる効果がある。
..
強盗公然わいせつ児童ポルノ法放火
固有値、固有ベクトルの計算をする時の固定条件の入れ方を教えて頂けないでしょうか？
ググっているのですがそれらしい文献を見つけられずに困っています。
詳しく書けと言われそうなので少し書いておきますが長いので面倒なら読み飛ばしてください
固有値、固有ベクトルの計算プログラムを作っているのですが、今までは固有ベクトルがぐちゃぐちゃだったのを
逆反復法で計算することでそれっぽい(と言ってもまだ怪しい)ものが出るようになったのですがoctaveで
プログラムしているのでとても遅いです。そこで比較的規模が大きくても良さそうなサブスペース法で組むことに
したのですがどうも小行列を作るとそれがコレスキー分解出来ないのです境界条件を設定していなく、剛体モード
を含んでいるのが原因ではないかなと思っています。その他にヤコビ法で小行列の固有値、固有ベクトルを求めると
サブスペース法の繰り返し計算が収束しない様で困っています。エラーが出ない行列を手作業で与えてやると
ビルトインのeigの結果と比較して符号が違う固有ベクトルがでていました。ハウスホルダー＋QR法＋逆反復法で
出した結果はeigと比べてパターンは違いますが符号が違う固有ベクトルが出ていますがサブスペース法は収束します
ここらもエスパーレスしてくれると児童買春・児童ポルノ法違反有難いです
剛体モードが出るっていうのは、つまり、逆行列が存在しない っていう意味ですかい？
(（1,1）
（2,2))
みたいな。
7月スタートのTBS・ドラマNEO枠『天魔さんがゆく』（毎週月曜 深0： 28）
で、『33分探偵』シリーズ（フジテレビ系）で組んだ福田雄一監督と全く新しいホラーコメディーに挑む。
2009年3月～4月放送の『帰ってこさせられた33分探偵』（全4話、フジテ
レビ系）を最後にドラマから遠ざかっていた剛。今回、「親友の福田監督か
らお声がけいただき、また面白い物が作れる、どんなものが飛び出してくる
のかとワクワクしています。福田雄一と堂本剛のワールドをぜひ楽しみにし
てください」と張り切っている。
演出による完全オリジナルドラマ。剛演じる主人公・天魔は、先祖代々
幽霊とコンタクトを取る能力を持っているが、その力に反して大の怖がりで、
「幽霊退治引き受けます」と掲げた会社の2代目という設定。米映画
『ゴースト・バスターズ』（1984年）のようだが、退治する意気地がない
天魔は、幽霊の言い分を聞き、この世に残している悔いを晴らしてあげたり、
優しく諭してあげたりしながらさまざまな幽霊退治の依頼を解決していく。
天魔の会社の女性社員に川口春奈、先代からの番頭でもある皆川猿時、
天魔を追いかけ回す近所のお姉さんに芹那、天魔と夢の中で対話する亡き
父親に佐藤二朗容疑者、天魔たちの周りを嗅ぎまわる刑事に森崎博之などが出演。
依頼人役や幽霊役のゲストにも注目だ。
できるだけ均質な（1要素の体積や辺の長さなどが近い）メッシュを作成したいのですが
どのようにすればよいでしょうか？
現在はTetmshでの自動生成を用いていますが，酷すぎて困っています．
よろしくお願いいたします．
ttp://http://www.hoshusokuhou.com/archives/29932855.html 先にメッシュシードを植える
本人たちがじゃれあっているつもりだったのでしょうが飛行機を待っていたお客さんは
ＪＡＬの服を着た人たちが乱闘している喧嘩していると勘違いした
それをｙｏｕｔｕｂｅにのってしまってＪＡＬは大恥
ＪＡＬグランドサービスはそんな会社です
俺たちの課長が代わって千代っていう馬鹿課長になったんさ
前の課長はまだ部下 したの人間の言い分も聞いて器のでかさがあったけど
その後きた千代って言うのがとんでもなくアホ
俺が周りの係長とかに苛められているという事実も今まで周りから不愉快な言いがかりや
暴言などを耐えていたという状態も知らず俺を悪い奴と決め付け絶対俺をやめさせると
いい 毎日のように粘着し ついには俺を会社から追い出すという暴挙にでた
「お前みたいな使えない奴を長くいさせるよりかは訴えられて賠償金を払ったほうが安くすむ」という暴言も吐かれた
それから「俺は今まで自分が使えないと思った奴を辞めさせてきた」と俺に豪語したこともあった
私はラバー（トイレ）を夢中で掃除をしていていました
そのときトントン トンと私の肩を叩く人がいたので振り向いてみるとその飛行機のクルーの一人が
目を吊り上げて私の前に立っていました
そして私に「あの人たち 口を動かしてばかりで仕事しない 注意してきてくれ」と私にいいました
私は急いで係長にそのことを報告しにいきました 班長はその横で口を緩ませてヘラヘラしているだけでした
他社からあいつら喋ってばかりで仕事をしないといわれても恥じもなく 適当な仕事を他社に見られても平気な
ＪＡＬグランドサービス
このことでその飛行機会社からは信用を失い二度とＪＡＬ系列に仕事を任せないと
決断してくれることを切に願うばかりである
あと係長が機内においてあった古新聞を一人で別の仕事をするとき車の中に持ちこんでいた
これ二つともやってはいけないこと
飛行機の中の物は外に勝手に持ち出したらいけないのにもしかして係長になったら許されるとおもってんのかな？
フィロソフィとかつくってもなにも役に立たないね
こんなことばかりやってるからＪＡＬグランドサービスは屑会社とか馬鹿会社とか言われるんだよなｗ
それから仕事中ガム噛んでいたりｗ
正社員の質の悪さっていったらそりゃもう
機内のものは持ち出し禁止なんだけどさ 新聞を持ち出しで車の中で読んでいるクズ係長もいるしｗ
なんでもありｗ あいつらの都合の悪いことは隠蔽したがるけど 俺はあのクズ会社辞めたからカキコするわｗ
つまり俺はＪＡＬグランドサービスの敵になったわけｗ
密告しまくりよｗ
さらに次の清掃場所にもっていって汚れたタオルを使えっていってきた
タオルからはなんともいえない臭いが漂い色もケチャップやなんやらがこびりついていて汚かった
お客さまのためにを毎日言っている会社なのになにこれ？ 矛盾してるよね？
青タオルはウンコがこびりついていて汚い それをキレイなタオルといいラバーのほかギャレー
床のゲロ掃除などさまざまなところで活用している
バケツにそのタオルを濡らすとき気持ちが悪い バケツもすごく汚いよね
バケツ洗っているところみたことがないし
ＪＡＬはこんな会社です
白タオルは生乾きで蒸れた悪臭がしている それでお客様がつかうテーブルを拭いているんだから
あのテーブルは菌をひろげているようなもんだ
ラバー（トイレの便器）やトイレの床を拭く尿やウンコつきの青タオルを洗う
バケツといっしょにいれてるから臭くなるんやｗ
それでお客様がつかうテーブルを拭いているのだからＪＡＬの衛生面は最悪
こいつらは所詮お客様の姿を金としかおもっていないからそういう掃除しかできないんや
口ではお客様のためとか言ってるけどなｗ
それを正せる上司が一人もいないというのが特徴の会社
業務妨害で訴えられちゃうよ。
解析実行して.datファイルがきちんと生成されるのに（内容は変位と応力とひずみだった）
cgxのコマンドに「data set」がなくて結果が読み込めないっぽい？
コマンドというかcgx作業画面の空白をクリックして出てくる機能の中にdatasetがない
結果データはdatファイルを読み込むわけじゃなくてfbdファイル内にそのまま入ってるのかなと思ってfdbをメモ帳で開いてみたけど結果は入ってない
ファイルの書式があればいいんだけどマニュアルにそこまで書いてあるのかな←いまここ
やっぱり導入時が一番苦労するなあ
fdbと書いている時点で怪しいかも(大分弄っていないので記憶が曖昧だけどfdbは違うんじゃないのかな)
もう見つけているかもしれないけど「有限要素法のフリーのソフトウエア CalculiX を試す」でググると良いかも
うちはLinux環境だからWindows特有の問題に巻き込まれているなら分からんよ
○fbd
うお、間違えた
もういい寝る
Workbenchですか？ Classic（APDL)ですか？
ちゃんと契約してるベンダーに聞けよ。
いやいや、ある面に熱伝達係数を設定した上で、その面の一部に接触面がある時の取扱いについては、
マニュアルに載っていないので、Quick-Aレベルだな。
classicです。
コンタクトウィザードの該当面の接触プロパティの
Thermal Contact Conductance が、デフォルトだと０になって断熱しているので、
この値をテキトーに入れてみてください
モデリング→作成→接触ペアのところでいいですか？
あと相互に熱伝導をしている状態を再現したいのですがたとえば物体A、Bが接触してるとして
コンタクトAターゲットBとコンタクトBターゲットAのように二回接触定義をすればいいんでしょうか？
２回やっちゃダメです。
互いに熱伝導っていうのはアタリマエで、片方に熱伝導するなんて非現実的な設定は、元々設定出来まセン。
了解しました。そうなるとThermal Contact Conductanceのテキトーな値というのはどう設定すべきなんでしょうか？
両者の熱伝導率の平均でいいんですかね？
違います。
これ以上の事はノウハウなので、サポートに聞いてください。
ttp://uni.2ch.net/test/read.cgi/rikei/1321761766/575 ↑ ↑ ↑ ↑ ↑ ↑
僕は働かないくせにお金ほしいクズ野郎です。
288
≫256
大丈夫。借金まみれなのにシステムトレードで不労所得目指してるクズ野郎いるよ。
ttp://http://manta.blog.jp/ 「糞汁は経済板から消えるべき」
「産業歴史博物館設立実行委員会（エンジン博物館）」
※スレタイに「北斗の軍」と入ってるスレ立てもこのクソ野郎の仕業
(´-ω-`)y-~（禁煙マスター） ← 史上最悪の中年荒らし （トリップを変えて別人になりすます技も覚えて来たようですｗ）
こいつのズリネタ
ttp://http://ameblo.jp/garo-farann/image-11515076726-12507503444.html 得意技：情報漏洩、ナマポ受給、糞ブログコピペ、ケムトレイル、人工地震、かまってちゃん、驚異の前言撤回ｗ、都合が悪いとすぐ話題脱線
ナマポで助けてもらい ttp://toro.5ch.net/test/read.cgi/occult/1354090242/78">
ttp://toro.2ch.net/test/read.cgi/occult/1354090242/78 日本人滅びろ発言
ttp://kohada.2ch.net/test/read.cgi/eco/1342954026/993 情報漏洩をきっかけに
ttp://kohada.2ch.net/test/read.cgi/eco/1326591504/646 ttp://kohada.2ch.net/test/read.cgi/eco/1342954026/926 別れの挨拶１！
ttp://kohada.2ch.net/test/read.cgi/eco/1356935230/329 別れの挨拶２！
ttp://toro.2ch.net/test/read.cgi/occult/1354090242/303 別れの挨拶３！
ttp://kohada.2ch.net/test/read.cgi/eco/1356855992/54 そして告発の逆恨み！
ttp://kohada.2ch.net/test/read.cgi/eco/1342954026/996 精神異常をきたした ttp://kohada.5ch.net/test/read.cgi/eco/1374582405/318">
ttp://kohada.2ch.net/test/read.cgi/eco/1374582405/318 超低学力！
ttp://kohada.2ch.net/test/read.cgi/eco/1369571885/587 さらにロリコンも暴露！
ttp://kohada.2ch.net/test/read.cgi/eco/1376273900/580-629 何を言ってるのか？
ttp://kohada.2ch.net/test/read.cgi/eco/1375448326/597 無駄にスレを作ってる張本人ｗ
ttp://toro.2ch.net/test/read.cgi/occult/1349421194/256 これはひどい
ttp://kohada.2ch.net/test/read.cgi/eco/1372265052/512 有限要素法っていうのでしょうか？
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。ひんがら目気色悪すぎこっち見んな死ね。
使う方としては最大手の一級品と比較してああしろこうしろと言ってしまうが。
パナのCAEの達人みたいな人がキレキレのトークで絶賛しながら滅多切りしてた。
ああいう人が上手く叱咤激励してくれるとソフトも育つ。
ansysもユーザー事例発表で東芝の人が不満点をぶちまけて早く直せと発破をかけてたな。今年のユーザー会は明日だが、さて
2chねらーに荒らされている。足立区に住んでいるそうだ。
ttp://http://inumenken.blog.jp/archives/7002197.html 2chねらーに荒らされている。足立区に住んでいるそうだ。
ttp://http://inumenken.blog.jp/archives/7002197.html LS dynaでS字のパイプに異方性材料を適用したいのよ。
内圧がかかる方向には高いヤング率を設定して曲げ方向には低いヤング率にする。
普通の異方性材料のカードでやるとひとつの座標系でしかベクトルを設定できないからもとから曲がりくねってるパイプにはうまく設定できない。
分かりにくくてごめん
それだと座標系を山ほど作らないといけないのでは？
それでいいんじゃね？
◎ +ニュース
・ 2NN
・ 2chTimes
◎ +ニュース新着
・ 2NN新着
・ Headline BBY
◎ +ニュース他
・ Desktop2ch
・ 記者別一覧
◎ 全板
・ 全板縦断勢いランキング
・ スレッドランキング総合ランキング
◎ 実況込み
・ 2勢
・ READ2CH
・ i-ikioi
※ 要サイト名検索
ttp://https://www.youtube.com/watch?v=dk3Kro3jZ0A&feature=youtu.be ttp://https://www.youtube.com/watch?v=YS_nZxxVvrk&feature=youtu.be ヽo,´-'─ 、 ♪
r, "~~~~"ヽ
i. ,'ﾉﾚノﾚ!ﾚ〉 ☆ 日本のカクブソウは絶対に必須です ☆
__ '!从.ﾟ ヮﾟﾉﾙ 総務省の『憲法改正国民投票法』のURLです。
ゝﾝ〈(つY_i(つ
ttp://http://www.soumu.go.jp/senkyo/kokumin_touhyou/index.html ｀,.く,§_,_,ゝ,
~i_ﾝｲﾉ
ttp://http://www.youtube.com/watch?v=HJKHJEJgxgo ユーザーの方もユーザーでない方もぜひご参加お待ちしております。
Nobuyuki Umetani
ttp://http://www.youtube.com/channel/UCqObO-PlZgCcz0z54GsE1YA 今Abaqusで荷重の代わりに強制変位を与えて板の弾性座屈解析を行っているんだけれども
この場合Abaqusが出力する固有値ってどんな意味があるんでしょう？
調べてもなかなか答えが見付からず・・・
ttp://http://civilyarou.web.fc2.com/WANtaroHP_F90_html5/f90_FDA/TeX_EGV_B.pdf についての調べる場合には､どんな論文集・本等を勉強すればよいか、どなたか教えていただけませんか。
遅くなりましたがありがとうございます
これを読みつつなんとか解決することができました
/:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::://ヽ:::::::::::::::|
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 | |
ttp://http://www.youtube.com/watch?v=z2qK2lhk9O0 ☆ 自民党、グッジョブですわ。 ☆
ttp://http://www.soumu.go.jp/senkyo/kokumin_touhyou/index.html ☆ 日本国民の皆様方、2016年7月の『第24回 参議院選挙』で、改憲の参議院議員が
3分の2以上を超えると日本国憲法の改正です。皆様方、必ず投票に自ら足を運んでください。
そして、私たちの日本国憲法を絶対に改正しましょう。☆
のxをt=0から順に計算していくにはどうすればいいの？
ttp://fluid.nuae.nagoya-u.ac.jp/archives/textbook/past/nensei/Nensei2001Chapt2Print.pdf 一からプログラム組んでる勢なんだけど死にたい．
これなんかどう？
＜解析塾秘伝＞有限要素法のつくり方！ＦＥＭプログラミングの手順とノウハウ
ttp://https://pub.nikkan.co.jp/books/detail/00002756 All together now～ All together now～ みんなで人間革命（Yeah～！）
ttp://http://goo.gl/5YOYqX ttp://http://goo.gl/7Sastj ttp://http://goo.gl/1ondZM 後藤達俊 バックドロップ43連発
ttp://http://goo.gl/HkEPgn ・沢山の人が、偏った意見を一貫して支持する
偏った意見でも、集団の中でその意見が信じられていれば、自分の考え方は間違っているのか、等と思わせる手法
・不利な質問をさせなくしたり、不利な質問には答えない、スルーする
誰にも質問や反論をさせないことにより、誰もが皆、疑いなど無いんだと信じ込ませる手法
偏った思想や考え方に染まっていたり、常識が通じない人間は、頭が悪いフリをしているカルト工作員の可能性が高い
靖 国 参 拝、皇 族、国 旗 国 歌、神 社 神 道を嫌う カ ル ト
10人に一人は カ ル ト か 外 国 人
「ガ ス ラ イ テ ィ ン グ」 で 検 索 を ！
まずは関数解析から
有限要素法って、いわゆる変分法のことだよ。
いろいろな解析対象を細かくメッシュで区切って、力学的に解析する手法。
一度、Ｒｕｂｙで作ってみたいなぁ。
ttp://http://ads00.com/affiliate/link.php?id=N0002971&adwares=A0000055 ttp://https://youtu.be/quIHgwuF6r4 ttp://https://youtu.be/quIHgwuF6r4 ttp://https://www.youtube.com/watch?v=quIHgwuF6r4&sns=em ttp://https://www.youtube.com/watch?v=quIHgwuF6r4&sns=em ttp://https://www.youtube.com/watch?v=quIHgwuF6r4&sns=em ttp://https://www.youtube.com/watch?v=quIHgwuF6r4&sns=em FORTRAN
ﾌﾞﾝﾃﾞｽ見れない
CLEL見れない
音楽聞けねえちきしょう
同和のクズ共死ねクソ共がざまあみろ気違い共
ほれ気違い共もっともっとドア閉めろ通れ
それしか能のない能無し共がざまあみろ地獄に落ちろ
悔しいか、ざまあみろくたばれクソ同和
お前らの恐ろしさをもっと見せてみろ。そんなんじゃなんともねえぞ
袋とじ見たぞ
悔しいか、ざまあみろくたばれクソ同和
嫌がらせがエスカレートするのが楽しみでしょうがない
今それだけが楽しみだ
俺の生き甲斐藁
それだけ怒ってるってことだもんな藁
ttp://https://goo.gl/Nximd0 これ本当？？
普通にショックだね。。
ttp://https://www.youtube.com/watch?v=CtiacppR5dk 9:27人工衛星（確実な部分）
ttp://https://www.youtube.com/watch?v=-Ls8O7jjK1A まぁ、それでもいいんだけどねぇ…
下手糞が組んだら何ができるものやら…
（どっかのFortran90/95による有限要素法うんたら（丸善）を見ながら）
昔の矢川宮崎本とか平井本とかのほうが優しいかもしれん。
平井本は入手不可能だが。
出来れば分野か、解く方程式書いたほうがわかりやすいよ。
ttp://http://www.data-max.co.jp/280113_ymh_02/ ↑ 自民改正案の真の狙いは言論の自由を奪うこと！
ttp://http://blog.goo.ne.jp/kimito39/e/ec37220f64a8e1d6ed732dd0ab95cbf0 ↑超危険な緊急事態条項で人権無視の内閣独裁に！
ttp://https://www.youtube.com/watch?v=h9x2n5CKhn8 ↑ 自民党は 国民に基本的人権はないと断言！
ttp://http://xn--nyqy26a13k.jp/archives/31687 ↑ 小池都知事も安倍と同じく 憲法改正で 人権
無視の大日本帝国憲法に戻す民主主義破壊論者！
ttp://http://www.mdsweb.jp/doc/1488/1488_03f.html ↑”9条自衛隊明記”は 9条無効化だった！
ttp://http://blog.goo.ne.jp/ngc2497/e/8899f65988fe0f35496934dc972e2489 ↑ ネトウヨ= 安倍サポーター工作員はネットで国民を騙す。
ttp://https://dot.asahi.com/aera/2016071100108.html?page=3 ttp://http://blog.goo.ne.jp/kimito39/e/c0dd73d58121b6446cf4165c96ebb674 ↑ 安倍自民を操るカルト右翼「日本会議」は国民主権否定。
国民投票や選挙では自民党、維新、小池新党に絶対に入れるな。
こういうのを仕事にしたかったのかもしれない
縮めようとすると伸ばそうとする力がはたらきます。
このような元に戻ろうとする力を弾性力といいます。
作用⋅反作用の法則にのっとった反作用の力です。
ばねの伸び、あるいは縮みとして目に見える形で現れます。物理では珍しいことです。
伸び縮みの長さは弾性力の大きさにきれいに比例します。
（といいますか、高校物理ではそのようなばねについて考えます）。
ばね秤ばかりは、ばねの伸びの長さによって物の重さを表示する装置です。
ばねの伸び縮みの長さを x [m] 、弾性力の大きさを F [N] とすると、
F = k x
という関係があります。フックの法則といいます。
k はばね定数と呼ばれる比例定数です＊。ばねの伸びにくさを表します。グラフの傾きです。k が大きいと、伸ばすのに、大きい F が必要になります。
すべては、ここから始まる。
しかし、バネの弾性力は目に見える
だから、物理の基礎を勉強するにはバネがいい
同じばねを３つ用意して、ばねが１つだけの場合と、ばねを２つ直列につなげた場合とで比較してみます。
どちらも20cm伸ばすとすると、フックの法則から、１つの場合の弾性力は F1=0.2×k [N] 。２つの場合は各10cmだけ伸びることになるから１つ当たりの弾性力は F2=0.1×k [N] 。
そうしますと、２つつなげた場合の弾性力は2倍の 2×F2=0.2×k [N] でしょうか？
途中が F2[N] ならどこもかしこも F2[N] です。
ばねを伸ばして静止した状態というのは力がつり合った状態です。
ばねの各微小部分同士が同じ力で引っ張り合ってるので静止しているのです。
ミクロな視点でいえば、ばねを構成する原子たちがお互いを F2[N] で引っ張り合ってつり合って静止しているのです。
端っこでも途中でもどの部分においても各微小部分同士は同じ力で引っ張り合ってつり合って静止しています。
というわけで２つつなげた場合の弾性力は 2×F2[N] ではなくて F2=0.1×k [N] です。
ばねが１つのときの F1=0.2×k [N] の半分です。
２つつなげた方が伸ばす力がいらない、すなわち伸ばしやすくなるのです。
このことは、ばねというものは、同じ材質、同じ形状なら長い方が伸ばしやすい、つまりばね定数が小さいということを意味します。
重い荷物を持つ場合は２人で持てば半分の力で済みます。
２人対１人で綱引きをすれば２人の方が勝ちます。
つまり同じばねを２つ並列につなげると、１つのときに比べて倍の力が必要です。
ばねは伸びにくくなるのです。
２つで１つのばねと見なしたときのばね定数は2倍になります。
並列の場合は単純に足すだけです。
ばねは伸びにくくなります。
これはばねが太くなったことにたとえられます。
ばねは太い方が伸びにくいです。
身の回りにある物体すべてに成り立ちます。
別な言い方をすれば、「す べてのものはバネである」と言っていいのです。
考えてみるとこれはすごいことです。
非常に幅広く成り立つ普遍的な法則といえるわけです。
なぜそんなに普遍的に成り立つのかというと、物質の構成要素である原子や分子の間に働く力に「フックの法則」が成り立っているから と考えられます。
ただし、その中身はマトリクスの形式で表されていて一見複雑ですが、やっていることはバネの公式と何ら変わりがありません。
[K]｛u}＝{f}
これが、有限要素法で最終的に解く式になります。
変な括弧が付いている他はバネの公式と同一です。
[ ]はマトリクスであることを表し、{ }はベクトルであることを表しています。
要素一つについての関係を表す場合には要素剛性マトリクス、解析モデル全体について重ね合わせたものを全体剛性マトリクスと呼びます。
{u}は変位ベクトルと呼び、それぞれの節点の変位になります。
これは拘束条件の入力にも用いますが、解析で求める値そのものになり、コンピュータはこれを求めるためにCPU時間のほとんどを使用します。
{f}は荷重ベクトルで荷重条件を設定した場合、各節点に対応した荷重値がここに入力されます。
この辺の詳しいところは境界条件の設定の項で詳しく説明します。
[K]｛u}＝{f}
何を求めたいかといえば、変位ベクトル｛u}を計算して求めたい。
{f}は荷重ベクトル
この２つが分かれば、変位ベクトル{u}も出てくる
伸び縮みの大きさはバネ定数に反比例し、バネに加えた荷重に比例します。
身の回りのモノはすべて、バネ、つまり弾性体と考えることができます。
変形は外力の大きさに比例し物体の剛性に反比例します。
物体の変形が外力および剛性と比例関係にある場合を線形(linear)といいます。
線形的な静的挙動を表す物体を有限要素法で表すと、[K]{u}={F}という行列方程式の問題に変換できます。
荷重ベクトルは物体にかかる荷重の大きさ
そして、求めたい変位ベクトルは、変形の大きさということになる
境界条件の設定で決まる。
あとは、剛性マトリクスが決まればいい。
ひずみが分かれば、応力も出る
だから、変位を求めることが第一目標
それが基本
中学の理科のレベルだけど、ここが一番大事
↑これが、フックの法則の行列表現
行列とベクトルが使われるからには、一応、高校の数学のレベルだ
この式をもとにして、変位ベクトル｛u}を求めるのが第一目標
そのために必要なのは、節点の剛性マトリクス
それがなんで行列になるのかについて、中学生に説明するのは難しい
節点に生じる変位の方向には、水平・垂直・回転の三つがある
つまり、１つの節点につき、3つの自由度がある。
たとえば、ピン支点では、鉛直方向と水平方向に固定され、回転だけ自由になっている
つまり、節点がピン支点になっている場合は、自由度が1つしかない
こんな具合にして、自由度の数を数える
ここでいう自由度の数、それこそが、求めるべき未知数の数ということになる
[K]｛u}＝{f}
のベクトル ｛u} と {f} のカッコ内にある要素数が、それぞれ8個になる
｛u}ｔ ＝ ｛u1 u2 u3 u4 u5 u6 u7 u8}
｛f}ｔ ＝ ｛f1 f2 f3 f4 f5 f6 f7 f8}
（つづく）
この場合、剛性マトリックス[K］は、8行8列の行列になる。
行列を中学生にもわかるように説明するなら、
「多くの数値をまとめて書く、便利な表記法」
といったところか。
ここで、荷重ベクトルを設定する。
２番目のところ、ｆ２で、荷重Pがかかっているとする。
それ以外のところ、ｆ1および、ｆ3～ｆ8では、荷重がかかっていない。
そうすると、荷重ベクトルはどうなるか。
荷重ベクトル ｛f} ＝ ｛0 P 0 0 0 0 0 0}
ということになる。こんな風に、行列やベクトルの中に0が多くないと、有限要素法らしくない（笑）
（転置記号は省略）
変位法とは「変位をもとめる方法」です。これはいたって明快です。
いや、そうはいっても「応力」だって出ますよね、という反論があるかもしれません。
しかし、変位が分かれば応力は自動的にもとまることになってますから、
応力は変位の後についてくるオマケみたいなものです。
変位法の実質的な計算は「変位をもとめる」ことに費やされている、
ということをまず頭に入れておいてください。
で、その「変位」とはなんのことかというと、これは「節点の変位」です。したがって、変位法で取り上げるべき最初の話題は「その節点とはいったい何なのか？」になります。
せっ‐てん【節点】
１ 構造物の骨組み部材の接合点。自由に回転する滑節点、角度の変わらない剛節点がある。
機械工学ならまだしも、土木工学でやる「地盤解析」では、節点の意義が本当にわかりづらい。
単に、人為的にメッシュを切った、縦線と横線の交点という感じしかしない。
まあ、節点とはそういうものなんだけど。
でも、最初の入口段階では、骨組みの接合点で考えたい。
《植物の節の意》
１ 活動・組織などの中心点。集合点。
２ グラフ理論の点と線分からなる図形における、点（頂点、節点）。
路線図の駅、電気回路の素子、ワールドワイドウェブのウェブページに相当する。
３ コンピューターネットワークを構成する個々の要素または中継点。
コンピューター、ハブ、ルーター、さらにはネットワークに接続されたプリンターなどが相当する。
の違いを教えて頂けないでしょうか？
構造物の形が元の状態から変わらなければ、ひずみは生じません。
逆に、形が元の状態から変わると、必ずひずみが生じます。
例として、x-y面内にある正方形板を想像してみましょう。
この正方形板が、力を受けて、xｙ面内で移動したとします。
すると、正方形板内の各点には、すべて変位が生じます。
もうひとつ、正方形板が、力を受けて、ｚ軸回りに回転したとしましょう。
すると、正方形板内の各点には、すべてｚ軸まわりの回転が生じるとともに、ｘｙ面内の変位も生じます。
以上の、変位したり回転したりした正方形板を、元の位置や向きに戻して、元の形状に重ねて見た時、ぴったりと重なれば、その正方形板には、ひずみは生じていません。要するに、”変形”していないのです。
もし、正方形板が力を受けて菱形状などに変形したとすれば、元の形状とはぴったりと重なりません。要するに、その正方形板には、ひずみ（＝変形）が生じています。
要するに、ひずみとは、ある構造物の形状が、元の形から変わった、その変化量を指すのです。
力を受けて移動した後の形状が、移動前の形状とぴったり重ならなければ、ひずみが生じていることになるわけです。
ある構造物に変位や回転が生じているからといって、その構造物にひずみ（＝変形）が生じているとは限りません。
フランス語では”ひずみ”に対応する言葉が、日本語の”変形”に対応するdeformationとなっていて、混乱が起きにくくなっています。ドイツ語など他の欧州言語でも同様です。
ひずみは、形状が元の状態から変わった”変化量”を指すので、定義は、ある２点間の距離の変化量で表すことができます。
要するに、数式で表せば、
ひずみ＝（移動後の距離－移動前の距離）／移動前の距離
などとなります。
ひずみの定義の仕方によって、分子や分母が微妙に変わってきますが、
ひずみ＝長さ／長さ
になることには変わりがありません。
その結果、ひずみは無次元の量になります。
無次元の量ということは、呼び名がないということにもなって、これでは不便です。
実は、角度や上記の回転も同様で、SI単位では無次元量です。しかし、これでは角度の数値を呼ぶときに不便なので、”ラジアン”などという便宜上の単位を設定しているのです。
最も一般的なのが、ひずみの大きさが通常は１よりも非常に小さいことに着目して
0.001 なら、0.1%
など、％で表す流儀です。
もうひとつ、主に測定関係者の間で使用されている流儀があります。
金属材料の通常の使用状態でのひずみは、10^-6のオーダーで表現すると適当な大きさで表すことができるようになることから、10^-6を”マイクロストレイン”などと呼んで、
0.001 なら、1000マイクロ（ストレイン）
などと呼ぶのです。
時計の針ですが、時々刻々と回転を続けます。
要するに、時計の針は、回転と変位を生じ続けています。
言いかえれば、時計の針は、向きと位置が変わり続けます。
しかし、時計の針の形自身は変わるわけではありません。
要するに、時計の針は、時刻の変化によって変位（と回転）は生じますが、ひずみは生じません。
お求めやすい価格にて販売しておりますので、是非ご覧ください！
店舗ホームページ
ttp://https://goodprice2018.shop-pro.jp/ 👀Rock54: Caution(BBR-MD5:fc5433912aa55592f73f2dda4d43bdf8)
みんなを殺した
理由なんて簡単さ
そこに弱いひとたちがいたから
女達は犯され
老人と子供は燃やされた
悪いひとたちはその土地に
家を建てて子供を生んだ
そして街ができ
悪いひとたちの子孫は増え続けた
朝鮮進駐軍 関東大震災 日本人10万人大虐殺
ttp://https://youtu.be/iBIA45CrE30 ttp://https://youtu.be/D0vgxFC04JQ ttp://https://www.youtube.com/watch?v=sYsrzIjKJBc ttp://https://www.youtube.com/watch?v=SiHp41uWo1I ttp://https://www.youtube.com/watch?v=zYBCTRryFP8 ttp://http://www.youtube.com/watch?v=AP4JfKlSdFg ttp://https://ja.m.wikipedia.org/wiki/%E7%B7%B4%E9%A6%AC%E5%8C%BA%E6%8A%95%E7%A5%A8%E6%89%80%E8%A5%B2%E6%92%83%E4%BA%8B%E4%BB%B6 なくなるの？
《超悪質！盗聴盗撮・つきまとい嫌がらせ犯罪首謀者の実名と住所／死ねっ!! 悪魔井口・千明っ!!》
【要注意!! 盗聴盗撮・つきまとい嫌がらせ犯罪工作員】
◎井口・千明（東京都葛飾区青戸６－２３－１６）
※盗聴盗撮・嫌がらせつきまとい犯罪者のリーダー的存在／犯罪組織の一員で様々な犯罪行為に手を染めている
低学歴で醜いほどの学歴コンプレックスの塊／超変態で食糞愛好家である／醜悪で不気味な顔つきが特徴的である
【超悪質！盗聴盗撮・嫌がらせつきまとい犯罪者の実名と住所／井口・千明の子分たち】
①宇野壽倫（東京都葛飾区青戸６－２３－２１ハイツニュー青戸２０２）
※宇野壽倫は過去に生活保護を不正に受給していた犯罪者です／どんどん警察や役所に通報・密告してやってください
②色川高志（東京都葛飾区青戸６－２３－２１ハイツニュー青戸１０３）
※色川高志は現在まさに、生活保護を不正に受給している犯罪者です／どんどん警察や役所に通報・密告してやってください
【通報先】
◎葛飾区福祉事務所（西生活課）
〒１２４－８５５５
東京都葛飾区立石５－１３－１
℡０３－３６９５－１１１１
③清水（東京都葛飾区青戸６－２３－１９）
※低学歴脱糞老女：清水婆婆 ☆☆低学歴脱糞老女・清水婆婆は高学歴家系を一方的に憎悪している☆☆
清水婆婆はコンプレックスの塊でとにかく底意地が悪い／醜悪な形相で嫌がらせを楽しんでいるまさに悪魔のような老婆である
④高添・沼田（東京都葛飾区青戸６－２６－６）
※犯罪首謀者井口・千明の子分／いつも逆らえずに言いなりになっている金魚のフン／親子孫一族そろって低能
⑤高橋（東京都葛飾区青戸６－２３－２３）
⑥長木義明（東京都葛飾区青戸６－２３－２０）
⑦若林豆腐店店主（東京都葛飾区青戸２－９－１４）
⑧肉の津南青戸店店主（東京都葛飾区青戸６－３５ー２
