問題タブ [finite-element-analysis]

割れた岩の立方体のブロックがあります。質問は：

• 上側から下側、または左側から右側への流体の流れをシミュレートする方法は?
• FEA (FEM,...) は唯一の実用的なソリューションですか?
• 上記の質問の最も単純な条件でそうである場合、つまり、流れはフラクチャを介してのみ発生する可能性があります。マトリックスと流体の間に相互作用はありません。など FEA を使って簡単にシミュレーションするにはどうすればよいですか?
• これは、FEA の専門知識を持つ人が数分で実行できる実用的なものですか? 適切なメッシュが既に生成されているとします。
• そうでない場合、そのような単純なケースを解決できるように迅速に開始するには、何をお勧めしますか?
• 同様の問題（フローモデリング）を経験した人はいますか？ある場合、何を使用し、どのように仕事を遂行しましたか?

FEniCS、OpenFoamなど、多くの FEM パッケージが利用可能であることを認識していることに注意してください。

有限要素デルファイライブラリを見つけることができる場所を誰かが教えてもらえますか？

できればオープンソース。

有限要素法(Strang and Fix, 1973)を使用して偏微分方程式を解くことができる R 関数はありますか?

ウィキペディアの有限要素ソフトウェア パッケージのリストに基づいて、メソッドは Fortran、C で実装されています。Fortran と C の両方を R 関数で呼び出すことができますが、CRAN にはそのような関数は見当たりません。

ストランとフィックス。有限要素法の分析。プレンティス・ホール、1973年。

メモリを節約するために、剛性マトリックスをスパース形式で FORTRAN に格納しようとしています。つまり、非ゼロ要素 (irows、icols、A) の 3 つのベクトルを使用しています。これらの配列のサイズを確認したら、次のステップはそれらに値を挿入することです。したがって、ガウス ポイントを使用しています。つまり、ガウス ポイントごとにローカル剛性マトリックスを見つけて、このローカル剛性マトリックスをグローバル (irows、icols、A) マトリックスに挿入します。

この挿入の主な問題は、新しい値がグローバル配列に存在するかどうかを毎回確認する必要があるため、値が存在する場合は新しい値を古い値に追加し、存在しない場合は末尾に追加することです。つまり、配列全体を検索して、値が存在するかどうかを確認する必要があります。これらの配列 (irows、icols、A) のサイズが大きい場合、この検索は計算コストが非常に高くなります。

グローバル剛性マトリックスの各ガウスポイントにローカル剛性マトリックスを挿入するより良い方法を提案できる人はいますか。

行列を作成するためだけに Fortran95 で Sparse Blas を使用したいのですが、ポイント エントリ構造を使用しています。コマンドを使用してマトリックスを作成した後

duscr_begin(n,n,a,istat) を呼び出す

ここで、a は n 行 n 列の行列のハンドルです。値を挿入した後、ハンドルを使用して最終的な行列を確認するにはどうすればよいですか? 他の操作に行列を使用したいので、行列を 3 つのベクトル (スパース) 形式 (row_index、Col_index、Value) で表示したいと考えています。

この Sparse Blas の詳細は第 3 章に記載されており、ここで見ることができます。

http://www.netlib.org/blas/blast-forum/

u4未知数、f1、f2、f3、f5およびf6特定のFEA問題について、この行列システムを解く必要があります。入力を使用して数値を含む行列を生成しましたが、これらの既知の行列を解く必要があります。

• コードは次のように動作する必要があります。

  0 + 0 + 0 + 150*1000*u4 + 0 + 0 = −1732

したがって、u4=-0.01154

• したがって、他の未知数は次のように見つけることができます。

  0 + 0 + 0 + (43:4)(−0:0115) 1000 + 0 + 0 = f1

したがってf1=-500 、など....

私は初心者で、numpy と python OOP の基本に精通しています。このためのコードを設定するのを手伝ってください。

有限要素解析アルゴリズムを実装しようとしています。を解いK u = fて変位uを求め、 でひずみを計算しu、次に応力を計算します。最後に、応力を使用してフォン ミーゼス応力を計算し、これを視覚化します。結果から、ひずみが四面体間で連続していないことがわかりました。

要素として 10 ノードの四面体を使用するため、変位はすべての要素で 2 次多項式になります。変位は、四面体間で連続するように強制する必要があります。そして、変位の一次導関数である歪みは、すべての四面体の内部で連続している必要があります。しかし、よくわかりません。これは、四面体間の境界面全体に当てはまりますか?

元の質問はこちらですが、sfepyと有限要素法を知っている人はここにいると思います。

さて、私はベクトル重力場をモデル化したい、つまり（弱定式化で）

ここで、G はターゲット ベクトル フィールド、rho は密度値です。これは、特定のマテリアルの定数パラメータです。

私の問題定義ファイル:

結果：

レンガに一様に向けられたフィールドを期待していましたが、そこから何かが発散していました（それは大きな問題ではなく、単なるサインの問題だと思います）。また、左下隅の値を増やすことに満足していません。

モデルを修正するにはどうすればよいですか?

前もって感謝します。