船体の永久磁気推定方法
【課題】脱磁処理を実施する前に、船体に磁場を印加したときの前記船体の永久磁気を推定することが可能な、船体の永久磁気推定方法を提供する。
【解決手段】脱磁対象の船体と同じ鋼材からなる試験片を作成し(ＳＴ１)、試験片への印加磁場を変化させて磁化特性を測定し(ＳＴ２)、測定した複数パターンの磁化特性から、鋼材のプライザッハ分布図を作成する(ＳＴ３)。公知の数値計算法に対応した船体数値計算モデルを作成し(ＳＴ４)、各要素の印加磁場を計算し(ＳＴ５)、計算した印加磁場に対応するプライザッハ図表を各要素について作成し(ＳＴ６)、プライザッハ分布図とプライザッハ図表とに基づいて各要素の永久磁気を算出し(ＳＴ７)、各要素の永久磁気推定値に基づいて数値計算法により船体周辺の永久磁場を推定し(ＳＴ８)、磁気処理後の船体周辺の永久磁場推定値を出力する(ＳＴ９)。
【解決手段】脱磁対象の船体と同じ鋼材からなる試験片を作成し(ＳＴ１)、試験片への印加磁場を変化させて磁化特性を測定し(ＳＴ２)、測定した複数パターンの磁化特性から、鋼材のプライザッハ分布図を作成する(ＳＴ３)。公知の数値計算法に対応した船体数値計算モデルを作成し(ＳＴ４)、各要素の印加磁場を計算し(ＳＴ５)、計算した印加磁場に対応するプライザッハ図表を各要素について作成し(ＳＴ６)、プライザッハ分布図とプライザッハ図表とに基づいて各要素の永久磁気を算出し(ＳＴ７)、各要素の永久磁気推定値に基づいて数値計算法により船体周辺の永久磁場を推定し(ＳＴ８)、磁気処理後の船体周辺の永久磁場推定値を出力する(ＳＴ９)。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、船体に磁場を印加したときの船体の永久磁気を推定する、船体の永久磁気推定方法に関する。
【背景技術】
【０００２】
船舶は多種多様な鋼材で構成されており、木製やＦＲＰ製等の小型船舶を除けばほとんど磁性体で占められている。船体磁気は、静磁界と動磁界で構成される。ただし、動磁界は静磁界に比べて非常に小さく、問題になることは少ない。静磁界は、誘導磁気と永久磁気に区別できる。誘導磁気は、船体が地磁気により磁化されることにより発生する。永久磁気は、船舶の建造時の溶接、裁断等により着磁するとともに、船体が船台に長期的に同じ方向に置かれると地磁気によって着磁する。さらに、就役後も波浪や地磁気を受けることによって永久磁気は変化する。したがって、就役直後は船体は特に大きな永久磁気を帯びており、時間の経過によってそれはさらに大きくなる。これに対処するためには、定期的に船体の永久磁気を脱磁する必要がある。
【０００３】
従来、鋼鉄によって構成された船体が保持する永久磁気を消去するために、船体に幾巻きかの胴巻輪線を装着して、それらに電流の極性を正負に変えながら、その大きさを漸減して通電する(つまり磁場を印加する)脱磁処理を実施している。そして、脱磁処理後に複数個の磁気検出器で船体周辺の永久磁場を測定し、脱磁処理の効果を評価している。
【０００４】
なお、積分方程式法又は有限要素法といった数値解析手法では、船舶計算モデルを作成して、船舶に地磁気が印加されることによる誘導磁気を計算することが可能である。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、以上の技術では、脱磁処理を実施した後でなくては複数個の磁気検出器による船体の永久磁気を把握することができないため、船体の形状及び材料に対応した最適な脱磁処理を実施することはできない。
【０００６】
本発明はこうした状況を認識してなされたものであり、その目的は、脱磁処理を実施する前に、船体に磁場を印加したときの前記船体の永久磁気を推定することが可能な、船体の永久磁気推定方法を提供することにある。
【課題を解決するための手段】
【０００７】
本発明のある態様は、船体の永久磁気推定方法である。この方法は、
船体に磁場を印加したときの前記船体の永久磁気を推定する、船体の永久磁気推定方法であって、
船体の磁性材を下記式で定義される磁化状態関数ｍ(Ｈa,Ｈb)
【数１】
によって磁化状態が表される磁気双極子であるヒステロンの集合体と仮定した場合の前記ヒステロンの分布数を前記Ｈa及び前記Ｈbの関数として表したプライザッハ分布関数Ｐ(Ｈa,Ｈb)を、前記Ｈa及び前記Ｈbの所定範囲で取得する取得ステップと、
前記船体の数値計算モデルの各要素について、所定の印加磁場に対応した前記磁化状態関数ｍ(Ｈa,Ｈb)を、前記Ｈa及び前記Ｈbの所定範囲で導出する導出ステップと、
取得した前記プライザッハ分布関数Ｐ(Ｈa,Ｈb)と、導出した前記磁化状態関数ｍ(Ｈa,Ｈb)とに基づいて、各要素の永久磁気を推定する推定ステップとを有する。
【０００８】
ある態様の方法において、前記プライザッハ分布関数Ｐ(Ｈa,Ｈb)は、前記船体と同じ磁性材からなる試験片に磁場を印加して取得した前記試験片の複数の磁化特性に基づいて算出したものであるとよい。
【０００９】
ある態様の方法において、前記数値計算モデルは、積分方程式法用モデルであってもよい。
【００１０】
あるいは、前記数値計算モデルは、有限要素法用モデルであってもよい。
【００１１】
ある態様の方法において、前記導出ステップは、印加磁場の極性の正負を変えながら前記印加磁場の大きさを漸減させた場合の各々の印加磁場に対応した前記磁化状態関数ｍ(Ｈa,Ｈb)を導出するものであるとよい。
【００１２】
ある態様の方法において、前記推定ステップは、下記式で表される永久磁気推定値Ｍ
【数２】
を計算するものであるとよい。
【００１３】
なお、以上の構成要素の任意の組合せ、本発明の表現を装置やシステム、プログラムなどの間で変換したものもまた、本発明の態様として有効である。
【発明の効果】
【００１４】
本発明によれば、脱磁処理を実施する前に、船体に磁場を印加したときの前記船体の永久磁気を推定することが可能となる。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施の形態に係る船体の永久磁気推定方法のフローチャート。
【図２】船体の脱磁処理の模式的説明図。
【図３】磁気双極子であるヒステロンが無数に集まって構成するメジャーループの説明図。
【図４】磁性体のプライザッハ図表。
【図５】一例としてＨs＝５としたときのプライザッハ図表とヒステリシス曲線の関係を示す説明図。
【図６】磁性体の例示的なプライザッハ分布図。
【図７】立体的に示したプライザッハ分布図。
【図８】磁性体の複数パターンの磁化特性図。
【図９】(Ａ)は、印加磁場をＨs→−Ｈs→１０と変化させたときの磁性体のプライザッハ図表及び磁化特性図。(Ｂ)は、印加磁場をＨs→２→１０と変化させたときの磁性体のプライザッハ図表及び磁化特性図。
【図１０】実施の形態で作成する試験片の例示的な形状説明図。
【図１１】船体の例示的な有限要素法用モデル。
【発明を実施するための形態】
【００１６】
以下、図面を参照しながら本発明の好適な実施の形態を詳述する。なお、各図面に示される同一または同等の構成要素、部材、処理等には同一の符号を付し、適宜重複した説明は省略する。また、実施の形態は発明を限定するものではなく例示であり、実施の形態に記述されるすべての特徴やその組み合わせは必ずしも発明の本質的なものであるとは限らない。
【００１７】
本実施の形態は、船体に幾巻きかの胴巻輪線を装着して通電する(磁場を印加する)脱磁処理に先だって、前記船体に磁場を印加したときの前記船体の永久磁気を推定するものである。説明の順序は下記のとおりとする。
１．実際の船体の脱磁処理
２．プライザッハ図表(永久磁気の推定に用いる)
３．プライザッハ分布図(永久磁気の推定に用いる)
４．船体の永久磁気推定方法の流れ
【００１８】
(１．実際の船体の脱磁処理)
図２は、船体の脱磁処理の模式的説明図である。本図に示されるように、脱磁処理の際には、船体１を係留用ワイヤ２とドルフィン３で固定した後に、船体１の外周に幾巻きかの胴巻輪線４を装着する。そして、磁気処理用電源５から磁気処理電流が通電され、船体の永久磁気を変化させる。船体１の下方の海底に設置される多数の磁気センサ６により、磁気処理電流通電後における船体周辺の磁場を計測し、コンピュータを有する測定室７に船体周辺の磁気波形が表示され、次の磁気処理電流が通電される。磁気処理電流値は、プラスマイナス交互に漸減しながら通電し、最終的に残留した永久磁場は磁気センサ６によって把握する。従来、磁気処理電流値は過去の実績を基に決めていた。本実施の形態では、磁気処理電流値を適切に定めるために、船体の脱磁処理に先だって、前記船体に磁場を印加したときの前記船体の永久磁気を推定するものである。
【００１９】
(２．プライザッハ図表)
プライザッハ図表は、本実施の形態で船体の永久磁気を推定するのに用いる。一般に、磁性体の磁化特性を表すＢ−Ｈ曲線(磁束密度Ｂと印加磁界Ｈの関係)は、飽和磁場±Ｈsを印加した時のメジャーループと呼ばれる曲線で示されている。しかし、船舶の磁気処理において実際に印加する磁界は、船体に装着した胴巻輪線(図２参照)への初期の巨大な通電では飽和磁界となるものの、通電量は次第に漸減し、その磁気特性はマイナーループをたどることとなる。そのため、磁気処理中及び磁気処理後の船体の永久磁気や周辺磁場を推定計算するには、船体を成す鋼材のマイナーループを把握する必要がある。メジャーループやある指定した磁場を印加した時のマイナーループについては、あらかじめ、船体を成す鋼材と同じ材質の磁性体からなる試験片を作成し、ヒステリシス特性を計測することが可能である。しかし、任意に磁場をかけたときのマイナーループを全て事前に計測しておくことはできない。
【００２０】
そこで、任意のマイナーループを推定するために、図３に示すように、ヒステロンと呼ばれる磁気双極子が無数に集まりメジャーループを構成していると仮定する。つまり、磁性体をヒステロンの集合体と仮定する。ヒステロンは、印加磁界がＨaとＨb(−Ｈs＜Ｈb＜Ｈa＜Ｈs)を境界に磁化のプラスとマイナスが反転する仕組みを備えている。磁性体には、ＨaとＨbの値が異なる様々なヒステロンが存在している。
【００２１】
ここで、あるＨaとＨbを持つヒステロンの磁化状態を表す関数(以下「磁化状態関数」と表記)ｍ(Ｈa,Ｈb)を下記式
【数３】
と定義すると、プライザッハ図表は、図４に示すように、横軸をＨa、縦軸をＨbにとる磁化状態関数ｍ(Ｈa,Ｈb)の分布図となる。完全に消磁された状態であれば、プラスとマイナスのヒステロンは同じ面積となる。
【００２２】
図５は、一例としてＨs＝５としたときのプライザッハ図表とヒステリシス曲線の関係を示す説明図である。完全消磁状態の状態Ｉからプラスの飽和状態である状態IIに変化した場合、０＜Ｈa＜５の範囲がプラスに変わる。次に状態IIから状態IIIに変化した場合、−３＜Ｈb＜５の範囲がマイナスに変わる。さらに状態IIIから状態IVに変化した場合、−３＜Ｈa＜−１の範囲がプラスに変わる。このようにプライザッハ図表は、磁性体のヒステリシス特性に関係なく、印加磁場のみで一意的に求めることができる。
【００２３】
(３．プライザッハ分布図)
プライザッハ分布図は、上述のプライザッハ図表とともに、本実施の形態で船体の永久磁気を推定するのに用いる。あるＨaとＨbをもったヒステロンがどのくらいの大きさで分布しているかを表す関数をプライザッハ分布関数Ｐ(Ｈa,Ｈb)と定義すると、プライザッハ分布図は、横軸をＨa、縦軸をＨbにとるプライザッハ分布関数Ｐ(Ｈa,Ｈb)の分布図となる。図６は、磁性体の例示的なプライザッハ分布図である。色の濃い領域ほど多くのヒステロンが存在することを示す。図７は、立体的に示したプライザッハ分布図である。プライザッハ分布図は、磁性体の材料固有のものであり、磁性体の複数個の磁化特性の計測データ(図８参照)により算出される。
【００２４】
プライザッハ分布関数Ｐ(Ｈa,Ｈb)の導出について説明する。図９(Ａ)は、印加磁場をＨs→−Ｈs→１０と変化させたときの磁性体のプライザッハ図表及び磁化特性図である。同図(Ｂ)は、印加磁場をＨs→２→１０と変化させたときの磁性体のプライザッハ図表及び磁化特性図である。なお、これらの図では便宜上、飽和磁界Ｈsを「１１」、−Ｈsを「１」としている。また、プライザッハ図表における横軸Ｉ及び縦軸Ｊは、ヒステロンのＨa及びＨbにそれぞれ対応する。
【００２５】
図９(Ａ),(Ｂ)を対比すると、１０≦Ｉ＜１１、１≦Ｊ＜２の範囲内のヒステロンの磁化方向は反転(プラスとマイナスが逆転)しているのに対し、他の範囲のヒステロンの磁化方向は変化していない。したがって、図９(Ａ),(Ｂ)の状態間における磁化の差ΔＭすなわちＭ(１０,２)−Ｍ(１０,１)は、１０≦Ｉ＜１１,１≦Ｊ＜２の範囲内のヒステロンの磁化方向が反転したことによるものであり、下記式
【数４】
と表される。なお、Ｍ(Ｉ,Ｊ)は、印加磁場をＨs→Ｊ→Ｉと変化させたときの磁性体の永久磁気である。これを任意の点(Ｉ,Ｊ)について一般化すると、下記式
【数５】
となる。したがって、任意の点(Ｉ,Ｊ)についてＰ(Ｉ,Ｊ)を一般化すれば、下記式
【数６】
となる。磁性体の複数個の磁化特性の計測データ(図８参照)を基に各点のＰ(Ｉ,Ｊ)を求めることで、図６及び図７に示すようなプライザッハ分布図を作成することができる。
【００２６】
(４．船体の永久磁気推定方法の流れ)
図１は、本発明の実施の形態に係る船体の永久磁気推定方法のフローチャートである。本方法では、船体と同じ鋼材からなる試験片に磁場を印加して前記試験片の複数の磁化特性を取得し、これに基づいて船体の鋼材のプライザッハ分布図を導出する(ＳＴ１〜ＳＴ３)。他方、船体の数値計算モデルを作成し、導出したプライザッハ分布図を基に数値計算モデルの各要素の永久磁気を推定し、さらに船体周辺の永久磁場を推定する(ＳＴ４〜ＳＴ９)。各ステップは、基本的に、コンピュータとソフトウェアの協働によって実現される。以下、詳細に説明する。
【００２７】
試験片の作成(ＳＴ１)… 脱磁対象の船体と同じ鋼材(磁性材)からなる試験片を作成する。図１０は、作成する試験片の例示的な形状説明図であり、(Ａ)は正面図、(Ｂ)は右側面図である。試験片は、例えば、内径３０±３ｍｍ、外径５０±５ｍｍ、厚さ１±０.５ｍｍ、質量約１０ｇのリング形状とする。
【００２８】
磁化特性の測定(ＳＴ２)… ＳＴ１で作成した試験片に印加コイル及び受信コイルをそれぞれ例えば５０回巻き、印加コイルへの印加電流値を変化させて(試験片への印加磁場を変化させて)、プライザッハ分布図を作成するために必要となる基本的な磁化特性を測定する。例えば、最初に試験片を成す鋼材の飽和磁界Ｈs(−３０×１０5ｎＴ)を印加した後、２.５×１０5ｎＴ単位で２５×１０5ｎＴから−２５×１０5ｎＴまで磁場を印加し、再度飽和磁界Ｈsを印加する。これにより、図８に示されるような、複数パターン(この例では２１個)の磁化特性を測定する。
【００２９】
プライザッハ分布図の作成(ＳＴ３)… ＳＴ２で測定した複数パターンの磁化特性から、メジャーループ内全域の傾きを取得して、試験片を成す鋼材のプライザッハ分布関数Ｐ(Ｈa,Ｈb)を導出し、プライザッハ分布図を作成する。なお、Ｐ(Ｈa,Ｈb)の導出手順は上述のとおりである。Ｐ(Ｈa,Ｈb)は、例えば、−Ｈs≦Ｈa≦Ｈs,−Ｈs≦Ｈb≦Ｈsの範囲内で導出する。
【００３０】
船体の数値計算モデル作成(ＳＴ４)… 一方、有限要素法や積分方程式法といった公知の数値計算法に対応した船体数値計算モデルを作成する。図１１に船体の有限要素法用モデルを示す。
【００３１】
各要素における印加磁場計算(ＳＴ５)… ＳＴ４で作成した船体の数値計算モデルの各要素(船体のメッシュ)における印加磁場を計算する。なお、ここで計算する印加磁場は、実際の脱磁処理における、船体に装着した胴巻輪線(図２参照)への通電量に対応するものである。
【００３２】
各要素におけるプライザッハ図表の作成(ＳＴ６)… 船体の数値計算モデルの各要素について、ＳＴ５で計算した印加磁場に対応する磁化状態関数ｍ(Ｈa,Ｈb)を導出し、プライザッハ図表を作成する。ｍ(Ｈa,Ｈb)は、例えば、−Ｈs≦Ｈa≦Ｈs,−Ｈs≦Ｈb≦Ｈsの範囲内で導出する。
【００３３】
各要素における永久磁気推定(ＳＴ７)… ＳＴ３で作成したプライザッハ分布図と、ＳＴ６で作成したプライザッハ図表とに基づいて、船体の数値計算モデルの各要素における永久磁気を算出する。具体的には、下記式
【数７】
で表される永久磁気推定値Ｍを計算する。
【００３４】
船体周辺の永久磁場推定(ＳＴ８)… ＳＴ７で計算した各要素の永久磁気推定値Ｍに基づいて、有限要素法や積分方程式法といった公知の数値計算法により、船体周辺の永久磁場を推定する。
【００３５】
ＳＴ５〜ＳＴ８は、実際の脱磁処理において胴巻輪線(図２参照)に電流の極性を正負に変えながらその大きさを漸減して通電する場合の各通電量に対して順番に行う。例えば、通電量＋１０Ａ,−８Ａ,＋６Ａ,−４Ａ,＋３Ａ,−２Ａ,＋１Ａ,−０.５Ａ,＋０.５Ａの各々に対して順番に、ＳＴ５〜ＳＴ８を繰り返す。ここで、必要に応じて各通電量に対応した船体周辺の永久磁場推定値(すなわち途中経過)を出力する。例えば、＋１０Ａを通電した後に０Ａとしたときの永久磁場推定値、−８Ａを通電した後に０Ａとしたときの永久磁場推定値、・・・、−０.５Ａを通電した後に０Ａとしたときの永久磁場推定値を出力する。
【００３６】
磁気処理後の船体周辺の永久磁場出力(ＳＴ９)… 最終的に通電量を０Ａ(印加磁場なし)とした場合すなわち漸減した最後の通電量(例えば＋０.５Ａ)の後に０Ａとした場合の船体周辺の永久磁場推定値を、磁気処理後の船体周辺の永久磁場推定値として出力する。
【００３７】
本実施の形態によれば、下記の効果を奏することができる。
【００３８】
(1) 磁性材をヒステロンと呼ばれる磁気双極子の集合体と仮定し、脱磁対象の船体と同じ鋼材(磁性材)からなる試験片を用いて前記船体の鋼材のプライザッハ分布関数Ｐ(Ｈa,Ｈb)を導出し、実際の脱磁処理における胴巻輪線への通電によって印加される磁場に対応する磁化状態関数ｍ(Ｈa,Ｈb)を船体の数値計算モデルの各要素について導出し、導出したプライザッハ分布関数Ｐ(Ｈa,Ｈb)及び磁化状態関数ｍ(Ｈa,Ｈb)に基づいて各要素における永久磁気を推定計算するので、実際の脱磁処理を行う前に脱磁処理中及び脱磁処理後の船体の永久磁気及び周辺の永久磁場を推定することが可能となる。
【００３９】
(2) 上記により、脱磁処理を行う前にあらかじめ、脱磁処理後の船体の永久磁気及び周辺の永久磁場を最小にするように磁気処理電流値や胴巻輪線の配置箇所の最適化を図ることができる。
【００４０】
以上、実施の形態を例に本発明を説明したが、実施の形態の各構成要素や各処理プロセスには請求項に記載の範囲で種々の変形が可能であることは当業者に理解されるところである。
【符号の説明】
【００４１】
１ 船体
２ 係留用ワイヤ
３ ドルフィン
４ 胴巻輪線
５ 磁気処理用電源
６ 磁気センサ
７ 測定室
【技術分野】
【０００１】
本発明は、船体に磁場を印加したときの船体の永久磁気を推定する、船体の永久磁気推定方法に関する。
【背景技術】
【０００２】
船舶は多種多様な鋼材で構成されており、木製やＦＲＰ製等の小型船舶を除けばほとんど磁性体で占められている。船体磁気は、静磁界と動磁界で構成される。ただし、動磁界は静磁界に比べて非常に小さく、問題になることは少ない。静磁界は、誘導磁気と永久磁気に区別できる。誘導磁気は、船体が地磁気により磁化されることにより発生する。永久磁気は、船舶の建造時の溶接、裁断等により着磁するとともに、船体が船台に長期的に同じ方向に置かれると地磁気によって着磁する。さらに、就役後も波浪や地磁気を受けることによって永久磁気は変化する。したがって、就役直後は船体は特に大きな永久磁気を帯びており、時間の経過によってそれはさらに大きくなる。これに対処するためには、定期的に船体の永久磁気を脱磁する必要がある。
【０００３】
従来、鋼鉄によって構成された船体が保持する永久磁気を消去するために、船体に幾巻きかの胴巻輪線を装着して、それらに電流の極性を正負に変えながら、その大きさを漸減して通電する(つまり磁場を印加する)脱磁処理を実施している。そして、脱磁処理後に複数個の磁気検出器で船体周辺の永久磁場を測定し、脱磁処理の効果を評価している。
【０００４】
なお、積分方程式法又は有限要素法といった数値解析手法では、船舶計算モデルを作成して、船舶に地磁気が印加されることによる誘導磁気を計算することが可能である。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、以上の技術では、脱磁処理を実施した後でなくては複数個の磁気検出器による船体の永久磁気を把握することができないため、船体の形状及び材料に対応した最適な脱磁処理を実施することはできない。
【０００６】
本発明はこうした状況を認識してなされたものであり、その目的は、脱磁処理を実施する前に、船体に磁場を印加したときの前記船体の永久磁気を推定することが可能な、船体の永久磁気推定方法を提供することにある。
【課題を解決するための手段】
【０００７】
本発明のある態様は、船体の永久磁気推定方法である。この方法は、
船体に磁場を印加したときの前記船体の永久磁気を推定する、船体の永久磁気推定方法であって、
船体の磁性材を下記式で定義される磁化状態関数ｍ(Ｈa,Ｈb)
【数１】
によって磁化状態が表される磁気双極子であるヒステロンの集合体と仮定した場合の前記ヒステロンの分布数を前記Ｈa及び前記Ｈbの関数として表したプライザッハ分布関数Ｐ(Ｈa,Ｈb)を、前記Ｈa及び前記Ｈbの所定範囲で取得する取得ステップと、
前記船体の数値計算モデルの各要素について、所定の印加磁場に対応した前記磁化状態関数ｍ(Ｈa,Ｈb)を、前記Ｈa及び前記Ｈbの所定範囲で導出する導出ステップと、
取得した前記プライザッハ分布関数Ｐ(Ｈa,Ｈb)と、導出した前記磁化状態関数ｍ(Ｈa,Ｈb)とに基づいて、各要素の永久磁気を推定する推定ステップとを有する。
【０００８】
ある態様の方法において、前記プライザッハ分布関数Ｐ(Ｈa,Ｈb)は、前記船体と同じ磁性材からなる試験片に磁場を印加して取得した前記試験片の複数の磁化特性に基づいて算出したものであるとよい。
【０００９】
ある態様の方法において、前記数値計算モデルは、積分方程式法用モデルであってもよい。
【００１０】
あるいは、前記数値計算モデルは、有限要素法用モデルであってもよい。
【００１１】
ある態様の方法において、前記導出ステップは、印加磁場の極性の正負を変えながら前記印加磁場の大きさを漸減させた場合の各々の印加磁場に対応した前記磁化状態関数ｍ(Ｈa,Ｈb)を導出するものであるとよい。
【００１２】
ある態様の方法において、前記推定ステップは、下記式で表される永久磁気推定値Ｍ
【数２】
を計算するものであるとよい。
【００１３】
なお、以上の構成要素の任意の組合せ、本発明の表現を装置やシステム、プログラムなどの間で変換したものもまた、本発明の態様として有効である。
【発明の効果】
【００１４】
本発明によれば、脱磁処理を実施する前に、船体に磁場を印加したときの前記船体の永久磁気を推定することが可能となる。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施の形態に係る船体の永久磁気推定方法のフローチャート。
【図２】船体の脱磁処理の模式的説明図。
【図３】磁気双極子であるヒステロンが無数に集まって構成するメジャーループの説明図。
【図４】磁性体のプライザッハ図表。
【図５】一例としてＨs＝５としたときのプライザッハ図表とヒステリシス曲線の関係を示す説明図。
【図６】磁性体の例示的なプライザッハ分布図。
【図７】立体的に示したプライザッハ分布図。
【図８】磁性体の複数パターンの磁化特性図。
【図９】(Ａ)は、印加磁場をＨs→−Ｈs→１０と変化させたときの磁性体のプライザッハ図表及び磁化特性図。(Ｂ)は、印加磁場をＨs→２→１０と変化させたときの磁性体のプライザッハ図表及び磁化特性図。
【図１０】実施の形態で作成する試験片の例示的な形状説明図。
【図１１】船体の例示的な有限要素法用モデル。
【発明を実施するための形態】
【００１６】
以下、図面を参照しながら本発明の好適な実施の形態を詳述する。なお、各図面に示される同一または同等の構成要素、部材、処理等には同一の符号を付し、適宜重複した説明は省略する。また、実施の形態は発明を限定するものではなく例示であり、実施の形態に記述されるすべての特徴やその組み合わせは必ずしも発明の本質的なものであるとは限らない。
【００１７】
本実施の形態は、船体に幾巻きかの胴巻輪線を装着して通電する(磁場を印加する)脱磁処理に先だって、前記船体に磁場を印加したときの前記船体の永久磁気を推定するものである。説明の順序は下記のとおりとする。
１．実際の船体の脱磁処理
２．プライザッハ図表(永久磁気の推定に用いる)
３．プライザッハ分布図(永久磁気の推定に用いる)
４．船体の永久磁気推定方法の流れ
【００１８】
(１．実際の船体の脱磁処理)
図２は、船体の脱磁処理の模式的説明図である。本図に示されるように、脱磁処理の際には、船体１を係留用ワイヤ２とドルフィン３で固定した後に、船体１の外周に幾巻きかの胴巻輪線４を装着する。そして、磁気処理用電源５から磁気処理電流が通電され、船体の永久磁気を変化させる。船体１の下方の海底に設置される多数の磁気センサ６により、磁気処理電流通電後における船体周辺の磁場を計測し、コンピュータを有する測定室７に船体周辺の磁気波形が表示され、次の磁気処理電流が通電される。磁気処理電流値は、プラスマイナス交互に漸減しながら通電し、最終的に残留した永久磁場は磁気センサ６によって把握する。従来、磁気処理電流値は過去の実績を基に決めていた。本実施の形態では、磁気処理電流値を適切に定めるために、船体の脱磁処理に先だって、前記船体に磁場を印加したときの前記船体の永久磁気を推定するものである。
【００１９】
(２．プライザッハ図表)
プライザッハ図表は、本実施の形態で船体の永久磁気を推定するのに用いる。一般に、磁性体の磁化特性を表すＢ−Ｈ曲線(磁束密度Ｂと印加磁界Ｈの関係)は、飽和磁場±Ｈsを印加した時のメジャーループと呼ばれる曲線で示されている。しかし、船舶の磁気処理において実際に印加する磁界は、船体に装着した胴巻輪線(図２参照)への初期の巨大な通電では飽和磁界となるものの、通電量は次第に漸減し、その磁気特性はマイナーループをたどることとなる。そのため、磁気処理中及び磁気処理後の船体の永久磁気や周辺磁場を推定計算するには、船体を成す鋼材のマイナーループを把握する必要がある。メジャーループやある指定した磁場を印加した時のマイナーループについては、あらかじめ、船体を成す鋼材と同じ材質の磁性体からなる試験片を作成し、ヒステリシス特性を計測することが可能である。しかし、任意に磁場をかけたときのマイナーループを全て事前に計測しておくことはできない。
【００２０】
そこで、任意のマイナーループを推定するために、図３に示すように、ヒステロンと呼ばれる磁気双極子が無数に集まりメジャーループを構成していると仮定する。つまり、磁性体をヒステロンの集合体と仮定する。ヒステロンは、印加磁界がＨaとＨb(−Ｈs＜Ｈb＜Ｈa＜Ｈs)を境界に磁化のプラスとマイナスが反転する仕組みを備えている。磁性体には、ＨaとＨbの値が異なる様々なヒステロンが存在している。
【００２１】
ここで、あるＨaとＨbを持つヒステロンの磁化状態を表す関数(以下「磁化状態関数」と表記)ｍ(Ｈa,Ｈb)を下記式
【数３】
と定義すると、プライザッハ図表は、図４に示すように、横軸をＨa、縦軸をＨbにとる磁化状態関数ｍ(Ｈa,Ｈb)の分布図となる。完全に消磁された状態であれば、プラスとマイナスのヒステロンは同じ面積となる。
【００２２】
図５は、一例としてＨs＝５としたときのプライザッハ図表とヒステリシス曲線の関係を示す説明図である。完全消磁状態の状態Ｉからプラスの飽和状態である状態IIに変化した場合、０＜Ｈa＜５の範囲がプラスに変わる。次に状態IIから状態IIIに変化した場合、−３＜Ｈb＜５の範囲がマイナスに変わる。さらに状態IIIから状態IVに変化した場合、−３＜Ｈa＜−１の範囲がプラスに変わる。このようにプライザッハ図表は、磁性体のヒステリシス特性に関係なく、印加磁場のみで一意的に求めることができる。
【００２３】
(３．プライザッハ分布図)
プライザッハ分布図は、上述のプライザッハ図表とともに、本実施の形態で船体の永久磁気を推定するのに用いる。あるＨaとＨbをもったヒステロンがどのくらいの大きさで分布しているかを表す関数をプライザッハ分布関数Ｐ(Ｈa,Ｈb)と定義すると、プライザッハ分布図は、横軸をＨa、縦軸をＨbにとるプライザッハ分布関数Ｐ(Ｈa,Ｈb)の分布図となる。図６は、磁性体の例示的なプライザッハ分布図である。色の濃い領域ほど多くのヒステロンが存在することを示す。図７は、立体的に示したプライザッハ分布図である。プライザッハ分布図は、磁性体の材料固有のものであり、磁性体の複数個の磁化特性の計測データ(図８参照)により算出される。
【００２４】
プライザッハ分布関数Ｐ(Ｈa,Ｈb)の導出について説明する。図９(Ａ)は、印加磁場をＨs→−Ｈs→１０と変化させたときの磁性体のプライザッハ図表及び磁化特性図である。同図(Ｂ)は、印加磁場をＨs→２→１０と変化させたときの磁性体のプライザッハ図表及び磁化特性図である。なお、これらの図では便宜上、飽和磁界Ｈsを「１１」、−Ｈsを「１」としている。また、プライザッハ図表における横軸Ｉ及び縦軸Ｊは、ヒステロンのＨa及びＨbにそれぞれ対応する。
【００２５】
図９(Ａ),(Ｂ)を対比すると、１０≦Ｉ＜１１、１≦Ｊ＜２の範囲内のヒステロンの磁化方向は反転(プラスとマイナスが逆転)しているのに対し、他の範囲のヒステロンの磁化方向は変化していない。したがって、図９(Ａ),(Ｂ)の状態間における磁化の差ΔＭすなわちＭ(１０,２)−Ｍ(１０,１)は、１０≦Ｉ＜１１,１≦Ｊ＜２の範囲内のヒステロンの磁化方向が反転したことによるものであり、下記式
【数４】
と表される。なお、Ｍ(Ｉ,Ｊ)は、印加磁場をＨs→Ｊ→Ｉと変化させたときの磁性体の永久磁気である。これを任意の点(Ｉ,Ｊ)について一般化すると、下記式
【数５】
となる。したがって、任意の点(Ｉ,Ｊ)についてＰ(Ｉ,Ｊ)を一般化すれば、下記式
【数６】
となる。磁性体の複数個の磁化特性の計測データ(図８参照)を基に各点のＰ(Ｉ,Ｊ)を求めることで、図６及び図７に示すようなプライザッハ分布図を作成することができる。
【００２６】
(４．船体の永久磁気推定方法の流れ)
図１は、本発明の実施の形態に係る船体の永久磁気推定方法のフローチャートである。本方法では、船体と同じ鋼材からなる試験片に磁場を印加して前記試験片の複数の磁化特性を取得し、これに基づいて船体の鋼材のプライザッハ分布図を導出する(ＳＴ１〜ＳＴ３)。他方、船体の数値計算モデルを作成し、導出したプライザッハ分布図を基に数値計算モデルの各要素の永久磁気を推定し、さらに船体周辺の永久磁場を推定する(ＳＴ４〜ＳＴ９)。各ステップは、基本的に、コンピュータとソフトウェアの協働によって実現される。以下、詳細に説明する。
【００２７】
試験片の作成(ＳＴ１)… 脱磁対象の船体と同じ鋼材(磁性材)からなる試験片を作成する。図１０は、作成する試験片の例示的な形状説明図であり、(Ａ)は正面図、(Ｂ)は右側面図である。試験片は、例えば、内径３０±３ｍｍ、外径５０±５ｍｍ、厚さ１±０.５ｍｍ、質量約１０ｇのリング形状とする。
【００２８】
磁化特性の測定(ＳＴ２)… ＳＴ１で作成した試験片に印加コイル及び受信コイルをそれぞれ例えば５０回巻き、印加コイルへの印加電流値を変化させて(試験片への印加磁場を変化させて)、プライザッハ分布図を作成するために必要となる基本的な磁化特性を測定する。例えば、最初に試験片を成す鋼材の飽和磁界Ｈs(−３０×１０5ｎＴ)を印加した後、２.５×１０5ｎＴ単位で２５×１０5ｎＴから−２５×１０5ｎＴまで磁場を印加し、再度飽和磁界Ｈsを印加する。これにより、図８に示されるような、複数パターン(この例では２１個)の磁化特性を測定する。
【００２９】
プライザッハ分布図の作成(ＳＴ３)… ＳＴ２で測定した複数パターンの磁化特性から、メジャーループ内全域の傾きを取得して、試験片を成す鋼材のプライザッハ分布関数Ｐ(Ｈa,Ｈb)を導出し、プライザッハ分布図を作成する。なお、Ｐ(Ｈa,Ｈb)の導出手順は上述のとおりである。Ｐ(Ｈa,Ｈb)は、例えば、−Ｈs≦Ｈa≦Ｈs,−Ｈs≦Ｈb≦Ｈsの範囲内で導出する。
【００３０】
船体の数値計算モデル作成(ＳＴ４)… 一方、有限要素法や積分方程式法といった公知の数値計算法に対応した船体数値計算モデルを作成する。図１１に船体の有限要素法用モデルを示す。
【００３１】
各要素における印加磁場計算(ＳＴ５)… ＳＴ４で作成した船体の数値計算モデルの各要素(船体のメッシュ)における印加磁場を計算する。なお、ここで計算する印加磁場は、実際の脱磁処理における、船体に装着した胴巻輪線(図２参照)への通電量に対応するものである。
【００３２】
各要素におけるプライザッハ図表の作成(ＳＴ６)… 船体の数値計算モデルの各要素について、ＳＴ５で計算した印加磁場に対応する磁化状態関数ｍ(Ｈa,Ｈb)を導出し、プライザッハ図表を作成する。ｍ(Ｈa,Ｈb)は、例えば、−Ｈs≦Ｈa≦Ｈs,−Ｈs≦Ｈb≦Ｈsの範囲内で導出する。
【００３３】
各要素における永久磁気推定(ＳＴ７)… ＳＴ３で作成したプライザッハ分布図と、ＳＴ６で作成したプライザッハ図表とに基づいて、船体の数値計算モデルの各要素における永久磁気を算出する。具体的には、下記式
【数７】
で表される永久磁気推定値Ｍを計算する。
【００３４】
船体周辺の永久磁場推定(ＳＴ８)… ＳＴ７で計算した各要素の永久磁気推定値Ｍに基づいて、有限要素法や積分方程式法といった公知の数値計算法により、船体周辺の永久磁場を推定する。
【００３５】
ＳＴ５〜ＳＴ８は、実際の脱磁処理において胴巻輪線(図２参照)に電流の極性を正負に変えながらその大きさを漸減して通電する場合の各通電量に対して順番に行う。例えば、通電量＋１０Ａ,−８Ａ,＋６Ａ,−４Ａ,＋３Ａ,−２Ａ,＋１Ａ,−０.５Ａ,＋０.５Ａの各々に対して順番に、ＳＴ５〜ＳＴ８を繰り返す。ここで、必要に応じて各通電量に対応した船体周辺の永久磁場推定値(すなわち途中経過)を出力する。例えば、＋１０Ａを通電した後に０Ａとしたときの永久磁場推定値、−８Ａを通電した後に０Ａとしたときの永久磁場推定値、・・・、−０.５Ａを通電した後に０Ａとしたときの永久磁場推定値を出力する。
【００３６】
磁気処理後の船体周辺の永久磁場出力(ＳＴ９)… 最終的に通電量を０Ａ(印加磁場なし)とした場合すなわち漸減した最後の通電量(例えば＋０.５Ａ)の後に０Ａとした場合の船体周辺の永久磁場推定値を、磁気処理後の船体周辺の永久磁場推定値として出力する。
【００３７】
本実施の形態によれば、下記の効果を奏することができる。
【００３８】
(1) 磁性材をヒステロンと呼ばれる磁気双極子の集合体と仮定し、脱磁対象の船体と同じ鋼材(磁性材)からなる試験片を用いて前記船体の鋼材のプライザッハ分布関数Ｐ(Ｈa,Ｈb)を導出し、実際の脱磁処理における胴巻輪線への通電によって印加される磁場に対応する磁化状態関数ｍ(Ｈa,Ｈb)を船体の数値計算モデルの各要素について導出し、導出したプライザッハ分布関数Ｐ(Ｈa,Ｈb)及び磁化状態関数ｍ(Ｈa,Ｈb)に基づいて各要素における永久磁気を推定計算するので、実際の脱磁処理を行う前に脱磁処理中及び脱磁処理後の船体の永久磁気及び周辺の永久磁場を推定することが可能となる。
【００３９】
(2) 上記により、脱磁処理を行う前にあらかじめ、脱磁処理後の船体の永久磁気及び周辺の永久磁場を最小にするように磁気処理電流値や胴巻輪線の配置箇所の最適化を図ることができる。
【００４０】
以上、実施の形態を例に本発明を説明したが、実施の形態の各構成要素や各処理プロセスには請求項に記載の範囲で種々の変形が可能であることは当業者に理解されるところである。
【符号の説明】
【００４１】
１ 船体
２ 係留用ワイヤ
３ ドルフィン
４ 胴巻輪線
５ 磁気処理用電源
６ 磁気センサ
７ 測定室
【特許請求の範囲】
【請求項１】
船体に磁場を印加したときの前記船体の永久磁気を推定する、船体の永久磁気推定方法であって、
船体の磁性材を下記式で定義される磁化状態関数ｍ(Ｈa,Ｈb)
【数８】
によって磁化状態が表される磁気双極子であるヒステロンの集合体と仮定した場合の前記ヒステロンの分布数を前記Ｈa及び前記Ｈbの関数として表したプライザッハ分布関数Ｐ(Ｈa,Ｈb)を、前記Ｈa及び前記Ｈbの所定範囲で取得する取得ステップと、
前記船体の数値計算モデルの各要素について、所定の印加磁場に対応した前記磁化状態関数ｍ(Ｈa,Ｈb)を、前記Ｈa及び前記Ｈbの所定範囲で導出する導出ステップと、
取得した前記プライザッハ分布関数Ｐ(Ｈa,Ｈb)と、導出した前記磁化状態関数ｍ(Ｈa,Ｈb)とに基づいて、各要素の永久磁気を推定する推定ステップとを有する、船体の永久磁気推定方法。
【請求項２】
請求項１に記載の方法において、前記プライザッハ分布関数Ｐ(Ｈa,Ｈb)は、前記船体と同じ磁性材からなる試験片に磁場を印加して取得した前記試験片の複数の磁化特性に基づいて算出したものである、船体の永久磁気推定方法。
【請求項３】
請求項１又は２に記載の方法において、前記数値計算モデルが積分方程式法用モデルである、船体の永久磁気推定方法。
【請求項４】
請求項１又は２に記載の方法において、前記数値計算モデルが有限要素法用モデルである、船体の永久磁気推定方法。
【請求項５】
請求項１から４のいずれかに記載の方法において、前記導出ステップは、印加磁場の極性の正負を変えながら前記印加磁場の大きさを漸減させた場合の各々の印加磁場に対応した前記磁化状態関数ｍ(Ｈa,Ｈb)を導出するものである、船体の永久磁気推定方法。
【請求項６】
請求項１から５のいずれかに記載の方法において、前記推定ステップは、下記式で表される永久磁気推定値Ｍ
【数９】
を計算するものである、船体の永久磁気推定方法。
【請求項１】
船体に磁場を印加したときの前記船体の永久磁気を推定する、船体の永久磁気推定方法であって、
船体の磁性材を下記式で定義される磁化状態関数ｍ(Ｈa,Ｈb)
【数８】
によって磁化状態が表される磁気双極子であるヒステロンの集合体と仮定した場合の前記ヒステロンの分布数を前記Ｈa及び前記Ｈbの関数として表したプライザッハ分布関数Ｐ(Ｈa,Ｈb)を、前記Ｈa及び前記Ｈbの所定範囲で取得する取得ステップと、
前記船体の数値計算モデルの各要素について、所定の印加磁場に対応した前記磁化状態関数ｍ(Ｈa,Ｈb)を、前記Ｈa及び前記Ｈbの所定範囲で導出する導出ステップと、
取得した前記プライザッハ分布関数Ｐ(Ｈa,Ｈb)と、導出した前記磁化状態関数ｍ(Ｈa,Ｈb)とに基づいて、各要素の永久磁気を推定する推定ステップとを有する、船体の永久磁気推定方法。
【請求項２】
請求項１に記載の方法において、前記プライザッハ分布関数Ｐ(Ｈa,Ｈb)は、前記船体と同じ磁性材からなる試験片に磁場を印加して取得した前記試験片の複数の磁化特性に基づいて算出したものである、船体の永久磁気推定方法。
【請求項３】
請求項１又は２に記載の方法において、前記数値計算モデルが積分方程式法用モデルである、船体の永久磁気推定方法。
【請求項４】
請求項１又は２に記載の方法において、前記数値計算モデルが有限要素法用モデルである、船体の永久磁気推定方法。
【請求項５】
請求項１から４のいずれかに記載の方法において、前記導出ステップは、印加磁場の極性の正負を変えながら前記印加磁場の大きさを漸減させた場合の各々の印加磁場に対応した前記磁化状態関数ｍ(Ｈa,Ｈb)を導出するものである、船体の永久磁気推定方法。
【請求項６】
請求項１から５のいずれかに記載の方法において、前記推定ステップは、下記式で表される永久磁気推定値Ｍ
【数９】
を計算するものである、船体の永久磁気推定方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【公開番号】特開２０１１−６５６１２（Ｐ２０１１−６５６１２Ａ）
【公開日】平成２３年３月３１日（２０１１．３．３１）
【国際特許分類】
【出願番号】特願２００９−２１８２１２（Ｐ２００９−２１８２１２）
【出願日】平成２１年９月１９日（２００９．９．１９）
【出願人】（３９００１４３０６）防衛省技術研究本部長 (169)
【Ｆターム（参考）】
【公開日】平成２３年３月３１日（２０１１．３．３１）
【国際特許分類】
【出願日】平成２１年９月１９日（２００９．９．１９）
【出願人】（３９００１４３０６）防衛省技術研究本部長 (169)
【Ｆターム（参考）】
[ Back to top ]