JP2011058888A - Apparatus and method for predicting crack initiation life - Google Patents
Apparatus and method for predicting crack initiation life Download PDFInfo
- Publication number
- JP2011058888A JP2011058888A JP2009207309A JP2009207309A JP2011058888A JP 2011058888 A JP2011058888 A JP 2011058888A JP 2009207309 A JP2009207309 A JP 2009207309A JP 2009207309 A JP2009207309 A JP 2009207309A JP 2011058888 A JP2011058888 A JP 2011058888A
- Authority
- JP
- Japan
- Prior art keywords
- load
- analysis
- element size
- crack
- life
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000000977 initiatory effect Effects 0.000 title claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 102
- 229910000679 solder Inorganic materials 0.000 claims abstract description 61
- 238000012937 correction Methods 0.000 claims abstract description 32
- 239000000758 substrate Substances 0.000 claims description 11
- 238000005259 measurement Methods 0.000 abstract description 17
- 238000012360 testing method Methods 0.000 description 18
- 238000001816 cooling Methods 0.000 description 14
- 239000000463 material Substances 0.000 description 5
- 230000035939 shock Effects 0.000 description 4
- 230000001186 cumulative effect Effects 0.000 description 3
- 238000010438 heat treatment Methods 0.000 description 3
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 241000352457 Shivajiella indica Species 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013001 point bending Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000005476 soldering Methods 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Electric Connection Of Electric Components To Printed Circuits (AREA)
Abstract
Description
本発明は、温度変化による熱負荷および荷重負荷が作用するはんだ接合部の亀裂発生寿命を予測する亀裂発生寿命予測装置および亀裂発生寿命予測方法に関する。 The present invention relates to a crack generation life prediction apparatus and a crack generation life prediction method for predicting a crack generation life of a solder joint to which a thermal load and a load load are applied due to a temperature change.
従来より、例えば電子部品をはんだ付けした表面実装型の基板では、電子部品のはんだ接合部に温度変化による熱負荷が繰り返し作用することによって亀裂が発生し、この亀裂が進展すると導通不良に至ることがある。このため、はんだ接合部における亀裂発生寿命を、寿命予測式を用いて予測することが行われている（例えば、特許文献１参照）。特許文献１の寿命予測では、はんだ接合部の解析モデルに熱負荷となる冷熱サイクルを繰り返し付与することによって生じるひずみ振幅を有限要素法（ＦＥＭ）を用いた解析により求め、このひずみ振幅から寿命予測式に基づいてはんだ接合部の亀裂発生寿命を算出している。 Conventionally, for example, in a surface-mount type substrate soldered with an electronic component, a crack is generated by repeatedly applying a thermal load due to a temperature change at a solder joint portion of the electronic component, and when this crack progresses, a conduction failure is caused. There is. For this reason, the crack generation life in a solder joint part is predicted using a life prediction formula (see, for example, Patent Document 1). In the life prediction of Patent Document 1, a strain amplitude generated by repeatedly applying a thermal cycle as a thermal load to an analysis model of a solder joint is obtained by analysis using a finite element method (FEM), and life prediction is performed from this strain amplitude. The crack generation life of the solder joint is calculated based on the equation.
しかしながら、上記寿命予測式は、はんだ接合部に熱負荷が作用する場合のみを考慮したものであるため、この熱負荷に加えて荷重負荷が作用する場合には、亀裂発生寿命を正確に予測することができないという問題がある。
そこで、はんだ接合部の解析モデルに、荷重負荷を作用させた状態で冷熱サイクルを付与することによって、はんだ接合部に生じる平均ひずみを有限要素法を用いた解析により求め、この平均ひずみから所定の低下定数に基づいて算出した寿命の低下率によって、上記寿命予測式を修正することが提案されている（例えば、非特許文献１参照）。低下定数は、亀裂発生寿命を試験により実測して求めた実測寿命と、亀裂発生寿命を解析により求めた解析寿命とに基づいて求められる。
However, since the above life prediction formula considers only the case where a thermal load acts on the solder joint, when the load load acts in addition to this thermal load, the crack generation life is accurately predicted. There is a problem that can not be.
Therefore, by applying a thermal cycle with a load applied to the analysis model of the solder joint, the average strain generated in the solder joint is obtained by analysis using the finite element method, and a predetermined strain is determined from this average strain. It has been proposed to correct the lifetime prediction formula based on the lifetime decrease rate calculated based on the decrease constant (see, for example, Non-Patent Document 1). The decrease constant is obtained based on the measured life obtained by actually measuring the crack generation life by the test and the analytical life obtained by analyzing the crack generation life.
ところで、上記平均ひずみを有限要素法を用いた解析によって求める際に、解析領域を複数の要素にメッシュ分割することが行われるが、その要素サイズは、一般的にはんだ接合部の解析領域の大きさに応じて異なるサイズに設定される。
例えば、電子部品が小型の場合は、はんだ接合部における解析領域が小さいため、解析時の要素サイズは、平均ひずみを高い解析精度で求めることができるように十数μｍと小さく設定される。一方、大型の電子部品、例えば車載コネクタの場合は、そのはんだ接合部における解析領域が大きいため、高い解析精度を求めて要素サイズを小さく設定すると、膨大な解析時間を必要とする。このため、大型の電子部品では、解析時間の短縮化を図るため、要素サイズを約５０μｍと大きく設定している。
By the way, when the average strain is obtained by analysis using the finite element method, the analysis area is divided into a plurality of elements, and the element size is generally larger than the analysis area of the solder joint. Depending on the size, different sizes are set.
For example, when the electronic component is small, since the analysis region in the solder joint is small, the element size at the time of analysis is set to a small value of several tens of micrometers so that the average strain can be obtained with high analysis accuracy. On the other hand, in the case of a large-sized electronic component, for example, an in-vehicle connector, the analysis area in the solder joint is large. Therefore, if an element size is set small for high analysis accuracy, a huge amount of analysis time is required. For this reason, in a large electronic component, the element size is set to a large value of about 50 μm in order to shorten the analysis time.
このように、要素サイズは解析領域の大きさによって異なるサイズに設定されるため、設定された要素サイズ毎にこの要素サイズに応じた適切な低下定数を求める必要がある。しかも、低下定数は、上述のように、実測と解析とにより求められるものであるため、低下定数を求めるために膨大な時間を必要とする。
したがって、この低下定数に基づく低下率の算出、ひいてはこの低下率に基づいて算出される亀裂発生寿命の予測に膨大な時間を必要とするという問題があった。
Thus, since the element size is set to a different size depending on the size of the analysis region, it is necessary to obtain an appropriate reduction constant corresponding to the element size for each set element size. Moreover, since the reduction constant is obtained by actual measurement and analysis as described above, it takes a long time to obtain the reduction constant.
Therefore, there is a problem that enormous time is required for calculation of the reduction rate based on this reduction constant, and in turn, the prediction of the crack generation life calculated based on this reduction rate.
本発明は前記問題点に鑑みてなされたものであり、荷重負荷に起因するひずみを有限要素法を用いて解析する際の要素サイズに応じて、亀裂発生寿命を短時間で予測することができる亀裂発生寿命予測装置および亀裂発生寿命予測方法を提供することを目的としている。 The present invention has been made in view of the above problems, and it is possible to predict the crack generation life in a short time according to the element size when analyzing strain caused by load using the finite element method. An object of the present invention is to provide a crack generation life prediction apparatus and a crack generation life prediction method.
本発明の亀裂発生寿命予測装置は、温度変化による熱負荷が作用するとともに荷重負荷が作用するはんだ接合部の亀裂発生寿命を予測する亀裂発生寿命予測装置であって、所定の第一要素サイズでメッシュ分割した前記はんだ接合部の解析モデルに、前記荷重負荷を作用させた状態で冷熱サイクルを付与することによって、前記熱負荷に起因するひずみ振幅、および当該荷重負荷に起因する平均ひずみを有限要素法を用いた解析により求める解析部と、前記第一要素サイズと異なるサイズの第二要素サイズにおいて平均ひずみから低下率を算出するための関係式を用いて、前記解析部で求めた平均ひずみから、前記荷重負荷に起因する前記亀裂発生寿命の低下率を算出する第一算出部と、前記解析部で求めた平均ひずみに対して、前記第一要素サイズと異なるサイズの第二要素サイズに基づく低下定数、および前記第一要素サイズに応じた所定の補正パラメータを作用させて、前記荷重負荷に起因する前記亀裂発生寿命の低下率を算出する第一算出部と、前記解析部で求めたひずみ振幅、および前記第一算出部で算出した低下率に基づいて、はんだ接合部の亀裂発生寿命を算出する第二算出部と、を備え、前記第一算出部は、前記関係式を用いて算出された値を、前記第一要素サイズに応じて所定の補正パラメータにより補正して低下率を求めることを特徴とする。 The crack generation life prediction apparatus of the present invention is a crack generation life prediction apparatus that predicts the crack generation life of a solder joint where a thermal load due to a temperature change acts and a load load acts, and has a predetermined first element size. By applying a cooling cycle with the load applied to the analysis model of the solder joint part divided into meshes, the strain amplitude resulting from the thermal load and the average strain resulting from the load are finite elements. From the average strain obtained by the analysis unit using the relational expression for calculating the decrease rate from the average strain in the second element size different from the first element size, and the analysis unit obtained by analysis using the method A first calculation unit for calculating a decrease rate of the crack initiation life due to the load load, and an average strain obtained by the analysis unit, A reduction constant based on a second element size of a size different from the size and a predetermined correction parameter corresponding to the first element size are applied to calculate a reduction rate of the crack initiation life due to the load load A calculation unit; and a second calculation unit that calculates a crack initiation life of the solder joint based on the strain amplitude obtained by the analysis unit and the decrease rate calculated by the first calculation unit, The calculation unit is characterized in that a value calculated using the relational expression is corrected with a predetermined correction parameter in accordance with the first element size to obtain a reduction rate.
また、本発明の亀裂発生寿命予測方法は、温度変化による熱負荷が作用するとともに荷重負荷が作用するはんだ接合部の亀裂発生寿命を予測する亀裂発生寿命予測方法であって、所定の第一要素サイズでメッシュ分割した前記はんだ接合部の解析モデルに、前記荷重負荷を作用させた状態で冷熱サイクルを付与することによって、前記熱負荷に起因するひずみ振幅、および当該荷重負荷に起因する平均ひずみを有限要素法を用いた解析により求める解析工程と、前記第一要素サイズと異なるサイズの第二要素サイズにおいて平均ひずみから低下率を算出するための関係式を用いて、前記解析部で求めた平均ひずみから、前記荷重負荷に起因する前記亀裂発生寿命の低下率を算出する第一算出工程と、前記解析工程で求めたひずみ振幅、および前記第一算出工程で算出した低下率に基づいて、はんだ接合部の亀裂発生寿命を算出する第二算出工程と、を備え、前記第一算出部は、前記関係式を用いて算出された値を、前記第一要素サイズに応じて所定の補正パラメータにより補正して低下率を求めることを特徴とする。 The crack generation life prediction method of the present invention is a crack generation life prediction method for predicting the crack generation life of a solder joint where a thermal load due to a temperature change acts and a load load acts. By applying a cooling cycle with the load applied to the analysis model of the solder joint divided into meshes by size, the strain amplitude caused by the heat load and the average strain caused by the load load are calculated. Using the relational expression for calculating the decrease rate from the average strain in the second element size different from the first element size, the analysis step obtained by the analysis using the finite element method, the average obtained by the analysis unit A first calculation step of calculating a decrease rate of the crack initiation life resulting from the load load from strain, a strain amplitude obtained in the analysis step, and A second calculation step of calculating a crack generation life of the solder joint based on the decrease rate calculated in the first calculation step, wherein the first calculation unit is a value calculated using the relational expression Is reduced by a predetermined correction parameter in accordance with the first element size to obtain a reduction rate.
本発明によれば、荷重負荷に起因する平均ひずみから低下率を算出する際に、第一要素サイズと異なるサイズの第二要素サイズにおいて平均ひずみから低下率を算出するための関係式により算出された値を、第一要素サイズに応じた補正パラメータにより補正することにより、低下率を算出することができる。このため、第二要素サイズに応じた低下定数を予め実測および解析により求めておけば、この低下定数に補正パラメータにより補正するだけで、第一要素サイズに応じた低下定数を実測および解析により求めることなく低下率を算出することができる。したがって、解析時の要素サイズに応じて、はんだ接合部における亀裂発生寿命を短時間で予測することができる。 According to the present invention, when calculating the reduction rate from the average strain caused by the load, it is calculated by a relational expression for calculating the reduction rate from the average strain in the second element size different from the first element size. The rate of decrease can be calculated by correcting the obtained value with a correction parameter corresponding to the first element size. For this reason, if the reduction constant corresponding to the second element size is obtained in advance by actual measurement and analysis, the reduction constant corresponding to the first element size is obtained by actual measurement and analysis only by correcting the reduction constant by the correction parameter. The reduction rate can be calculated without any problem. Therefore, the crack generation life in the solder joint can be predicted in a short time according to the element size at the time of analysis.
また、前記第二要素サイズが、前記第一要素サイズよりも小さいことが好ましい。この場合、第一要素サイズを大きく設定することができるので、この第一要素サイズにより平均ひずみを解析する際の解析時間を短縮することができる。したがって、亀裂発生寿命をより短時間で予測することができる。 The second element size is preferably smaller than the first element size. In this case, since the first element size can be set large, the analysis time for analyzing the average strain can be shortened by the first element size. Therefore, the crack generation life can be predicted in a shorter time.
また、前記第二算出部が、Ｍａｎｓｏｎ−Ｃｏｆｆｉｎ則に基づく下記式により前記亀裂発生寿命を算出することが好ましい。
Ｎ＝１０００×（Δεi／０．０１）-1.24×（１−ｋ×αεa）
ここで、Ｎは亀裂発生寿命、Δεiはひずみ振幅、ｋは補正パラメータ、αは低下定数、εaは平均ひずみ、ｋ×αεaは低下率（ただし、０≦ｋ×αεa≦１）である。
この場合は、Ｍａｎｓｏｎ−Ｃｏｆｆｉｎ則に基づく上記式により亀裂発生寿命を算出するようにしたので、はんだ接合部における亀裂発生寿命を的確に予測することができる。
Moreover, it is preferable that said 2nd calculation part calculates the said crack generation lifetime by the following formula based on Manson-Coffin rule.
N = 1000 × (Δε i /0.01 ) −1.24 × (1−k × αε a )
Here, N is the crack life, Δε i is the strain amplitude, k is the correction parameter, α is the reduction constant, ε a is the average strain, and k × αε a is the reduction rate (where 0 ≦ k × αε a ≦ 1) It is.
In this case, since the crack generation life is calculated by the above formula based on the Manson-Coffin rule, the crack generation life in the solder joint can be accurately predicted.
また、前記はんだ接合部が、基板上に実装された車載コネクタのはんだ接合部であることが好ましい。
この場合、車載コネクタのはんだ接合部において、解析時の要素サイズに応じて、熱負荷および荷重負荷が作用する場合のはんだ接合部における亀裂発生寿命を短時間で予測することができる。
Moreover, it is preferable that the said solder joint part is a solder joint part of the vehicle-mounted connector mounted on the board | substrate.
In this case, in the solder joint portion of the in-vehicle connector, the crack generation life in the solder joint portion when a thermal load and a load load act can be predicted in a short time according to the element size at the time of analysis.
本発明によれば、第一要素サイズに応じた低下定数を実測および解析により求めることなく低下率を算出することができるので、解析時の要素サイズに応じて、はんだ接合部における亀裂発生寿命を短時間で予測することができる。 According to the present invention, since the reduction rate can be calculated without obtaining a reduction constant according to the first element size by actual measurement and analysis, the crack generation life in the solder joint can be reduced according to the element size at the time of analysis. Can be predicted in a short time.
以下、本発明の実施の形態を図面に基づいて説明する。
図１は、本発明の一実施形態に係る亀裂発生寿命予測装置により寿命予測されるはんだ接合部を有する車載コネクタの概略平面図であり、図２は、図１におけるＡ矢視図である。
図１および図２において、車載コネクタ１は、ハウジング２と、その後面（図１の下側）に接続された複数のリード端子３と、左右両側面に固定された略Ｌ字型の固定ペグ４とを備えている。車載コネクタ１は、このリード端子３の先端部および固定ペグ４の水平部４ａを電子制御ユニット（ＥＣＵ）の基板５にそれぞれはんだ付けされることにより、この基板５上に実装されている。
Hereinafter, embodiments of the present invention will be described with reference to the drawings.
FIG. 1 is a schematic plan view of an in-vehicle connector having a solder joint portion whose life is predicted by a crack generation life prediction apparatus according to an embodiment of the present invention, and FIG. 2 is a view as seen from the direction of arrow A in FIG.
1 and 2, an in-vehicle connector 1 includes a housing 2, a plurality of lead terminals 3 connected to the rear surface (the lower side in FIG. 1), and a substantially L-shaped fixed peg fixed to both left and right side surfaces. 4 is provided. The in-vehicle connector 1 is mounted on the substrate 5 by soldering the tip portion of the lead terminal 3 and the horizontal portion 4a of the fixed peg 4 to the substrate 5 of the electronic control unit (ECU).
ハウジング２の前面（図１の上側）には複数の雄側端子（図示せず）が配設されており、この雄側端子に相手側のハウジング６の雌側端子（図示せず）が嵌合されている。この相手側のハウジング６にはワイヤハーネス７の一端が接続されている。 A plurality of male terminals (not shown) are disposed on the front surface of the housing 2 (upper side in FIG. 1), and female terminals (not shown) of the mating housing 6 are fitted into the male terminals. Are combined. One end of a wire harness 7 is connected to the counterpart housing 6.
車載コネクタ１および基板５は、車内の雰囲気温度の変化が激しいため、この雰囲気温度の変化によって熱膨張したり熱収縮したりする。その際、車載コネクタ１と基板５とでは材料が異なるため、両者の熱膨張率の差によってはんだ接合部８（図２参照）には熱負荷によるせん断力が作用する。
また、ワイヤハーネス７は、車内での配線スペースが限られているため、無理に引っ張られた状態で車載コネクタ１に接続されていることが多い。この場合、はんだ接合部８には、ワイヤハーネス７の引張荷重による荷重負荷が一方向（図１の矢符Ｂ方向）に作用することになる。
Since the in-vehicle connector 1 and the substrate 5 have a drastic change in the ambient temperature in the vehicle, the in-vehicle connector 1 and the substrate 5 are thermally expanded or contracted by the change in the ambient temperature. At that time, since the in-vehicle connector 1 and the substrate 5 are made of different materials, a shear force due to a thermal load acts on the solder joint portion 8 (see FIG. 2) due to a difference in thermal expansion coefficient between them.
Further, since the wiring harness 7 has a limited wiring space in the vehicle, the wire harness 7 is often connected to the in-vehicle connector 1 while being forcibly pulled. In this case, a load load due to the tensile load of the wire harness 7 acts on the solder joint portion 8 in one direction (the arrow B direction in FIG. 1).
したがって、車載コネクタ１のはんだ接合部８には、上記熱負荷が繰り返し作用するとともに上記荷重負荷が常に作用することにより亀裂が発生し、この亀裂が進展する。最悪の場合は、はんだ接合部８が剥離して基板５から車載コネクタ１が外れることになる。
本実施形態では、亀裂発生寿命予測装置１０により、複合負荷（熱負荷および荷重負荷）が作用する車載コネクタ１のはんだ接合部８における亀裂発生寿命を予測するものである。
Therefore, a crack is generated in the solder joint portion 8 of the in-vehicle connector 1 when the thermal load repeatedly acts and the load load always acts, and the crack progresses. In the worst case, the solder joint 8 is peeled off and the in-vehicle connector 1 is detached from the substrate 5.
In this embodiment, the crack occurrence life prediction apparatus 10 predicts the crack occurrence life in the solder joint portion 8 of the in-vehicle connector 1 on which the composite load (thermal load and load load) acts.
図３は、亀裂発生寿命予測装置１０の全体構成を示すブロック図であり、図４は、亀裂発生寿命予測装置１０により亀裂発生寿命Ｎを予測する方法を示すフローチャートである。
図３および図４において、亀裂発生寿命予測装置１０は、上記熱負荷に起因するひずみ振幅Δεiおよび上記荷重負荷に起因する平均ひずみεa50をそれぞれ有限要素法を用いた解析により求める解析部１１（解析工程）と、低下率β50を算出する第一算出部１２（第一算出工程）と、はんだ接合部８の亀裂発生寿命Ｎを算出する第二算出部１３（第二算出工程）とによって構成されている。
FIG. 3 is a block diagram showing the overall configuration of the crack occurrence life prediction apparatus 10, and FIG. 4 is a flowchart showing a method for predicting the crack occurrence life N by the crack occurrence life prediction apparatus 10.
3 and 4, the crack generation life prediction apparatus 10 includes an analysis unit 11 that obtains a strain amplitude Δε i caused by the thermal load and an average strain ε a50 caused by the load load by analysis using a finite element method. (Analysis step), a first calculation unit 12 (first calculation step) for calculating the decrease rate β 50 , and a second calculation unit 13 (second calculation step) for calculating the crack occurrence life N of the solder joint 8 It is constituted by.
解析部１１（解析工程）では、亀裂のないはんだ接合部８をモデル化した解析モデルを用いて、はんだ接合部８におけるひずみ振幅Δεiおよび平均ひずみεa50を求める（ステップＳ１）。解析モデルは、例えば１／２モデルにするとともに、一次の四面体要素でメッシュ分割する。メッシュ分割時の要素サイズ（第一要素サイズ）は、５０μｍとし、解析精度よりも解析時間を優先した値に設定している。 The analysis unit 11 (analysis process) obtains the strain amplitude Δε i and the average strain ε a50 in the solder joint 8 using an analysis model obtained by modeling the solder joint 8 having no crack (step S1). The analysis model is, for example, a ½ model and is divided into meshes by primary tetrahedral elements. The element size (first element size) at the time of mesh division is set to 50 μm, and the analysis time is set to have a priority over the analysis accuracy.
また、ワイヤハーネス７は梁要素でモデル化し、その先端に車載コネクタ１の垂直上方への荷重負荷を加えるように点荷重を付与する。さらに、はんだ材料の特性としては、塑性特性およびクリープ特性を与える。温度条件としては、例えば−４０〜１２０℃（各３０分）／サイクルとする冷熱サイクルを付与する。 The wire harness 7 is modeled with a beam element, and a point load is applied to the tip of the wire harness 7 so as to apply a load load vertically above the in-vehicle connector 1. Furthermore, as the characteristics of the solder material, plastic characteristics and creep characteristics are given. As a temperature condition, for example, a cooling / heating cycle of −40 to 120 ° C. (each 30 minutes) / cycle is applied.
上記条件により解析を行うと、冷熱サイクル中のはんだ接合部８に生じる応力およびひずみは、図５に示すように、最初の１サイクルでＯ1点（原点）からＡ点を経由してＢ点に至り、次の１サイクルでＢ点からＡ点を経由してＢ点に戻る。そして、その後の冷熱サイクル中における応力およびひずみは、Ａ点とＢ点との間を一定振幅で変動する。したがって、この冷熱サイクル中におけるひずみ振幅Δεiは、冷熱サイクルを最初に２回繰り返した時点で求めることができる。 When the analysis is performed under the above conditions, the stress and strain generated in the solder joint 8 during the cooling / heating cycle are as shown in FIG. 5 in the first cycle from point O 1 (origin) to point B via point A. In the next cycle, the point B returns to the point B via the point A. The stress and strain during the subsequent cooling cycle fluctuate between point A and point B with a constant amplitude. Therefore, the strain amplitude Δε i during the cooling cycle can be obtained when the cooling cycle is first repeated twice.
一方、冷熱サイクル中には上記荷重負荷も作用するため、はんだ接合部８にはこの荷重負荷に起因してクリープ変形が生じる。このため、冷熱サイクル中のはんだ接合部８に生じるひずみは、クリープ変形により時間経過とともに徐々に増大する。このひずみの増大により、図５に示すように、初期のひずみ振幅Δεiの中心となるＯ1点は、所定時間経過後にはＯ3点まで移動し、このＯ3点を中心としたひずみ振幅Δεiとなる。
このように、荷重負荷に起因するひずみは、時間経過とともに変化するものであるため、冷熱サイクルが所定回数（例えば３０回）経過したときの累積ひずみから、１サイクル当たりのひずみ（図５のＯ2点におけるひずみ）を算出し、このひずみを荷重負荷に起因する平均ひずみεa50として用いる。
On the other hand, since the above load is also applied during the cooling cycle, creep deformation occurs in the solder joint 8 due to this load. For this reason, the distortion which arises in the solder joint part 8 in a cooling-heat cycle increases gradually with time progress by creep deformation. Due to this increase in strain, as shown in FIG. 5, the O 1 point that is the center of the initial strain amplitude Δε i moves to the O 3 point after a predetermined time has elapsed, and the strain amplitude centered on this O 3 point. Δε i .
As described above, since the strain caused by the load changes with time, the strain per cycle (O in FIG. 5) is calculated from the accumulated strain when a predetermined number of times (for example, 30 times) of the thermal cycle. (Strain at two points) is calculated, and this strain is used as the average strain ε a50 due to the load.
第一算出部１２（第一算出工程）では、解析部１１で求めた平均ひずみεa50から低下率β50を算出する（ステップＳ２）。低下率βは、荷重負荷に起因する亀裂発生寿命Ｎの低下率を示すものである。すなわち、低下率βは、熱負荷および荷重負荷に起因する亀裂発生寿命Ｎを試験により実測して求めた実測寿命と、熱負荷のみに起因する亀裂発生寿命を解析により求めた解析寿命との間に生じる誤差の程度を示すものであり、上記実測に基づく低下率βの算出式は、下記式（１）となる。
β＝１−実測寿命／解析寿命 ・・・（１）
In the first calculation unit 12 (first calculation step), the decrease rate β 50 is calculated from the average strain ε a50 obtained by the analysis unit 11 (step S2). The decrease rate β indicates the decrease rate of the crack generation life N caused by the load. That is, the reduction rate β is between the measured life obtained by actually measuring the thermal load and the crack initiation life N caused by the load load by the test, and the analytical life obtained by analyzing the crack initiation life caused only by the thermal load. The calculation formula of the decrease rate β based on the actual measurement is the following formula (1).
β = 1−actual life / analysis life (1)
図６は、三点曲げ試験等により求めた実測寿命、および要素サイズ（第二要素サイズ）を１２．５μｍとして解析により求めた解析寿命から上記式（１）に基づいて低下率β12.5を算出し、この低下率β12.5と第二要素サイズで解析された平均ひずみεa12.5との関係を示すグラフである。低下率β12.5と平均ひずみεa12.5とは、図６のグラフから、平均ひずみεa12.5の増大にともなって低下率β12.5が直線的に増大する比例関係にあることがわかる。したがって、その比例定数を低下定数α12.5とすると、第二要素サイズにおいて平均ひずみεa12.5から低下率β12.5を算出するための関係式は、β12.5＝αεa12.5となる。この関係式を使用すると、第一要素サイズに応じた低下率β50を、下記式（２）により、第二要素サイズに基づく平均ひずみεa12.5を用いて算出することができる。
β50＝ｋ×αεa12.5 ・・・（２）
ここで、ｋは補正パラメータであり、β50は０≦β50≦１である。
FIG. 6 shows the decrease rate β 12.5 based on the above equation (1) from the measured life obtained by the three-point bending test and the analysis life obtained by analysis with the element size (second element size) being 12.5 μm. FIG. 6 is a graph showing the relationship between the decrease rate β 12.5 and the average strain ε a12.5 analyzed with the second element size. It can be seen from the graph of FIG. 6 that the decrease rate β 12.5 and the average strain ε a12.5 are in a proportional relationship in which the decrease rate β 12.5 increases linearly as the average strain ε a12.5 increases. Therefore, if the proportionality constant is the reduction constant α 12.5 , the relational expression for calculating the reduction rate β 12.5 from the average strain ε a12.5 in the second element size is β 12.5 = αε a12.5 . When this relational expression is used, the reduction rate β 50 corresponding to the first element size can be calculated using the average strain ε a12.5 based on the second element size according to the following expression (2).
β 50 = k × αε a12.5 (2)
Here, k is a correction parameter, and β 50 is 0 ≦ β 50 ≦ 1.
補正パラメータｋは、第二要素サイズに基づく低下定数α12.5を、第一要素サイズに適応した低下定数に補正するためのパラメータであり、本実施形態ではｋ＝２として設定されている。
このように、平均ひずみεa50の解析時における第一要素サイズに応じて、予め実測および解析により求められた第二要素サイズに基づく低下定数α12.5に補正パラメータｋを乗算することにより、第一要素サイズに応じた低下率β50を算出することができる。
The correction parameter k is a parameter for correcting the decrease constant α 12.5 based on the second element size to a decrease constant adapted to the first element size, and is set as k = 2 in the present embodiment.
In this way, by multiplying the reduction parameter α 12.5 based on the second element size obtained in advance by actual measurement and analysis by the correction parameter k according to the first element size at the time of analysis of the average strain ε a50 , The reduction rate β 50 according to the element size can be calculated.
第二算出部１３（第二算出工程）では、解析部１１で求めたひずみ振幅Δεi、および第一算出部１２で算出した低下率β50に基づいて、はんだ接合部８の亀裂発生寿命Ｎを算出する（ステップＳ３）。亀裂発生寿命Ｎは、はんだ接合部８に亀裂が発生するまでに繰り返される冷熱サイクル数で表され、下記式（３）に基づいて算出される。
Ｎ＝１０００×（Δεi／０．０１）-1.24×（１−β50） ・・・（３）
ここで、式（３）中の１０００×（Δεi／０．０１）-1.24は、熱負荷が作用している場合におけるＭａｎｓｏｎ−Ｃｏｆｆｉｎ則に基づく亀裂発生寿命の予測式である。
In the second calculation unit 13 (second calculation step), based on the strain amplitude Δε i obtained by the analysis unit 11 and the decrease rate β 50 calculated by the first calculation unit 12, the crack occurrence life N of the solder joint portion 8. Is calculated (step S3). The crack generation life N is represented by the number of cooling cycles repeated until a crack is generated in the solder joint 8 and is calculated based on the following formula (3).
N = 1000 × (Δε i /0.01)-1.24×(1-β 50 ) (3)
Here, 1000 × (Δε i /0.01)−1.24 in the formula (3) is a prediction formula of the crack generation life based on the Manson-Coffin rule when a thermal load is applied.
このように、熱負荷（ひずみ振幅Δεi）に起因する予測式と、荷重負荷（平均ひずみεa50）に起因する低下率β50とによって、亀裂発生寿命Ｎを算出するようにしたので、複合負荷が作用する場合の亀裂発生寿命Ｎを、試験を行うことなく、有限要素法による解析に基づいて算出することができる。 As described above, the crack initiation life N is calculated based on the prediction formula resulting from the thermal load (strain amplitude Δε i ) and the decrease rate β 50 resulting from the load (average strain ε a50 ). The crack generation life N when a load is applied can be calculated based on analysis by a finite element method without performing a test.
以上、本発明の実施形態に係る亀裂発生寿命予測装置１０によれば、車載コネクタ１におけるはんだ接合部８の解析モデルを第一要素サイズ（５０μｍ）でメッシュ分割して解析することにより、荷重負荷に起因する平均ひずみεa50を求め、この平均ひずみεa50に対して、第一要素サイズと異なるサイズであってかつ予め実測および解析により求められた第二要素サイズ（１２．５μｍ）に基づく低下定数εa12.5、および第一要素サイズに応じた所定の補正パラメータｋを作用させて、亀裂発生寿命Ｎの低下率β50を算出することができる。このため、第二要素サイズに応じた低下定数α12.5を予め実測および解析により求めておけば、この低下定数α12.5に補正パラメータｋを作用させるだけで、第一要素サイズに応じた低下定数α50を実測および解析により求めることなく低下率β50を算出することができる。したがって、解析時の要素サイズに応じて、はんだ接合部８における亀裂発生寿命Ｎを短時間で予測することができる。 As described above, according to the crack occurrence life prediction apparatus 10 according to the embodiment of the present invention, by analyzing the analysis model of the solder joint portion 8 in the in-vehicle connector 1 by dividing the mesh with the first element size (50 μm), obtain an average strain epsilon a50 due to decrease the relative average strain epsilon a50, based on the second element size determined by a size different from the first element size and measuring beforehand and analysis (12.5 .mu.m) The reduction rate β 50 of the crack initiation life N can be calculated by applying the constant ε a12.5 and a predetermined correction parameter k corresponding to the first element size. For this reason, if the reduction constant α 12.5 corresponding to the second element size is obtained in advance by actual measurement and analysis, the reduction constant α corresponding to the first element size can be obtained simply by applying the correction parameter k to the reduction constant α 12.5. The decrease rate β 50 can be calculated without obtaining 50 by actual measurement and analysis. Therefore, the crack generation life N in the solder joint 8 can be predicted in a short time according to the element size at the time of analysis.
また、第一要素サイズを第二要素サイズよりも大きく設定することができるので、この第一要素サイズにより平均ひずみεa50を解析する際の解析時間を短縮することができる。したがって、亀裂発生寿命をより短時間で予測することができる。
また、亀裂発生寿命Ｎを、熱負荷に起因するひずみ振幅Δεiを用いたＭａｎｓｏｎ−Ｃｏｆｆｉｎ則に基づく上記式（３）により算出するようにしたので、はんだ接合部８における亀裂発生寿命Ｎを的確に予測することができる。
In addition, since the first element size can be set larger than the second element size, the analysis time for analyzing the average strain ε a50 can be shortened by the first element size. Therefore, the crack generation life can be predicted in a shorter time.
In addition, since the crack initiation life N is calculated by the above equation (3) based on the Manson-Coffin rule using the strain amplitude Δε i caused by the thermal load, the crack initiation life N in the solder joint 8 is accurately determined. Can be predicted.
次に、本発明を実施例に基づいて、さらに詳細に説明するが、本発明は、かかる実施例のみに限定されるものではない。
この実施例では、亀裂発生寿命予測装置１０により算出した亀裂発生寿命Ｎの予測精度を、亀裂発生後の亀裂進展について試験および有限要素法を用いた解析（以下、ＦＥＭ解析という。）を行うことにより比較検証した。
Next, the present invention will be described in more detail based on examples, but the present invention is not limited to such examples.
In this embodiment, the prediction accuracy of the crack generation life N calculated by the crack generation life prediction apparatus 10 is subjected to a test and analysis using a finite element method (hereinafter referred to as FEM analysis) for crack propagation after the crack is generated. And verified by comparison.
試験は、サーマルショック試験（冷熱衝撃試験）により行った。試験に使用した車載コネクタ１のハウジング２の材質は、はんだ接合部８における亀裂の発生および進展を把握し易いように、線膨張係数の大きい材質とした。また、はんだ材料は、Ｓｎ３．０Ａｇ０．５Ｃｕを使用した。
試験では、車載コネクタ１が下側となるように基板５を治具に固定し、サーマルショック試験槽内に配置した。温度条件としては、実使用環境よりも厳しい−４０〜１２０℃（各３０分）／サイクルとする冷熱サイクルを付与する。また、試験条件は、ワイヤハーネス７に錘を吊り下げてはんだ接合部８に剥離方向へ荷重負荷を作用させる負荷条件と、ワイヤハーネス７に錘を吊り下げない無負荷条件との２パターンとし、各条件でそれぞれ試験を行った。なお、サイクル数は各試験条件ともに、３０００サイクルまで試験を行い、はんだ接合部８に生じる亀裂の進展長さを調査した。
The test was conducted by a thermal shock test (cold thermal shock test). The material of the housing 2 of the in-vehicle connector 1 used for the test was a material having a large linear expansion coefficient so that the occurrence and progress of cracks in the solder joints 8 can be easily grasped. Moreover, Sn3.0Ag0.5Cu was used for the solder material.
In the test, the board | substrate 5 was fixed to the jig | tool so that the vehicle-mounted connector 1 might become a lower side, and it has arrange | positioned in the thermal shock test tank. As the temperature condition, a cooling / heating cycle that is -40 to 120 ° C. (each 30 minutes) / cycle, which is severer than the actual use environment, is applied. The test conditions are two patterns: a load condition in which a weight is suspended from the wire harness 7 and a load is applied to the solder joint 8 in the peeling direction, and a no-load condition in which no weight is suspended from the wire harness 7. Each test was performed under each condition. The number of cycles was tested up to 3000 cycles for each test condition, and the length of cracks that occurred in the solder joint 8 was investigated.
ＦＥＭ解析は、図７に示すように、亀裂進展解析により亀裂進展寿命を予測するフローチャートに沿って行った。解析工程（ステップＳ１０１）、第一算出工程（ステップＳ１０２）および第二算出工程（ステップＳ１０３）は、上記実施形態における解析工程、第一算出工程および第二算出工程と同一の工程である。すなわち、解析工程で求めたひずみ振幅Δεiと平均ひずみεa50、および第一算出工程で第二要素サイズに基づく低下定数α12.5および第一要素サイズに応じた補正パラメータｋを乗算して算出した低下率β50に基づいて、第二算出工程で亀裂発生寿命Ｎを算出する。なお、補正パラメータｋは、上記実施形態と同様に、ｋ＝２として設定した。
亀裂発生寿命Ｎの算出後は、累積線形被害則による削除要素決定を行う（ステップＳ１０４）。すなわち、累積線形被害則に基づいてはんだ接合部８の亀裂による剥離部分の要素を決定し、この要素を削除することによって剥離部分を解析対象から除外する。
As shown in FIG. 7, the FEM analysis was performed according to a flowchart for predicting the crack propagation life by crack propagation analysis. The analysis step (step S101), the first calculation step (step S102), and the second calculation step (step S103) are the same steps as the analysis step, the first calculation step, and the second calculation step in the above embodiment. That is, it was calculated by multiplying the strain amplitude Δε i and the average strain ε a50 obtained in the analysis step, the reduction constant α 12.5 based on the second element size in the first calculation step, and the correction parameter k corresponding to the first element size. Based on the decrease rate β 50 , the crack generation life N is calculated in the second calculation step. The correction parameter k is set as k = 2 as in the above embodiment.
After the crack generation life N is calculated, the deletion element is determined based on the cumulative linear damage law (step S104). That is, based on the cumulative linear damage law, an element of a peeling portion due to a crack in the solder joint portion 8 is determined, and the peeling portion is excluded from the analysis target by deleting this element.
続いて、計算終了するか否かを判断する（ステップＳ１０５）。冷熱サイクル数が３０００サイクルに到達していない場合は、ステップ１０１に戻り、当該ステップＳ１０１〜１０４の工程を繰り返す。冷熱サイクル数が３０００サイクルに到達すると、亀裂進展寿命を算出し（ステップ１０６）、亀裂進展解析を終了する。
なお、ＦＥＭ解析は、市販の構造解析ソフト（ＡＢＡＱＵＳ／ＳＴＡＮＤＡＲＤ ｖｅｒ６．７）を使用し、累積線形被害則による削除要素決定は、ｐｙｔｈｏｎプログラムにより行った。また、ＦＥＭ解析は、上記実施形態と同一の解析条件とし、サーマルショック試験の負荷条件に対応する解析と、同試験の無負荷条件に対応する解析とを行った。さらに、負荷条件に対応する解析では、第一算出工程で補正パラメータｋにより補正する場合と、補正パラメータｋを補正しない場合との２種類の解析を行った。
Subsequently, it is determined whether or not to end the calculation (step S105). When the number of cooling cycles has not reached 3000 cycles, the process returns to step 101 and the processes of steps S101 to S104 are repeated. When the number of cooling cycles reaches 3000, the crack propagation life is calculated (step 106), and the crack propagation analysis is terminated.
For the FEM analysis, commercially available structural analysis software (ABAQUS / STANDARD ver 6.7) was used, and deletion element determination based on the cumulative linear damage law was performed by the Python program. The FEM analysis was performed under the same analysis conditions as in the above embodiment, and an analysis corresponding to the load condition of the thermal shock test and an analysis corresponding to the no-load condition of the test were performed. Furthermore, in the analysis corresponding to the load condition, two types of analysis were performed, that is, when the correction parameter k is corrected in the first calculation step and when the correction parameter k is not corrected.
図８〜図１０は、上記試験による実測結果、および上記ＦＥＭ解析による解析結果を示すグラフである。各グラフ中の「ＭＡＸ」は実測結果の最大値、「ＭＩＮ」は実測結果の最小値、「ＡＶＥ」は実測結果の平均値をそれぞれ意味する。
図８は、負荷条件かつ第一算出工程で補正パラメータｋにより補正した場合における亀裂進展長さ（μｍ）と冷熱サイクル数との関係を示すグラフである。図８に示された結果から、ＦＥＭ解析による解析結果は、サイクル数の増加にともなって亀裂進展長さが増大しており、その変化量も試験による実測結果と近似していることから、高い解析精度を有していることがわかった。
8 to 10 are graphs showing the actual measurement result by the above test and the analysis result by the FEM analysis. In each graph, “MAX” means the maximum value of the actual measurement results, “MIN” means the minimum value of the actual measurement results, and “AVE” means the average value of the actual measurement results.
FIG. 8 is a graph showing the relationship between the crack growth length (μm) and the number of thermal cycles when the load parameter is corrected by the correction parameter k in the first calculation step. From the results shown in FIG. 8, the analysis result by the FEM analysis is high because the crack growth length increases with an increase in the number of cycles, and the amount of change approximates the actual measurement result by the test. It was found to have analysis accuracy.
図９は、無負荷条件かつ第一算出工程で補正パラメータｋにより補正した場合における亀裂進展長さ（μｍ）と冷熱サイクル数との関係を示すグラフである。図８および図９のグラフを比較すると、図８の負荷条件による解析結果の亀裂進展長さは、図９の無負荷条件による解析結果の亀裂進展長さよりも大きくなるという結果となった。したがって、ＦＥＭ解析により、実測結果と同様に複合負荷に起因して亀裂進展寿命が短くなることを再現できることがわかった。 FIG. 9 is a graph showing the relationship between the crack growth length (μm) and the number of cooling cycles when the load is corrected with the correction parameter k in the first calculation step. Comparing the graphs of FIG. 8 and FIG. 9, the crack growth length of the analysis result under the load condition of FIG. 8 is larger than the crack growth length of the analysis result under the no-load condition of FIG. 9. Therefore, it was found by FEM analysis that the crack growth life can be shortened due to the combined load as in the actual measurement result.
図１０は、負荷条件かつ第一算出工程で補正パラメータｋにより補正しなかった場合における亀裂進展長さ（μｍ）と冷熱サイクル数との関係を示すグラフである。図８および図１０のグラフを比較すると、両者は同じ負荷条件により解析したにも関わらず、図１０の解析結果による亀裂進展長さのほうが短くなるという結果となった。また、図９および図１０のグラフを比較すると、図１０の解析結果による亀裂進展長さは、負荷条件により解析しているにも関わらず、図９の無負荷条件により解析した亀裂進展長さとほぼ近似する結果となった。これらの結果から、負荷条件において補正パラメータｋにより補正しなかった場合の解析結果は、無負荷条件の解析結果と近似する結果となり、負荷条件の実測結果を再現できないことがわかった。換言すると、負荷条件において補正パラメータｋにより補正すれば、その解析結果が実測結果を再現できることがわかった。 FIG. 10 is a graph showing the relationship between the crack growth length (μm) and the number of cooling cycles when the load condition is not corrected by the correction parameter k in the first calculation step. Comparing the graphs of FIG. 8 and FIG. 10, although both were analyzed under the same load condition, the crack propagation length according to the analysis result of FIG. 10 was shorter. Further, comparing the graphs of FIG. 9 and FIG. 10, the crack growth length according to the analysis result of FIG. 10 is the same as the crack growth length analyzed under the no-load condition of FIG. The result was almost approximate. From these results, it was found that the analysis result when the correction parameter k was not corrected under the load condition approximated the analysis result under the no-load condition, and the actual measurement result under the load condition could not be reproduced. In other words, it was found that the analysis result can reproduce the actual measurement result if the correction is made with the correction parameter k under the load condition.
なお、今回開示された実施の形態はすべての点で例示であって制限的なものでないと考えられるべきである。本発明の範囲は、上記した意味ではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味、および範囲内でのすべての変更が含まれることが意図される。
例えば、上記実施形態では、上記式（２）において低下率βを算出する際に、補正パラメータｋを乗算しているが、これに限定されず、補正パラメータｋを加算、減算もしくは除算するものであってもよい。
また、平均ひずみから低下率を算出するための関係式において、比例定数である低下定数ｋを用いているが、平均ひずみεaと低下率βとが曲線的に変化する関係の場合には、この関係に基づく値を用いればよい。その際、補正パラメータｋは、一定値ではなく、関数により設定されるものであってもよい。
さらに、第一要素サイズを５０μｍ、第二要素サイズを１２．５μｍとしているが、第一要素サイズを１２．５μｍ、第二要素サイズを５０μｍとしてもよい。
また、第一要素サイズおよび第二要素サイズは、上記１２．５μｍおよび５０μｍに限定されるものではなく、第一要素サイズと第二要素サイズとが異なるサイズに設定されていればよい。
また、上記式（３）においてＭａｎｓｏｎ−Ｃｏｆｆｉｎ則に基づく予測式を用いているが、その他の予測式を用いることも可能である。
また、亀裂発生寿命予測装置１０は、車載コネクタ１のはんだ接合部８における亀裂発生寿命Ｎを予測するものであるが、その他の車載電子部品、または車載以外の電子部品のはんだ接合部における亀裂発生寿命を予測するものであってもよい。
The embodiment disclosed this time should be considered as illustrative in all points and not restrictive. The scope of the present invention is defined by the terms of the claims, rather than the meanings described above, and is intended to include any modifications within the scope and meaning equivalent to the terms of the claims.
For example, in the above embodiment, the correction parameter k is multiplied when calculating the reduction rate β in the above equation (2). However, the present invention is not limited to this, and the correction parameter k is added, subtracted, or divided. There may be.
In the relational expression for calculating the reduction rate from the average strain, the reduction constant k, which is a proportional constant, is used. In the case where the average strain ε a and the reduction rate β change in a curve, A value based on this relationship may be used. At this time, the correction parameter k may be set by a function instead of a constant value.
Furthermore, although the first element size is 50 μm and the second element size is 12.5 μm, the first element size may be 12.5 μm and the second element size may be 50 μm.
Further, the first element size and the second element size are not limited to the above 12.5 μm and 50 μm, and the first element size and the second element size may be set to different sizes.
Moreover, although the prediction formula based on the Manson-Coffin rule is used in the above formula (3), other prediction formulas can also be used.
Further, the crack generation life prediction device 10 predicts the crack generation life N in the solder joint portion 8 of the in-vehicle connector 1, but crack generation in the solder joint portion of other in-vehicle electronic components or electronic components other than the in-vehicle devices. The life may be predicted.
１ 車載コネクタ
５ 基板
８ はんだ接合部
１０ 亀裂発生寿命予測装置
１１ 解析部
１２ 第一算出部
１３ 第二算出部
ｋ 補正パラメータ
Ｎ 亀裂発生寿命
α12.5 第二の要素サイズに基づく低下定数
β50 低下率
Δεi ひずみ振幅
εa50 平均ひずみ
Ｓ１ 解析工程
Ｓ２ 第一算出工程
Ｓ３ 第二算出工程
DESCRIPTION OF SYMBOLS 1 In-vehicle connector 5 Board | substrate 8 Solder joint 10 Crack generation life prediction apparatus 11 Analysis part 12 First calculation part 13 Second calculation part k Correction parameter N Crack generation life α 12.5 Decrease constant based on second element size β 50 Reduction rate Δε i Strain amplitude ε a50 Average strain S1 Analysis step S2 First calculation step S3 Second calculation step
Claims (5)
所定の第一要素サイズでメッシュ分割した前記はんだ接合部の解析モデルに、前記荷重負荷を作用させた状態で冷熱サイクルを付与することによって、前記熱負荷に起因するひずみ振幅、および当該荷重負荷に起因する平均ひずみを有限要素法を用いた解析により求める解析部と、
前記第一要素サイズと異なるサイズの第二要素サイズにおいて平均ひずみから低下率を算出するための関係式を用いて、前記解析部で求めた平均ひずみから、前記荷重負荷に起因する前記亀裂発生寿命の低下率を算出する第一算出部と、
前記解析部で求めたひずみ振幅、および前記第一算出部で算出した低下率に基づいて、はんだ接合部の亀裂発生寿命を算出する第二算出部と、を備え、
前記第一算出部は、前記関係式を用いて算出された値を、前記第一要素サイズに応じて所定の補正パラメータにより補正して低下率を求める
ことを特徴とする亀裂発生寿命予測装置。 A crack generation life prediction device that predicts the crack generation life of a solder joint at which a thermal load due to a temperature change acts and a load load acts,
By applying a thermal cycle with the load applied to the analysis model of the solder joint divided into meshes with a predetermined first element size, the strain amplitude caused by the thermal load and the load load are applied. An analysis unit for determining the resulting average strain by analysis using the finite element method;
Using the relational expression for calculating the decrease rate from the average strain in the second element size having a size different from the first element size, the crack initiation life resulting from the load load from the average strain obtained by the analysis unit A first calculation unit for calculating a decrease rate of
Based on the strain amplitude obtained by the analysis unit and the decrease rate calculated by the first calculation unit, a second calculation unit that calculates the crack initiation life of the solder joint,
The crack generation life prediction apparatus characterized in that the first calculation unit corrects a value calculated using the relational expression with a predetermined correction parameter according to the first element size to obtain a decrease rate.
Ｎ＝１０００×（Δεi／０．０１）-1.24×（１−ｋ×αεa）
ここで、Ｎは亀裂発生寿命、Δεiはひずみ振幅、ｋは補正パラメータ、αは低下定数、εaは平均ひずみ、ｋ×αεaは低下率（ただし、０≦ｋ×αεa≦１）である。 The crack generation life prediction apparatus according to claim 1 or 2, wherein the second calculation unit calculates the crack generation life by the following formula based on the Manson-Coffin rule.
N = 1000 × (Δε i /0.01 ) −1.24 × (1−k × αε a )
Here, N is the crack life, Δε i is the strain amplitude, k is the correction parameter, α is the reduction constant, ε a is the average strain, and k × αε a is the reduction rate (where 0 ≦ k × αε a ≦ 1) It is.
所定の第一要素サイズでメッシュ分割した前記はんだ接合部の解析モデルに、前記荷重負荷を作用させた状態で冷熱サイクルを付与することによって、前記熱負荷に起因するひずみ振幅、および当該荷重負荷に起因する平均ひずみを有限要素法を用いた解析により求める解析工程と、
前記第一要素サイズと異なるサイズの第二要素サイズにおいて平均ひずみから低下率を算出するための関係式を用いて、前記解析部で求めた平均ひずみから、前記荷重負荷に起因する前記亀裂発生寿命の低下率を算出する第一算出工程と、
前記解析工程で求めたひずみ振幅、および前記第一算出工程で算出した低下率に基づいて、はんだ接合部の亀裂発生寿命を算出する第二算出工程と、を備え、
前記第一算出部は、前記関係式を用いて算出された値を、前記第一要素サイズに応じて所定の補正パラメータにより補正して低下率を求める
ことを特徴とする亀裂発生寿命予測方法。 A crack generation life prediction method for predicting a crack generation life of a soldered joint where a heat load acts due to a temperature change and a load load acts,
By applying a thermal cycle with the load applied to the analysis model of the solder joint divided into meshes with a predetermined first element size, the strain amplitude caused by the thermal load and the load load are applied. An analysis process for determining the resulting average strain by analysis using the finite element method;
Using the relational expression for calculating the decrease rate from the average strain in the second element size having a size different from the first element size, the crack initiation life resulting from the load load from the average strain obtained by the analysis unit A first calculation step of calculating a decrease rate of
Based on the strain amplitude obtained in the analysis step and the decrease rate calculated in the first calculation step, a second calculation step for calculating the crack generation life of the solder joint, and
The crack generation life prediction method, wherein the first calculation unit obtains a decrease rate by correcting a value calculated using the relational expression with a predetermined correction parameter according to the first element size.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2009207309A JP5356164B2 (en) | 2009-09-08 | 2009-09-08 | Crack generation life prediction apparatus and crack generation life prediction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2009207309A JP5356164B2 (en) | 2009-09-08 | 2009-09-08 | Crack generation life prediction apparatus and crack generation life prediction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2011058888A true JP2011058888A (en) | 2011-03-24 |
| JP5356164B2 JP5356164B2 (en) | 2013-12-04 |
Family
ID=43946702
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2009207309A Active JP5356164B2 (en) | 2009-09-08 | 2009-09-08 | Crack generation life prediction apparatus and crack generation life prediction method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5356164B2 (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109932251A (en) * | 2017-12-15 | 2019-06-25 | 天津大学 | Creep induction period prediction method considering load-independent constraint parameters under plastic condition |
| CN109932242A (en) * | 2017-12-15 | 2019-06-25 | 天津大学 | Creep induction period prediction method considering load-independent constraint parameters under steady-state creep condition |
| CN110967267A (en) * | 2019-11-25 | 2020-04-07 | 中国民用航空飞行学院 | Test method for judging fatigue crack initiation life |
| CN112179765A (en) * | 2020-06-09 | 2021-01-05 | 西北工业大学 | Fatigue life prediction method based on CT test piece |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10915679B2 (en) | 2015-02-26 | 2021-02-09 | Nippon Steel Corporation | Break prediction method, break prediction device, program, recording medium, and break discernment standard calculation method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006071406A (en) * | 2004-09-01 | 2006-03-16 | Toyota Motor Corp | Method and apparatus for estimating fatigue life of solder joint |
| JP2008002869A (en) * | 2006-06-21 | 2008-01-10 | Denso Corp | Method for estimating life of solder |
| JP2009174936A (en) * | 2008-01-23 | 2009-08-06 | Fujitsu Ltd | Crack growth evaluation apparatus and method |
-
2009
- 2009-09-08 JP JP2009207309A patent/JP5356164B2/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006071406A (en) * | 2004-09-01 | 2006-03-16 | Toyota Motor Corp | Method and apparatus for estimating fatigue life of solder joint |
| JP2008002869A (en) * | 2006-06-21 | 2008-01-10 | Denso Corp | Method for estimating life of solder |
| JP2009174936A (en) * | 2008-01-23 | 2009-08-06 | Fujitsu Ltd | Crack growth evaluation apparatus and method |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109932251A (en) * | 2017-12-15 | 2019-06-25 | 天津大学 | Creep induction period prediction method considering load-independent constraint parameters under plastic condition |
| CN109932242A (en) * | 2017-12-15 | 2019-06-25 | 天津大学 | Creep induction period prediction method considering load-independent constraint parameters under steady-state creep condition |
| CN109932251B (en) * | 2017-12-15 | 2021-04-09 | 天津大学 | Creep induction period prediction method considering load-independent constraint parameters under plastic condition |
| CN109932242B (en) * | 2017-12-15 | 2021-08-24 | 天津大学 | Creep induction period prediction method considering load-independent constraint parameters under steady-state creep condition |
| CN110967267A (en) * | 2019-11-25 | 2020-04-07 | 中国民用航空飞行学院 | Test method for judging fatigue crack initiation life |
| CN112179765A (en) * | 2020-06-09 | 2021-01-05 | 西北工业大学 | Fatigue life prediction method based on CT test piece |
Also Published As
| Publication number | Publication date |
|---|---|
| JP5356164B2 (en) | 2013-12-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Motalab et al. | Determination of Anand constants for SAC solders using stress-strain or creep data | |
| JP5356164B2 (en) | Crack generation life prediction apparatus and crack generation life prediction method | |
| KR101526313B1 (en) | Method for predicting fatigue life | |
| JP4855947B2 (en) | Crack growth evaluation device, crack growth evaluation method, and crack growth evaluation program | |
| Esfandyari et al. | Simulation, optimization and experimental verification of the over–pressure reflow soldering process | |
| JP2013058657A (en) | Electronic device life estimation method and electronic device design method using the same | |
| JP2015094717A (en) | Thermal fatigue life prediction device, thermal fatigue life prediction method, and program | |
| KR101550878B1 (en) | Method for fatigue life prediction of SMT solder joints | |
| JP5581860B2 (en) | Thermal fatigue life diagnosis method for solder joints | |
| JP4071703B2 (en) | Thermal fatigue life diagnosis method and apparatus for solder joints | |
| JP2791174B2 (en) | Electronic component solder connection life evaluation method | |
| JP2004237304A (en) | Method for estimating life time of solder joining | |
| JP2005026250A (en) | Method and device for predicting crack initiating timing in soldered junction | |
| Benabou et al. | Continuum damage approach for fatigue life prediction of viscoplastic solder joints | |
| Grams et al. | A geometry-independent lifetime modelling method for aluminum heavy wire bond joints | |
| JP4626577B2 (en) | Solder life prediction method | |
| Höhne et al. | Improved damage modeling for solder joints under combined vibration and temperature cycling loading | |
| Gomez et al. | Damage mechanics modeling of concurrent thermal and vibration loading on electronics packaging | |
| Libot et al. | Experimental SAC305 Shear Stress–Strain Hysteresis Loop Construction Using Hall's One-Dimensional Model Based on Strain Gages Measurements | |
| JP4179248B2 (en) | Method and apparatus for estimating fatigue life of solder joints | |
| JP4702140B2 (en) | Method for evaluating the life of solder joints | |
| JP2007255926A (en) | Joint evaluation method and joint evaluation device | |
| JP2004085397A (en) | Life estimation method for solder junction part | |
| Johansson et al. | A computational method for evaluating the damage in a solder joint of an electronic package subjected to thermal loads | |
| Lee et al. | Probabilistic Fatigue Life Estimation of Film Capacitor Assemblies Under Thermal Cycling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120802 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20120802 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130614 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130625 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130718 |
|
| 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: 20130820 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130828 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5356164 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |