問題タブ [finite-element-analysis]

暗黙の有限差分法を使用して、木製シリンダー内の熱伝導をモデル化しようとしています。円筒形と球形に使用している一般的な熱方程式は次のとおりです。

ここで、p は形状係数で、円柱の場合は p = 1、球の場合は p = 2 です。境界条件には、表面での対流が含まれます。モデルの詳細については、以下の Matlab コードのコメントを参照してください。

メインの m ファイルは次のとおりです。

解く連立方程式を作成する関数 m ファイル funcACbar は次のとおりです。

最後に、動的反応を処理する関数 funcY は次のとおりです。

上記のコードを実行すると、木製シリンダーの中心と表面の温度プロファイルが得られます。

このプロットからわかるように、何らかの理由で中心温度と表面温度が 2 分マークで急速に収束しますが、これは正しくありません。

これを修正する方法、または問題を解決するためのより効率的な方法を作成する方法について何か提案はありますか?

連成非線形積分偏微分方程式の膨大なセットがあります。方程式を単純化し、少なくとも半解析的に解こうと長い間試みた結果、私には効率的な数値的方法しか残されていないという結論に達しました。有限要素は、弱い形式の解を与えるガラーキン法に基づいているため、最も従順なように思われるため、最終的に方程式を解くことができることを大いに期待しています。しかし同時に、私はこの分野にまったく慣れていないため、コードをすべてゼロから作成します。

有限要素解析を効率的に実行できる Python ライブラリは既に利用可能ですか?

また、FEniCS/Dolphin のようなソフトウェアが積分微分方程式を解くことにも興味がありましたか?

私たちはエンジニアの学生で、有限要素を使用してシミュレーションを計算する必要があります。私たちのプロジェクトは、オブジェクトの端にあるすべてのノードを見つけて、それらを定義された線と比較することです。Python の経験はなく、a について読んでいます。今週、運がなかったので、今この投稿を書こうとしています。生の日付は次のようになります

必要なことは、3 列目から 6 列目に 2 回または 3 回出現する数字を検索し、それらを別のテキスト ドキュメントに出力することです。その後、同じテキスト ドキュメントでこれらのノードの座標を見つける必要があります。これらは次のようになります。

すべてのエッジ座標が得られたので、シミュレーションのエッジを必要なエッジと比較できます。

Python でのプログラミングを学びたいので、最終的な解決策は必要ありません。また、試験でこれを説明できなければなりません。適切な文献や、再生できる基本的なスクリプトを誰かが教えてくれることを願っています。このフォーラムで私たちのニーズに合ったものが見つからなかったので、

前もって感謝します

ABAQUS と組み合わせた Python コードに matplotlib を使用している方がいるかどうか疑問に思っています。Python 2.6を搭載したABAQUS 6.12を使用しています。また、Windows 64 ビットを使用しています。

matplotlib をダウンロードしましたが、matplotlib を ABAQUS-Python スクリプト インターフェイスにインポートすると、次のエラーが発生します。

ImportError: matplotlib には dateutil が必要です

誰か助けてくれませんか？これを修正するにはどうすればよいですか?

どんな助けでも大歓迎です！ありがとう

注: FEA ソフトウェアに慣れていない方のために説明すると、Abaqus Scripting Interface は Python オブジェクト指向プログラミング言語の拡張であり、ABAQUS/CAE に付属しています。

有限要素を使用して、2d 長方形のポアソン方程式の解をコーディングしています。コードを単純化するために、基底関数へのハンドルを配列に格納し、これらの基底関数をループして、行列と右辺を作成します。これに関する問題は、非常に粗いグリッドの場合でも、法外に遅いことです。9x9 グリッド (ディリクレ BC を使用すると、解決するノードが 49 あります) の場合、約 20 秒かかります。プロファイルを使用して、約半分の時間が基底関数へのアクセス (実行ではない) に費やされていることに気付きました。

プロファイラーによるmatrix_assembly>@(x,y)bilinearBasisFunction(x,y,xc(k-1),xc(k),xc(k+1),yc(j-1),yc(j),yc(j+1)) (156800 calls, 11.558 sec)と、自己時間 (双一次基底コードを実行していない) は 9 秒を超えています。なぜこれがとても遅いのかについてのアイデアはありますか?

コードの一部を次に示します。必要に応じてさらに投稿できます。

どの方法を試しても、Cgal-bindigs のインストールが停止し、次の問題が発生します。

すべての前提条件をインストールしたと確信しています。エラーがわかりません:S. 役立つ場合は、「make VERBOSE=1」の出力を貼り付けます: http://pastebin.com/V1e9knAh 生物医学の調査と有限要素解析の基礎となる四面体化に使用する必要があるのを助けてください。