1, はじめに 電磁導波路の伝達問題のシミュレータを実装したので、今回は電子レンジの電磁波シミュレーションを行います。 2. モデル 2.45GHzのTE10モードを導波管を通して給電したときの、電子レンジ空洞内の電界分布を計算しています。 電子レンジ内には寒…

1. はじめに 2次四面体辺要素を構成し、3次元導波路伝達問題を解きます。 2. 定式化 pdfにまとめました。 2nd Order Tetrahedral Edge Element http://starlightparade.usamimi.info/ivyfem/doc/2ndOrderTetrahedralEdgeElement.pdf?p=0 面内パラメータは隣…

1. はじめに 辺要素FEMを用いてマイクロストリップ線路(micro-strip line)のエアーブリッジ(air bridge)の伝達問題を計算しています。 2. 定式化 過去の記事を参照してください。 3. 計算結果 文献と一致しています。 4. まとめ マイクロストリップ線路のエ…

1. はじめに 四面体辺要素(tetrahedral edge element)を使ってマイクロストリップ線路(micro-strip line)のステップ不連続(step discontinuity)およびテーパー不連続(taper discontinuity)の反射・透過特性を計算する。 2. 定式化 次の記事を参照ください。 …

1. はじめに 3次元の電磁導波路(electro-magnetic waveguide)の伝達問題（不連続問題）を解いています。 伝達問題とは回路に入力したときの反射、透過応答を計算する問題です。回路素子は一様導波路でなく不連続な形状になるため不連続問題ともいわれます。 …

1. はじめに 富岳のようなことがしたい。 咳による飛沫のシミュレーションが活発に行われています。 cough droplet finite elementで検索して文献を調べるといくつか定式化されているものを見つけることができました。 そのどれもが乱流を仮定して速度場を計…

1. はじめに Finite Calculus (FIC) Formulationは、無限小領域に対して成立する支配方程式を有限領域に対応するように再定式化したものです。 支配方程式に新たな項が追加されます。これは安定化数値解法の一種です。 2. FICの定式化 下記の下線を引いた部…

1. はじめに これまで弾性波プレート導波路の伝達問題を周波数領域FEMで解いてきました。 今回時間領域FEM(Time Domain Finite Element Method)でLamb wave弾性波導波路の伝達問題を解いています。 解析端の境界条件にはPML(Perfectly Matched Layers)を用い…

1. はじめに 前にSH wave弾性波導波路の伝達問題を計算しましたが、そのときポート境界条件は固有モード展開を用いました。 本記事ではポート境界条件にPMLを装荷した場合と、ABC(吸収境界条件, Absorbing Boundary Conditions)を用いた場合の２つの定式化を…

1. はじめに これまでにLamb waveの弾性波プレート導波路の伝達問題(transmission problem, 不連続問題 discontinuity problem, 散乱問題 scattering problem)について半無限長導波路を表現するために固有モード展開を用いる方法とPML(Perfectly Matched Lay…

1. はじめに 前に弾性波導波路(elastic waveguides)のLamb wave伝達問題(transmission problem, 不連続問題 discontinuity problem、散乱問題 scattering problem)をFEMで解きましたが、ポート境界は固有モード展開を用いて表現しました。 本記事では、ポー…

1. はじめに 前記事で弾性波プレート導波路を伝搬するLamb waveの伝達問題を定式化、計算しました。 本記事ではSH wave(Shear Horizontal wave)の伝達問題を定式化、計算します。 2. SH waveの伝達問題の定式化 pdfにまとめました。 Frequency Domain FEM Fo…

1. はじめに 弾性波プレート導波路を伝搬するLamb waveの基本モードをポート1から入射したとき、不連続部により反射する、またはポート2に透過するのをFEMで計算しています。 基本モードは対称な分布のモードのため、中央面(mid-plane)を境に半分にした領域…

1. はじめに 変位電流が無視できるような周波数帯は準定常場とよばれ、渦電流場問題で使われます。 支配方程式は、ベクトルポテンシャルAとスカラーポテンシャルφで表したものが使われることが多いようです(A-φ法)。ポテンシャルにはなんらかのゲージを課す…

1. はじめに 前記事ではMITC3シェル要素を用いた幾何学的非線形な変形の計算を行いました。 本記事では非圧縮性超弾性体(Mooney-Rivlinモデル)をMITC3シェル要素で解いています。 前記事ではSt.Venant-Kirchhoffモデルを用いましたが、厚さは変化しないもの…

1. はじめに これまで薄板の曲げの計算にDKT(discrete Kirchhoff Theory)要素、Mindlin薄板要素を用いました。これらはKirchhoff-Love板理論やReisner-Mindlin板理論といった特別な仮定のもとで定式化したものです。 一方3次元のソリッド要素、例えば三角柱…

1. はじめに 前記事ではDiscrete Kirchhoff板理論による三角形要素(DKT要素)を使いましたが、今回はMindlin板理論による三角形要素を使って薄板の曲げを計算します。 2. 定式化 pdfにまとめました。 Reissner-Mindlin Plate Bending Triagular Elements http…

1. はじめに 薄板の曲げの理論にはKirchhoffの板理論があります。ただ、この理論に基づく2次元要素に有効なものはないとのことです。 Kirchhoffの板理論を一般化したReissner-Mindlinの板理論を基礎とし、2次三角形要素の6節点上でのみKirchhoffの仮定を満た…

1. はじめに 幾何学的非線形梁をTotal Lagrange法を用いて計算します。 Euler-Bernoulli梁の場合は、変位と回転角に関係式が存在するため、変位２つ(uX、uY)と回転角(θ)は自由度について非線形な補間関数を構成する必要があります。これをfield-consistentな…

1. はじめに 幾何学的非線形な梁をCo-rotational法を用いて計算しています。 Co-rotational法は変位を剛体運動と微小変形に分離して考えます。 Euler-Bernoulli Shallow Arch梁とTimoshenko梁について定式化しています。 2. Co-rotational Formulation pdfに…

1. はじめに フレーム要素を下記のような梁＋トラス構造の解析に使いましたが、間違っていることに気付いたので記します。 2. 複数のビームが接続する節点の連続性 ３つのビームが1つの節点に接続されている場合、X方向変位u、Y方向変位vは連続としていいの…

1. はじめに 前記事で1次元構造要素のトラス、ビーム、フレーム要素によるFEM変位計算を行いましたが、そこで用いたビーム、フレームはEuler-Bernoulli梁をもとにしていました。 本記事ではTimoshenko梁のビーム、フレーム要素を用いて計算しています。 2.Eu…

1. はじめに ビーム要素とフレーム要素を実装したので、両端支持梁の固有振動を計算しました。 比較するために線形弾性体、Saint Venant超弾性体でも固有振動を計算しています。 2.両端支持梁 棒の断面は、 b=0.2 m h = (1/4)b 材料定数は、 E=169 GPa ν=0.2…

1. はじめに 円柱や角柱のような棒をつなげて構造物を作るとき、柱の断面が長手方向に比べて十分小さい場合は、1次元の梁理論が適用できます。 弾性体を三角形要素で分割して曲げなどを計算しましたが、梁要素を使えば1次元として扱うことができ、線要素で分…

1. はじめに 前に遮蔽形導波路のベクトル波固有モードをedge/nodal FEMを用いて計算しました。 今回は遮蔽ボックスを取り払った場合、つまり開放形(Open-type, unbounded)の導波路をFEMを用いて計算しています。 開放形導波路の場合は、無限領域を内部の有限…

1. はじめに Navier-Stokes方程式は速度ベクトルvと静水圧pを解く方程式でしたが、対流を表すvorticityベクトル ωと流線を表すstream functionベクトル ψを用いた定式化もできます。 2次元問題の場合、vorticity、stream functionは１方向しか成分を持たない…

1. はじめに Navier-Stokes方程式は速度ベクトルvと静水圧pを解く方程式でしたが、対流を表すvorticityベクトル ωと流線を表すstream functionベクトル ψを用いた定式化もできます。 2次元問題の場合、vorticity、stream functionは１方向しか成分を持たない…

流体力学の支配方程式、Navier-Stokes方程式を導く過程をpdfにまとめました。 以下の順に説明しています。 3 時間導関数 4 Eular の膨張公式 5 Reynolds の輸送定理 6 連続の式 7 Cauchy の第一運動法則 8 Cauchy の第二運動法則 9 非圧縮性の式 10 流体の構…

1. はじめに 以前、Navier-Stokes方程式 + div vの2式をガラーキン法で定式化し、Newmark β法とRK4(4次Runge Kutta)で解きました。またNewmark β法については粘性係数が小さい場合に対応するためSUPGの定式化も行いました。 今回、Navier-Stokes方程式と、Pr…

1. はじめに 前記事で導波路ベクトル波(full-wave)固有モードをedge/nodal 三角形要素を用いて計算しましたが、用いたedge/nodal elementsは基本要素(1次要素)でした。 今回2次のedge/nodal elementsを用いてベクトル波導波路の固有モードを計算しています。…