JP6034909B2 - 異方性導電物質の電磁界解析方法 - Google Patents
異方性導電物質の電磁界解析方法 Download PDFInfo
- Publication number
- JP6034909B2 JP6034909B2 JP2015062483A JP2015062483A JP6034909B2 JP 6034909 B2 JP6034909 B2 JP 6034909B2 JP 2015062483 A JP2015062483 A JP 2015062483A JP 2015062483 A JP2015062483 A JP 2015062483A JP 6034909 B2 JP6034909 B2 JP 6034909B2
- Authority
- JP
- Japan
- Prior art keywords
- lattice
- electromagnetic field
- degree
- calculation
- along
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims description 61
- 230000005672 electromagnetic field Effects 0.000 title claims description 41
- 239000004020 conductor Substances 0.000 title claims description 17
- 230000005684 electric field Effects 0.000 claims description 44
- 230000002093 peripheral effect Effects 0.000 claims description 23
- 239000000835 fiber Substances 0.000 claims description 22
- 238000000034 method Methods 0.000 claims description 12
- 229920002430 Fibre-reinforced plastic Polymers 0.000 claims description 3
- 239000011151 fibre-reinforced plastic Substances 0.000 claims description 3
- 239000011347 resin Substances 0.000 claims description 3
- 229920005989 resin Polymers 0.000 claims description 3
- 239000002131 composite material Substances 0.000 description 15
- 238000006243 chemical reaction Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 9
- 238000010521 absorption reaction Methods 0.000 description 5
- 239000002648 laminated material Substances 0.000 description 5
- 239000004918 carbon fiber reinforced polymer Substances 0.000 description 4
- 239000000463 material Substances 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 229920000049 Carbon (fiber) Polymers 0.000 description 1
- 240000006829 Ficus sundaica Species 0.000 description 1
- 239000004917 carbon fiber Substances 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000010030 laminating Methods 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0807—Measuring electromagnetic field characteristics characterised by the application
- G01R29/0814—Field measurements related to measuring influence on or from apparatus, components or humans, e.g. in ESD, EMI, EMC, EMP testing, measuring radiation leakage; detecting presence of micro- or radiowave emitters; dosimetry; testing shielding; measurements related to lightning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Electromagnetism (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Measuring Magnetic Variables (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Complex Calculations (AREA)
Description
本発明は、異方性導電率を有する異方性導電物質の電磁界解析方法に関する。
航空機の機体設計においては、雷撃によるスパーク発生に起因する燃料引火などを防止すべく、落雷対策が必要とされている。この落雷対策を講じるうえでは、雷撃時における機体中の電流分布を明らかにすることが重要であり、そのための高精度な電磁界解析手法が望まれている。
この種の電磁界解析では、ＦＤＴＤ法（Finite-difference time-domain method；時間領域差分法）が広く用いられている（例えば、特許文献１参照）。ＦＤＴＤ法は、直交する３軸の電磁界成分（電界と磁界）を、Ｙｅｅ格子と呼ばれる立方体の計算格子によって計算する手法である。
ところで、航空機の機体構造材においては、比強度等に優れる炭素繊維強化プラスチック（ＣＦＲＰ）などの複合材料の適用が進んでいる。なかでも、例えば図９に示すように、繊維方向が互いに４５°異なる複数の繊維層を積層させた擬似等方積層材が広く用いられてきている。
このような複合材料の電磁界解析では、各繊維が導電性を有するため、繊維方向に沿った導電率の設定が必要となる。つまり、０°方向や９０°方向だけでなく、４５°方向に沿った異方性導電率も模擬する必要がある。
このような複合材料の電磁界解析では、各繊維が導電性を有するため、繊維方向に沿った導電率の設定が必要となる。つまり、０°方向や９０°方向だけでなく、４５°方向に沿った異方性導電率も模擬する必要がある。
しかしながら、Ｙｅｅ格子を単純に用いたＦＤＴＤ法では、直交３軸に沿った導電率しか設定できないため、０°（９０°）方向に沿ってＹｅｅ格子を配列させた場合に、４５°方向に沿った異方性導電率を設定することができない。そのため、この異方性導電率を有する異方性導電物質の電磁界解析を精度良く行うことは困難であった。
本発明は、上記課題を鑑みてなされたもので、異方性導電率を有する異方性導電物質の電磁界解析を高精度に行うことを目的とする。
上記目的を達成するために、請求項１に記載の発明は、所定の第一の方向に沿った導電率が他の方向のものとは異なる一の層と、前記第一の方向とは異なる第二の方向に沿った導電率が他の方向のものとは異なる他の層とが積層されてなる異方性導電物質の電磁界解析方法であって、
コンピュータが、
前記第一の方向に沿った辺を有する第一の計算格子により前記一の層を格子分割するとともに、前記第二の方向に沿った辺を有する第二の計算格子により前記他の層を格子分割し、
前記第一の格子における電界及び磁界のうちの一方の電磁界成分の分布から前記第二の格子における当該一方の電磁界成分を補間計算し、
前記第二の格子における電界及び磁界のうちの他方の電磁界成分の分布から前記第一の格子における当該他方の電磁界成分を補間計算することを特徴とする。
コンピュータが、
前記第一の方向に沿った辺を有する第一の計算格子により前記一の層を格子分割するとともに、前記第二の方向に沿った辺を有する第二の計算格子により前記他の層を格子分割し、
前記第一の格子における電界及び磁界のうちの一方の電磁界成分の分布から前記第二の格子における当該一方の電磁界成分を補間計算し、
前記第二の格子における電界及び磁界のうちの他方の電磁界成分の分布から前記第一の格子における当該他方の電磁界成分を補間計算することを特徴とする。
請求項２に記載の発明は、請求項１に記載の異方性導電物質の電磁界解析方法において、
コンピュータが、
解析空間の周縁境界を滑らかな面状に構成しつつ、当該解析空間の周縁部に前記第一の格子及び前記第二の格子のいずれか一方のみを配列することを特徴とする。
コンピュータが、
解析空間の周縁境界を滑らかな面状に構成しつつ、当該解析空間の周縁部に前記第一の格子及び前記第二の格子のいずれか一方のみを配列することを特徴とする。
請求項３に記載の発明は、請求項１または２に記載の異方性導電物質の電磁界解析方法において、
前記一の層及び前記他の層は、いずれもｘｙ平面内に配置されるとともに、ｚ軸方向に互いに積層され、
前記第一の格子は、ｘｙｚの直交３軸に沿った立方体状の計算格子であり、
前記第二の格子は、前記第一の格子に内接する大きさのものであって、ｘｙ平面内で±４５°方向に沿った２辺とｚ軸に沿った１辺とからなる直方体状の計算格子であることを特徴とする。
前記一の層及び前記他の層は、いずれもｘｙ平面内に配置されるとともに、ｚ軸方向に互いに積層され、
前記第一の格子は、ｘｙｚの直交３軸に沿った立方体状の計算格子であり、
前記第二の格子は、前記第一の格子に内接する大きさのものであって、ｘｙ平面内で±４５°方向に沿った２辺とｚ軸に沿った１辺とからなる直方体状の計算格子であることを特徴とする。
請求項４に記載の発明は、請求項１〜３のいずれか一項に記載の異方性導電物質の電磁界解析方法において、
前記異方性導電物質は、前記第一の方向及び前記第二の方向に沿った導電性の繊維で樹脂が強化された繊維強化プラスチックであることを特徴とする。
前記異方性導電物質は、前記第一の方向及び前記第二の方向に沿った導電性の繊維で樹脂が強化された繊維強化プラスチックであることを特徴とする。
請求項１に記載の発明によれば、第一の方向に沿った辺を有する第一の計算格子により一の層が格子分割されるとともに、第一の方向とは異なる第二の方向に沿った辺を有する第二の計算格子により他の層が格子分割される。これにより、一の層における第一の方向に沿った導電率と、他の層における第二の方向に沿った導電率とを個別に設定することができる。
また、第一の格子における電界及び磁界のうちの一方の電磁界成分の分布から第二の格子における当該一方の電磁界成分が補間計算され、第二の格子における電界及び磁界のうちの他方の電磁界成分の分布から第一の格子における当該他方の電磁界成分が補間計算される。これにより、第一の格子と第二の格子との間での電磁波伝搬計算をシームレスに実行することができる。
したがって、異方性導電率を有する一の層と他の層とが積層されてなる異方性導電物質の電磁界解析を、高精度に行うことができる。
また、第一の格子における電界及び磁界のうちの一方の電磁界成分の分布から第二の格子における当該一方の電磁界成分が補間計算され、第二の格子における電界及び磁界のうちの他方の電磁界成分の分布から第一の格子における当該他方の電磁界成分が補間計算される。これにより、第一の格子と第二の格子との間での電磁波伝搬計算をシームレスに実行することができる。
したがって、異方性導電率を有する一の層と他の層とが積層されてなる異方性導電物質の電磁界解析を、高精度に行うことができる。
請求項２に記載の発明によれば、解析空間の周縁境界が滑らかな面状に構成されつつ、当該解析空間の周縁部に第一の格子及び第二の格子のいずれか一方のみが配列される。
これにより、計算実行時における周縁境界での発散や反射の発生を抑え、有限空間である解析空間を無限空間として取り扱う吸収境界の計算を成立させることができる。
これにより、計算実行時における周縁境界での発散や反射の発生を抑え、有限空間である解析空間を無限空間として取り扱う吸収境界の計算を成立させることができる。
以下、本発明の実施形態について、図面を参照して説明する。
まず、本実施形態における異方性導電物質の電磁界解析方法（以下、単に「電磁界解析方法」という）の解析モデルについて説明する。
図１は、本実施形態における電磁界解析方法の解析モデルを示す図であり、図２は、計算格子である後述の９０度格子Ｍ１及び４５度格子Ｍ２を示す図であり、図３は、後述の解析空間Ｓの周縁部での格子分割を説明するための図である。
図１は、本実施形態における電磁界解析方法の解析モデルを示す図であり、図２は、計算格子である後述の９０度格子Ｍ１及び４５度格子Ｍ２を示す図であり、図３は、後述の解析空間Ｓの周縁部での格子分割を説明するための図である。
本実施形態における電磁界解析方法は、異方性導電率を有する物質（本実施形態では複合材料ＣＭ）における電流分布等の電磁的特性を、ＦＤＴＤ法（Finite-difference time-domain method；時間領域差分法）を用いた電磁界解析ソフトにより算出するものである。
その解析対象である複合材料ＣＭは、本実施形態においては、樹脂を導電性の炭素繊維で強化した炭素繊維強化プラスチック（ＣＦＲＰ）である。より詳しくは、図１に示すように、複合材料ＣＭは、繊維方向が互いに４５°異なる４層の繊維層を２セット（すなわち合計８層）積層させた略平板状の擬似等方積層材であり、本実施形態では、８つの繊維層（第１層Ｌ１〜第８層Ｌ８）が解析空間Ｓにおけるｘｙ平面内に配置されつつｚ方向に互いに積層されている。
本実施形態における電磁界解析方法では、複合材料ＣＭとその周囲の十分に広い空間とを含む解析空間Ｓ全体がモデル化される。
具体的には、図２（ａ），（ｂ）に示すように、９０度格子Ｍ１と４５度格子Ｍ２との２種類の計算格子を用いて解析空間Ｓ内が格子分割される。
このうち、９０度格子Ｍ１は、各辺がｘｙｚの直交３軸に沿った立方体状の計算格子であり、つまり、当該９０度格子Ｍ１では、ｘｙ平面内の２辺が、ｘ軸を基準とする０°（ｘ軸）及び９０°（ｙ軸）の２方向に沿っている。
一方、４５度格子Ｍ２は、９０度格子Ｍ１に内接する大きさのものであって、ｘｙ平面内で±４５°方向に沿った２辺とｚ軸に沿った１辺とからなる直方体状の計算格子である。
具体的には、図２（ａ），（ｂ）に示すように、９０度格子Ｍ１と４５度格子Ｍ２との２種類の計算格子を用いて解析空間Ｓ内が格子分割される。
このうち、９０度格子Ｍ１は、各辺がｘｙｚの直交３軸に沿った立方体状の計算格子であり、つまり、当該９０度格子Ｍ１では、ｘｙ平面内の２辺が、ｘ軸を基準とする０°（ｘ軸）及び９０°（ｙ軸）の２方向に沿っている。
一方、４５度格子Ｍ２は、９０度格子Ｍ１に内接する大きさのものであって、ｘｙ平面内で±４５°方向に沿った２辺とｚ軸に沿った１辺とからなる直方体状の計算格子である。
この格子分割では、複合材料ＣＭの各層Ｌ１〜Ｌ８に対し、その繊維方向に沿って導電率を設定できるように、当該繊維方向に沿った辺を有する計算格子で格子分割する。つまり、複合材料ＣＭのうち、繊維方向が０°または９０°である第２，４，５，７層Ｌ２，Ｌ４，Ｌ５，Ｌ７を９０度格子Ｍ１で格子分割し、繊維方向が４５°または−４５°である第１，３，６，８層Ｌ１，Ｌ３，Ｌ６，Ｌ８を４５度格子Ｍ２で格子分割する。
またこの格子分割では、解析空間Ｓのうち、複合材料ＣＭの周囲の空間を、複合材料ＣＭと連なる計算格子で格子分割する。つまり、例えば複合材料ＣＭの第１層Ｌ１とｘｙ平面内で連なる周囲の空間は、当該第１層Ｌ１と同じ４５度格子Ｍ２で格子分割される。
但し、解析空間Ｓの周縁部では、図３に示すように、ｘｙｚいずれの方向においても９０度格子Ｍ１が所定列数（本実施形態では３列）だけ配列されている。このとき、９０度格子Ｍ１と４５度格子Ｍ２との境界（以下、「格子境界」という）Ｂ１では、４５度格子Ｍ２の各節点が９０度格子Ｍ１の各辺の中点に位置するように、互いの一部が重複している。
但し、解析空間Ｓの周縁部では、図３に示すように、ｘｙｚいずれの方向においても９０度格子Ｍ１が所定列数（本実施形態では３列）だけ配列されている。このとき、９０度格子Ｍ１と４５度格子Ｍ２との境界（以下、「格子境界」という）Ｂ１では、４５度格子Ｍ２の各節点が９０度格子Ｍ１の各辺の中点に位置するように、互いの一部が重複している。
そのため、解析空間Ｓの周縁部では９０度格子Ｍ１が一様に配列され、当該解析空間Ｓの周縁部の境界（以下、「周縁境界」という）Ｂ２が、ｘｙｚの各軸に沿った滑らかな面状に構成される。これにより、周縁境界Ｂ２における吸収境界の計算が成立する，すなわち、有限空間である解析空間Ｓをあたかも無限空間であるかのように取り扱うことができる。
ｘｙ平面における解析空間Ｓの周縁部に４５度格子Ｍ２を単純に配列してしまうと、周縁境界Ｂ２が滑らかな面状に構成されず（ギザギザになる）、周縁境界Ｂ２に電磁波が垂直に入射しなくなるため、計算実行時に当該周縁境界Ｂ２で発散や反射が発生してしまう。そこで、本実施形態では、４５度格子Ｍ２で格子分割された領域の周囲に９０度格子Ｍ１を配列させて周縁境界Ｂ２を滑らかな面状に構成することで、このような発散や反射の発生を抑え、吸収境界の計算を成立させている。
ｘｙ平面における解析空間Ｓの周縁部に４５度格子Ｍ２を単純に配列してしまうと、周縁境界Ｂ２が滑らかな面状に構成されず（ギザギザになる）、周縁境界Ｂ２に電磁波が垂直に入射しなくなるため、計算実行時に当該周縁境界Ｂ２で発散や反射が発生してしまう。そこで、本実施形態では、４５度格子Ｍ２で格子分割された領域の周囲に９０度格子Ｍ１を配列させて周縁境界Ｂ２を滑らかな面状に構成することで、このような発散や反射の発生を抑え、吸収境界の計算を成立させている。
続いて、９０度格子Ｍ１と４５度格子Ｍ２との間での電界及び磁界の変換計算方法について説明する。
図４（ａ），（ｂ）は、この変換計算を説明するための概念図であり、図５（ａ）は、電界の変換計算を説明するための図であり、図５（ｂ）は、磁界の変換計算を説明するための図である。
図４（ａ），（ｂ）は、この変換計算を説明するための概念図であり、図５（ａ）は、電界の変換計算を説明するための図であり、図５（ｂ）は、磁界の変換計算を説明するための図である。
図４（ａ）に示すように、本実施形態では、９０度格子Ｍ１から４５度格子Ｍ２へ磁界が不可逆的に伝搬するものとして、９０度格子Ｍ１での磁界分布から、４５度格子Ｍ２での磁界が補間計算される。また、各格子Ｍ１，Ｍ２では、電界及び磁界が可逆的に伝搬するものとする。
一方、４５度格子Ｍ２から９０度格子Ｍ１へは、図４（ｂ）に示すように、電界が不可逆的に伝搬するものとして、４５度格子Ｍ２での電界分布から、９０度格子Ｍ１での電界が補間計算される。また、各格子Ｍ１，Ｍ２では、電界及び磁界が可逆的に伝搬するものとする。
一方、４５度格子Ｍ２から９０度格子Ｍ１へは、図４（ｂ）に示すように、電界が不可逆的に伝搬するものとして、４５度格子Ｍ２での電界分布から、９０度格子Ｍ１での電界が補間計算される。また、各格子Ｍ１，Ｍ２では、電界及び磁界が可逆的に伝搬するものとする。
より詳しくは、図５（ａ）に示すように、格子間での電界の変換計算は、以下のように行われる。
９０度格子Ｍ１と４５度格子Ｍ２とが重複した格子境界Ｂ１では、４５度格子Ｍ２の電界が９０度格子Ｍ１の磁界と４５度格子Ｍ２の磁界から計算される（ア部）とともに、９０度格子Ｍ１の電界が４５度格子Ｍ２の電界から補間計算される（イ部）。
但し、ｚ軸に沿った４５度格子Ｍ２の電界や、ｘ軸に沿った９０度格子Ｍ１の電界は、各格子の磁界から計算される（ウ，エ部）。また、格子境界Ｂ１以外の部分における電界も、各格子の磁界から計算される（オ部）。
９０度格子Ｍ１と４５度格子Ｍ２とが重複した格子境界Ｂ１では、４５度格子Ｍ２の電界が９０度格子Ｍ１の磁界と４５度格子Ｍ２の磁界から計算される（ア部）とともに、９０度格子Ｍ１の電界が４５度格子Ｍ２の電界から補間計算される（イ部）。
但し、ｚ軸に沿った４５度格子Ｍ２の電界や、ｘ軸に沿った９０度格子Ｍ１の電界は、各格子の磁界から計算される（ウ，エ部）。また、格子境界Ｂ１以外の部分における電界も、各格子の磁界から計算される（オ部）。
また、格子間での磁界の変換計算は、図５（ｂ）に示すように、以下のように行われる。
９０度格子Ｍ１と４５度格子Ｍ２とが重複した格子境界Ｂ１では、４５度格子Ｍ２の磁界が９０度格子Ｍ１の磁界から補間計算される（カ部）。
但し、ｚ軸に沿った４５度格子Ｍ２の磁界や、９０度格子Ｍ１の磁界は、各格子の磁界から計算される（キ，ク部）。このうち、ｚ軸に沿った９０度格子Ｍ１の磁界は、４５度格子Ｍ２の磁界としても利用される（ケ部）。また、格子境界Ｂ１以外の部分における磁界は、各格子の電界から計算される（コ部）。
９０度格子Ｍ１と４５度格子Ｍ２とが重複した格子境界Ｂ１では、４５度格子Ｍ２の磁界が９０度格子Ｍ１の磁界から補間計算される（カ部）。
但し、ｚ軸に沿った４５度格子Ｍ２の磁界や、９０度格子Ｍ１の磁界は、各格子の磁界から計算される（キ，ク部）。このうち、ｚ軸に沿った９０度格子Ｍ１の磁界は、４５度格子Ｍ２の磁界としても利用される（ケ部）。また、格子境界Ｂ１以外の部分における磁界は、各格子の電界から計算される（コ部）。
図６は、ｚ軸方向の電界Ｅｚの格子間に亘る伝搬態様を例示したイメージ図である。
本実施形態における電磁界解析方法では、９０度格子Ｍ１での電界Ｅｚの電磁波は、図６に示すように４５度格子Ｍ２へ伝搬する。
このとき、格子境界Ｂ１では、ｚ軸に沿った４５度格子Ｍ２の電界が、補間計算で求められた４５度格子Ｍ２の磁界から計算される（サ部）。また、４５°方向に沿った４５度格子Ｍ２の磁界は、９０度格子Ｍ１の磁界から補間計算され、４５度格子Ｍ２の電界からは計算されない（シ部）。また、ｚ軸に沿った９０度格子Ｍ１の電界は、４５度格子Ｍ２の電界から補間計算される（ス部）。
本実施形態における電磁界解析方法では、９０度格子Ｍ１での電界Ｅｚの電磁波は、図６に示すように４５度格子Ｍ２へ伝搬する。
このとき、格子境界Ｂ１では、ｚ軸に沿った４５度格子Ｍ２の電界が、補間計算で求められた４５度格子Ｍ２の磁界から計算される（サ部）。また、４５°方向に沿った４５度格子Ｍ２の磁界は、９０度格子Ｍ１の磁界から補間計算され、４５度格子Ｍ２の電界からは計算されない（シ部）。また、ｚ軸に沿った９０度格子Ｍ１の電界は、４５度格子Ｍ２の電界から補間計算される（ス部）。
図７は、格子間での補間計算を説明するための図である。
９０度格子Ｍ１の磁界分布から４５度格子Ｍ２の磁界を求める補間計算では、９０度格子Ｍ１内で磁界が直線的に分布するものと仮定して、以下の式１〜式４により線形補間を行う。
ここで、図７（ａ）に示すように、Ｈｐ，Ｈｑは、４５度の座標系をｐｑｒの直交３軸とする４５度格子Ｍ２の磁界であり、Ｈｘ，Ｈｙは、９０度格子Ｍ１の磁界である。
９０度格子Ｍ１の磁界分布から４５度格子Ｍ２の磁界を求める補間計算では、９０度格子Ｍ１内で磁界が直線的に分布するものと仮定して、以下の式１〜式４により線形補間を行う。
また、４５度格子Ｍ２の電界分布から９０度格子Ｍ１の電界を求める補間計算では、４５度格子Ｍ２内で電界が直線的に分布するものと仮定して、以下の式５及び式６により線形補間を行う。
ここで、図７（ｂ）に示すように、Ｅｙ，Ｅｚは、９０度格子Ｍ１の電界であり、Ｅｐ，Ｅｑ，Ｅｒは、４５度の座標系をｐｑｒの直交３軸とする４５度格子Ｍ２の電界である。
続いて、格子境界Ｂ１を跨ぐ電磁波の伝搬態様について、解析例を示して説明する。
この解析例では、中心部を４５度格子Ｍ２で格子分割し、その周囲を９０度格子Ｍ１で格子分割した解析モデルの中心に電磁波パルスを印加し、この電磁波が格子境界Ｂ１を跨いで伝搬する様子を確認した。
この解析例では、中心部を４５度格子Ｍ２で格子分割し、その周囲を９０度格子Ｍ１で格子分割した解析モデルの中心に電磁波パルスを印加し、この電磁波が格子境界Ｂ１を跨いで伝搬する様子を確認した。
図８は、この解析例の結果を示す図であって、電磁波パルス印加から１０ｎｓ後，５０ｎｓ後，１００ｎｓ後におけるｘｙ平面での電界強度を示す図である。
この図に示すように、本解析例では、印加された電磁波パルスによる電磁波が、二点鎖線で示す格子境界Ｂ１を跨いで４５度格子Ｍ２から９０度格子Ｍ１へ滑らかに伝搬していることが分かる。
また、解析空間の周縁境界Ｂ２においても、発散や反射が発生することなく、好適に吸収境界の計算を実行できていることが分かる。
この図に示すように、本解析例では、印加された電磁波パルスによる電磁波が、二点鎖線で示す格子境界Ｂ１を跨いで４５度格子Ｍ２から９０度格子Ｍ１へ滑らかに伝搬していることが分かる。
また、解析空間の周縁境界Ｂ２においても、発散や反射が発生することなく、好適に吸収境界の計算を実行できていることが分かる。
以上のように、本実施形態によれば、ｘｙ平面内で０°方向及び９０°方向に沿った辺を有する９０度格子Ｍ１により、複合材料ＣＭのうちの第２，４，５，７層が格子分割されるとともに、ｘｙ平面内で±４５°方向に沿った辺を有する４５度格子Ｍ２により、残りの第１，３，６，８層が格子分割される。これにより、第２，４，５，７層における０°方向及び９０°方向に沿った導電率と、第１，３，６，８層における±４５°方向に沿った導電率とを個別に設定することができる。
また、９０度格子Ｍ１における磁界分布から４５度格子Ｍ２における磁界が補間計算され、４５度格子Ｍ２における電界分布から９０度格子Ｍ１における電界が補間計算される。これにより、９０度格子Ｍ１と４５度格子Ｍ２との間での電磁波伝搬計算をシームレスに実行することができる。
したがって、異方性導電率を有する第２，４，５，７層と第１，３，６，８層とが積層されてなる複合材料ＣＭの電磁界解析を、高精度に行うことができる。
また、９０度格子Ｍ１における磁界分布から４５度格子Ｍ２における磁界が補間計算され、４５度格子Ｍ２における電界分布から９０度格子Ｍ１における電界が補間計算される。これにより、９０度格子Ｍ１と４５度格子Ｍ２との間での電磁波伝搬計算をシームレスに実行することができる。
したがって、異方性導電率を有する第２，４，５，７層と第１，３，６，８層とが積層されてなる複合材料ＣＭの電磁界解析を、高精度に行うことができる。
また、解析空間Ｓの周縁境界Ｂ２が滑らかな面状に構成されつつ、当該解析空間Ｓの周縁部に９０度格子Ｍ１のみが配列される。
これにより、計算実行時における周縁境界Ｂ２での発散や反射の発生を抑え、有限空間である解析空間Ｓを無限空間として取り扱う吸収境界の計算を成立させることができる。
これにより、計算実行時における周縁境界Ｂ２での発散や反射の発生を抑え、有限空間である解析空間Ｓを無限空間として取り扱う吸収境界の計算を成立させることができる。
なお、本発明を適用可能な実施形態は、上述した実施形態に限定されることなく、本発明の趣旨を逸脱しない範囲で適宜変更可能である。
例えば、上記実施形態では、９０度格子Ｍ１における磁界分布から４５度格子Ｍ２における磁界が補間計算され、４５度格子Ｍ２における電界分布から９０度格子Ｍ１における電界が補間計算されることとしたが、これらの補間計算における電界と磁界を入れ替えてもよい。つまり、９０度格子Ｍ１における電界分布から４５度格子Ｍ２における電界が補間計算され、４５度格子Ｍ２における磁界分布から９０度格子Ｍ１における磁界が補間計算されることとしてもよい。
また、解析空間Ｓの周縁部には、９０度格子Ｍ１及び４５度格子Ｍ２のいずれか一方のみが配列されればよく、周縁境界Ｂ２が滑らかな面状に構成可能であれば、当該周縁部には、９０度格子Ｍ１でなく、４５度格子Ｍ２のみが配列されることとしてもよい。
また、格子間での電界及び磁界の補間計算では、各格子内で電界及び磁界が直線的に分布するものと仮定して線形補間することとしたが、各格子内での分布状況に応じて線形補間以外の補間方法を用いてもよい。
また、本発明に係る異方性導電物質は、所定の第一の方向に沿った導電率が他の方向のものとは異なる一の層と、第一の方向とは異なる第二の方向に沿った導電率が他の方向のものとは異なる他の層とが積層されたものであれば、特に限定されず、擬似等方積層材でなくともよいし、繊維強化プラスチックでなくともよい。
また、本実施形態では、解析モデルとして、図１に示すような繊維方向が互いに４５°異なる４層の繊維層を２セット（すなわち合計８層）積層させた略平板状の擬似等方積層材を用いたが、本発明の解析モデルとしては、積層数が何層であってもよく、また繊維方向は０°、９０°、４５°、−４５°がどのような順番で積層されていてもよい。
さらに当該解析モデルは、互いに異なる繊維方向の複数の繊維層が積層されたものであればよく、複数の繊維層の繊維方向が４５°以外の角度で異なっていてもよい。但し、この場合に、繊維方向に沿った辺を有する計算格子を用いることや、繊維方向間の角度に対応した計算式を用いて格子間での電磁界成分の補間計算を行うことなどは勿論である。
さらに当該解析モデルは、互いに異なる繊維方向の複数の繊維層が積層されたものであればよく、複数の繊維層の繊維方向が４５°以外の角度で異なっていてもよい。但し、この場合に、繊維方向に沿った辺を有する計算格子を用いることや、繊維方向間の角度に対応した計算式を用いて格子間での電磁界成分の補間計算を行うことなどは勿論である。
ＣＭ 複合材料
Ｌ１〜Ｌ８ 第１層〜第８層
Ｍ１ ９０度格子
Ｍ２ ４５度格子
Ｂ１ 格子境界
Ｂ２ 周縁境界
Ｓ 解析空間
Ｌ１〜Ｌ８ 第１層〜第８層
Ｍ１ ９０度格子
Ｍ２ ４５度格子
Ｂ１ 格子境界
Ｂ２ 周縁境界
Ｓ 解析空間
Claims (4)
- 所定の第一の方向に沿った導電率が他の方向のものとは異なる一の層と、前記第一の方向とは異なる第二の方向に沿った導電率が他の方向のものとは異なる他の層とが積層されてなる異方性導電物質の電磁界解析方法であって、
コンピュータが、
前記第一の方向に沿った辺を有する第一の計算格子により前記一の層を格子分割するとともに、前記第二の方向に沿った辺を有する第二の計算格子により前記他の層を格子分割し、
前記第一の格子における電界及び磁界のうちの一方の電磁界成分の分布から前記第二の格子における当該一方の電磁界成分を補間計算し、
前記第二の格子における電界及び磁界のうちの他方の電磁界成分の分布から前記第一の格子における当該他方の電磁界成分を補間計算することを特徴とする異方性導電物質の電磁界解析方法。 - コンピュータが、
解析空間の周縁境界を滑らかな面状に構成しつつ、当該解析空間の周縁部に前記第一の格子及び前記第二の格子のいずれか一方のみを配列することを特徴とする請求項１に記載の異方性導電物質の電磁界解析方法。 - 前記一の層及び前記他の層は、いずれもｘｙ平面内に配置されるとともに、ｚ軸方向に互いに積層され、
前記第一の格子は、ｘｙｚの直交３軸に沿った立方体状の計算格子であり、
前記第二の格子は、前記第一の格子に内接する大きさのものであって、ｘｙ平面内で±４５°方向に沿った２辺とｚ軸に沿った１辺とからなる直方体状の計算格子であることを特徴とする請求項１または２に記載の異方性導電物質の電磁界解析方法。 - 前記異方性導電物質は、前記第一の方向及び前記第二の方向に沿った導電性の繊維で樹脂が強化された繊維強化プラスチックであることを特徴とする請求項１〜３のいずれか一項に記載の異方性導電物質の電磁界解析方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2015062483A JP6034909B2 (ja) | 2015-03-25 | 2015-03-25 | 異方性導電物質の電磁界解析方法 |
| US15/054,027 US9805149B2 (en) | 2015-03-25 | 2016-02-25 | Electromagnetic field analysis method for anisotropic conductive material |
| EP16158513.8A EP3073278B1 (en) | 2015-03-25 | 2016-03-03 | Electromagnetic field analysis method for anisotropic conductive material |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2015062483A JP6034909B2 (ja) | 2015-03-25 | 2015-03-25 | 異方性導電物質の電磁界解析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2016181225A JP2016181225A (ja) | 2016-10-13 |
| JP6034909B2 true JP6034909B2 (ja) | 2016-11-30 |
Family
ID=55484863
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2015062483A Active JP6034909B2 (ja) | 2015-03-25 | 2015-03-25 | 異方性導電物質の電磁界解析方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US9805149B2 (ja) |
| EP (1) | EP3073278B1 (ja) |
| JP (1) | JP6034909B2 (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107748834B (zh) * | 2017-11-22 | 2018-08-14 | 中南大学 | 一种计算起伏观测面磁场的快速、高精度数值模拟方法 |
| CN111079326B (zh) * | 2019-11-20 | 2023-04-28 | 杭州电子科技大学 | 二维各向异性网格单元度量张量场光滑化方法 |
| CN111090958B (zh) * | 2019-12-30 | 2021-11-30 | 西北工业大学 | 一种基于亚网格技术的电磁波时域高效数值混合方法 |
| CN111381114B (zh) * | 2020-04-01 | 2020-12-18 | 华中科技大学 | 一种用混合场域等效无穷远边界传导介质的方法及*** |
| CN112507647B (zh) * | 2020-12-15 | 2023-07-21 | 重庆邮电大学 | 空间电磁场作用分叉线的电磁耦合时域建模分析方法 |
| CN112836415B (zh) * | 2021-02-26 | 2023-12-05 | 英特工程仿真技术（大连）有限公司 | 一种电磁场非匹配棱边元的插值方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AUPN499795A0 (en) * | 1995-08-25 | 1995-09-14 | University Of Queensland, The | Computer simulation of magnetic resonance spectra |
| JP3415530B2 (ja) | 1999-12-24 | 2003-06-09 | 日本電信電話株式会社 | 電磁界シミュレーション方法 |
| JP3987377B2 (ja) * | 2002-05-21 | 2007-10-10 | オリンパス株式会社 | Ｆｄｔｄ電磁場シミュレーションをコンピュータに実行させるためのコンピュータプログラムが記録された記録媒体及びｆｄｔｄ電磁場シミュレーション方法 |
| JP4319843B2 (ja) * | 2003-02-05 | 2009-08-26 | 株式会社エヌ・ティ・ティ・ドコモ | 電磁界特性解析方法 |
| JP2008275403A (ja) * | 2007-04-27 | 2008-11-13 | Hitachi Ltd | 磁気異方性磁界解析法 |
| GB0723224D0 (en) | 2007-11-27 | 2008-01-09 | Fujitsu Ltd | A multi-level gris embedding process with uniform and non-uniform grid refinement for multigurid-fdtd electromagnetic solver |
| JP2010262338A (ja) | 2009-04-30 | 2010-11-18 | Fujitsu Ltd | 電磁流分布処理装置、電磁流分布処理方法、及び電磁流分布処理プログラム |
-
2015
- 2015-03-25 JP JP2015062483A patent/JP6034909B2/ja active Active
-
2016
- 2016-02-25 US US15/054,027 patent/US9805149B2/en active Active
- 2016-03-03 EP EP16158513.8A patent/EP3073278B1/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2016181225A (ja) | 2016-10-13 |
| US9805149B2 (en) | 2017-10-31 |
| EP3073278A1 (en) | 2016-09-28 |
| EP3073278B1 (en) | 2023-11-22 |
| US20160283620A1 (en) | 2016-09-29 |
| EP3073278C0 (en) | 2023-11-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6034909B2 (ja) | 異方性導電物質の電磁界解析方法 | |
| Matt et al. | Macro-fiber composite piezoelectric rosettes for acoustic source location in complex structures | |
| Man et al. | Semi-analytical analysis for piezoelectric plate using the scaled boundary finite-element method | |
| US20030204343A1 (en) | Electromagnetic field analysis method based on FDTD method, medium representation method in electromagnetic field analysis, simulation device, and storage medium | |
| Ge et al. | Accurate modeling of PZT-induced Lamb wave propagation in structures by using a novel spectral finite element method | |
| Huang et al. | GO/PO and PTD with virtual divergence factor for fast analysis of scattering from concave complex targets | |
| D’Ottavio et al. | Assessment of plate theories for free-edge effects | |
| Zhang et al. | Performance of a massively parallel higher-order method of moments code using thousands of CPUs and its applications | |
| Carrera et al. | Free vibration analysis of thin-walled cylinders reinforced with longitudinal and transversal stiffeners | |
| Shi et al. | Spectral element method for band-structure calculations of 3D phononic crystals | |
| Li et al. | 3D FEM simulation of SAW resonators using hierarchical cascading technique and general purpose graphic processing unit | |
| Chen et al. | A semianalytical spectral element method for the analysis of 3-D layered structures | |
| Senjanović et al. | Analytical solution for free vibrations of a moderately thick rectangular plate | |
| Bathe | The finite element method with “overlapping finite elements” | |
| Shang et al. | 8-node unsymmetric distortion-immune element based on Airy stress solutions for plane orthotropic problems | |
| Cao et al. | Design for additive manufacturing: Optimization of piping network in compact system with enhanced path-finding approach | |
| JP6317776B2 (ja) | 航空機機体構造材の耐雷試験方法 | |
| Codecasa et al. | Discrete geometric formulation of admittance boundary conditions for frequency domain problems over tetrahedral dual grids | |
| CN109214021B (zh) | 一种电磁仿真中斜口面激励源的生成方法 | |
| Topal et al. | Frequency optimization of laminated composite skew sandwich plates | |
| JP2018028768A (ja) | 異方性導電物質の電磁界解析方法 | |
| Jafari et al. | Vibration optimization of skew composite plates using the Rayleigh-Ritz & response surface methods | |
| Li et al. | Finite element method for analysis of band structures of phononic crystal slabs with archimedean-like tilings | |
| Ansarizadeh | Calculation of the shielding effectiveness of carbon-fiber composite structures | |
| Cai | Evaluation of layered multiple-scattering method for antiplane shear wave scattering from gratings |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160706 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20161004 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20161028 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6034909 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |