Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7310558B2 - Prediction device, prediction program, generation device and generation program - Google Patents
[go: Go Back, main page]

JP7310558B2 - Prediction device, prediction program, generation device and generation program - Google Patents

Prediction device, prediction program, generation device and generation program Download PDF

Info

Publication number
JP7310558B2
JP7310558B2 JP2019204724A JP2019204724A JP7310558B2 JP 7310558 B2 JP7310558 B2 JP 7310558B2 JP 2019204724 A JP2019204724 A JP 2019204724A JP 2019204724 A JP2019204724 A JP 2019204724A JP 7310558 B2 JP7310558 B2 JP 7310558B2
Authority
JP
Japan
Prior art keywords
substrate
test
feature amount
value
amount
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
JP2019204724A
Other languages
Japanese (ja)
Other versions
JP2021076521A (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.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
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 Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2019204724A priority Critical patent/JP7310558B2/en
Publication of JP2021076521A publication Critical patent/JP2021076521A/en
Application granted granted Critical
Publication of JP7310558B2 publication Critical patent/JP7310558B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)

Description

本発明は、予測装置、予測プログラム、生成装置及び生成プログラムに関する。 The present invention relates to a prediction device, a prediction program, a generation device and a generation program.

電子機器の高度化に伴い、機器内に搭載される半導体に対しても高い品質が求められている。近年、多品種少量生産によって様々な種類の半導体が作られている中で、品質を保証する信頼性試験はますます重要となっている。一方で、半導体の集積化が進展すると、信頼性試験の工数も増加する。 With the advancement of electronic equipment, high quality is required for the semiconductors mounted in the equipment. In recent years, various types of semiconductors have been produced by high-mix low-volume production, and reliability tests for quality assurance have become more and more important. On the other hand, as the integration of semiconductors progresses, the man-hours for reliability tests also increase.

従来、半導体の信頼性試験を効率化するための方法として、試験の実績データから半導体の基板の品質を予測する方法が提案されている。例えば、従来の方法では、はんだ接合部の寿命サイクル数等の品質に関わる情報と、接合部断面の硬さとの関係式を実績データから求めることで、基板の品質の予測が行われている。 Conventionally, a method of predicting the quality of a semiconductor substrate from actual test data has been proposed as a method for improving the efficiency of a semiconductor reliability test. For example, in the conventional method, the quality of the board is predicted by obtaining the relational expression between the information related to the quality such as the number of life cycles of the solder joint and the hardness of the cross section of the joint from actual data.

特開2013-32931号公報JP 2013-32931 A 特開2002-310888号公報Japanese Patent Application Laid-Open No. 2002-310888 特開平3-128431号公報JP-A-3-128431

しかしながら、上記の方法では、基板の品質の予測精度が低下する場合があるという問題がある。例えば、従来の方法では、サイズや材料等の設計値の変更で基板の特性が大きく変化するような場合に、変更後の品質の予測精度が低下してしまう恐れがある。 However, the above method has a problem that the accuracy of predicting the quality of the substrate may be lowered. For example, in the conventional method, when the design values such as the size and material are changed, and the characteristics of the substrate are greatly changed, there is a possibility that the post-change quality prediction accuracy is lowered.

1つの側面では、基板の品質の予測精度の低下を抑止することを目的とする。 An object of one aspect is to suppress deterioration in accuracy of predicting the quality of a substrate.

1つの態様において、予測装置は、抽出部、生成部及び算出部を有する。抽出部は、基板を変形させる試験における基板の変形量の変化を、周波数領域で表現したときの、試験の複数の進行度合いのそれぞれに対応する所定の周波数の範囲における振幅の積算値の合計値を、基板の特徴量として抽出する。生成部は、特徴量と、基板の所定の設計値との関係を表すモデルを生成する。算出部は、生成部によって生成されたモデルを用いて、基板の設計値から特徴量の予測値を算出する。 In one aspect, a predictor comprises an extractor, a generator and a calculator. The extraction unit obtains the total value of integrated values of amplitudes in a predetermined frequency range corresponding to each of a plurality of degrees of progress of the test, when changes in the amount of deformation of the board in the test of deforming the board are expressed in the frequency domain. is extracted as the feature quantity of the substrate. The generator generates a model representing the relationship between the feature amount and the predetermined design value of the substrate. The calculator uses the model generated by the generator to calculate the predicted value of the feature amount from the design value of the board.

1つの側面では、基板の品質の予測精度の低下を抑止することができる。 In one aspect, it is possible to prevent deterioration in accuracy of predicting the quality of the board.

図1は、実施例1に係る予測装置の構成例を示す図である。FIG. 1 is a diagram illustrating a configuration example of a prediction device according to a first embodiment; 図2は、変形量を説明するための図である。FIG. 2 is a diagram for explaining the amount of deformation. 図3は、設計情報の一例を示す図である。FIG. 3 is a diagram showing an example of design information. 図4は、ヒステリシス曲線を説明するための図である。FIG. 4 is a diagram for explaining a hysteresis curve. 図5は、変形量の傾きを説明するための図である。FIG. 5 is a diagram for explaining the slope of the deformation amount. 図6は、特徴量の抽出方法を説明するための図である。FIG. 6 is a diagram for explaining a feature amount extraction method. 図7は、モデルの生成方法を説明するための図である。FIG. 7 is a diagram for explaining a model generation method. 図8は、予測値の算出方法を説明するための図である。FIG. 8 is a diagram for explaining a method of calculating predicted values. 図9は、予測値の例を示す図である。FIG. 9 is a diagram showing an example of predicted values. 図10は、予測処理の流れを示すフローチャートである。FIG. 10 is a flowchart showing the flow of prediction processing. 図11は、特徴量を抽出する処理の流れを示すフローチャートである。FIG. 11 is a flow chart showing the flow of processing for extracting feature amounts. 図12は、ハードウェア構成例を説明する図である。FIG. 12 is a diagram illustrating a hardware configuration example.

以下に、本発明に係る予測装置、予測プログラム、生成装置及び生成プログラムの実施例を図面に基づいて詳細に説明する。なお、この実施例により本発明が限定されるものではない。また、各実施例は、矛盾のない範囲内で適宜組み合わせることができる。 Exemplary embodiments of a prediction device, a prediction program, a generation device, and a generation program according to the present invention will be described below in detail with reference to the drawings. It should be noted that the present invention is not limited by these examples. Moreover, each embodiment can be appropriately combined within a range without contradiction.

[機能構成]
図1を用いて、実施例に係る予測装置の構成を説明する。図1は、実施例1に係る予測装置の構成例を示す図である。図1に示すように、予測装置10は、インタフェース部11、記憶部12及び制御部13を有する。
[Function configuration]
The configuration of the prediction device according to the embodiment will be described with reference to FIG. FIG. 1 is a diagram illustrating a configuration example of a prediction device according to a first embodiment; As shown in FIG. 1 , the prediction device 10 has an interface section 11 , a storage section 12 and a control section 13 .

インタフェース部11は、データの入出力、及び他の装置との間でのデータの通信を行うためのインタフェースである。また、例えば、インタフェース部11は、NIC(Network Interface Card)であり、インターネットを介してデータの通信を行ってもよい。 The interface unit 11 is an interface for data input/output and data communication with other devices. Further, for example, the interface unit 11 is a NIC (Network Interface Card), and may perform data communication via the Internet.

記憶部12は、データや制御部13が実行するプログラム等を記憶する記憶装置の一例であり、例えばハードディスクやメモリ等である。記憶部12は、実績情報121、設計情報122及びモデル情報123を記憶する。 The storage unit 12 is an example of a storage device that stores data, programs executed by the control unit 13, and the like, such as a hard disk or memory. The storage unit 12 stores performance information 121 , design information 122 and model information 123 .

実績情報121は、半導体の基板の温度サイクル試験(T/C試験)の結果である。T/C試験は、試験対象の基板の温度を繰り返し変化させる試験である。実績情報121には、後述するヒステリシス曲線を作成するために必要な情報が含まれる。例えば、実績情報121には、試験の繰り返し回数及び温度ごとの変形量等が含まれる。また、図2に示すように、各温度における変形量は、基板の所定の位置と基準となる位置との距離で表される。図2は、変形量を説明するための図である。 The performance information 121 is the result of a temperature cycle test (T/C test) of a semiconductor substrate. A T/C test is a test in which the temperature of the substrate under test is changed repeatedly. The performance information 121 includes information necessary for creating a hysteresis curve, which will be described later. For example, the performance information 121 includes the number of repetitions of the test, the amount of deformation for each temperature, and the like. Further, as shown in FIG. 2, the amount of deformation at each temperature is represented by the distance between a predetermined position of the substrate and a reference position. FIG. 2 is a diagram for explaining the amount of deformation.

設計情報122は、基板ごとの1つ以上の設計値及び品質である。図3は、設計情報の一例を示す図である。図3に示すように、設計情報122には、基板の設計値として、厚さ、層数、線膨張係数及びヤング率が含まれる。また、設計情報122には、基板の品質として、故障するまでのT/C試験の繰り返し回数が含まれる。以降、故障するまでのT/C試験の繰り返し回数を故障T/C回数と呼ぶ場合がある。 Design information 122 is one or more design values and qualities for each board. FIG. 3 is a diagram showing an example of design information. As shown in FIG. 3, the design information 122 includes the thickness, number of layers, coefficient of linear expansion, and Young's modulus as design values of the substrate. The design information 122 also includes the number of T/C test repetitions until a failure occurs as the board quality. Hereinafter, the number of times the T/C test is repeated until a failure occurs may be referred to as the number of failure T/Cs.

また、品質は、故障T/C回数のような連続値で表されてもよいし、「High」、「Medium」、「Low」のような離散値のラベルで表されてもよい。なお、ラベルは、故障T/C回数に応じて決定されてもよい。例えば、ラベルは、0<故障T/C回数<400である場合は「Low」、400≦故障T/C回数<600である場合は「Medium」、600≦故障T/C回数である場合は「High」のように決定されてもよい。 Also, the quality may be represented by a continuous value such as the number of failure T/Cs, or may be represented by discrete value labels such as "High", "Medium", and "Low". Note that the label may be determined according to the number of failure T/Cs. For example, the label is "Low" if 0 < number of failure T/Cs < 400, "Medium" if 400 ≤ number of failure T/Cs < 600, and It may be determined as "High".

図3の例では、基板αの設計値(i)(厚さ)は1mmであり、設計値(ii)(層数)は5層であり、設計値(iii)(線膨張係数)は20ppm/℃であり、設計値(iv)(ヤング率)は20GPaであり、故障T/C回数は300回である。また、基板αの品質は、ラベル「Low」で表されてもよい。 In the example of FIG. 3, the design value (i) (thickness) of the substrate α is 1 mm, the design value (ii) (number of layers) is 5 layers, and the design value (iii) (linear expansion coefficient) is 20 ppm. /°C, the design value (iv) (Young's modulus) is 20 GPa, and the number of failure T/C is 300 times. Also, the quality of the substrate α may be represented by the label “Low”.

モデル情報123は、生成されたモデルのパラメータ等である。例えば、モデル情報123は、後に説明する応答曲面を作成するためのモデル式等である。 The model information 123 is parameters and the like of the generated model. For example, the model information 123 is a model formula or the like for creating a response surface to be described later.

制御部13は、例えば、CPU(Central Processing Unit)、MPU(Micro Processing Unit)、GPU(Graphics Processing Unit)等によって、内部の記憶装置に記憶されているプログラムがRAMを作業領域として実行されることにより実現される。また、制御部13は、例えば、ASIC(Application Specific Integrated Circuit)やFPGA(Field Programmable Gate Array)等の集積回路により実現されるようにしてもよい。制御部13は、抽出部131、選択部132、生成部133及び算出部134を有する。 The control unit 13 executes a program stored in an internal storage device using, for example, a CPU (Central Processing Unit), MPU (Micro Processing Unit), GPU (Graphics Processing Unit), etc., as a work area. It is realized by Also, the control unit 13 may be implemented by an integrated circuit such as an ASIC (Application Specific Integrated Circuit) or an FPGA (Field Programmable Gate Array). The control unit 13 has an extraction unit 131 , a selection unit 132 , a generation unit 133 and a calculation unit 134 .

抽出部131は、基板を変形させる試験における基板の変形量の変化を、周波数領域で表現したときの、試験の複数の進行度合いのそれぞれに対応する所定の周波数の範囲における振幅の積算値の合計値を、基板の特徴量として抽出する。 The extraction unit 131 obtains the sum of amplitude integrated values in a predetermined frequency range corresponding to each of a plurality of degrees of progress of the test, when the change in the amount of deformation of the board in the test for deforming the board is expressed in the frequency domain. A value is extracted as a feature quantity of the substrate.

本実施例では、まず、抽出部131は、T/C試験における温度変化に対する基板の変形量の変化を計算する。なお、T/C試験は、温度を所定の範囲で繰り返し増減させることで基板を変形させる試験ということができる。 In this embodiment, first, the extraction unit 131 calculates changes in the amount of deformation of the substrate with respect to temperature changes in the T/C test. The T/C test can be said to be a test in which the substrate is deformed by repeatedly increasing and decreasing the temperature within a predetermined range.

具体的には、抽出部131は、基板の温度を40℃から140℃の間で繰り返し変化させたときの変形量を表すヒステリシス曲線を各店で微分することで、変形量の変化を得る。なお、抽出部131は、実績情報121からヒステリシス曲線を取得する。 Specifically, the extraction unit 131 obtains the change in the amount of deformation by differentiating the hysteresis curve representing the amount of deformation when the temperature of the substrate is repeatedly changed between 40° C. and 140° C. at each point. Note that the extraction unit 131 acquires the hysteresis curve from the performance information 121 .

図4は、ヒステリシス曲線を説明するための図である。ヒステリシス曲線は、T/C回数ごとの、昇温時の変形量を表す曲線と降温時の変形量を表す曲線とを組み合わせた曲線である。例えば、基板の温度は、昇温時には温度40℃から140℃に推移し、降温時には140℃から40℃に推移する。 FIG. 4 is a diagram for explaining a hysteresis curve. The hysteresis curve is a curve obtained by combining a curve representing the amount of deformation during temperature increase and a curve representing the amount of deformation during temperature decrease for each T/C number. For example, the substrate temperature changes from 40° C. to 140° C. when the temperature rises, and changes from 140° C. to 40° C. when the temperature drops.

図4の曲線201cは、T/C回数が1回の場合の昇温時の基板の変形量を表している。また、曲線202cは、T/C回数が1回の場合の降温時の基板の変形量を表している。また、曲線231cは、T/C回数が600回の場合の昇温時の基板の変形量を表している。また、曲線232cは、T/C回数が600回の場合の降温時の基板の変形量を表している。 A curve 201c in FIG. 4 represents the amount of deformation of the substrate when the temperature is raised when the number of times of T/C is one. A curve 202c represents the deformation amount of the substrate when the temperature is lowered when the number of T/C is one. A curve 231c represents the deformation amount of the substrate when the temperature is raised when the T/C number is 600 times. A curve 232c represents the deformation amount of the substrate when the temperature is lowered when the number of times of T/C is 600 times.

本実施例では、基板の塑性変形に対応する部分の変形量を基に特徴量を抽出する。ここで、弾性変形及び塑性変形について説明する。弾性変形は、荷重により変形した物体が、荷重を取り去ると元の形に戻る変形である。一方、塑性変形は、荷重により変形した物体が、荷重を取り去っても元の形に戻らない変形である。なお、T/C試験における温度変化は荷重の一例である。 In this embodiment, the feature amount is extracted based on the deformation amount of the portion corresponding to the plastic deformation of the substrate. Here, elastic deformation and plastic deformation will be described. Elastic deformation is deformation in which an object deformed by a load returns to its original shape when the load is removed. On the other hand, plastic deformation is deformation in which an object deformed by a load does not return to its original shape even after the load is removed. Note that the temperature change in the T/C test is an example of the load.

ここで、基板は高温状態で成形される。また、基板の故障に繋がる亀裂や破断は、降温時の塑性変形の際に生じる。これは、ストレスフリーである高温状態を起点として、降温時に熱応力が生じ始めるためである。 Here, the substrate is molded in hot conditions. In addition, cracks and breaks that lead to failure of the substrate occur during plastic deformation during cooling. This is because thermal stress begins to occur when the temperature is lowered, starting from a stress-free high temperature state.

弾性変形は温度変化に対して線形であるが、塑性変形では温度変化に対して変形量が低下する。このことを利用して、図5に示すように、抽出部131は、降温時のヒステリシス曲線の微分値、すなわち傾きが閾値以下である範囲の変形量を基に特徴量を抽出する。図5は、変形量の傾きを説明するための図である。 Elastic deformation is linear with respect to temperature change, but plastic deformation decreases with respect to temperature change. Utilizing this fact, as shown in FIG. 5, the extracting unit 131 extracts a feature amount based on the differential value of the hysteresis curve when the temperature is lowered, that is, the amount of deformation in the range where the slope is equal to or less than the threshold. FIG. 5 is a diagram for explaining the slope of the deformation amount.

次に、抽出部131は、塑性変形の範囲の変形量のデータを、FFT(Fast Fourier Transform)等の変換手法により、周波数に対する振幅に変換する。言い換えると、抽出部131は、温度領域で表現されていたデータを周波数領域のスペクトルに変換する。 Next, the extraction unit 131 converts the deformation amount data in the range of plastic deformation into amplitude with respect to frequency by a conversion method such as FFT (Fast Fourier Transform). In other words, the extraction unit 131 transforms the data expressed in the temperature domain into a spectrum in the frequency domain.

T/C試験を繰り返すと、熱応力で基板内に表れる不良箇所の規模が変化することにより、塑性変形モードが変化する。不良個所の規模は、図6に示すような規模感によって表現される。図6は、特徴量の抽出方法を説明するための図である。 When the T/C test is repeated, the plastic deformation mode changes due to the change in the size of the defects appearing in the substrate due to thermal stress. The scale of the defective portion is represented by a sense of scale as shown in FIG. FIG. 6 is a diagram for explaining a feature amount extraction method.

図6に示すように、ここでは、T/C回数が200回の場合、基板の四隅の辺の1/100以下程度の範囲である点を目安とする規模感の不良が生じるものとする。また、T/C回数が400回の場合、基板の面の1/100以下程度の範囲である辺を目安とする規模感の不良が生じるものとする。また、T/C回数が600回の場合、基板面積の1/100以上の範囲を目安とする規模感の不良が生じるものとする。 As shown in FIG. 6, here, when the number of times of T/C is 200, it is assumed that a defect of a sense of scale occurs with a range of about 1/100 or less of the four corner sides of the substrate as a guideline. Also, when the number of T/Cs is 400, it is assumed that a defect of a scale occurs with a side that is in the range of about 1/100 or less of the surface of the substrate as a guideline. Also, when the number of times of T/C is 600, it is assumed that a defective scale occurs in a range of 1/100 or more of the substrate area as a guideline.

このように、T/C試験を繰り返すことで、離散的に発生する点状の小規模な亀裂が増加して線状に連続的につながり、さらに面状に規模が拡がっていく。このとき、塑性変形の変形量のデータには、点状の不良に応じて高周波数の変動が加わり、辺状の不良に応じて中間周波数の変動が加わり、面状の不良に応じて低周波数の変動が加わる。 In this way, by repeating the T/C test, the number of discretely generated small-scale cracks in the form of dots increases, connects continuously in a linear pattern, and further expands in size in a plane. At this time, the deformation amount data of plastic deformation is subjected to high-frequency fluctuations corresponding to point-like defects, medium-frequency fluctuations corresponding to side-like defects, and low-frequency fluctuations corresponding to surface-like defects. change is added.

図6に示すように、加わる変動の周波数に対応して、塑性変形モードが、高周波数、中間周波数、低周波数と変化する。さらに、変形量のデータを周波数領域で表した場合、変動として加わった周波数に相当する領域に振幅のピークが現れる。 As shown in FIG. 6, the plastic deformation mode changes from high frequency, intermediate frequency, and low frequency, corresponding to the frequency of the applied fluctuation. Furthermore, when the deformation amount data is represented in the frequency domain, an amplitude peak appears in the domain corresponding to the frequency applied as the variation.

なお、基板に現れる不良の規模感は、設計の違いや個体差によって異なることが考えられる。そのため、図6を用いて説明した規模感は目安に過ぎず、高、中間、低周波のそれぞれに対するT/C回数は、規模感に関係なくあらゆる基板について固定的に決められたものであってもよい。 It is conceivable that the scale of defects appearing on the board may differ depending on differences in design and individual differences. Therefore, the sense of scale explained using FIG. 6 is only a guideline, and the number of T/Cs for each of high, intermediate, and low frequencies is fixedly determined for all substrates regardless of sense of scale. good too.

抽出部131は、試験の複数の繰り返し回数のそれぞれに対応する周波数の範囲における振幅の積算値の合計値を、基板の特徴量として抽出する。本実施例では、抽出部131は、T/C回数200回に対応する1~3[10-2Hz]の周波数の範囲の振幅の第1の積算値を計算する。また、抽出部131は、T/C回数400回に対応する4~8[10-3Hz]の周波数の範囲の振幅の第2の積算値を計算する。また、抽出部131は、T/C回数600回に対応する4~8[10-4Hz]の周波数の範囲の振幅の第3の積算値を計算する。そして、抽出部131は、第1の積算値と第2の積算値と第3の積算値の合計値を計算し、計算した合計値を特徴量として抽出する。 The extraction unit 131 extracts the total value of the amplitude integrated values in the frequency range corresponding to each of the plurality of repetition times of the test as the feature amount of the board. In this embodiment, the extraction unit 131 calculates the first integrated value of the amplitude in the frequency range of 1 to 3 [10 -2 Hz] corresponding to 200 T/C times. Also, the extraction unit 131 calculates a second integrated value of the amplitude in the frequency range of 4 to 8 [10 −3 Hz] corresponding to 400 T/C times. Also, the extraction unit 131 calculates the third integrated value of the amplitude in the frequency range of 4 to 8 [10 −4 Hz] corresponding to 600 T/C times. Then, the extraction unit 131 calculates the total value of the first integrated value, the second integrated value, and the third integrated value, and extracts the calculated total value as a feature amount.

抽出部131は、周波数の振幅を周波数で積分した値を積算値として得ることができる。図6に示すように、不良個所の規模感の変化に応じて、周波数領域に新たな振幅のピークが表れ、振幅の積分値も増加する。このため、抽出部131は、振幅の積分値を基に、基板の熱疲労の蓄積を表す特徴量を計算することができる。 The extraction unit 131 can obtain a value obtained by integrating the amplitude of the frequency with respect to the frequency as an integrated value. As shown in FIG. 6, a new amplitude peak appears in the frequency domain and the integral value of the amplitude increases as the scale of the defective portion changes. Therefore, the extracting unit 131 can calculate the feature amount representing the accumulation of thermal fatigue of the substrate based on the integrated value of the amplitude.

なお、T/C試験は、基板を変形させる試験、及び温度を所定の範囲で繰り返し増減させることで基板を変形させる試験の一例である。ここでの試験は、T/C試験に限定されず、例えば、圧力及び電圧等を加えて物体を変形させる試験であってもよい。 Note that the T/C test is an example of a test in which the substrate is deformed and a test in which the substrate is deformed by repeatedly increasing and decreasing the temperature within a predetermined range. The test here is not limited to the T/C test, and may be, for example, a test in which an object is deformed by applying pressure, voltage, or the like.

選択部132は、基板の設計値と基板が故障に至るまでの試験の繰り返し回数との関係を表すモデルを基に、設計値の中から、基板が故障に至るまでの試験の繰り返し回数に対する寄与率が高い順に所定の数の設計値を選択する。 The selection unit 132 selects, from among the design values, the contribution to the number of test repetitions until the board fails, based on a model representing the relationship between the board design value and the number of test repetitions until the board fails. A predetermined number of design values are selected in descending order of rate.

本実施例では、選択部132は、図3に示す各設計値を説明変数とし、品質を目的変数とする教師ありの学習モデルにおいて、寄与率が大きい高い順に2つの設計値を選択する。また、学習モデルは、モデルベース特徴量選択手法であるランダムフォレストであってもよい。例えば、選択部132は、品質の予測に用いる設計値を入れ替えながら予測精度を評価し、予測精度が大きく変わった場合は入れ替えた設計値の寄与率を上げて、そうでない場合は寄与率を下げることで、各設計値の寄与率を計算する。 In this embodiment, the selection unit 132 selects two design values in descending order of contribution rate in a supervised learning model with each design value shown in FIG. 3 as an explanatory variable and quality as an objective variable. Also, the learning model may be a random forest, which is a model-based feature selection technique. For example, the selection unit 132 evaluates the prediction accuracy while replacing the design values used for quality prediction, and increases the contribution rate of the replaced design value when the prediction accuracy has changed significantly, and otherwise decreases the contribution rate. By doing so, the contribution rate of each design value is calculated.

生成部133は、特徴量と基板の所定の設計値との関係を表すモデルを生成する。具体的には、生成部133は、抽出部131によって抽出された特徴量と、選択部132によって選択された設計値との関係を表すモデルを生成する。本実施例では、生成部133は、複数の繰り返し回数ごとに、設計値を入力とし、特徴量を応答とする応答曲面を表すモデルを生成する。 The generation unit 133 generates a model representing the relationship between the feature amount and the predetermined design value of the board. Specifically, the generation unit 133 generates a model representing the relationship between the feature amount extracted by the extraction unit 131 and the design value selected by the selection unit 132 . In this embodiment, the generation unit 133 generates a model representing a response surface having a design value as an input and a feature amount as a response for each of a plurality of iterations.

図7を用いて、生成部133によるモデルの生成方法を説明する。図7は、モデルの生成方法を説明するための図である。ここで、選択部132は、設計値(i)及び設計値(iii)を選択したものとする。まず、生成部133は、各基板の設計値に対するT/C回数ごとの特徴量を、モデルの空間にプロットする。ここでは、モデルの空間は3次元であり、x軸が設計値(i)、y軸が設計値(iii)、z軸が特徴量に対応しているものとする。 A model generation method by the generation unit 133 will be described with reference to FIG. FIG. 7 is a diagram for explaining a model generation method. Here, it is assumed that the selection unit 132 selects the design value (i) and the design value (iii). First, the generation unit 133 plots the feature amount for each T/C number with respect to the design value of each board in the space of the model. Here, it is assumed that the model space is three-dimensional, with the x-axis corresponding to the design value (i), the y-axis corresponding to the design value (iii), and the z-axis corresponding to the feature amount.

図7の例では、xα1は、基板αの設計値に対する、T/C回数が200回の場合の特徴量である。また、xα2は、基板αの設計値に対する、T/C回数が400回の場合の特徴量である。また、xα3は、基板αの設計値に対する、T/C回数が600回の場合の特徴量である。また、xβ1は、基板βの設計値に対する、T/C回数が200回の場合の特徴量である。また、xβ2は、基板βの設計値に対する、T/C回数が400回の場合の特徴量である。また、xβ3は、基板βの設計値に対する、T/C回数が600回の場合の特徴量である。 In the example of FIG. 7, x α1 is the feature quantity when the number of times of T/C is 200 times with respect to the design value of the substrate α. Also, x α2 is a feature amount when the number of times of T/C is 400 times with respect to the design value of the substrate α. Also, x α3 is a feature amount when the number of times of T/C is 600 times with respect to the design value of the substrate α. Also, x β1 is a feature amount when the number of T/Cs is 200 times with respect to the design value of the substrate β. Also, x β2 is a feature amount when the T/C count is 400 times with respect to the design value of the substrate β. Also, x β3 is a feature amount when the T/C count is 600 times with respect to the design value of the substrate β.

さらに、生成部133は、各基板のT/C回数が故障T/C回数である場合の特徴量をプロットする。なお、故障T/C回数は、基板が完全に破損したときのT/C回数に限られず、半導体への使用に耐えられない程度に電気的な特性が失われた時点のT/C回数等であってもよい。また、基板が完全に破損した場合、有効な特徴量を得ることが困難になるため、破損する直前のT/C試験で得られた特徴量が、故障T/C回数における特徴量であってもよい。 Furthermore, the generation unit 133 plots the feature amount when the number of T/Cs of each substrate is the number of failure T/Cs. The number of failure T/Cs is not limited to the number of T/Cs when the substrate is completely damaged, but the number of T/Cs when the electrical characteristics are lost to the extent that it cannot be used for semiconductors. may be In addition, if the substrate is completely damaged, it becomes difficult to obtain effective feature values. good too.

そして、生成部133は、プロットの結果から応答曲面210S、220S、230S及び290Sを生成する。生成部133は、プロットしたデータに対してガウス分布を重ね合わせて曲面を作成するガウシアンプロセスにより応答曲面を生成することができる。以降、これらの応答曲面を疲労蓄積モデルと呼ぶ場合がある。 Then, the generator 133 generates response surfaces 210S, 220S, 230S and 290S from the plotted results. The generation unit 133 can generate a response surface by a Gaussian process that superimposes a Gaussian distribution on the plotted data to create a surface. Henceforth, these response surfaces may be called a fatigue accumulation model.

応答曲面210Sは、T/C回数が200回の場合の特徴量に対応している。また、応答曲面220Sは、T/C回数が400回の場合の特徴量に対応している。また、応答曲面230Sは、T/C回数が600回の場合の特徴量に対応している。また、応答曲面290Sは、T/C回数が故障T/C回数の場合の特徴量に対応している。 The response surface 210S corresponds to the feature amount when the T/C count is 200 times. Also, the response surface 220S corresponds to the feature amount when the T/C count is 400 times. Also, the response surface 230S corresponds to the feature amount when the T/C count is 600 times. Also, the response curved surface 290S corresponds to the feature amount when the number of T/Cs is the number of failure T/Cs.

算出部134は、生成部133によって生成された疲労蓄積モデルを用いて、基板の設計値から特徴量の予測値を算出する。算出部134は、応答曲線モデルを基に、設計値が既知であり、特徴量が未知である対象基板の、T/C回数と特徴量の関係を表す回帰式を求め、当該回帰式を使って予測値を算出する。 The calculation unit 134 uses the fatigue accumulation model generated by the generation unit 133 to calculate the predicted value of the feature amount from the board design values. Based on the response curve model, the calculation unit 134 obtains a regression equation representing the relationship between the number of times of T/C and the feature value of the target substrate whose design value is known and whose feature value is unknown, and uses the regression equation. to calculate the predicted value.

図8を用いて、予測値の算出方法を説明する。図8は、予測値の算出方法を説明するための図である。算出部134は、まず、対象基板の設計値を疲労蓄積モデルに入力し、疲労蓄積式を得る。ここで、疲労蓄積式は、T/C回数に対する特徴量を表す式である。 A method of calculating the predicted value will be described with reference to FIG. FIG. 8 is a diagram for explaining a method of calculating predicted values. The calculator 134 first inputs the design values of the target substrate into the fatigue accumulation model to obtain the fatigue accumulation formula. Here, the fatigue accumulation formula is a formula representing a feature amount with respect to the number of times of T/C.

疲労蓄積式からは200回、400回等の特定のT/C回数についての特徴量しか得られない。そのため、算出部134は、疲労蓄積式に基づき、最小二乗法を用いた回帰分析により回帰式を求める。回帰式は、例えばy=ax2+bx+c(ただし、yは特徴量、xはT/C回数、a、b、cは係数)のような二次式であってもよい。そして、算出部134は、回帰式から任意のT/C回数における特徴量及び故障T/C回数の予測値を得ることができる。なお、図8のx1はT/C回数が200回の場合の特徴量である。また、x2はT/C回数が400回の場合の特徴量である。また、x3はT/C回数が600回の場合の特徴量である。また、x9はT/C回数が故障T/C回数の場合の特徴量である。 From the fatigue accumulation formula, we can only obtain feature values for specific T/C times such as 200 times and 400 times. Therefore, the calculation unit 134 obtains a regression formula by regression analysis using the least squares method based on the fatigue accumulation formula. The regression formula may be a quadratic formula such as y=ax 2 +bx+c (where y is a feature amount, x is the number of T/Cs, and a, b, and c are coefficients). Then, the calculation unit 134 can obtain the feature quantity and the predicted value of the number of failure T/Cs at an arbitrary number of T/Cs from the regression equation. Note that x1 in FIG. 8 is the feature amount when the T/C count is 200 times. Also, x 2 is the feature amount when the T/C count is 400 times. Also, x 3 is a feature amount when the T/C count is 600 times. Also, x9 is a feature quantity when the number of T/Cs is the number of failure T/Cs.

さらに、算出部134は、基板の設計値、及び基板に対して所定の回数だけ試験を行って得られた特徴量を基に、モデルを用いて所定の回数より大きい回数における特徴量を算出することができる。例えば、算出部134は、対象基板について、予測される故障T/C回数の1/10の回数だけ試験を行って得られた特徴量を用いてもよい。 Further, the calculation unit 134 calculates the feature amount for a number of times greater than a predetermined number of times using a model, based on the design value of the board and the feature amount obtained by performing the test a predetermined number of times on the board. be able to. For example, the calculation unit 134 may use the feature amount obtained by performing the test on the target substrate for the number of times of 1/10 of the predicted number of failure T/Cs.

図8の例では、予測されるT/C故障回数は750回なので、算出部134は、75回のT/C試験で得られた特徴量を実測値として回帰式と同じ平面上にプロットする。これにより、算出部134は、回帰性の妥当性を確認することや、回帰式の微調整を行うことができる。例えば、算出部134は、プロットした特徴量と対応する回帰式の値との誤差が閾値以下であれば、当該回帰式を採用し、故障T/C回数やT/C試験を終了する回数を求めるようにしてもよい。また、算出部134は、最小二乗法等を用いて、プロットした特徴量に対してさらに回帰式を最適化するようにしてもよい。 In the example of FIG. 8, the predicted number of T/C failures is 750, so the calculation unit 134 plots the feature values obtained from the 75 T/C tests as actual measurements on the same plane as the regression equation. . As a result, the calculation unit 134 can confirm the validity of the regression and finely adjust the regression equation. For example, if the error between the plotted feature quantity and the value of the corresponding regression formula is equal to or less than a threshold, the calculation unit 134 adopts the regression formula, You can ask for it. Further, the calculation unit 134 may further optimize the regression equation for the plotted feature amount using the method of least squares or the like.

また、特徴量はT/C回数に対して急激に増加する場合がある。このため、図9に示すように、回帰式は対数等の所定のスケールで表されてもよい。その場合も同様に、算出部134は、特定の特徴量に対するT/C回数の予測値を算出することができる。 Also, the feature amount may increase rapidly with the number of T/Cs. Therefore, as shown in FIG. 9, the regression equation may be represented by a predetermined scale such as logarithm. In that case, similarly, the calculation unit 134 can calculate the predicted value of the number of times of T/C for a specific feature amount.

[処理の流れ]
図10を用いて、予測装置10による予測処理の流れを説明する。図10は、予測処理の流れを示すフローチャートである。図10に示すように、まず、予測装置10は、T/C試験の実績情報121から、疲労を表す特徴量を抽出する(ステップS11)。
[Process flow]
The flow of prediction processing by the prediction device 10 will be described with reference to FIG. FIG. 10 is a flowchart showing the flow of prediction processing. As shown in FIG. 10, first, the prediction device 10 extracts a feature quantity representing fatigue from the T/C test performance information 121 (step S11).

次に、予測装置10は、基板の設計値のうち、品質に寄与率の高い設計値を選択する(ステップS12)。そして、予測装置10は、選択した設計値に対する特徴量の応答曲面を生成する(ステップS13)。ステップS11からS13までの処理は、学習工程の処理である。学習工程の完了後には、運用工程の処理が行われる。 Next, the prediction device 10 selects a design value with a high contribution rate to quality from the design values of the board (step S12). Then, the prediction device 10 generates a response surface of the feature quantity for the selected design value (step S13). The processing from steps S11 to S13 is the processing of the learning process. After the learning process is completed, the operation process is processed.

運用工程において、まず予測装置10は、予測対象基板の設計値を基に、応答曲面から予測対象基板の疲労蓄積式を抽出する(ステップS21)。そして、予測装置10は、疲労蓄積式を用いて、故障T/C回数の予測値を算出する(ステップS22)。また、予測装置10は、特定の特徴量に至るT/C回数の予測値を算出してもよい。 In the operation process, the prediction device 10 first extracts the fatigue accumulation formula of the prediction target board from the response surface based on the design values of the prediction target board (step S21). Then, the prediction device 10 uses the fatigue accumulation formula to calculate the predicted value of the number of failure T/Cs (step S22). Also, the prediction device 10 may calculate a predicted value of the number of times of T/C leading to a specific feature amount.

図11を用いて、疲労を表す特徴量を抽出する処理(図10のステップS11)の詳細を説明する。図11は、特徴量を抽出する処理の流れを示すフローチャートである。図11に示すように、まず、予測装置10は、実績情報121に基づくヒステリシス曲線を各点で微分して傾きを算出する(ステップS111)。次に、予測装置10は、ヒステリシス曲線の傾きが閾値以下の範囲を塑性変形データとして抽出する(ステップS112)。 The details of the process of extracting the feature quantity representing fatigue (step S11 in FIG. 10) will be described with reference to FIG. FIG. 11 is a flow chart showing the flow of processing for extracting feature amounts. As shown in FIG. 11, first, the prediction device 10 differentiates the hysteresis curve based on the performance information 121 at each point to calculate the slope (step S111). Next, the prediction device 10 extracts a range in which the slope of the hysteresis curve is equal to or less than the threshold as plastic deformation data (step S112).

ここで、予測装置10は、塑性変形データをFFTにより周波数領域のデータに変換する(ステップS113)。そして、予測装置10は、周波数領域のデータから、高、中間、低周波数領域の振幅の積算値を算出する(ステップS114)。さらに、予測装置10は、高、中間、低周波数領域の振幅の積算値の合計値を特徴量として算出する(ステップS115)。 Here, the prediction device 10 transforms the plastic deformation data into frequency domain data by FFT (step S113). Then, the prediction device 10 calculates integrated values of the amplitudes in the high, intermediate, and low frequency regions from the frequency domain data (step S114). Furthermore, the prediction device 10 calculates the total value of the integrated values of the amplitudes in the high, intermediate, and low frequency regions as a feature amount (step S115).

[効果]
上述したように、抽出部131は、基板を変形させる試験における基板の変形量の変化を、周波数領域で表現したときの、試験の複数の進行度合いのそれぞれに対応する所定の周波数の範囲における振幅の積算値の合計値を、基板の特徴量として抽出する。また、生成部133は、特徴量と、基板の所定の設計値との関係を表すモデルを生成する。また、算出部134は、生成部133によって生成されたモデルを用いて、基板の設計値から特徴量の予測値を算出する。
[effect]
As described above, the extraction unit 131 extracts amplitudes in a predetermined frequency range corresponding to each of a plurality of degrees of progress of the test when the change in the amount of deformation of the board in the test for deforming the board is expressed in the frequency domain. is extracted as the feature amount of the substrate. The generation unit 133 also generates a model representing the relationship between the feature amount and the predetermined design value of the board. The calculation unit 134 also uses the model generated by the generation unit 133 to calculate the predicted value of the feature amount from the design value of the board.

このように、予測装置10は、基板の設計値だけでなく、疲労を表す特徴量を使って故障に至るまでの試験の回数等を予測することができる。このため、本実施例によれば、設計値が変更された場合であっても、基板の品質の予測精度の低下を抑止することができる。 In this manner, the prediction device 10 can predict the number of tests until failure, etc., using not only the design value of the board but also the characteristic quantity representing fatigue. Therefore, according to the present embodiment, even when the design values are changed, it is possible to prevent deterioration of the prediction accuracy of the board quality.

さらに、本実施例によれば、精度良く品質を予測することができるようになるため、品質試験にようする工数を短縮することができる。また、T/C試験の実績情報を蓄積していくことで応答曲面の精度が向上し、その結果、品質の予測精度が向上する。 Furthermore, according to this embodiment, the quality can be predicted with high accuracy, so that the man-hour required for the quality test can be reduced. In addition, by accumulating T/C test performance information, the accuracy of the response surface improves, and as a result, the accuracy of quality prediction improves.

また、抽出部131は、温度を所定の範囲で繰り返し増減させることで基板を変形させる試験における、温度変化に対する基板の変形量の変化を、周波数に対する振幅に変換する。さらに、抽出部131は、変換したデータから、試験の複数の繰り返し回数のそれぞれに対応する周波数の範囲における振幅の積算値の合計値を、基板の特徴量として抽出する。 In addition, the extraction unit 131 converts the change in the amount of deformation of the substrate with respect to temperature change in a test in which the substrate is deformed by repeatedly increasing and decreasing the temperature within a predetermined range, into amplitude with respect to frequency. Furthermore, the extracting unit 131 extracts, from the converted data, the total value of integrated values of the amplitudes in the frequency range corresponding to each of the plurality of repetitions of the test, as the characteristic amount of the board.

このように、本実施例はT/C試験に適用することができる。また、本実施例では、各周波数の範囲における振幅の蓄積値を合計した1次元の特徴量が得られるため、モデル生成のための十分な情報を確保しつつ、計算量を抑えることが可能になる。 Thus, this embodiment can be applied to the T/C test. In addition, in this embodiment, a one-dimensional feature amount is obtained by summing the accumulated amplitude values in each frequency range, so it is possible to reduce the amount of calculation while securing sufficient information for model generation. Become.

また、選択部132は、基板の設計値と基板が故障に至るまでの試験の繰り返し回数との関係を表すモデルを基に、設計値の中から、基板が故障に至るまでの試験の繰り返し回数に対する寄与率が高い順に所定の数の設計値を選択する。また、生成部133は、特徴量と、選択部132によって選択された設計値との関係を表すモデルを生成する。このように、寄与率の高い設計値を選択することで、モデルの精度を向上させつつ、計算量を抑えることが可能になる。 Further, the selection unit 132 selects the number of test repetitions until failure of the board from among the design values based on the model representing the relationship between the design value of the board and the number of repetitions of the test until failure of the board. A predetermined number of design values are selected in descending order of the contribution rate to . The generation unit 133 also generates a model representing the relationship between the feature quantity and the design value selected by the selection unit 132 . In this way, by selecting a design value with a high contribution rate, it is possible to reduce the amount of calculation while improving the accuracy of the model.

また、生成部133は、複数の繰り返し回数ごとに、設計値を入力とし、特徴量を応答とする応答曲面を表すモデルを生成する。これにより、設計値が判明していれば、予測値を算出するための計算式を容易に得ることが可能になる。 In addition, the generation unit 133 generates a model representing a response surface having design values as inputs and characteristic amounts as responses for each of a plurality of iterations. As a result, if the design value is known, it becomes possible to easily obtain the calculation formula for calculating the predicted value.

また、算出部134は、基板の設計値、及び基板に対して所定の回数だけ試験を行って得られた特徴量を基に、モデルを用いて所定の回数より大きい回数における特徴量を算出する。これにより、少ない試験回数で予測値の算出精度をさらに向上させることが可能になる。 Further, the calculation unit 134 uses a model to calculate the feature amount for the number of times greater than the predetermined number of times, based on the design value of the board and the feature amount obtained by performing the test a predetermined number of times on the board. . This makes it possible to further improve the calculation accuracy of the predicted value with a small number of tests.

なお、本実施例では、予測装置10が学習工程及び運用工程の両方の処理を行うものとして説明した。一方で、学習工程と運用工程は異なる装置によって行われてもよい。その場合、例えば、生成装置が学習工程の処理を実行し、予測装置が運用工程の処理を実行する。 In addition, in the present embodiment, the prediction device 10 has been described as performing both the learning process and the operation process. On the other hand, the learning process and the operational process may be performed by different devices. In that case, for example, the generation device executes the processing of the learning process, and the prediction device executes the processing of the operation process.

また、予測装置10によって生成されるモデルは応答曲面に限られない。例えば、予測装置10は、1つ以上の設計値及びT/C回数を説明変数とし、特徴量又は故障T/C回数を目的変数とする多変数の回帰モデルを生成してもよい。 Also, the model generated by the prediction device 10 is not limited to a response surface. For example, the prediction device 10 may generate a multi-variable regression model with one or more design values and the number of T/Cs as explanatory variables, and a feature amount or the number of failure T/Cs as an objective variable.

[システム]
上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に変更することができる。また、実施例で説明した具体例、分布、数値等は、あくまで一例であり、任意に変更することができる。
[system]
Information including processing procedures, control procedures, specific names, and various data and parameters shown in the above documents and drawings can be arbitrarily changed unless otherwise specified. Further, the specific examples, distributions, numerical values, etc. described in the examples are only examples, and can be arbitrarily changed.

また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散や統合の具体的形態は図示のものに限られない。つまり、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。さらに、各装置にて行われる各処理機能は、その全部又は任意の一部が、CPU及び当該CPUにて解析実行されるプログラムにて実現され、あるいは、ワイヤードロジックによるハードウェアとして実現され得る。 Also, each component of each device illustrated is functionally conceptual, and does not necessarily need to be physically configured as illustrated. That is, the specific forms of distribution and integration of each device are not limited to those shown in the drawings. That is, all or part of them can be functionally or physically distributed and integrated in arbitrary units according to various loads and usage conditions. Furthermore, all or any part of each processing function performed by each device can be implemented by a CPU and a program analyzed and executed by the CPU, or implemented as hardware based on wired logic.

[ハードウェア]
図12は、ハードウェア構成例を説明する図である。図12に示すように、予測装置10は、通信インタフェース10a、HDD(Hard Disk Drive)10b、メモリ10c、プロセッサ10dを有する。また、図12に示した各部は、バス等で相互に接続される。
[hardware]
FIG. 12 is a diagram illustrating a hardware configuration example. As shown in FIG. 12, the prediction device 10 has a communication interface 10a, a HDD (Hard Disk Drive) 10b, a memory 10c, and a processor 10d. 12 are interconnected by a bus or the like.

通信インタフェース10aは、ネットワークインタフェースカード等であり、他のサーバとの通信を行う。HDD10bは、図1に示した機能を動作させるプログラムやDBを記憶する。 The communication interface 10a is a network interface card or the like, and communicates with other servers. The HDD 10b stores programs and DBs that operate the functions shown in FIG.

プロセッサ10dは、図1に示した各処理部と同様の処理を実行するプログラムをHDD10b等から読み出してメモリ10cに展開することで、図1等で説明した各機能を実行するプロセスを動作させるハードウェア回路である。すなわち、このプロセスは、予測装置10が有する各処理部と同様の機能を実行する。具体的には、プロセッサ10dは、抽出部131、選択部132、生成部133及び算出部134と同様の機能を有するプログラムをHDD10b等から読み出す。そして、プロセッサ10dは、抽出部131、選択部132、生成部133及び算出部134等と同様の処理を実行するプロセスを実行する。 The processor 10d reads from the HDD 10b or the like a program that executes the same processing as each processing unit shown in FIG. It is a wear circuit. That is, this process performs the same function as each processing unit of the prediction device 10 . Specifically, the processor 10d reads a program having functions similar to those of the extraction unit 131, the selection unit 132, the generation unit 133, and the calculation unit 134 from the HDD 10b or the like. Then, the processor 10d executes processes for executing the same processes as those of the extraction unit 131, the selection unit 132, the generation unit 133, the calculation unit 134, and the like.

このように予測装置10は、プログラムを読み出して実行することで学習類方法を実行する情報処理装置として動作する。また、予測装置10は、媒体読取装置によって記録媒体から上記プログラムを読み出し、読み出された上記プログラムを実行することで上記した実施例と同様の機能を実現することもできる。なお、この他の実施例でいうプログラムは、予測装置10によって実行されることに限定されるものではない。例えば、他のコンピュータ又はサーバがプログラムを実行する場合や、これらが協働してプログラムを実行するような場合にも、本発明を同様に適用することができる。 In this manner, the prediction device 10 operates as an information processing device that executes a learning method by reading and executing a program. Also, the prediction device 10 can read the program from the recording medium by the medium reading device, and execute the read program to realize the same function as the embodiment described above. It should be noted that the programs referred to in other embodiments are not limited to being executed by the prediction device 10 . For example, the present invention can be applied in the same way when another computer or server executes the program, or when they cooperate to execute the program.

このプログラムは、インターネット等のネットワークを介して配布することができる。また、このプログラムは、ハードディスク、フレキシブルディスク(FD)、CD-ROM、MO(Magneto-Optical disk)、DVD(Digital Versatile Disc)等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することができる。 This program can be distributed via a network such as the Internet. In addition, this program is recorded on a computer-readable recording medium such as a hard disk, flexible disk (FD), CD-ROM, MO (Magneto-Optical disk), DVD (Digital Versatile Disc), etc., and can be downloaded from the recording medium by a computer. It can be executed by being read.

10 予測装置
11 インタフェース部
12 記憶部
13 制御部
121 実績情報
122 設計情報
123 モデル情報
131 抽出部
132 選択部
133 生成部
134 算出部
201c、202c、231c、232c 曲線
210S、220S、230S、290S 応答曲面
10 prediction device 11 interface unit 12 storage unit 13 control unit 121 performance information 122 design information 123 model information 131 extraction unit 132 selection unit 133 generation unit 134 calculation unit 201c, 202c, 231c, 232c curves 210S, 220S, 230S, 290S response surface

Claims (8)

基板を変形させる試験における前記基板の変形量の変化を、周波数領域で表現したときの、前記試験の複数の進行度合いのそれぞれに対応する所定の周波数の範囲における振幅の積算値の合計値を、前記基板の特徴量として抽出する抽出部と、
前記特徴量と、前記基板の所定の設計値との関係を表すモデルを生成する生成部と、
前記生成部によって生成されたモデルを用いて、基板の設計値から特徴量の予測値を算出する算出部と、
を有することを特徴とする予測装置。
When the change in the amount of deformation of the substrate in the test for deforming the substrate is expressed in the frequency domain, the total value of the integrated values of the amplitudes in a predetermined frequency range corresponding to each of the plurality of degrees of progress of the test, an extraction unit that extracts as a feature amount of the substrate;
a generation unit that generates a model representing a relationship between the feature amount and a predetermined design value of the substrate;
a calculation unit that calculates a predicted value of the feature amount from the design value of the board using the model generated by the generation unit;
A prediction device characterized by comprising:
前記抽出部は、温度を所定の範囲で繰り返し増減させることで基板を変形させる試験における、温度変化に対する前記基板の変形量の変化を、周波数に対する振幅に変換したデータから、前記試験の複数の繰り返し回数のそれぞれに対応する周波数の範囲における振幅の積算値の合計値を、前記基板の特徴量として抽出することを特徴とする請求項1に記載の予測装置。 The extraction unit extracts a change in the amount of deformation of the substrate with respect to temperature change in a test in which the substrate is deformed by repeatedly increasing and decreasing the temperature within a predetermined range, from data converted into amplitude with respect to frequency, a plurality of repetitions of the test. 2. The prediction apparatus according to claim 1, wherein a total value of integrated values of amplitudes in a frequency range corresponding to each number of times is extracted as the feature amount of the substrate. 前記基板の設計値と前記基板が故障に至るまでの前記試験の繰り返し回数との関係を表すモデルを基に、前記設計値の中から、前記基板が故障に至るまでの前記試験の繰り返し回数に対する寄与率が高い順に所定の数の設計値を選択する選択部をさらに有し、
前記生成部は、前記特徴量と、前記選択部によって選択された設計値との関係を表すモデルを生成することを特徴とする請求項2に記載の予測装置。
Based on the model representing the relationship between the design value of the board and the number of repetitions of the test until the board fails, the number of repetitions of the test until the board fails is selected from the design values further comprising a selection unit that selects a predetermined number of design values in descending order of contribution rate;
3. The prediction device according to claim 2, wherein the generation unit generates a model representing a relationship between the feature amount and the design value selected by the selection unit.
前記生成部は、前記複数の繰り返し回数ごとに、前記設計値を入力とし、前記特徴量を応答とする応答曲面を表すモデルを生成することを特徴とする請求項2に記載の予測装置。 3. The prediction apparatus according to claim 2, wherein the generation unit receives the design value as an input and generates a model representing a response surface having the feature value as a response for each of the plurality of iterations. 前記算出部は、前記基板の設計値、及び前記基板に対して所定の回数だけ試験を行って得られた前記特徴量を基に、前記モデルを用いて前記所定の回数より大きい回数における特徴量を算出することを特徴とする請求項1に記載の予測装置。 The calculation unit uses the model based on the design value of the board and the feature amount obtained by performing the test a predetermined number of times on the board. 2. The prediction device according to claim 1, which calculates 基板を変形させる試験における前記基板の変形量の変化を、周波数領域で表現したときの、前記試験の複数の進行度合いのそれぞれに対応する所定の周波数の範囲における振幅の積算値の合計値を、前記基板の特徴量として抽出し、
前記特徴量と、前記基板の所定の設計値との関係を表すモデルを生成し、
前記生成する処理によって生成されたモデルを用いて、基板の設計値から特徴量の予測値を算出する
処理をコンピュータに実行させることを特徴とする予測プログラム。
When the change in the amount of deformation of the substrate in the test for deforming the substrate is expressed in the frequency domain, the total value of the integrated values of the amplitudes in a predetermined frequency range corresponding to each of the plurality of degrees of progress of the test, extracted as a feature quantity of the substrate,
generating a model representing a relationship between the feature amount and a predetermined design value of the substrate;
A prediction program for causing a computer to execute a process of calculating a predicted value of a feature amount from a design value of a substrate using the model generated by the generating process.
基板を変形させる試験における前記基板の変形量の変化を、周波数領域で表現したときの、前記試験の複数の進行度合いのそれぞれに対応する所定の周波数の範囲における振幅の積算値の合計値を、前記基板の特徴量として抽出する抽出部と、
前記特徴量と、前記基板の所定の設計値との関係を表すモデルを生成する生成部と、
を有することを特徴とする生成装置。
When the change in the amount of deformation of the substrate in the test for deforming the substrate is expressed in the frequency domain, the total value of the integrated values of the amplitudes in a predetermined frequency range corresponding to each of the plurality of degrees of progress of the test, an extraction unit that extracts as a feature amount of the substrate;
a generation unit that generates a model representing a relationship between the feature amount and a predetermined design value of the substrate;
A generating device comprising:
基板を変形させる試験における前記基板の変形量の変化を、周波数領域で表現したときの、前記試験の複数の進行度合いのそれぞれに対応する所定の周波数の範囲における振幅の積算値の合計値を、前記基板の特徴量として抽出し、
前記特徴量と、前記基板の所定の設計値との関係を表すモデルを生成する
処理をコンピュータに実行させることを特徴とする生成プログラム。
When the change in the amount of deformation of the substrate in the test for deforming the substrate is expressed in the frequency domain, the total value of the integrated values of the amplitudes in a predetermined frequency range corresponding to each of the plurality of degrees of progress of the test, extracted as a feature quantity of the substrate,
A generating program that causes a computer to execute a process of generating a model representing a relationship between the feature amount and a predetermined design value of the substrate.
JP2019204724A 2019-11-12 2019-11-12 Prediction device, prediction program, generation device and generation program Active JP7310558B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019204724A JP7310558B2 (en) 2019-11-12 2019-11-12 Prediction device, prediction program, generation device and generation program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019204724A JP7310558B2 (en) 2019-11-12 2019-11-12 Prediction device, prediction program, generation device and generation program

Publications (2)

Publication Number Publication Date
JP2021076521A JP2021076521A (en) 2021-05-20
JP7310558B2 true JP7310558B2 (en) 2023-07-19

Family

ID=75898637

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019204724A Active JP7310558B2 (en) 2019-11-12 2019-11-12 Prediction device, prediction program, generation device and generation program

Country Status (1)

Country Link
JP (1) JP7310558B2 (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006084248A (en) 2004-09-15 2006-03-30 Furukawa Electric Co Ltd:The Reliability evaluation method for fatigue behavior of soldered joints
JP2008241432A (en) 2007-03-27 2008-10-09 Toshiba Corp Load calculation device and load calculation method
JP2015232578A (en) 2015-09-09 2015-12-24 株式会社東芝 Electronic component and measuring method
JP2016080677A (en) 2014-10-09 2016-05-16 日本精工株式会社 Fatigue life evaluation method, fatigue life evaluation system, and fatigue life evaluation program
JP2019121326A (en) 2018-01-11 2019-07-22 富士通株式会社 Distortion calculation program, distortion calculation method, and information processing device for calculating distortion

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006084248A (en) 2004-09-15 2006-03-30 Furukawa Electric Co Ltd:The Reliability evaluation method for fatigue behavior of soldered joints
JP2008241432A (en) 2007-03-27 2008-10-09 Toshiba Corp Load calculation device and load calculation method
JP2016080677A (en) 2014-10-09 2016-05-16 日本精工株式会社 Fatigue life evaluation method, fatigue life evaluation system, and fatigue life evaluation program
JP2015232578A (en) 2015-09-09 2015-12-24 株式会社東芝 Electronic component and measuring method
JP2019121326A (en) 2018-01-11 2019-07-22 富士通株式会社 Distortion calculation program, distortion calculation method, and information processing device for calculating distortion

Also Published As

Publication number Publication date
JP2021076521A (en) 2021-05-20

Similar Documents

Publication Publication Date Title
Singh et al. Bearing remaining useful life estimation using an adaptive data-driven model based on health state change point identification and K-means clustering
JP5011830B2 (en) DATA PROCESSING METHOD, DATA PROCESSING PROGRAM, RECORDING MEDIUM CONTAINING THE PROGRAM, AND DATA PROCESSING DEVICE
Ma et al. Locally linear embedding on Grassmann manifold for performance degradation assessment of bearings
CN101520652B (en) A method for evaluating the service reliability of NC equipment
Shen et al. Recognition of rolling bearing fault patterns and sizes based on two-layer support vector regression machines
CN111275235B (en) System and method for optimizing a manufacturing process based on inspection of a component
Lu et al. Prognosis of bearing degradation using gradient variable forgetting factor RLS combined with time series model
JP7707124B2 (en) Information processing device, information processing method, and program
Behzad et al. Prognostics of rolling element bearings with the combination of paris law and reliability method
Ruan et al. Improvement by Monte Carlo for trajectory similarity-based RUL prediction
Bondoc et al. Employing LIVE digital twin in prognostic and health management: Identifying location of the sensors
Dorst et al. Influence of synchronization within a sensor network on machine learning results
Kim et al. A robust health indicator for rotating machinery under time-varying operating conditions
Bernstein et al. Reliability prediction for microelectronics
JP7310558B2 (en) Prediction device, prediction program, generation device and generation program
Lu et al. Bearing health assessment based on chaotic characteristics
An et al. Improved MCMC method for parameter estimation based on marginal probability density function
Chen et al. A multi-source data fusion driven power field effect transistor health state assessment and remaining useful life prediction method
Jeon et al. Extremely rare anomaly detection pipeline in semiconductor bonding process with digital twin-driven data augmentation method
CN109829245B (en) Bearing fault early warning method and device
CN119990051B (en) A flip chip FC packaging technology method
Huynh et al. Degradation modeling and RUL estimation of deteriorating systems in s-plane
Thukral et al. Ai-enabled board level vibration testing: Unveiling the physics of degradation
Sun et al. Reliability Analysis with Multiple Dependent Features from a Vibration‐Based Accelerated Degradation Test
Thoppil et al. An integrated learning algorithm for vibration feature selection and remaining useful life estimation of lathe spindle unit

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220708

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230510

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

R150 Certificate of patent or registration of utility model

Ref document number: 7310558

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150