JP6317776B2 - 航空機機体構造材の耐雷試験方法 - Google Patents
航空機機体構造材の耐雷試験方法 Download PDFInfo
- Publication number
- JP6317776B2 JP6317776B2 JP2016069876A JP2016069876A JP6317776B2 JP 6317776 B2 JP6317776 B2 JP 6317776B2 JP 2016069876 A JP2016069876 A JP 2016069876A JP 2016069876 A JP2016069876 A JP 2016069876A JP 6317776 B2 JP6317776 B2 JP 6317776B2
- Authority
- JP
- Japan
- Prior art keywords
- electric field
- analysis
- field component
- along
- conductivity
- 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
- 239000000463 material Substances 0.000 title claims description 18
- 238000010998 test method Methods 0.000 title claims description 12
- 238000004458 analytical method Methods 0.000 claims description 74
- 230000005684 electric field Effects 0.000 claims description 59
- 238000000034 method Methods 0.000 claims description 26
- 230000005672 electromagnetic field Effects 0.000 claims description 23
- 239000000835 fiber Substances 0.000 claims description 19
- 229920002430 Fibre-reinforced plastic Polymers 0.000 claims description 15
- 239000011151 fibre-reinforced plastic Substances 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 6
- 239000004020 conductor Substances 0.000 description 6
- 239000004918 carbon fiber reinforced polymer Substances 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002250 progressing effect Effects 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
Classifications
-
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/60—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating electrostatic variables, e.g. electrographic flaw testing
-
- 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
- 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
- G01R29/0842—Measurements related to lightning, e.g. measuring electric disturbances, warning systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/005—Testing of electric installations on transport means
- G01R31/008—Testing of electric installations on transport means on air- or spacecraft, railway rolling stock or sea-going vessels
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/28—Fuselage, exterior or interior
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Shielding Devices Or Components To Electric Or Magnetic Fields (AREA)
Description
この解析手法では、直交3軸の座標系で定義された導電率をz軸回りに回転させたときに現れる導電率テンソルの非対角要素によって、異方性導電率を表現している。そして、導電性テンソルをマクスウェルの式に適用して差分法を用いることで、各軸についての離散化式を得る。
Yee格子の定義では離散化された各軸の電界が格子上の異なる位置に存在するが、得られた離散化式では、異なる位置の電界(Yee格子では定義されていない位置の電界)が混在しており、この式のままでは計算することができない(後述の実施形態で詳細に説明する)。そこで、従来の解析手法では、これらの異なる位置の電界を等しいものと仮定して解法を試みている。
繊維方向に沿った第一の導電率と、他の方向に沿った前記第一の導電率とは異なる第二の導電率とを有する繊維強化プラスチックを備える航空機機体構造材の前記繊維方向の両端に電流を印加する航空機耐雷試験と、
互いに直交する第一の辺と第二の辺を有する複数の解析格子を、前記繊維方向及び前記他の方向のうちの少なくとも一方が、前記複数の解析格子の前記第一の辺及び前記第二の辺のいずれとも異なる方向に沿うように適用した前記繊維強化プラスチックの解析モデルを用い、コンピュータが、前記第一の辺上に位置するとともに前記第二の辺に沿った第一の電界成分を、当該電界成分の周囲の複数の第二の辺上に位置するとともに第二の辺に沿った複数の第二の電界成分に基づいて算出する電磁界解析と、
を備え、
前記コンピュータは、前記複数の第二の電界成分から前記第一の電界成分を算出するときに、反復法を適用して、前記第一の電界成分の変化率が所定の閾値以下となるまで繰り返し計算を行うことを特徴とする。
前記コンピュータは、前記第一の電界成分を前記複数の第二の電界成分の平均値として算出することを特徴とする。
前記コンピュータは、前記第一の電界成分が前記繊維強化プラスチックの周縁に沿った第一の辺上に位置する場合、前記複数の第二の電界成分のうち、前記繊維強化プラスチックに含まれるものに基づいて、当該第一の電界成分を算出することを特徴とする。
前記解析格子は、直交3軸に沿った立方体状の解析格子であり、
前記複数の第二の電界成分は、前記第一の電界成分が位置する第一の辺の両端の2つの節点で当該第一の辺と連なる4つの第二の辺上に個別に位置する4つの電界成分であることを特徴とする。
前記電磁界解析は、FDTD法を用いたものであり、
前記解析格子は、直交3軸に沿った立方体状の解析格子であり、
前記コンピュータは、
直交3軸の座標系で定義された導電率を前記繊維方向及び前記他の方向のいずれとも直交する軸回りに回転させたときの導電率テンソルをアンペールの式に適用した後に、差分法を用いて離散化させ、
当該離散化させた離散化式に反復法を適用して繰り返し計算を行うときに、当該離散化式に含まれる未定義の電界成分を前記第一の電界成分として前記複数の第二の電界成分から算出することを特徴とする。
つまり、第二の辺上に定義されたものであって第一の辺上には定義されていない第一の電界成分が、その周辺の複数の第二の辺上に位置する複数の第二の電界成分に基づいて補間計算される。
したがって、異なる位置の電界を単純に等しいと定義していた従来に比べ、異方性導電率を有する繊維強化プラスチックの電磁界解析を高精度に行うことができ、ひいては航空機機体構造材の耐雷試験を効率的かつ高精度に行うことができる。
以下、この解析手法について詳細に説明する。
そこで、電界の位置を「|i,j,k 」と表記し、差分法を用いて式(6)を離散化すると、以下の式(7)〜(9)のようになる。
この問題に対し、従来の解析手法では、これらの異なる位置の電界を、以下の式(10)のように等しいものと仮定したうえで、式(7)と式(8)を連立させている。
具体的には、以下の式(11),(12)を用いて、未定義Ex及び未定義Eyを算出する。
例えば、図1(a)において、図中の右側の解析格子が解析対象に含まれ、左側の解析格子がその周囲空間のものであった場合、未定義Exの算出にはi+1/2の2つの定義Exのみを使用する。
未定義Eyについても同様である。
しかし、本実施形態では、より効率的な方法として、式(7),(8)を連立させずに、ガウスザイデル法に代表される反復法を適用し、Ex、Eyの変化率が十分に小さくなる(所定の閾値以下となる)まで繰り返し計算する。そして、この繰り返し計算の中で、式(11),(12)を用いて、未定義Ex,Eyを複数の定義Ex,Eyで適宜補間する。
これにより、異方性導電率物質の電磁的特性を解析することができる。
続いて、本実施形態の電磁界解析方法を用いた解析例について説明する。
本解析例では、本実施形態の電磁界解析方法によるものを、上述した従来の解析手法によるものと比較し、それぞれの妥当性を検証した。
比較を行うための解析モデルは、図2(a)に示すような一方向材平板とした。この平板の大きさは300mm×300mm×10mmとし、メッシュ数は21×21×1とした。
格子上の一本の直線で定義された金属棒を解析対象と吸収境界に接続することで当該金属棒を電極と定義し、印加電極とリターン電極を、繊維方向と平行な平板の対角線上の両端に接続した。
印加電流は、SAE ARP-5416で規定される航空機耐雷試験用Component A波形で、ピーク値を3kAとした。
図2(b),(c)は、本解析例の結果を示す図であって、印加電流がピーク値となる印加後6.4μsにおける電流密度を示す図であり、(b)が従来の解析手法によるもの、(c)が本実施形態の電磁界解析方法によるものである。
なお、解析結果の表示には、解析格子が形成するセル毎に電流密度分布をカラーグラデーションにより塗り潰した。また、各セルでの電流の流れる方向を表示するため、セル内部には式(12),(13)で求めたJx、Jyが形成するベクトル方向を矢印で示した。
これに対し、本実施形態の電磁界解析方法では、図2(c)に示すように、両電極間を結ぶように平板の対角線上に強い電流分布が計算されており、繊維方向への異方性導電率に応じた電流分布を再現できていると言える。
以上のように、本実施形態の電磁界解析方法によれば、解析格子のy軸に沿った辺上に位置する未定義Exが、当該未定義Exの周囲のx軸に沿った複数の辺上に位置する複数の定義Exに基づいて算出され、未定義Eyについても同様に、複数の定義Eyに基づいて算出される。
つまり、未定義Ex,Eyがその周辺の複数の定義Ex,Eyに基づいて補間計算される。
したがって、異なる位置の電界を単純に等しいと定義していた従来に比べ、異方性導電率を有する異方性導電物質の電磁界解析を高精度に行うことができる。
なお、本発明を適用可能な実施形態は、上述した実施形態に限定されることなく、本発明の趣旨を逸脱しない範囲で適宜変更可能である。
Claims (5)
- 繊維方向に沿った第一の導電率と、他の方向に沿った前記第一の導電率とは異なる第二の導電率とを有する繊維強化プラスチックを備える航空機機体構造材の前記繊維方向の両端に電流を印加する航空機耐雷試験と、
互いに直交する第一の辺と第二の辺を有する複数の解析格子を、前記繊維方向及び前記他の方向のうちの少なくとも一方が、前記複数の解析格子の前記第一の辺及び前記第二の辺のいずれとも異なる方向に沿うように適用した前記繊維強化プラスチックの解析モデルを用い、コンピュータが、前記第一の辺上に位置するとともに前記第二の辺に沿った第一の電界成分を、当該電界成分の周囲の複数の第二の辺上に位置するとともに第二の辺に沿った複数の第二の電界成分に基づいて算出する電磁界解析と、
を備え、
前記コンピュータは、前記複数の第二の電界成分から前記第一の電界成分を算出するときに、反復法を適用して、前記第一の電界成分の変化率が所定の閾値以下となるまで繰り返し計算を行うことを特徴とする航空機機体構造材の耐雷試験方法。 - 前記コンピュータは、前記第一の電界成分を前記複数の第二の電界成分の平均値として算出することを特徴とする請求項1に記載の航空機機体構造材の耐雷試験方法。
- 前記コンピュータは、前記第一の電界成分が前記繊維強化プラスチックの周縁に沿った第一の辺上に位置する場合、前記複数の第二の電界成分のうち、前記繊維強化プラスチックに含まれるものに基づいて、当該第一の電界成分を算出することを特徴とする請求項1または2に記載の航空機機体構造材の耐雷試験方法。
- 前記解析格子は、直交3軸に沿った立方体状の解析格子であり、
前記複数の第二の電界成分は、前記第一の電界成分が位置する第一の辺の両端の2つの節点で当該第一の辺と連なる4つの第二の辺上に個別に位置する4つの電界成分であることを特徴とする請求項1〜3のいずれか一項に記載の航空機機体構造材の耐雷試験方法。 - 前記電磁界解析は、FDTD法を用いたものであり、
前記解析格子は、直交3軸に沿った立方体状の解析格子であり、
前記コンピュータは、
直交3軸の座標系で定義された導電率を前記繊維方向及び前記他の方向のいずれとも直交する軸回りに回転させたときの導電率テンソルをアンペールの式に適用した後に、差分法を用いて離散化させ、
当該離散化させた離散化式に反復法を適用して繰り返し計算を行うときに、当該離散化式に含まれる未定義の電界成分を前記第一の電界成分として前記複数の第二の電界成分から算出することを特徴とする請求項1〜4のいずれか一項に記載の航空機機体構造材の耐雷試験方法。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016069876A JP6317776B2 (ja) | 2016-03-31 | 2016-03-31 | 航空機機体構造材の耐雷試験方法 |
US15/417,085 US10119936B2 (en) | 2016-03-31 | 2017-01-26 | Electromagnetic field analysis method for anisotropic conductive material |
EP17157394.2A EP3226148B1 (en) | 2016-03-31 | 2017-02-22 | Electromagnetic field analysis method for anisotropic conductive material |
ES17157394T ES2710619T3 (es) | 2016-03-31 | 2017-02-22 | Métodos de análisis de campos electromagnéticos para materiales conductores anisotrópicos |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016069876A JP6317776B2 (ja) | 2016-03-31 | 2016-03-31 | 航空機機体構造材の耐雷試験方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2017182516A JP2017182516A (ja) | 2017-10-05 |
JP6317776B2 true JP6317776B2 (ja) | 2018-04-25 |
Family
ID=58191272
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016069876A Active JP6317776B2 (ja) | 2016-03-31 | 2016-03-31 | 航空機機体構造材の耐雷試験方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US10119936B2 (ja) |
EP (1) | EP3226148B1 (ja) |
JP (1) | JP6317776B2 (ja) |
ES (1) | ES2710619T3 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110119586B (zh) * | 2019-05-21 | 2023-02-28 | 中煤科工集团西安研究院有限公司 | 轴向电导率各向异性瞬变电磁三分量三维fdtd正演方法 |
JP2021110997A (ja) * | 2020-01-07 | 2021-08-02 | 住友重機械工業株式会社 | シミュレーション方法、シミュレーション装置、及びプログラム |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3415530B2 (ja) | 1999-12-24 | 2003-06-09 | 日本電信電話株式会社 | 電磁界シミュレーション方法 |
JP4040269B2 (ja) * | 2001-06-22 | 2008-01-30 | 日本電信電話株式会社 | 二次元電磁界シミュレーション方法 |
JP3630120B2 (ja) * | 2001-07-16 | 2005-03-16 | Tdk株式会社 | 磁気変換デバイスの特性解析方法及びプログラム |
JP3824561B2 (ja) * | 2002-07-05 | 2006-09-20 | 独立行政法人科学技術振興機構 | 電磁環境解析方法、電磁環境解析プログラム及びそのプログラムを記録した記録媒体 |
US7248041B2 (en) * | 2003-07-28 | 2007-07-24 | Cummins, Inc. | Device and method for measuring transient magnetic performance |
JP4984586B2 (ja) * | 2006-03-17 | 2012-07-25 | 富士通株式会社 | 電磁界シミュレータ、及び電磁界シミュレートプログラム |
JP4262762B2 (ja) | 2006-12-11 | 2009-05-13 | シャープ株式会社 | 電磁界解析プログラム |
JP2010262338A (ja) * | 2009-04-30 | 2010-11-18 | Fujitsu Ltd | 電磁流分布処理装置、電磁流分布処理方法、及び電磁流分布処理プログラム |
JP5149321B2 (ja) | 2010-03-24 | 2013-02-20 | 株式会社東芝 | 電磁場シミュレーション方法、電磁場シミュレーション装置、半導体装置の製造方法 |
KR102160483B1 (ko) | 2012-05-23 | 2020-09-28 | 앤시스 루머리컬 아이피, 엘엘씨 | 비등방성 매체를 시뮬레이션 하기 위해서 좌표계를 변환하기 위한 장치 및 방법 |
-
2016
- 2016-03-31 JP JP2016069876A patent/JP6317776B2/ja active Active
-
2017
- 2017-01-26 US US15/417,085 patent/US10119936B2/en active Active
- 2017-02-22 EP EP17157394.2A patent/EP3226148B1/en active Active
- 2017-02-22 ES ES17157394T patent/ES2710619T3/es active Active
Also Published As
Publication number | Publication date |
---|---|
EP3226148A1 (en) | 2017-10-04 |
JP2017182516A (ja) | 2017-10-05 |
EP3226148B1 (en) | 2018-11-21 |
ES2710619T3 (es) | 2019-04-26 |
US20170284964A1 (en) | 2017-10-05 |
US10119936B2 (en) | 2018-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huynh et al. | The extended finite element method for fracture in composite materials | |
US9805149B2 (en) | Electromagnetic field analysis method for anisotropic conductive material | |
CN103714216B (zh) | 基于三维微观晶须构成的复合材料力学性能评估方法 | |
JP6317776B2 (ja) | 航空機機体構造材の耐雷試験方法 | |
Fletcher et al. | Fast modelling of field evaporation in atom probe tomography using level set methods | |
US9753072B2 (en) | Magnetic coupling for electrical conductivity assessment | |
CN115618577A (zh) | 一种低电导率复合材料飞行器的雷电效应仿真评估方法 | |
Gaynutdinov et al. | Study of impact lightning at of electromagnetic and thermal stability of fuselage of unmanned aerial vehicle | |
Deluzet et al. | A two field iterated asymptotic-preserving method for highly anisotropic elliptic equations | |
Aubertin et al. | Periodic and anti-periodic boundary conditions with the Lagrange multipliers in the FEM | |
Jalali et al. | Accuracy assessment of finite volume discretizations of convective fluxes on unstructured meshes | |
Albanese et al. | Accuracy assessment of numerical tracing of three-dimensional magnetic field lines in tokamaks with analytical invariants | |
Wang | On the current source implementation for the ADI-FDTD method | |
Weber et al. | Validation of numerical simulation approach for lightning transient analysis of a transport category aircraft | |
Grayver et al. | Comparison of iterative and direct solvers for 3D CSEM modeling | |
Zhang et al. | Three-dimensional magnetostatic analysis using structured mesh and nodal elements | |
JP4597559B2 (ja) | 電磁界解析装置、電磁界解析方法および電磁界解析プログラム | |
Elakkad et al. | A mixed finite element method for elasticity problem | |
Poggie | DC glow discharges: a computational study for flow control applications | |
SATAKE et al. | Development of electromagnetic simulation supporting lightning protection design of mitsubishi regional jet | |
JP2004004054A (ja) | Fdtd法を用いた電磁界解析方法、電磁界解析における媒質表現方法、シミュレーション装置、及びプログラム | |
Borzi et al. | Solution of large complex BEM systems derived from high-resolution human models | |
Wang et al. | The study about cloud of points reconstruction with the framework of meshfree method for viscous flows | |
Mifune | Interface homogenization technique for electromagnetic finite element analysis including anisotropic medea | |
Albanese et al. | Numerical formulations for accurate magnetic field flow tracing in fusion tokamaks |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20171107 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20171227 |
|
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: 20180306 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20180330 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6317776 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
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 |