JP7299615B2 - Material analysis system, material analysis method, and material analysis program - Google Patents

Material analysis system, material analysis method, and material analysis program Download PDF

Info

Publication number
JP7299615B2
JP7299615B2 JP2019161138A JP2019161138A JP7299615B2 JP 7299615 B2 JP7299615 B2 JP 7299615B2 JP 2019161138 A JP2019161138 A JP 2019161138A JP 2019161138 A JP2019161138 A JP 2019161138A JP 7299615 B2 JP7299615 B2 JP 7299615B2
Authority
JP
Japan
Prior art keywords
transport
image
microstructure
model
phase
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
Application number
JP2019161138A
Other languages
Japanese (ja)
Other versions
JP2021039021A (en
Inventor
明保 山本
拓也 小原
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NATIONAL UNIVERSITY CORPORATION TOKYO UNIVERSITY OF AGRICULUTURE & TECHNOLOGY
Original Assignee
NATIONAL UNIVERSITY CORPORATION TOKYO UNIVERSITY OF AGRICULUTURE & TECHNOLOGY
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NATIONAL UNIVERSITY CORPORATION TOKYO UNIVERSITY OF AGRICULUTURE & TECHNOLOGY filed Critical NATIONAL UNIVERSITY CORPORATION TOKYO UNIVERSITY OF AGRICULUTURE & TECHNOLOGY
Priority to JP2019161138A priority Critical patent/JP7299615B2/en
Publication of JP2021039021A publication Critical patent/JP2021039021A/en
Application granted granted Critical
Publication of JP7299615B2 publication Critical patent/JP7299615B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Superconductor Devices And Manufacturing Methods Thereof (AREA)

Description

本開示の一側面は材料解析システム、材料解析方法、および材料解析プログラムに関する。 One aspect of the present disclosure relates to a material analysis system, material analysis method, and material analysis program.

材料の輸送特性を得るための手法が知られている。例えば、特許文献1には、電極内の活物質による電子伝導ネットワークを評価する方法が記載されている。この方法は、活物質層および集電体を含む断面の観察像と、該断面に存在する個々の活物質の導電性評価結果とを相関させることにより、電極内における電気的に活性な活物質の分布状態を測定するものである。 Techniques for obtaining the transport properties of materials are known. For example, Patent Literature 1 describes a method of evaluating an electron conduction network by an active material in an electrode. This method correlates an observed image of a cross section containing an active material layer and a current collector with the conductivity evaluation result of each active material present in the cross section to determine the electrically active active material in the electrode. It measures the distribution of

特開2014-81362号公報JP 2014-81362 A

特許文献1に記載の方法では、導電性評価結果を得るために活物質と集電体との間の導電性を測定する必要があり、したがって、その測定のための設備が必要である。そこで、材料の輸送特性を簡易に得るための手法が望まれている。 In the method described in Patent Document 1, it is necessary to measure the electrical conductivity between the active material and the current collector in order to obtain the electrical conductivity evaluation result, and therefore equipment for the measurement is required. Therefore, a technique for easily obtaining the transport properties of materials is desired.

本開示の一側面に係る材料解析システムは、少なくとも一つのプロセッサを備える。少なくとも一つのプロセッサは、材料の微細組織を写した材料画像を取得し、材料画像を解析することで観察対象の相の位置を識別し、相の位置に基づいて材料の微細組織に対応する輸送モデルを生成し、輸送モデルに基づいて輸送特性を算出し、輸送特性を示す輸送特性データを出力する。 A material analysis system according to one aspect of the present disclosure includes at least one processor. At least one processor acquires a material image depicting a microstructure of the material, analyzes the material image to identify the position of the observed phase, and determines transport corresponding to the material microstructure based on the phase position. A model is generated, transport characteristics are calculated based on the transport model, and transport characteristics data representing the transport characteristics are output.

このような側面においては、材料の微細組織を写した材料画像から、該微細組織に対応する輸送モデルが生成され、その輸送モデルに基づいて輸送特性が算出される。したがって、材料画像さえ準備すれば材料の輸送特性を簡易に得ることができる。 In this aspect, a transport model corresponding to the microstructure is generated from a material image representing the microstructure of the material, and the transport characteristics are calculated based on the transport model. Therefore, the transportation characteristics of the material can be easily obtained by preparing the material image.

本開示の一側面によれば、特別な設備を必要とすることなく材料の輸送特性を得ることができる。 According to one aspect of the present disclosure, material transport properties can be obtained without the need for special equipment.

実施形態に係る材料解析システムを構成するコンピュータのハードウェア構成の一例を示す図である。It is a figure which shows an example of the hardware constitutions of the computer which comprises the material analysis system which concerns on embodiment. 実施形態に係る材料解析システムの機能構成の一例を示す図である。It is a figure showing an example of functional composition of a material analysis system concerning an embodiment. 実施形態に係る材料解析システムの動作の一例を示すフローチャートである。4 is a flow chart showing an example of the operation of the material analysis system according to the embodiment; 輸送モデルを説明するための図である。It is a figure for demonstrating a transportation model. 輸送特性の具体的な計算手順の一例を示すフローチャートである。4 is a flow chart showing an example of a specific calculation procedure for transportation characteristics; 可視化された輸送特性データの一例を示す図である。It is a figure which shows an example of the visualized transport characteristic data. 図6の一部を拡大して示す図である。It is a figure which expands and shows a part of FIG. 実施形態に係る材料解析プログラムの構成の一例を示す図である。It is a figure which shows an example of a structure of the material-analysis program which concerns on embodiment. 実施例におけるサンプルAに関する画像を示す図である。FIG. 10 is a diagram showing an image of sample A in the example; 実施例におけるサンプルB1に関する画像を示す図である。FIG. 10 is a diagram showing an image regarding sample B1 in the example. 実施例におけるサンプルB2に関する画像を示す図である。FIG. 10 is a diagram showing an image regarding sample B2 in the example.

以下、添付図面を参照しながら本開示での実施形態を詳細に説明する。なお、図面の説明において同一または同等の要素には同一の符号を付し、重複する説明を省略する。 Hereinafter, embodiments of the present disclosure will be described in detail with reference to the accompanying drawings. In the description of the drawings, the same or equivalent elements are denoted by the same reference numerals, and overlapping descriptions are omitted.

[システムの構成]
実施形態に係る材料解析システム10は、材料の微細組織を写した画像を解析することで該材料の輸送特性を推定するコンピュータシステムである。材料とは、何らかの製品を作製する際に用いられる物質のことをいう。材料は限定されず、例えば超伝導体の材料、すなわち超伝導材料でもよいし、電池の材料(例えば電極の材料)でもよいし、熱伝導材料でもよいし、鉄筋コンクリート等の構造物材料でもよい。材料の微細組織とは、マイクロメートル、ナノメートルなどの小さい単位のスケールの分解能によって認識することができる微細な構造のことをいう。画像は、例えば、顕微鏡により得られる画像、すなわち顕微鏡画像のことをいう。顕微鏡の種類は限定されず、例えば、電子顕微鏡、走査型プローブ顕微鏡、X線顕微鏡、超音波顕微鏡、または光学顕微鏡でもよい。本開示では、材料の微細組織を写した画像を「材料画像」ともいう。
[System configuration]
A material analysis system 10 according to an embodiment is a computer system that estimates transport characteristics of a material by analyzing an image of the microstructure of the material. A material is a substance that is used in making some product. The material is not limited, and may be, for example, a superconductor material, that is, a superconducting material, a battery material (for example, an electrode material), a thermally conductive material, or a structural material such as reinforced concrete. The microstructure of a material refers to a fine structure that can be recognized with a resolution of a small unit scale such as micrometers or nanometers. An image refers to, for example, an image obtained by a microscope, that is, a microscopic image. The type of microscope is not limited, and may be, for example, an electron microscope, a scanning probe microscope, an X-ray microscope, an ultrasonic microscope, or an optical microscope. In the present disclosure, an image showing the microstructure of a material is also referred to as a "material image".

輸送特性とは、材料中における、輸送媒体を介した物理量の移動に関する特有の性質のことをいう。物理量および輸送媒体は限定されない。例えば、物理量は質量、電荷、エネルギ、熱、または運動量でもよく、これに対応して、輸送媒体は電子、イオン、流体、または格子振動でもよい。したがって、輸送特性は電気伝導度、イオン伝導度、熱伝導度などの様々な指標によって表現され得る。材料解析システム10は、輸送媒体の生成および消滅が生じない限り、任意の材料および現象について輸送特性を推定することができる。 Transport properties refer to the specific properties of a material regarding the movement of physical quantities through a transport medium. The physical quantity and transport medium are not limited. For example, the physical quantity may be mass, charge, energy, heat, or momentum, and correspondingly the transport medium may be electrons, ions, fluids, or lattice vibrations. Therefore, transport properties can be expressed by various indices such as electrical conductivity, ionic conductivity, and thermal conductivity. The material analysis system 10 is capable of estimating transport properties for arbitrary materials and phenomena as long as the creation and destruction of transport media do not occur.

材料解析システム10は1台のコンピュータで構成されてもよいし、複数台のコンピュータで構成されてもよい。複数台のコンピュータを用いる場合には、これらのコンピュータがインターネット、イントラネットなどの通信ネットワークを介して接続されることで、論理的に一つの材料解析システム10が構築される。 The material analysis system 10 may be composed of one computer, or may be composed of a plurality of computers. When using a plurality of computers, these computers are connected via a communication network such as the Internet or an intranet to logically construct one material analysis system 10 .

図1は材料解析システム10を構成するコンピュータ100のハードウェア構成の一例を示す図である。コンピュータ100は、オペレーティングシステム、アプリケーション・プログラムなどを実行するプロセッサ(例えばCPU)101と、ROMおよびRAMで構成される主記憶部102と、ハードディスク、フラッシュメモリなどで構成される補助記憶部103と、ネットワークカードまたは無線通信モジュールで構成される通信制御部104と、キーボード、マウスなどの入力装置105と、モニタなどの出力装置106とを備える。 FIG. 1 is a diagram showing an example of the hardware configuration of a computer 100 that constitutes the material analysis system 10. As shown in FIG. The computer 100 includes a processor (e.g., CPU) 101 that executes an operating system, application programs, etc., a main memory unit 102 composed of a ROM and a RAM, an auxiliary memory unit 103 composed of a hard disk, a flash memory, etc. It has a communication control unit 104 configured by a network card or a wireless communication module, an input device 105 such as a keyboard and a mouse, and an output device 106 such as a monitor.

材料解析システム10の各機能要素は、プロセッサ101または主記憶部102の上に所定のソフトウェア(例えば、後述する材料解析プログラムP1)を読み込ませてプロセッサ101にそのソフトウェアを実行させることで実現される。プロセッサ101はそのソフトウェアに従って、通信制御部104、入力装置105、または出力装置106を動作させ、主記憶部102または補助記憶部103におけるデータの読み出しおよび書き込みを行う。処理に必要なデータまたはデータベースは主記憶部102または補助記憶部103内に格納される。 Each functional element of the material analysis system 10 is realized by causing the processor 101 or the main storage unit 102 to load predetermined software (for example, a material analysis program P1 to be described later) and causing the processor 101 to execute the software. . The processor 101 operates the communication control unit 104, the input device 105, or the output device 106 according to the software, and reads and writes data in the main storage unit 102 or the auxiliary storage unit 103. FIG. Data or databases necessary for processing are stored in the main memory unit 102 or the auxiliary memory unit 103 .

図2は材料解析システム10の機能構成の一例を示す図である。材料解析システム10は機能的構成要素として画像解析部11、モデル生成部12、および特性算出部13を備える。画像解析部11は、材料画像を解析することで観察対象の相の位置を識別する機能要素である。観察対象とは、微細組織の中で、輸送特性を観察するために着目する物質のことをいう。相とは、化学的組成および物理的状態が均一であることによって他と区別される部分のことをいう。モデル生成部12はその相の位置に基づいて、材料の微細組織に対応する輸送モデルを生成する機能要素である。輸送モデルとは、微細組織おいて物理量が伝達される系を抽象的に表現するモデルのことをいう。すなわち、輸送モデルは、微細組織に対応する物理量に関するモデルである。特性算出部13はその輸送モデルに基づいて輸送特性を算出する機能要素である。 FIG. 2 is a diagram showing an example of the functional configuration of the material analysis system 10. As shown in FIG. The material analysis system 10 includes an image analysis section 11, a model generation section 12, and a property calculation section 13 as functional components. The image analysis unit 11 is a functional element that identifies the position of the phase to be observed by analyzing the material image. Observation target refers to a substance of interest for observing transport properties in a microstructure. A phase is a portion that is distinguished by uniformity in chemical composition and physical state. The model generator 12 is a functional element that generates a transport model corresponding to the microstructure of the material based on the position of the phase. A transport model is a model that abstractly expresses a system in which physical quantities are transmitted in a microstructure. That is, the transport model is a model relating to physical quantities corresponding to microstructures. The characteristic calculator 13 is a functional element that calculates the transport characteristic based on the transport model.

[システムの動作]
材料解析システム10の動作を説明するとともに本実施形態に係る材料解析方法について説明する。図3は材料解析システム10の動作の一例を処理フローS1として示すフローチャートである。
[System operation]
The operation of the material analysis system 10 will be described, and the material analysis method according to this embodiment will be described. FIG. 3 is a flow chart showing an example of the operation of the material analysis system 10 as a processing flow S1.

ステップS11では、画像解析部11が材料画像を取得する。材料画像の取得方法は限定されない。例えば、画像解析部11は顕微鏡から入力された材料画像を取得してもよい。あるいは、画像解析部11は所与のメモリまたはデータベースにアクセスしてその記憶装置から材料画像を読み出してもよい。あるいは、画像解析部11は他のコンピュータから送信された材料画像を受信してもよい。あるいは、画像解析部11はユーザ操作により入力された材料画像を取得してもよい。一例では、材料画像は256階調のグレースケール画像である。 In step S11, the image analysis unit 11 acquires a material image. The acquisition method of the material image is not limited. For example, the image analysis unit 11 may acquire material images input from a microscope. Alternatively, the image analysis unit 11 may access a given memory or database and read the material image from that storage device. Alternatively, the image analysis section 11 may receive a material image transmitted from another computer. Alternatively, the image analysis unit 11 may acquire a material image input by a user's operation. In one example, the material image is a 256-level grayscale image.

ステップS12では、画像解析部11がその材料画像を解析することで観察対象の相の位置を識別する。画像解析部11は材料画像に応じた任意の画像解析手法を用いて相の位置を識別する。画像解析手法は限定されない。一例として、材料画像が256階調のグレースケール画像であれば、画像解析部11は、画像の階調を下げる技術であるポスタリゼーションを用いて材料画像の個々の画素をいくつかの数のグループに分類することで、観察対象を含むいくつかの相の位置を識別してもよい。材料が複数の物質で構成される場合には、画像解析部11はその複数の物質のうちの一つを観察対象として処理して、その観察対象の相の位置を識別してもよい。 In step S12, the image analysis unit 11 analyzes the material image to identify the position of the phase to be observed. The image analysis unit 11 identifies the position of the phase using any image analysis method suitable for the material image. The image analysis method is not limited. As an example, if the material image is a 256-gradation grayscale image, the image analysis unit 11 divides the individual pixels of the material image into several groups using posterization, which is a technique for lowering the gradation of the image. Classification may identify the location of several phases containing the observed object. When the material is composed of a plurality of substances, the image analysis unit 11 may process one of the plurality of substances as an observation target and identify the position of the phase of the observation target.

一例では、画像解析部11は観察対象の相の位置と、観察対象ではない相の位置とを識別する。あるいは、画像解析部11は観察対象の相の位置を識別し、残りの部分を、観察対象ではない相として識別してもよい。観察対象ではない相の例として不純物相および空隙相が挙げられるが、これらに限定されない。画像解析部11は、それぞれの相を特定することで、個々の画素がどの相に属するかを認識できる。超伝導材料である場合には、一例では、画像解析部11は超伝導相を観察対象の相として処理してその位置を識別し、不純物相および空隙相のそれぞれについても位置を識別する。別の例として、画像解析部11は超伝導相の位置を識別し、残りの部分を超伝導相ではない相として識別してもよい。 In one example, the image analysis unit 11 identifies the positions of the phases to be observed and the positions of the phases not to be observed. Alternatively, the image analysis unit 11 may identify the position of the phase to be observed and identify the remaining portion as the phase not to be observed. Examples of phases that are not of interest include, but are not limited to, impurity phases and void phases. By specifying each phase, the image analysis unit 11 can recognize which phase each pixel belongs to. In the case of a superconducting material, in one example, the image analysis unit 11 treats the superconducting phase as a phase to be observed and identifies its position, and also identifies the positions of each of the impurity phase and the void phase. As another example, the image analysis unit 11 may identify the position of the superconducting phase and identify the remaining portion as a non-superconducting phase.

画像解析部11は、観察対象の相の位置に加えて、その相に対応する複数の画素のそれぞれについて結晶方位を特定してもよい。結晶方位とは、結晶の方向をミラー指数によって表現したものをいう。結晶方位の特定方法は限定されない。例えば、画像解析部11は、各画素の色に基づいて結晶方位を特定してもよい。この場合には、画像解析部11は、予め用意された、RGB値と結晶方位との対応関係を示すデータを用いる。画像解析部11は材料画像から個々の画素のRGB値を取得し、その対応関係を参照することで、それぞれの画素について結晶方位(θ、φ、ω)を特定する。あるいは、画像解析部11は、電子顕微鏡から得られるElectron BackScatter Diffraction(EBSD)データから直接に結晶方位を特定してもよい。 In addition to the position of the phase to be observed, the image analysis unit 11 may specify the crystal orientation for each of the plurality of pixels corresponding to the phase. "Crystal orientation" refers to the orientation of a crystal expressed by Miller indices. A method for specifying the crystal orientation is not limited. For example, the image analysis unit 11 may identify the crystal orientation based on the color of each pixel. In this case, the image analysis unit 11 uses previously prepared data indicating the correspondence between RGB values and crystal orientations. The image analysis unit 11 acquires the RGB values of each pixel from the material image, and refers to the corresponding relationship to identify the crystal orientation (θ, φ, ω) of each pixel. Alternatively, the image analysis unit 11 may specify the crystal orientation directly from Electron Backscatter Diffraction (EBSD) data obtained from an electron microscope.

ステップS13では、モデル生成部12が観察対象の相の位置に基づいて輸送モデルを生成する。一例では、モデル生成部12は、観察対象の相の位置に対応して閉サイトを設定し、観察対象ではない相の位置に対応して開サイトを設定することで、輸送モデルを生成する。閉サイトとは物理量が流れる領域のことをいい、開サイトとは物理量が流れない領域のことをいう。したがって、輸送モデルは個々の閉サイトの位置と個々の開サイトの位置とを示すモデルである。モデル生成部12は閉サイトおよび開サイトを、材料画像を構成する複数の画素のそれぞれに対応して設定する。材料が超伝導材料であれば、モデル生成部12は、電流の経路である超伝導相の位置に対応して閉サイトを設定し、超伝導相以外の相の位置(例えば、不純物相および空隙相)の位置に対応して開サイトを設定することで、輸送モデルを生成する。輸送モデルの表現方法は限定されない。例えば、モデル生成部12は基礎方程式、物性パラメータなどの様々な手法または定義を用いて輸送モデルを表現してよい。 In step S13, the model generation unit 12 generates a transport model based on the position of the observation target phase. In one example, the model generation unit 12 generates a transport model by setting a closed site corresponding to the position of the phase to be observed and setting an open site corresponding to the position of the phase not to be observed. A closed site means a region through which a physical quantity flows, and an open site means a region where no physical quantity flows. A transport model is therefore a model that indicates the position of each closed site and the position of each open site. The model generation unit 12 sets closed sites and open sites corresponding to each of a plurality of pixels forming the material image. If the material is a superconducting material, the model generator 12 sets a closed site corresponding to the position of the superconducting phase, which is the current path, and the position of the phase other than the superconducting phase (for example, impurity phases and voids). A transport model is generated by setting the open sites corresponding to the positions of the phases. The transportation model representation method is not limited. For example, the model generator 12 may express the transport model using various methods or definitions such as basic equations and physical parameters.

図4は輸送モデルを説明するための図である。図4の例(a)は、実際の微細組織に合わせて個々のサイトを3次元で表現したものであり、図4の例(b)は、材料画像に合わせてそれらのサイトを2次元で表現したものである。便宜のために、例(a),(b)のそれぞれには3軸の座標も示す。図4は、開サイト21、閉サイト22、および閉サイト23がこの順に1列に並び、y軸方向に沿って物理量29が閉サイト22から閉サイト23へと流れる例を示す。 FIG. 4 is a diagram for explaining the transportation model. Example (a) of FIG. 4 is a three-dimensional representation of individual sites matched to an actual microstructure, and example (b) of FIG. 4 is a two-dimensional representation of those sites matched to a material image. It is an expression. For convenience, 3-axis coordinates are also shown in each of examples (a) and (b). FIG. 4 shows an example in which an open site 21, a closed site 22, and a closed site 23 are arranged in this order, and a physical quantity 29 flows from the closed site 22 to the closed site 23 along the y-axis direction.

隣接するサイト間をボンドとして考えると、閉サイト同士を結ぶボンドの抵抗率を1と定義し、開サイトに接するボンドの抵抗率を無限大と定義することができる。閉サイト同士を結ぶボンドの抵抗率は1以外の値でもよい。モデル生成部12は複数のボンドのそれぞれの抵抗率を輸送モデルに含めてもよい。モデル生成部12は他の情報を輸送モデルに含めてもよく、例えば、観察対象の相に対応する結晶の最大輸送量を輸送モデルに含めてもよい。モデル生成部12は個々のサイトに対応する個々の結晶方位を輸送モデルに含めてもよく、この場合には、モデル生成部12はその結晶方位に依存する抵抗率を個々のサイトに設定する。モデル生成部12は、例えば、材料の異方性パラメータと、異方的ギンツブルグ・ランダウ方程式(ρ(θ、φ、ω),ρ(θ、φ、ω),ρ(θ、φ、ω))とに基づいてそれぞれのサイトの抵抗率を算出する。そして、モデル生成部12は、閉サイト同士を結ぶボンドの抵抗率を、隣接するサイトの抵抗率の平均値に設定する。モデル生成部12は、ボンドの抵抗率、結晶の最大輸送量、異方的ギンツブルグ・ランダウ方程式などの各種の参照データを所与の記憶部から読み出し、その参照データも用いて輸送モデルを生成することができる。 Considering adjacent sites as bonds, the resistivity of a bond connecting closed sites can be defined as 1, and the resistivity of a bond in contact with an open site can be defined as infinity. The resistivity of the bond connecting the closed sites may be a value other than one. The model generator 12 may include the resistivity of each of the multiple bonds in the transport model. The model generator 12 may include other information in the transport model, for example, the maximum transport amount of the crystal corresponding to the observed phase may be included in the transport model. The model generator 12 may include individual crystallographic orientations corresponding to individual sites in the transport model, in which case the model generator 12 sets resistivities for individual sites that depend on the crystallographic orientation. For example, the model generation unit 12 generates an anisotropic parameter of a material, an anisotropic Ginzburg-Landau equation (ρ x (θ, φ, ω), ρ y (θ, φ, ω), ρ z (θ, φ , ω)) and the resistivity of each site is calculated. Then, the model generator 12 sets the resistivity of the bond connecting the closed sites to the average resistivity of the adjacent sites. The model generation unit 12 reads out various reference data such as bond resistivity, crystal maximum transport amount, anisotropic Ginzburg-Landau equation, etc. from a given storage unit, and also uses the reference data to generate a transport model. be able to.

ステップS14では、特性算出部13がその輸送モデルに基づいて輸送特性を算出する。図5は輸送特性の具体的な計算手順の一例を示すフローチャートである。 In step S14, the characteristic calculator 13 calculates the transport characteristic based on the transport model. FIG. 5 is a flow chart showing an example of a specific calculation procedure for transportation characteristics.

ステップS141では、特性算出部13が、材料画像を構成する個々の画素を示す配列を定義する。これは、材料画像を構成する複数の画素に対応して複数の配列要素が設定されることを意味する。 In step S141, the property calculation unit 13 defines an array indicating individual pixels forming the material image. This means that a plurality of array elements are set corresponding to a plurality of pixels forming the material image.

ステップS142では、特性算出部13が輸送モデルの情報を配列に設定する。具体的には、特性算出部13はそれぞれの配列要素に、閉サイトか開サイトか否かを示す二値のサイト情報を設定する。或る一つの画素Pxに対応する配列要素がExであるとすると、特性算出部13は、画素Pxに対応する位置が閉サイトであれば、閉サイトを示す値を配列要素Exに設定し、その位置が開サイトであれば、開サイトを示す値を配列要素Exに設定する。輸送モデルが他の情報を含む場合には、特性算出部13はその情報の少なくとも一部を配列に設定してもよく、例えば結晶方位(θ、φ、ω)および抵抗率のうちの少なくとも一つを配列に格納してよい。 In step S142, the characteristic calculation unit 13 sets the information of the transport model in an array. Specifically, the characteristic calculation unit 13 sets binary site information indicating whether the site is a closed site or an open site for each array element. Assuming that the array element corresponding to one pixel Px is Ex, if the position corresponding to the pixel Px is a closed site, the characteristic calculation unit 13 sets a value indicating the closed site to the array element Ex, If the position is an open site, a value indicating the open site is set in the array element Ex. If the transport model includes other information, the property calculator 13 may set at least part of the information in the array, for example, at least one of crystal orientation (θ, φ, ω) and resistivity. can be stored in an array.

ステップS143では、特性算出部13がクラスタ解析を実行することでパーコレートクラスタを設定する。パーコレートクラスタとは、物理量が通過する可能性がある閉サイトの集合のことをいう。材料の微細組織の入力端から、該入力端の反対に位置する出力端への物理量の伝導特性が変化する現象(より具体的には、その伝導特性が飛躍的に上昇する現象)をパーコレーション転移という。 In step S143, the characteristic calculation unit 13 sets percolate clusters by performing cluster analysis. A percolate cluster is a collection of closed sites through which a physical quantity may pass. Percolation transition is a phenomenon in which the conduction characteristics of a physical quantity from the input end of the microstructure of a material to the output end located opposite to the input end change (more specifically, the phenomenon in which the conduction characteristics rise dramatically). It says.

ステップS144では、特性算出部13がそのパーコレートクラスタに基づいて各サイトでのポテンシャル量を計算する。ポテンシャル量とは、所与の基準位置に基づいて設定される物理量である。ポテンシャル量が具体的に何を示すかは物理量の種類に依存する。例えば、電気輸送特性を推定する場合には電位がポテンシャル量に相当し、したがって、特性算出部13はこのステップにおいて各サイトでの電位を計算する。特性算出部13は、輸送モデルの入力端に所与の物理量またはポテンシャル量(例えば電圧)を与えることを仮定する。そして、特性算出部13は複数のサイトのそれぞれについて、該サイトと、該サイトに隣接する1以上のサイトとの間の物理量の出入りについて、物理量の移動に関する一般的な基礎方程式を応用して、該サイトでのポテンシャル量を算出する。基礎方程式は限定されず、現象または物理量に応じて選択されてよい。例えば、電気に関する場合には、キルヒホッフの法則およびオームの法則の一方または双方が用いられてよい。キルヒホッフの法則は、一つのサイトに流入する物理量(例えば電流)の総和は、該サイトから流出する物理量(例えば電流)の総和に等しい、というものである。他の基礎方程式の例として、フーリエの法則、フィックの拡散の法則などが挙げられる。特性算出部13は、少なくともそれぞれの閉サイトについてポテンシャル量を算出する。 In step S144, the characteristic calculator 13 calculates the potential amount at each site based on the percolate cluster. A potential quantity is a physical quantity set based on a given reference position. What the potential quantity specifically indicates depends on the type of physical quantity. For example, when estimating the electrotransport property, the potential corresponds to the amount of potential, so the property calculator 13 calculates the potential at each site in this step. It is assumed that the characteristic calculator 13 gives a given physical quantity or potential quantity (for example, voltage) to the input terminal of the transport model. Then, for each of the plurality of sites, the characteristic calculation unit 13 applies a general basic equation regarding the movement of physical quantities with respect to the entry and exit of physical quantities between the site and one or more sites adjacent to the site, A potential amount at the site is calculated. The underlying equation is not limited and may be selected according to phenomena or physical quantities. For example, when it comes to electricity, one or both of Kirchhoff's law and Ohm's law may be used. Kirchhoff's law states that the sum of physical quantities (eg, current) flowing into one site is equal to the sum of physical quantities (eg, current) flowing out of the site. Examples of other fundamental equations include Fourier's law, Fick's law of diffusion, and the like. The characteristic calculation unit 13 calculates the potential amount at least for each closed site.

ステップS145では、特性算出部13が各サイトでのポテンシャル量に基づいて各断面での移動物理量を計算する。断面とは、材料の入力端から出力端に向かう方向と直交する断面のことをいう。移動物理量とは、或る断面からその隣の断面へと移動する物理量のことをいう。 In step S145, the characteristic calculator 13 calculates the movement physical quantity at each cross section based on the potential quantity at each site. A cross section refers to a cross section perpendicular to the direction from the input end of the material to the output end. A moving physical quantity is a physical quantity that moves from one cross section to the next cross section.

ステップS146では、特性算出部13が各断面の移動物理量が収束したか否かを判定する。断面を流れる物理量はそれぞれの断面で同じであるという保存則が成立する。物理量が電流であれば、その保存則は電流保存則である。各断面の移動物理量が収束するとは、個々の断面の間で移動物理量が一致するか、または所与の誤差範囲内に収まることをいう。移動物理量が収束しない場合には処理はステップS144に戻り(ステップS146においてNO)、特性算出部13は各サイトでのポテンシャル量と各断面での移動物理量とを再計算する。 In step S146, the characteristic calculation unit 13 determines whether or not the movement physical quantity of each cross section has converged. Conservation law holds that the physical quantity flowing through a cross section is the same for each cross section. If the physical quantity is electric current, its conservation law is the electric current conservation law. Convergence of the moving physical amounts of the respective cross sections means that the moving physical amounts of the individual cross sections match or fall within a given error range. If the movement physical quantity does not converge, the process returns to step S144 (NO in step S146), and the characteristic calculator 13 recalculates the potential quantity at each site and the movement physical quantity at each cross section.

移動物理量が収束した場合には処理はステップS147に進む(ステップS146においてYES)。ステップS147では、特性算出部13が物理量の伝導度を計算する。この伝導度はコネクティビティともいわれる。伝導度は、開サイトが一つも含まれない系における断面の移動物理量に対する、対象とする系における断面の移動物理量の比で表される。特性算出部13は各サイトでのポテンシャル量に基づいて、物理量を出力するサイトを出力サイトとして識別する。そして、特性算出部13は、開サイトが一つも含まれない系における出力端の断面の移動物理量に対する、対象とする系における出力端の断面の移動物理量の比を伝導度として得る。特性算出部13は、算出された比を所与のフィッティング式に代入することで、材料の立体形状に対応する伝導度を算出してもよい。 If the movement physical quantity converges, the process proceeds to step S147 (YES in step S146). In step S147, the characteristic calculator 13 calculates the conductivity of the physical quantity. This conductivity is also called connectivity. The conductivity is expressed by the ratio of the physical quantity of cross section movement in a system of interest to the physical quantity of cross section movement in a system containing no open sites. Based on the potential amount at each site, the characteristic calculator 13 identifies the site that outputs the physical quantity as the output site. Then, the characteristic calculator 13 obtains, as the conductivity, the ratio of the physical quantity of movement of the cross section of the output end in the target system to the physical quantity of movement of the cross section of the output end in the system containing no open sites. The characteristic calculator 13 may calculate the conductivity corresponding to the three-dimensional shape of the material by substituting the calculated ratio into a given fitting formula.

特性算出部13は各サイトでのポテンシャル量に基づいて様々な輸送特性を算出できる。例えば、特性算出部13は微細組織内の移動物理量の分布と、伝導度とのうちの少なくとも一つを輸送特性として求めることができる。移動物理量が電気であれば、特性算出部13は、各サイトでの電位に基づいて、電位分布、電流分布、および電気伝導度のうちの少なくとも一つを輸送特性(電気輸送特性)として算出してもよい。 The characteristic calculator 13 can calculate various transport characteristics based on the potential amount at each site. For example, the property calculation unit 13 can obtain at least one of the distribution of the moving physical quantity in the microstructure and the conductivity as the transport property. If the moving physical quantity is electricity, the property calculation unit 13 calculates at least one of the potential distribution, the current distribution, and the electrical conductivity as the transport property (electric transport property) based on the potential at each site. may

図3に戻って、ステップS15では、特性算出部13が輸送特性データを出力する。輸送特性データは輸送特性を示す電子データである。輸送特性データの出力方法は限定されない。例えば、特性算出部13は輸送特性データをモニタ上に出力してもよいし、メモリ、データベースなどの所与の記憶装置に輸送特性データを格納してもよいし、他のコンピュータに向けて輸送特性データを送信してもよい。輸送特性データの表現方法も限定されない。例えば、特性算出部13は輸送特性データを数字、記号、文字列、テキスト、画像(コンピュータグラフィックス)、またはこれらのうちの任意の2以上の手法の組合せによって表現してもよい。あるいは、特性算出部13は、コンピュータグラフィックスによって生成した輸送特性データの画像(輸送特性画像)を材料画像に重畳することで該輸送特性データを表現してもよい。輸送特性画像を材料画像に重畳することで、微細組織における局所的または大局的な品質を容易に把握することができる。例えば、超伝導材料において電流が集中している箇所を特定することができる。 Returning to FIG. 3, in step S15, the property calculator 13 outputs the transport property data. Transport property data is electronic data indicative of transport properties. The method of outputting the transportation characteristic data is not limited. For example, the characteristic calculation unit 13 may output the transport characteristic data on a monitor, store the transport characteristic data in a given storage device such as a memory or database, or transmit the transport characteristic data to another computer. Characteristic data may be transmitted. The method of expressing the transport characteristic data is also not limited. For example, the property calculator 13 may express the transport property data by numbers, symbols, character strings, texts, images (computer graphics), or a combination of any two or more of these methods. Alternatively, the property calculator 13 may express the transport property data by superimposing an image of the transport property data (transport property image) generated by computer graphics on the material image. By superimposing the transport property image on the material image, local or global quality in the microstructure can be easily grasped. For example, it is possible to identify locations in a superconducting material where current is concentrated.

図6は可視化された輸送特性データの一例を示す図である。この例では、材料は、超伝導多結晶材料の一例である二ホウ化マグネシウム(MgB)焼結体である。微細組織の左端から右端に向かって電流が流れると仮定して各サイトでの電流値を算出し、その計算結果を電流分布として可視化することで、図6のような画像が得られる。グレースケールで示す図6では分かりづらいが、カラー画像では個々のサイトでの電流量が多階調の色で示される。したがって、電流が集中する部分などの局所的な傾向と、微細組織の大局的な傾向とを容易に確認することができる。図7は、図6に示す領域30を拡大して示す図である。この領域における部分41は、電気が流れない領域、すなわち不純物相および空隙相である。部分42および部分43はいずれも電気が流れる領域であり、すなわち超伝導相(観察対象の相)である。部分42,43の色の違いは電流値の大きさを示し、この例では、部分43での電流値の方が部分42での電流値よりも高い。個々のサイト、すなわち個々の画素について輸送特性が得られるので、図6および図7に示すように輸送特性の局所的な分布を確認できる。例えば、どこで電流が集中しているかを容易に特定することが可能になる。 FIG. 6 is a diagram showing an example of visualized transport property data. In this example, the material is magnesium diboride (MgB 2 ) sintered body, which is an example of a superconducting polycrystalline material. Assuming that the current flows from the left end to the right end of the microstructure, the current value at each site is calculated, and by visualizing the calculation results as a current distribution, an image such as that shown in FIG. 6 is obtained. Although it is difficult to understand in FIG. 6, which is shown in grayscale, the current amount at each site is shown in multi-gradation colors in the color image. Therefore, it is possible to easily confirm a local tendency such as a portion where current concentrates and a global tendency of the microstructure. FIG. 7 is an enlarged view of the area 30 shown in FIG. Portions 41 in this region are regions where electricity does not flow, ie impurity and void phases. Both portion 42 and portion 43 are regions in which electricity flows, ie, the superconducting phase (the phase to be observed). The difference in color between portions 42 and 43 indicates the magnitude of the current value, with the current value at portion 43 being higher than the current value at portion 42 in this example. Since the transport properties are obtained for individual sites, ie individual pixels, the local distribution of the transport properties can be ascertained as shown in FIGS. For example, it becomes possible to easily identify where the current is concentrated.

[プログラム]
図8を参照しながら、材料解析システム10を実現するための材料解析プログラムP1を説明する。図8は実施形態に係る材料解析プログラムP1の構成の一例を示す図である。
[program]
A material analysis program P1 for realizing the material analysis system 10 will be described with reference to FIG. FIG. 8 is a diagram showing an example of the configuration of the material analysis program P1 according to the embodiment.

材料解析プログラムP1は、メインモジュールP10、画像解析モジュールP11、モデル生成モジュールP12、および特性算出モジュールP13を備える。メインモジュールP10は材料解析を統括的に制御する部分である。画像解析モジュールP11、モデル生成モジュールP12、および特性算出モジュールP13を実行することにより画像解析部11、モデル生成部12、および特性算出部13が実現する。 The material analysis program P1 includes a main module P10, an image analysis module P11, a model generation module P12, and a property calculation module P13. The main module P10 is a part that controls material analysis in an integrated manner. The image analysis unit 11, the model generation unit 12, and the characteristic calculation unit 13 are realized by executing the image analysis module P11, the model generation module P12, and the characteristic calculation module P13.

材料解析プログラムP1は、CD-ROM、DVD-ROM、半導体メモリなどの非一時的な記録媒体に固定的に記録された上で提供されてもよい。あるいは、材料解析プログラムP1は、搬送波に重畳されたデータ信号として通信ネットワークを介して提供されてもよい。 The material analysis program P1 may be provided after being fixedly recorded in a non-temporary recording medium such as a CD-ROM, a DVD-ROM, or a semiconductor memory. Alternatively, the material analysis program P1 may be provided via a communication network as a data signal superimposed on a carrier wave.

[効果]
以上説明したように、本開示の一側面に係る材料解析システムは、少なくとも一つのプロセッサを備える。少なくとも一つのプロセッサは、材料の微細組織を写した材料画像を取得し、材料画像を解析することで観察対象の相の位置を識別し、相の位置に基づいて材料の微細組織に対応する輸送モデルを生成し、輸送モデルに基づいて輸送特性を算出し、輸送特性を示す輸送特性データを出力する。
[effect]
As described above, the material analysis system according to one aspect of the present disclosure includes at least one processor. At least one processor acquires a material image depicting a microstructure of the material, analyzes the material image to identify the position of the observed phase, and determines transport corresponding to the material microstructure based on the phase position. A model is generated, transport characteristics are calculated based on the transport model, and transport characteristics data representing the transport characteristics are output.

本発明の一側面に係る材料解析方法は、少なくとも一つのプロセッサを備えるコンピュータシステムにより実行される。材料解析方法は、材料の微細組織を写した材料画像を取得するステップと、材料画像を解析することで観察対象の相の位置を識別するステップと、相の位置に基づいて材料の微細組織に対応する輸送モデルを生成するステップと、輸送モデルに基づいて輸送特性を算出するステップと、輸送特性を示す輸送特性データを出力するステップとを含む。 A material analysis method according to one aspect of the present invention is performed by a computer system comprising at least one processor. The material analysis method includes the steps of acquiring a material image showing the microstructure of the material, identifying the position of the phase to be observed by analyzing the material image, and identifying the microstructure of the material based on the position of the phase. generating a corresponding transport model; calculating transport characteristics based on the transport model; and outputting transport characteristic data indicative of the transport characteristics.

本発明の一側面に係る材料解析プログラムは、材料の微細組織を写した材料画像を取得するステップと、材料画像を解析することで観察対象の相の位置を識別するステップと、相の位置に基づいて材料の微細組織に対応する輸送モデルを生成するステップと、輸送モデルに基づいて輸送特性を算出するステップと、輸送特性を示す輸送特性データを出力するステップとをコンピュータシステムに実行させる。 A material analysis program according to one aspect of the present invention includes a step of acquiring a material image showing a microstructure of a material, a step of identifying the position of a phase to be observed by analyzing the material image, and a phase position. generate a transport model corresponding to the microstructure of the material based on the transport model; calculate transport properties based on the transport model; and output transport property data representing the transport properties.

このような側面においては、材料の微細組織を写した材料画像から、該微細組織に対応する輸送モデルが生成され、その輸送モデルに基づいて輸送特性が算出される。したがって、材料画像さえ準備すれば材料の輸送特性を簡易に得ることができる。材料の輸送特性を得るための実測作業を要することなく、材料画像に基づいて詳細な計算が可能になるので、材料解析のさらなる高度化および高速化が期待できる。輸送特性データは、材料の開発、材料の品質の判定、材料の欠陥の検出などの様々な目的のために用いることができる。 In this aspect, a transport model corresponding to the microstructure is generated from a material image representing the microstructure of the material, and the transport characteristics are calculated based on the transport model. Therefore, the transportation characteristics of the material can be easily obtained by preparing the material image. Since detailed calculations can be performed based on material images without the need for actual measurement work to obtain the transport properties of materials, further sophistication and speeding up of material analysis can be expected. Transport property data can be used for a variety of purposes, such as material development, material quality determination, and material defect detection.

他の側面に係る材料解析システムでは、輸送モデルが、微細組織に対応する物理量に関するモデルであってもよい。物理量に関する輸送モデルを生成することで、材料の輸送特性を詳細に算出することができる。 In a material analysis system according to another aspect, the transport model may be a model relating to physical quantities corresponding to microstructures. By generating a transport model for physical quantities, the transport properties of materials can be calculated in detail.

他の側面に係る材料解析システムでは、少なくとも一つのプロセッサが、相の位置に基づいて、材料画像の複数の画素のそれぞれに、物理量が流れる領域である閉サイト、または物理量が流れない領域である開サイトを設定し、1以上の閉サイトおよび1以上の開サイトのそれぞれの位置を示す輸送モデルを生成してもよい。個々の画素に対応して閉サイトおよび開サイトを設定することで精密な輸送モデルが生成されるので、詳細な輸送特性を算出することができる。 In a material analysis system according to another aspect, at least one processor is a closed site, which is a region through which a physical quantity flows, or a region in which a physical quantity does not flow, in each of a plurality of pixels of a material image based on the position of a phase. An open site may be set and a transportation model may be generated indicating the position of each of the one or more closed sites and the one or more open sites. By setting closed sites and open sites corresponding to individual pixels, a precise transport model is generated, so detailed transport properties can be calculated.

他の側面に係る材料解析システムでは、少なくとも一つのプロセッサが、輸送特性データの画像を材料画像に重畳することで該輸送特性データを表現してもよい。このような表現によって、微細組織における局所的または大局的な品質を容易に把握することができる。 In the material analysis system according to another aspect, at least one processor may represent the transport property data by superimposing an image of the transport property data on the material image. Such a representation makes it easy to grasp local or global quality in microstructures.

他の側面に係る材料解析システムでは、材料が超伝導材料であり、輸送特性が電気輸送特性であってもよい。この場合には、超伝導材料の電気輸送特性を簡易に得ることができる。 In a material analysis system according to another aspect, the material may be a superconducting material and the transport properties may be electrotransport properties. In this case, the electrotransport properties of superconducting materials can be easily obtained.

以下、実施例を具体的に説明するが、本開示はそれらに何ら限定されるものではない。 Examples will be specifically described below, but the present disclosure is not limited to them.

下記の参考文献で紹介されている、二ホウ化マグネシウム(MgB)焼結体の被研磨断面を走査電子顕微鏡(SEM)により観察し、その顕微鏡画像を材料画像として取得した。本実施例では、参考資料においてBulk-Sとして示されている材料をサンプルAとし、Bulk-Dとして示されている材料をサンプルB1およびB2として扱った。 A polished cross-section of a magnesium diboride (MgB 2 ) sintered body introduced in the following references was observed with a scanning electron microscope (SEM), and its microscopic image was acquired as a material image. In this example, the material indicated as Bulk-S in the reference material was treated as sample A, and the material indicated as Bulk-D was treated as samples B1 and B2.

(参考文献)
S. Mizutani et al., Understanding routesfor high connectivity in ex situ MgB2 by self-sintering, SuperconductorScience and Technology, Volume 27, Issue 4, 044012 (2014).
(References)
S. Mizutani et al., Understanding routes for high connectivity in ex situ MgB 2 by self-sintering, SuperconductorScience and Technology, Volume 27, Issue 4, 044012 (2014).

材料画像を解析することで、観察対象である二ホウ化マグネシウムの相(超伝導相)を分析した。そして、材料画像を超伝導相とそれ以外の相(不純物相および空隙相)とに二値化する処理を実行した。超伝導相の位置に対応させて閉サイトを設定し、不純物相および空隙相に対応させて開サイトを設定することで輸送モデルを生成した。 By analyzing the material image, the phase (superconducting phase) of magnesium diboride, which is the object of observation, was analyzed. Then, the material image was binarized into the superconducting phase and other phases (impurity phase and void phase). A transport model was generated by setting closed sites corresponding to the positions of the superconducting phase and setting open sites corresponding to the impurity and void phases.

輸送特性を計算するために、FORTRAN90というプログラム言語を用いて、特性算出部13に対応するプログラムを作成した。そして、そのプログラムを用いて個々のサイトでの電流と、材料の有効電気伝導度とを求め、計算結果をCSVファイル形式で出力した。 In order to calculate transport characteristics, a program corresponding to the characteristics calculator 13 was created using a program language called FORTRAN90. Then, using the program, the current at each site and the effective electrical conductivity of the material were obtained, and the calculation results were output in CSV file format.

図9はサンプルAの顕微鏡画像、可視化した輸送モデル、および可視化した輸送特性データを示す。図9と同様に、サンプルB1について図10に示し、サンプルB2について図11に示す。輸送モデルの画像は、閉サイト(超伝導相)を黒で示し、開サイト(不純物相および空隙相)を白で示す。輸送特性データの画像は電流分布を示し、実際にはカラー画像である。図9~図11では、電流の経路は、概して、明るいグレーで示される。これらの画像から、材料の微細組織内の絶縁成分の影響によってその経路が蛇行していることが分かる。実際のカラー画像では個々のサイトでの電流量が多階調の色で示され、図9~図11では、電流量が多い部分が白で表され、電流量が無いかまたは少ない部分が黒で表されている。したがって、電流が集中する部分などの局所的な傾向と、微細組織の大局的な傾向とを画像から容易に確認することができる。 FIG. 9 shows a microscopic image of sample A, a visualized transport model, and visualized transport property data. Similar to FIG. 9, the sample B1 is shown in FIG. 10 and the sample B2 is shown in FIG. Images of the transport model show closed sites (superconducting phase) in black and open sites (impurity and void phases) in white. The image of the transport property data shows the current distribution and is actually a color image. In FIGS. 9-11 the current paths are generally shown in light gray. From these images, it can be seen that the path is tortuous due to the influence of insulating components within the microstructure of the material. In an actual color image, the amount of current at each site is shown in multi-gradation colors. In FIGS. 9 to 11, portions with a large amount of current are represented in white, and portions with no or little amount of current are represented in black. is represented by Therefore, it is possible to easily confirm the local tendency such as the part where the current concentrates and the global tendency of the microstructure from the image.

サンプルA,B1,B2のそれぞれについての、有効電気伝導度の実測値と、本開示の手法による推定値とを表1に示す。推定値(2D)は画像解析から得られた結果そのものであり、画像という2次元(2D)空間における推定値である。推定値(3D)は、下記のフィッティング式(1)にその推定値(2D)を代入することで得られた値であり、材料の立体形状に対応する推定値である。この式(1)は、コンピュータで開サイトと閉サイトとをランダムに分布させて仮想的に発生させた2次元の系における有効電気伝導度の計算結果と、同様の手法で発生させた3次元の系における有効電気伝導度の計算結果とを最小二乗法で近似することで得られた関数である。式(1)において、K3Dは推定値(3D)であり、K2Dは推定値(2D)である。表1では伝導度を百分率で示す。
3D=0.7729K2D+0.2271 …(1)
Table 1 shows the measured effective electrical conductivity values and estimated values obtained by the method of the present disclosure for each of samples A, B1, and B2. The estimated value (2D) is the result obtained from the image analysis itself, and is the estimated value in the two-dimensional (2D) space of the image. The estimated value (3D) is a value obtained by substituting the estimated value (2D) into the following fitting equation (1), and is an estimated value corresponding to the three-dimensional shape of the material. This formula (1) is a calculation result of the effective electrical conductivity in a two-dimensional system virtually generated by randomly distributing open sites and closed sites on a computer, and a three-dimensional system generated by a similar method. It is a function obtained by approximating the calculation result of the effective electrical conductivity in the system by the least squares method. In equation (1), K 3D is the estimate (3D) and K 2D is the estimate (2D). Table 1 shows the conductivity as a percentage.
K3D = 0.7729K2D +0.2271 (1)

Figure 0007299615000001
Figure 0007299615000001

図9~図11および表1に示すように、微細組織の局所的な傾向と、大局的な傾向との双方を特定することができた。 As shown in FIGS. 9-11 and Table 1, both local and global trends in microstructure could be identified.

[変形例]
以上、本発明をその実施形態および実施例に基づいて詳細に説明した。しかし、本発明は上記の実施形態および実施例に限定されるものではない。本発明は、その要旨を逸脱しない範囲で様々な変形が可能である。
[Modification]
The present invention has been described in detail above based on its embodiments and examples. However, the invention is not limited to the above embodiments and examples. Various modifications are possible for the present invention without departing from the gist thereof.

本開示において、「少なくとも一つのプロセッサが、第1の処理を実行し、第2の処理を実行し、…第nの処理を実行する。」との表現、またはこれに対応する表現は、第1の処理から第nの処理までのn個の処理の実行主体(すなわちプロセッサ)が途中で変わる場合を含む概念を示す。すなわち、この表現は、n個の処理のすべてが同じプロセッサで実行される場合と、n個の処理においてプロセッサが任意の方針で変わる場合との双方を含む概念を示す。 In the present disclosure, the expression “at least one processor executes the first process, the second process, . . . The concept is shown including the case where the executing subject (that is, the processor) of n processes from process 1 to process n changes in the middle. That is, this expression shows a concept including both the case where all of the n processes are executed by the same processor and the case where the processors are changed according to an arbitrary policy in the n processes.

少なくとも一つのプロセッサにより実行される方法の処理手順は上記実施形態での例に限定されない。例えば、上述したステップ(処理)の一部が省略されてもよいし、別の順序で各ステップが実行されてもよい。また、上述したステップのうちの任意の2以上のステップが組み合わされてもよいし、ステップの一部が修正又は削除されてもよい。あるいは、上記の各ステップに加えて他のステップが実行されてもよい。 The processing steps of the method executed by at least one processor are not limited to the examples in the above embodiments. For example, some of the steps (processes) described above may be omitted, or the steps may be performed in a different order. Also, any two or more of the steps described above may be combined, and some of the steps may be modified or deleted. Alternatively, other steps may be performed in addition to the above steps.

コンピュータシステムまたはコンピュータ内で二つの数値の大小関係を比較する際には、「以上」および「よりも大きい」という二つの基準のどちらを用いてもよく、「以下」および「未満」という二つの基準のうちのどちらを用いてもよい。このような基準の選択は、二つの数値の大小関係を比較する処理についての技術的意義を変更するものではない。 When comparing two numerical values within a computer system or within a computer, either of the two criteria "greater than" and "greater than" may be used, and the two criteria "less than" and "less than" may be used. Either of the criteria may be used. Selection of such a criterion does not change the technical significance of the process of comparing two numerical values.

10…材料解析システム、11…画像解析部、12…モデル生成部、13…特性算出部、P1…材料解析プログラム、P10…メインモジュール、P11…画像解析モジュール、P12…モデル生成モジュール、P13…特性算出モジュール。 DESCRIPTION OF SYMBOLS 10... Material analysis system 11... Image analysis part 12... Model generation part 13... Property calculation part P1... Material analysis program P10... Main module P11... Image analysis module P12... Model generation module P13... Property calculation module.

Claims (7)

少なくとも一つのプロセッサを備え、
前記少なくとも一つのプロセッサが、
材料の微細組織を写した材料画像を取得し、
前記材料画像を解析することで観察対象の相の位置を識別し、
前記相の位置に基づいて前記材料の微細組織に対応する輸送モデルを生成し、
前記輸送モデルに基づいて、前記微細組織に対応する移動物理量の分布を含む輸送特性を算出し、
前記輸送特性を示す輸送特性データを出力する、
材料解析システム。
comprising at least one processor,
the at least one processor
Acquire a material image that shows the microstructure of the material,
identifying the position of the phase to be observed by analyzing the material image;
generating a transport model corresponding to the microstructure of the material based on the phase positions;
Based on the transport model , calculate the transport characteristics including the distribution of the moving physical quantity corresponding to the microstructure ,
outputting transport property data indicative of said transport properties;
Material analysis system.
前記輸送モデルが、前記微細組織に対応する物理量に関するモデルである、
請求項1に記載の材料解析システム。
wherein the transport model is a model relating to physical quantities corresponding to the microstructure;
The material analysis system according to claim 1.
前記少なくとも一つのプロセッサが、
前記相の位置に基づいて、前記材料画像の複数の画素のそれぞれに、物理量が流れる領域である閉サイト、または前記物理量が流れない領域である開サイトを設定し、
1以上の前記閉サイトおよび1以上の前記開サイトのそれぞれの位置を示す前記輸送モデルを生成する、
請求項1または2に記載の材料解析システム。
the at least one processor
Based on the position of the phase, in each of the plurality of pixels of the material image, set a closed site that is a region where the physical quantity flows, or an open site that is a region where the physical quantity does not flow;
generating the transport model indicating the position of each of the one or more closed sites and the one or more open sites;
The material analysis system according to claim 1 or 2.
前記少なくとも一つのプロセッサが、前記輸送特性データの画像を前記材料画像に重畳することで該輸送特性データを表現する、
請求項1~3のいずれか一項に記載の材料解析システム。
wherein the at least one processor represents the transport property data by overlaying an image of the transport property data on the material image;
The material analysis system according to any one of claims 1-3.
前記材料が超伝導材料であり、前記輸送特性が電気輸送特性である、
請求項1~4のいずれか一項に記載の材料解析システム。
said material is a superconducting material and said transport property is an electrotransport property;
The material analysis system according to any one of claims 1-4.
少なくとも一つのプロセッサを備えるコンピュータシステムにより実行される材料解析方法であって、
材料の微細組織を写した材料画像を取得するステップと、
前記材料画像を解析することで観察対象の相の位置を識別するステップと、
前記相の位置に基づいて前記材料の微細組織に対応する輸送モデルを生成するステップと、
前記輸送モデルに基づいて、前記微細組織に対応する移動物理量の分布を含む輸送特性を算出するステップと、
前記輸送特性を示す輸送特性データを出力するステップと
を含む材料解析方法。
A material analysis method performed by a computer system comprising at least one processor, comprising:
obtaining a material image that captures the microstructure of the material;
identifying the position of the observed phase by analyzing the material image;
generating a transport model corresponding to the microstructure of the material based on the phase positions;
a step of calculating a transport characteristic including a distribution of movement physical quantities corresponding to the microstructure based on the transport model;
and outputting transport property data indicating the transport property.
材料の微細組織を写した材料画像を取得するステップと、
前記材料画像を解析することで観察対象の相の位置を識別するステップと、
前記相の位置に基づいて前記材料の微細組織に対応する輸送モデルを生成するステップと、
前記輸送モデルに基づいて、前記微細組織に対応する移動物理量の分布を含む輸送特性を算出するステップと、
前記輸送特性を示す輸送特性データを出力するステップと
をコンピュータシステムに実行させる材料解析プログラム。
obtaining a material image that captures the microstructure of the material;
identifying the position of the observed phase by analyzing the material image;
generating a transport model corresponding to the microstructure of the material based on the phase positions;
a step of calculating a transport characteristic including a distribution of movement physical quantities corresponding to the microstructure based on the transport model;
A material analysis program for causing a computer system to execute a step of outputting transport property data indicating the transport property.
JP2019161138A 2019-09-04 2019-09-04 Material analysis system, material analysis method, and material analysis program Active JP7299615B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019161138A JP7299615B2 (en) 2019-09-04 2019-09-04 Material analysis system, material analysis method, and material analysis program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019161138A JP7299615B2 (en) 2019-09-04 2019-09-04 Material analysis system, material analysis method, and material analysis program

Publications (2)

Publication Number Publication Date
JP2021039021A JP2021039021A (en) 2021-03-11
JP7299615B2 true JP7299615B2 (en) 2023-06-28

Family

ID=74849197

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019161138A Active JP7299615B2 (en) 2019-09-04 2019-09-04 Material analysis system, material analysis method, and material analysis program

Country Status (1)

Country Link
JP (1) JP7299615B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023283781A1 (en) * 2021-07-12 2023-01-19 西门子股份公司 Method for generating microstructure data of material

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008058323A (en) 2002-11-22 2008-03-13 Showa Denko Kk Method of evaluating metal structure of r-t-b magnet alloy
JP2014081362A (en) 2012-09-28 2014-05-08 Sumika Chemical Analysis Service Ltd Evaluation method of electrode and evaluation device
JP2018190631A (en) 2017-05-09 2018-11-29 株式会社日立製作所 MgB2 SUPERCONDUCTING THIN FILM WIRE AND MANUFACTURING METHOD THEREFOR

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02153803A (en) * 1988-06-06 1990-06-13 Nippon Steel Corp Oxide superconductor bulk material and production thereof

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008058323A (en) 2002-11-22 2008-03-13 Showa Denko Kk Method of evaluating metal structure of r-t-b magnet alloy
JP2014081362A (en) 2012-09-28 2014-05-08 Sumika Chemical Analysis Service Ltd Evaluation method of electrode and evaluation device
JP2018190631A (en) 2017-05-09 2018-11-29 株式会社日立製作所 MgB2 SUPERCONDUCTING THIN FILM WIRE AND MANUFACTURING METHOD THEREFOR

Also Published As

Publication number Publication date
JP2021039021A (en) 2021-03-11

Similar Documents

Publication Publication Date Title
Spowart Automated serial sectioning for 3-D analysis of microstructures
Breen et al. Correlating atom probe crystallographic measurements with transmission Kikuchi diffraction data
Wargo et al. Resolving macro-and micro-porous layer interaction in polymer electrolyte fuel cells using focused ion beam and X-ray computed tomography
Strassburg et al. Reconstruction of electrostatic force microscopy images
da Rocha et al. Ultimate conductivity performance in metallic nanowire networks
Reed et al. One‐stop stereology: the estimation of 3D parameters using isotropic rulers
Lefebvre et al. Application of Delaunay tessellation for the characterization of solute-rich clusters in atom probe tomography
Larson et al. On the use of simulated field-evaporated specimen apex shapes in atom probe tomography data reconstruction
TWI463161B (en) Distributed analysis device
Sperry et al. Comparison of EBSD, DIC, AFM, and ECCI for active slip system identification in deformed Ti-7Al
Lucocq et al. Systems biology in 3D space–enter the morphome
Morelly et al. Three-dimensional visualization of conductive domains in battery electrodes with contrast-enhancing nanoparticles
Rolland et al. New atom probe tomography reconstruction algorithm for multilayered samples: Beyond the hemispherical constraint
JP7299615B2 (en) Material analysis system, material analysis method, and material analysis program
Kozikowski Extracting three-dimensional information from SEM images by means of photogrammetry
Wijaya et al. Multi-method characterization approach to facilitate a strategy to design mechanical and electrical properties of sintered copper
Korneev et al. A data-driven multiscale framework to estimate effective properties of lithium-ion batteries from microstructure images
Xiong et al. Multiscale modeling of interface-mediated mechanical, thermal, and mass transport in heterogeneous materials: Perspectives and applications
Bagherian et al. A new statistical descriptor for the physical characterization and 3D reconstruction of heterogeneous materials
Ghadban et al. 3D microscale modeling of NMC cathodes using multi-resolution FIB-SEM tomography
Holm et al. Three-dimensional materials science
Cresti et al. Theoretical imaging of current profiles in two-dimensional devices
Takahashi et al. Approach for three-dimensional observation of mesoscopic precipitates in alloys by coherent x-ray diffraction microscopy
JP2018109577A (en) Analyzer
Mouton et al. Calibration of atom probe tomography reconstructions through correlation with electron micrographs

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220629

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230331

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230404

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20230525

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: 20230606

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230609

R150 Certificate of patent or registration of utility model

Ref document number: 7299615

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150