Break prediction method, break prediction device, program, recording medium, and break discernment standard calculation method

Abstract

A break prediction method according to the present invention predicts a break of a joint portion of an object to be analyzed including a pair of members joined to each other by using a finite element method, and includes a first step of acquiring at least an element size of a base material portion, from among parameters set in an element model for the object to be analyzed; a second step of calculating, as a break discernment standard, a break limit moment defined by a function including the element size of the base material portion as a variable; and a third step of discerning whether the moment applied to the joint portion in a deformation analysis of the element model for the object to be analyzed exceeds the break limit moment, and outputting the result of the discernment as a break prediction result for the joint portion.

Claims

1. A break prediction method which predicts a break of a joint portion of an object to be analyzed including a pair of members joined to each other by using a finite element method, the method comprising: a first step of acquiring at least an element size of a base material portion, from among parameters set in an element model for the object to be analyzed; a second step of calculating, as a break discernment standard, a break limit moment based on at least the element size of the base material portion; a third step of discerning whether the moment applied to the joint portion in a deformation analysis of the element model for the object to be analyzed exceeds the break limit moment, and outputting a result of the discernment as a break prediction result for the joint portion, and a step of designing a structure of the joint portion of the object for preventing break at the time of collision based on the base material and the break prediction result for the joint portion, wherein in the first step, from among the parameters set in the element model, a tensile strength of the object to be analyzed and the element size of the base material portion are acquired; and wherein in the second step, the calculation of the break limit moment is changed depending on the tensile strength, wherein in a case where the tensile strength is equal to or greater than 980 MPa, in the second step, the break limit moment is calculated by using Equation (1),
Mf=Me.Math.F(Me,t,D,W,L,e)(1) where, Mf: break limit moment (N.Math.mm) Me: modified elastic bending moment (N.Math.mm) F (Me, t, D, W, L, e): correction terms t: thickness of the object to be analyzed (mm) D: nugget diameter in a case where a joint portion is a spot-welded portion (mm) W: effective width that a spot-welded portion receives a load (mm) L: arm length (mm) e: an element size of a base material portion (mm), and wherein in Equation (1), the modified elastic bending moment Me is defined by Equation (2), the correction terms F(Me, t, D, W, L, e) are defined by Equation (3),
Me=(el/L).Math.(E.Math.D.Math.t.sup.3)/12(2) where, el: total elongation of an object to be analyzed () E: Young's modulus of the object to be analyzed (MPa)
F(Me,t,D,W,L,e)=f(Me).Math.f(t).Math.f(D).Math.f(W).Math.f(L).Math.f(e)(3) where, f(Me): correction term of the modified elastic bending moment Me f(t): correction term of the thickness t f(D): correction term of the nugget diameter D f(W): correction term of the effective width W f(L): correction term of the arm length L f(e): correction term of the element size e of the base material portion.

2. The break prediction method according to claim 1, wherein f(Me), f(t), f(D), f(W), f(L), and f(e) are equations which are determined from multiple regression such that an error between a moment applied to a spot-welded portion, obtained by loading a breaking load which is confirmed by experiments to a test piece edge and the break limit moment Mf calculated using Equation (1) is minimized, in a finite element method model generated using an L-shaped joint type test piece.

3. The break prediction method according to claim 1, wherein in a case where the tensile strength is less than 980 MPa, in the second step, the break limit moment is calculated using Equation (4),
Mf=Mp.Math.F(Mp,t,D,W,el,e)(4) where, Mf: break limit moment (N.Math.mm) Mp: total plastic bending moment (N.Math.mm) F(Mp, t, D, W, el, e): correction terms t: thickness of the object to be analyzed (mm) D: nugget diameter in a case where a joint portion is a spot-welded portion (mm) W: effective width that a spot-welded portion receives a load (mm) el: total elongation of the object to be analyzed () e: an element size of a base material portion (mm).

4. The break prediction method according to claim 3, wherein in Equation (4), the total plastic bending moment Mp is defined by Equation (5), and the correction terms F(Mp, t, D, W, el, e) are defined by Equation (6),
Mp=(TS.Math.W.Math.t.sup.2)/4(5) where, TS: tensile strength of the object to be analyzed (MPa)
F(Mp,t,D,W,el,e)=f(Mp).Math.f(t).Math.f(D).Math.f(W).Math.f(el).Math.f(e)(6) where, f(Mp): correction term of the total plastic bending moment Mp f(t): correction term of the thickness t f(D): correction term of the nugget diameter D f(W): correction term of the effective width W f(el): correction term of total elongation el f(e): correction term of the element size e of the base material portion.

5. The break prediction method according to claim 4, wherein f(Mp), f(t), f(D), f(W), f(el), and f(e) are equations which are determined from multiple regression such that an error between a moment applied to a spot-welded portion, obtained by loading a breaking load which is confirmed by experiments to a test piece edge and the break limit moment Mf calculated using Equation (4) is minimized, in a finite element method model generated using an L-shaped joint type test piece.

6. The break prediction method according to claim 1, wherein in the third step, in a case where a relationship between the moment M1 applied to the joint portion in a deformation analysis of an element model of the object to be analyzed and the break limit moment Mf satisfies Equation (9), a result indicating that there is a break is output as the break prediction result,
M1/Mf1(9).

7. The break prediction method according to claim 1, wherein in the third step, in a case where an axial force applied to the joint portion in a deformation analysis of an element model of the object to be analyzed is a compressed axial force, a result indicating that there is no break is forcibly output as the break prediction result.

8. A break prediction device which predicts a break of a joint portion of an object to be analyzed including a pair of members joined to each other by using a finite element method, the break prediction device comprising: a parameter acquisition processing circuitry to acquire at least an element size of a base material portion, from among parameters set in an element model for the object to be analyzed; a storage that stores a break limit moment that is based on at least the element size of the base material portion; a break discernment standard calculation processing circuitry to calculate the break limit moment as a break discernment standard by using the element size of the base material portion which is acquired by the parameter acquisition processing circuitry; and a break discernment processing circuitry to discern whether the moment applied to the joint portion in a deformation analysis of the element model for the object to be analyzed exceeds the break limit moment, output a result of the discernment as a break prediction result for the joint portion, and a structure of the joint portion of the object for preventing break at the time of collision based on the base material and the break prediction result for the joint portion, wherein the parameter acquisition processing circuitry acquires a tensile strength of the object to be analyzed and the element size of the base material portion, from among the parameters set in the element model, wherein the storage stores a plurality of functions corresponding to the tensile strength, and wherein the break discernment standard calculation processing circuitry reads the function corresponding to the tensile strength acquired by the parameter acquisition processing circuitry from the storage and calculates the break limit moment, wherein in a case where the tensile strength is equal to or greater than 980 MPa, the break discernment standard calculation processing circuitry reads Equation (1) from the storage and calculates the break limit moment,
Mf=Me.Math.F(Me,t,D,W,L,e)(1) where, Mf: break limit moment (N.Math.mm) Me: modified elastic bending moment (N.Math.mm) F (Me, t, D, W, L, e): correction terms t: thickness of the object to be analyzed (mm) D: nugget diameter in a case where a joint portion is a spot-welded portion (mm) W: effective width that a spot-welded portion receives a load (mm) L: arm length (mm) e: an element size of a base material portion (mm), and wherein in Equation (1), the modified elastic bending moment Me is defined by Equation (2), the correction terms F(Me, t, D, W, L, e) are defined by Equation (3),
Me=(el/L).Math.(E.Math.D.Math.t.sup.3)/12(2) where, el: total elongation of an object to be analyzed () E: Young's modulus of the object to be analyzed (MPa)
F(Me,t,D,W,L,e)=f(Me).Math.f(t).Math.f(D).Math.f(W).Math.f(L).Math.f(e)(3) where, f(Me): correction term of the modified elastic bending moment Me f(t): correction term of the thickness t f(D): correction term of the nugget diameter D f(W): correction term of the effective width W f(L): correction term of the arm length L f(e): correction term of the element size e of the base material portion.

9. The break prediction device according to claim 8, wherein f(Me), f(t), f(D), f(W), f(L), and f(e) are equations which are determined from multiple regression such that an error between a moment applied to a spot-welded portion, obtained by loading a breaking load which is confirmed by experiments to a test piece edge and the break limit moment Mf calculated using Equation (1) is minimized, in a finite element method model generated using an L-shaped joint type test piece.

10. The break prediction device according to claim 8, wherein in a case where the tensile strength is less than 980 MPa, the break discernment standard calculation processing circuitry reads Equation (4) from the storage and calculates the break limit moment,
Mf=Mp.Math.F(Mp,t,D,W,el,e)(4) where, Mf: break limit moment (N.Math.mm) Mp: total plastic bending moment (N.Math.mm) F(Mp, t, D, W, el, e): correction terms t: thickness of the object to be analyzed (mm) D: nugget diameter in a case where a joint portion is a spot-welded portion (mm) W: effective width that a spot-welded portion receives a load (mm) el: total elongation of the object to be analyzed () e: an element size of a base material portion (mm).

11. The break prediction device according to claim 10, wherein in Equation (4), the total plastic bending moment Mp is defined by Equation (5), and the correction terms F(Mp, t, D, W, el, e) are defined by Equation (6),
Mp=(TS.Math.W.Math.t.sup.2)/4(5) where, TS: tensile strength of the object to be analyzed (MPa)
F(Mp,t,D,W,el,e)=f(Mp).Math.f(t).Math.f(D).Math.f(W).Math.f(el).Math.f(e)(6) where, f(Mp): correction term of the total plastic bending moment Mp f(t): correction term of the thickness t f(D): correction term of the nugget diameter D f(W): correction term of the effective width W f(el): correction term of total elongation el f(e): correction term of the element size e of the base material portion.

12. The break prediction device according to claim 11, wherein f(Mp), f(t), f(D), f(W), f(el), and f(e) are equations which are determined from multiple regression such that an error between a moment applied to a spot-welded portion, obtained by loading a breaking load which is confirmed by experiments to a test piece edge and the break limit moment Mf calculated using Equation (4) is minimized, in a finite element method model generated using an L-shaped joint type test piece.

13. The break prediction device according to claim 8, wherein in a case where a relationship between the moment M1 applied to the joint portion in a deformation analysis of an element model of the object to be analyzed and the break limit moment Mf satisfies Equation (9), the break discernment processing circuitry outputs a result indicating that there is a break as the break prediction result,
M1/Mf1(9).

14. The break prediction device according to claim 8, wherein in a case where an axial force applied to the joint portion in a deformation analysis of an element model of the object to be analyzed is a compressed axial force, the break discernment processing circuitry forcibly outputs a result indicating that there is no break as the break prediction result.

15. A non-transitory computer readable medium having a program causing a computer to execute a process of predicting a break of a joint portion of an object to be analyzed including a pair of members joined to each other by using a finite element method, the process comprising: a first process of acquiring at least an element size of a base material portion, from among parameters set in an element model for the object to be analyzed; a second process of calculating, as a break discernment standard, a break limit moment based on at least the element size of the base material portion; a third process of discerning whether the moment applied to the joint portion in a deformation analysis of the element model for the object to be analyzed exceeds the break limit moment, and outputting a result of the discernment as a break prediction result for the joint portion, and a process of designing a structure of the joint portion of the object for preventing break at the time of collision based on the base material and the break prediction result for the joint portion, wherein in the first process, the computer is caused to execute a process of acquiring a tensile strength of the object to be analyzed and the element size of the base material portion from among the parameters set in the element model, and wherein in the second process, the computer is caused to execute a process of changing the calculation of the break limit moment depending on the tensile strength, wherein in a case where the tensile strength is equal to or greater than 980 MPa, in the second process, the computer is caused to execute a process of calculating the break limit moment by using Equation (1),
Mf=Me.Math.F(Me,t,D,W,L,e)(1) where, Mf: break limit moment (N.Math.mm) Me: modified elastic bending moment (N.Math.mm) F (Me, t, D, W, L, e): correction terms t: thickness of the object to be analyzed (mm) D: nugget diameter in a case where a joint portion is a spot-welded portion (mm) W: effective width that a spot-welded portion receives a load (mm) L: arm length (mm) e: an element size of a base material portion (mm), and wherein in Equation (1), the modified elastic bending moment Me is defined by Equation (2), the correction terms F(Me, t, D, W, L, e) are defined by Equation (3),
Me=(el/L).Math.(E.Math.D.Math.t.sup.3)/12(2) where, el: total elongation of an object to be analyzed () E: Young's modulus of the object to be analyzed (MPa)
F(Me,t,D,W,L,e)=f(Me).Math.f(t).Math.f(D).Math.f(W).Math.f(L).Math.f(e)(3) where, f(Me): correction term of the modified elastic bending moment Me f(t): correction term of the thickness t f(D): correction term of the nugget diameter D f(W): correction term of the effective width W f(L): correction term of the arm length L f(e): correction term of the element size e of the base material portion.

16. The non-transitory computer readable medium having the program according to claim 15, wherein f(Me), f(t), f(D), f(W), f(L), and f(e) are equations which are determined from multiple regression such that an error between a moment applied to a spot-welded portion, obtained by loading a breaking load which is confirmed by experiments to a test piece edge and the break limit moment Mf calculated using Equation (1) is minimized, in a finite element method model generated using an L-shaped joint type test piece.

17. The non-transitory computer readable medium having the program according to claim 15, wherein in a case where the tensile strength is less than 980 MPa, in the second process, the computer is caused to execute a process of calculating the break limit moment by using Equation (4),
Mf=Mp.Math.F(Mp,t,D,W,el,e)(4) where, Mf: break limit moment (N.Math.mm) Mp: total plastic bending moment (N.Math.mm) F(Mp, t, D, W, el, e): correction terms t: thickness of the object to be analyzed (mm) D: nugget diameter in a case where a joint portion is a spot-welded portion (mm) W: effective width that a spot-welded portion receives a load (mm) el: total elongation of the object to be analyzed (c) e: an element size of a base material portion (mm).

18. The non-transitory computer readable medium having the program according to claim 17, wherein in Equation (4), the total plastic bending moment Mp is defined by Equation (5), and the correction terms F(Mp, t, D, W, el, e) are defined by Equation (6),
Mp=(TS.Math.W.Math.t.sup.2)/4(5) where, TS: tensile strength of the object to be analyzed (MPa)
F(Mp,t,D,W,el,e)=f(Mp).Math.f(t).Math.f(D).Math.f(W).Math.f(el).Math.f(e)(6) where, f(Mp): correction term of the total plastic bending moment Mp f(t): correction term of the thickness t f(D): correction term of the nugget diameter D f(W): correction term of the effective width W f(el): correction term of total elongation el f(e): correction term of the element size e of the base material portion.

19. The non-transitory computer readable medium having the program according to claim 18, wherein f(Mp), f(t), f(D), f(W), f(el), and f(e) are equations which are determined from multiple regression such that an error between a moment applied to a spot-welded portion, obtained by loading a breaking load which is confirmed by experiments to a test piece edge and the break limit moment Mf calculated using Equation (4) is minimized, in a finite element method model generated using an L-shaped joint type test piece.

20. The non-transitory computer readable medium having the program according to claim 15, wherein in the third process, in a case where a relationship between the moment M1 applied to the joint portion in a deformation analysis of an element model of the object to be analyzed and the break limit moment Mf satisfies Equation (9), the computer is caused to execute a process of outputting a result indicating that there is a break as the break prediction result,
M1/Mf1(9).

21. The non-transitory computer readable medium having the program according to claim 15, wherein in the third process, in a case where an axial force applied to the joint portion in a deformation analysis of an element model of the object to be analyzed is a compressed axial force, the computer is caused to execute a process of forcibly outputting a result indicating that there is no break as the break prediction result.

22. A break discernment standard calculation method which calculates a break discernment standard used when predicting a break of a joint portion of an object to be analyzed including a pair of members joined to each other by using a finite element method, the method comprising: a first step of acquiring at least an element size of a base material portion, from among parameters set in an element model for the object to be analyzed; a second step of calculating, as a break discernment standard, a break limit moment based on at least the element size of the base material portion; and a step of designing a structure of the joint portion of the object for preventing break at the time of collision based on the base material and the break prediction result for the joint portion, wherein in the first step, from among the parameters set in the element model, a tensile strength of the object to be analyzed and the element size of the base material portion are acquired; and wherein in the second step, the calculation of the break limit moment is changed depending on the tensile strength, wherein in a case where the tensile strength is equal to or greater than 980 MPa, in the second step, the break limit moment is calculated by using Equation (1),
Mf=Me.Math.F(Me,t,D,W,L,e)(1) where, Mf: break limit moment (N.Math.mm) Me: modified elastic bending moment (N.Math.mm) F (Me, t, D, W, L, e): correction terms t: thickness of the object to be analyzed (mm) D: nugget diameter in a case where a joint portion is a spot-welded portion (mm) W: effective width that a spot-welded portion receives a load (mm) L: arm length (mm) e: an element size of a base material portion (mm) and, wherein in Equation (1), the modified elastic bending moment Me is defined by Equation (2), the correction terms F(Me, t, D, W, L, e) are defined by Equation (3),
Me=(el/L).Math.(E.Math.D.Math.t.sup.3)/12(2) where, el: total elongation of an object to be analyzed () E: Young's modulus of the object to be analyzed (MPa)
F(Me,t,D,W,L,e)=f(Me).Math.f(t).Math.f(D).Math.f(W).Math.f(L).Math.f(e)(3) where, f(Me): correction term of the modified elastic bending moment Me f(t): correction term of the thickness t f(D): correction term of the nugget diameter D f(W): correction term of the effective width W f(L): correction term of the arm length L f(e): correction term of the element size e of the base material portion.

23. The break prediction method according to claim 22, wherein f(Me), f(t), f(D), f(W), f(L), and f(e) are equations which are determined from multiple regression such that an error between a moment applied to a spot-welded portion, obtained by loading a breaking load which is confirmed by experiments to a test piece edge and the break limit moment Mf calculated using Equation (1) is minimized, in a finite element method model generated using an L-shaped joint type test piece.

24. The break prediction method according to claim 22, wherein in a case where the tensile strength is less than 980 MPa, in the second step, the break limit moment is calculated using Equation (4),
Mf=Mp.Math.F(Mp,t,D,W,el,e)(4) where, Mf: break limit moment (N.Math.mm) Mp: total plastic bending moment (N.Math.mm) F(Mp, t, D, W, el, e): correction terms t: thickness of the object to be analyzed (mm) D: nugget diameter in a case where a joint portion is a spot-welded portion (mm) W: effective width that a spot-welded portion receives a load (mm) el: total elongation of the object to be analyzed () e: an element size of a base material portion (mm).

25. The break prediction method according to claim 24, wherein in Equation (4), the total plastic bending moment Mp is defined by Equation (5), and the correction terms F(Mp, t, D, W, el, e) are defined by Equation (6),
Mp=(TS.Math.W.Math.t.sup.2)/4(5) where, TS: tensile strength of the object to be analyzed (MPa)
F(Mp,t,D,W,el,e)=f(Mp).Math.f(t).Math.f(D).Math.f(W).Math.f(el).Math.f(e)(6) where, f(Mp): correction term of the total plastic bending moment Mp f(t): correction term of the thickness t f(D): correction term of the nugget diameter D f(W): correction term of the effective width W f(el): correction term of total elongation el f(e): correction term of the element size e of the base material portion.

26. The break prediction method according to claim 25, wherein f(Mp), f(t), f(D), f(W), f(el), and f(e) are equations which are determined from multiple regression such that an error between a moment applied to a spot-welded portion, obtained by loading a breaking load which is confirmed by experiments to a test piece edge and the break limit moment Mf calculated using Equation (4) is minimized, in a finite element method model generated using an L-shaped joint type test piece.

27. The break prediction method according to claim 22, wherein in the third step, in a case where a relationship between the moment M1 applied to the joint portion in a deformation analysis of an element model of the object to be analyzed and the break limit moment Mf satisfies Equation (9), a result indicating that there is a break is output as the break prediction result,
M1/Mf1(9).

28. The break prediction method according to claim 22, wherein in the third step, in a case where an axial force applied to the joint portion in a deformation analysis of an element model of the object to be analyzed is a compressed axial force, a result indicating that there is no break is forcibly output as the break prediction result.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a schematic diagram showing a schematic configuration of a break prediction device according to an embodiment of the present invention.

(2) FIG. 2 is a flowchart showing a break prediction method according to an embodiment of the present invention in a step order.

(3) FIG. 3A is a diagram showing a result of investigating a relationship between a predicted breaking load obtained by the break prediction method according to the embodiment of the present invention and an experimental breaking load obtained by experiment, in a case where the tensile strength of an L-shaped joint which is an object to be analyzed is equal to or greater than 980 MPa.

(4) FIG. 3B is a diagram showing a result of investigating a relationship between a predicted breaking load obtained by the break prediction method according to the embodiment of the present invention and an experimental breaking load obtained by experiment, in a case where the tensile strength of an L-shaped joint which is an object to be analyzed is less than 980 MPa.

(5) FIG. 4 is a schematic diagram showing an internal structure of a personal user terminal device.

EMBODIMENTS OF THE INVENTION

(6) Hereinafter, an embodiment of the present invention will be described in detail with reference to the drawings.

(7) In the present embodiment, a break (in particular, a break caused by the application of a moment to the joint portion) of a joint portion of an object to be analyzed including a pair of members joined to each other is predicted by using a finite element method. In the element model, the base material portion is a shell element or a solid element, and the joint portion is a beam element, a solid element, or a shell element. In the present embodiment, an L-shaped joint type test piece including a pair of L-shaped steel plates joined by spot welding is exemplified as the object to be analyzed, and a case of predicting a break of a spot-welded portion of the test piece will be described. In order to stably obtain high break prediction accuracy without depending on the element (mesh) size using FEM, the break prediction method according to the present embodiment includes a first step of acquiring at least an element size of a base material portion, from among parameters set in an element model of an L-shaped joint type test piece; a second step of calculating, as a break discernment standard, a break limit moment defined by a function including the element size of the base material portion as a variable; and a third step of discerning whether the moment applied to the spot-welded portion in a deformation analysis of the element model for the object to be analyzed exceeds the break limit moment, and outputting the result of the discernment as a break prediction result for the spot-welded portion.

(8) In order to increase the break prediction accuracy, it is preferable that in the first step, the tensile strength of the base material portion of the L-shaped joint type test piece is acquired together with the element size of the base material portion, among the parameters set in the element model of an L-shaped joint type test piece, and in the second step, the function used to calculate the break limit moment varies depending on the tensile strength.

(9) For example, in a case where the tensile strength of the L-shaped joint type test piece is 980 MPa or more (in a case where the base material of the L-shaped joint type test piece is ultra-high-tensile material), in the second step, the break limit moment is calculated by using Equation (1) as the function.
Mf=Me.Math.F(Me,t,D,W,L,e)(1)

(10) where, Mf: break limit moment (N mm) Me: modified elastic bending moment (N mm)

(11) F(Me, t, D, W, L, e): correction terms t: thickness of a base material of an L-shaped joint type test piece (mm) D: nugget diameter of a spot-welded portion (mm) W: effective width that a spot-welded portion receives a load (mm) L: arm length (mm) e: an element size of a base material portion (mm)

(12) In Equation (1), the modified elastic bending moment Me may be defined by Equation (2), and the correction terms F(Me, t, D, W, L, e) may be defined by Equation (3).
Me=(el/L).Math.(E.Math.D.Math.t.sup.3)/12(2)

(13) where, el: total elongation of the base material of the L-shaped joint type test piece () E: Young's modulus of a base material of the L-shaped joint type test piece (MPa)

(14) F(Me, t, D, W, L, e)
=F(Me).Math.f(t).Math.f(D).Math.f(W).Math.f(L).Math.f(e)(3)

(15) where, f(Me): correction term of the modified elastic bending moment Me

(16) f(t): correction term of the thickness t

(17) f(D): correction term of the nugget diameter D

(18) f(W): correction term of the effective width W

(19) f(L): correction term of the arm length L

(20) f(e): correction term of the element size e of the base material portion

(21) In a case where the tensile strength of the L-shaped joint type test piece is 980 MPa or more, the break limit moment is corrected based on the modified elastic bending moment, but the modified elastic bending moment will be described below.

(22) In a case where the tensile strength of the L-shaped joint type test piece is 980 MPa or more, plastic deformation is observed locally in the vicinity of the spot-welded portion, but it is found that the joint is almost in the elastic deformed state as a whole, and a spot-welded portion is broken. From this, it is found that in materials having a tensile strength of 980 MPa or more, as a result of examining the relationship between a break limit moment and an elastic bending moment, it becomes clear that there is a certain correlation. Although the general equation of an elastic bending moment is (1/).Math.(E.Math.W.Math.t.sup.3)/12, it is difficult to obtain the curvature (1/) of the base material of the spot-welded portion at the time of break occurrence as a parameter from an element model of an object to be analyzed. Therefore, in order to replace the curvature with the value of an obtainable parameter, as a result of investigating the relation of each parameter, it is found out that the product of the curvature and the arm length is proportional to the total elongation of base material. Therefore, the curvature (1/) is replaced by the value (el/L) obtained by dividing the total elongation (el) of the base material by the arm length (L). Further, in a case where the tensile strength is 980 MPa or more, since the effect of the test piece width on the break limit moment is small, and the distribution of strain in an FEM model of the L-shaped joint type test piece is concentrated in a spot-welded portion, the width W that receives a load matches the nugget diameter D, and W in the general equation of an elastic bending moment is replaced with the nugget diameter D. The modified elastic bending moment generated in this manner is corrected to include correction term of the element size, and thus the break limit moment is calculated.

(23) On the other hand, for example, in a case where the tensile strength of the L-shaped joint type test piece is less than 980 MPa (in a case where the the tensile strength of a base material (steel plate) of the L-shaped joint type test piece is 780 MPa or less, as strength classification), in the second step, the break limit moment is calculated using Equation (4) as the function.
Mf=Mp.Math.F(Mp,t,D,W,el,e)(4)

(24) where, Mf: break limit moment (N.Math.mm) Mp: total plastic bending moment (N.Math.mm)

(25) F(Mp, t, D, W, el, e): correction terms t: thickness of a base material of an L-shaped joint type test piece (mm) D: nugget diameter of a spot-welded portion (mm) W: effective width that a spot-welded portion receives a load (mm) el: total elongation of the base material of the L-shaped joint type test piece () e: an element size of a base material portion (mm)

(26) In Equation (4), the total plastic bending moment Mp is defined by Equation (5), and the correction terms F(Mp, t, D, W, el, e) are defined by Equation (6).
Mp=(TS.Math.W.Math.t.sup.2)/4(5)

(27) where, TS: tensile strength of the base material of the L-shaped joint type test piece (MPa)

(28) F(Mp, t, D, W, el, e)
=f(Mp).Math.f(t).Math.f(D).Math.f(W).Math.f(el).Math.f(e)(6)

(29) where, f(Mp): correction term of the total plastic bending moment Mp

(30) f(t): correction term of the thickness t

(31) f(D): correction term of the nugget diameter D

(32) f(W): correction term of the effective width W

(33) f(el): correction term of total elongation el

(34) f(e): correction term of the element size e of the base material portion

(35) In a case where the test piece having a tensile strength less than 980 MPa is an object to be analyzed, it has been found that there is a certain correlation between the break limit moment and the total plastic bending moment, and the break limit moment is calculated by correcting to include a correction term of an element size based on the total plastic bending moment.

(36) Hereinafter, a method of determining each correction term of Equation (3) and Equation (6) will be described.

(37) First, in an FEM model generated using an L-shaped joint type test piece, the breaking load (maximum load) confirmed in an experiment is loaded on the test piece edge (chuck part), and the moment applied to a spot-welded portion is acquired as the break limit moment. The acquisition of the break limit moment is performed by various kinds of steel, a thickness, an L-shaped joint shape, a nugget diameter, the element size of a base material portion, or the like.

(38) Subsequently, each correction term is determined from each multiple regression obtained such that an error between the break limit moment obtained as described above and the break limit moment calculated by Equation (3) or (6) is minimized. Specifically, each correction term is defined as follows. The equation form of the correction term is not particularly limited, and for example, a quadratic equation may be used instead of the linear equation.

(39) 1. In a case where a test piece having a tensile strength of 980 MPa or more is an object to be analyzed
f(Me)=(A1/Me)+A2
f(t)=B1.Math.t+B2
f(D)=C1.Math.D+C2
f(W)=D1.Math.W+D2
f(L)=E1.Math.L+E2
f(e)=F1.Math.e+F2(7)

(40) 2. In a case where a test piece having a tensile strength less than 980 MPa is an object to be analyzed
f(Mp)=(a1/Mp)+a2
f(t)=b1.Math.t+b2
f(D)=c1.Math.D+c2
f(W)=d1.Math.W+d2
f(el)=el.Math.(e1)+e2
f(e)=f1.Math.e+f2(8)

(41) As described above, respective coefficients A1, A2, B1, B2, C1, C2, D1, D2, F1 and F2 in Equation (7) and respective coefficients a1, a2, b1, b2, c1, c2, d1, d2, e1, e2, f1 and f2 of Equation (8) are determined such that an error between the break limit moment obtained in the FEM model and the break limit moment calculated by Equation (3) or (6) is minimized, and the break limit moment Mf of Equation (3) or (6) is obtained.

(42) In connection with the present embodiment, in an FEM model generated using an L-shaped joint type test piece, the breaking load (maximum load) confirmed in an experiment is loaded on the test piece edge (chuck part), the moment applied to a spot-welded portion is acquired, and the value is directly used as the break limit moment. However, there are numerous combinations of conditions (the kind of steel, the thickness, the shape, the nugget diameter, the element size of a base material portion, or the like) of an L-shaped joint for which the user actually tries to predict a break by simulation. It is impossible to perform experiment with all these combinations and obtain a breaking load (maximum load). Therefore, in the present embodiment, a break is predicted by using the above equation (3) or (6).

(43) FIG. 1 is a schematic diagram showing a schematic configuration of a break prediction device according to the present embodiment. FIG. 2 is a flowchart showing a break prediction method according to the present embodiment in a step order.

(44) As shown in FIG. 1, the break prediction device according to the present embodiment includes a parameter acquisition unit 1, a break discernment standard calculation unit 2, a break discernment unit 3, and a storage unit 4.

(45) The break prediction device according to the present embodiment can be realized by a computer such as a personal computer. The storage unit 4 is a nonvolatile storage device such as a flash memory, a hard disk, or a read only memory (ROM) provided in such a computer.

(46) The parameter acquisition unit 1, the break discernment standard calculation unit 2, and the break discernment unit 3 are functions realized when an arithmetic processing unit (not shown in FIG. 1) such as a central processing unit (CPU) provided in the computer operates according to the program stored in the storage unit 4.

(47) Here, the program is application software constructed with a computer-readable machine language to realize the above break prediction method according to the present embodiment by a computer. The program can be stored in the storage unit 4 by downloading the program from a portable recording medium such as a universal serial bus (USB) memory or a CD-ROM to the computer.

(48) For example, in a case of predicting break of a spot-welded portion in collision FEM analysis of a hat type member (L-shaped joint type members) joined by spot welding, LS-DYNA which is general purpose collision analysis software is used as a main routine program, and it is possible to link the program to LS-DYNA as a subroutine program of LS-DYNA. That is, the arithmetic processing unit of a computer (break prediction device) operates in accordance with the main routine program LS-DYNA to perform a collision deformation analysis process of a hat type member which is an object to be analyzed, and operates in accordance with the program which is a subroutine program to perform a break prediction process of determining whether or not a break occurs in the spot-welded portion while interlocking with the collision deformation analysis process.

(49) Therefore, not only the program but also LS-DYNA which is the main routine program are stored in the storage unit 4. In addition, LS-DYNA may be stored in a storage unit different from the storage unit 4. Further, other programs necessary for the operation of a computer, such as the operating system (OS) program, may be stored in the storage unit 4, or may be stored in a storage unit different from the storage unit 4.

(50) Further, each correction term of Equation (7) where each of coefficients A1, A2, B1, B2, C1, C2, D1, D2, F1 and F2 is determined, and each correction term of Equation (8) where each of coefficients a1, a2, b1, b2, c1, c2, d1, d2, e1, e2, f1 and f2 is determined, and Equations (1) to (6) are stored in the storage unit 4. These data items can be stored in the storage unit 4 by downloading these data items together with the program from a portable recording medium such as a USB memory or a CD-ROM to the computer.

(51) Hereinafter, with reference to FIG. 2, the break prediction method (functions of the parameter acquisition unit 1, the break discernment standard calculation unit 2, and the break discernment unit 3) realized by the arithmetic processing unit of the break prediction device (computer) operating in accordance with the program will be described.

(52) As shown in FIG. 2, the parameter acquisition unit 1 acquires a tensile strength TS, a Young's modulus E, an effective width W, a thickness t, a nugget diameter D, an arm length L, a total elongation el, and an element size e of the base material portion, from among parameters set in an element model of a hat type member (step S1: first process).

(53) As described above, the arithmetic processing unit of the break prediction device executes a collision deformation analysis process of a hat type member based on LS-DYNA and a break prediction process based on the program in parallel. In order to execute the collision deformation analysis process of a hat type member based on LS-DYNA, it is necessary to prepare an element model of the hat type member in advance, and therefore it is necessary to set various parameters. Therefore, the parameter acquisition unit 1 can easily acquire the tensile strength TS, the Young's modulus E, the effective width W, the thickness t, the nugget diameter D, the arm length L, the total elongation el, and the element size e of the base material portion, from among various parameters previously set to generate an element model of a hat type member.

(54) In addition, these parameters may be data input by the input device (not shown in FIG. 1) provided in the break prediction device, at the start of the execution of the break prediction process based on the program.

(55) Further, in step S1, the average element size around the spot-welded portion may be acquired as the element size e of the base material portion, by automatically reading various parameters of the element model of a hat type member which is the object to be analyzed from the input file for collision analysis, using input auxiliary software, and searching for the elements of a base material portion connected to the spot-welded portion.

(56) Next, the break discernment standard calculation unit 2 calculates the break limit moment Mf as a break discernment standard (step S2: second process).

(57) Specifically, in a case where the tensile strength TS acquired in step S1 is 980 MPa or more, the break discernment standard calculation unit 2 reads Equations (1) to (3) and each correction term of Equation (7) where each of coefficients A1, A2, B1, B2, C1, C2, D1, D2, F1 and F2 is determined, from the storage unit 4, and calculates the break limit moment Mf, by substituting the Young's modulus E, the effective width W, the thickness t, the nugget diameter D, the arm length L, the total elongation el, and the element size e of the base material portion into each equation, among the parameters acquired in step S1.

(58) On the other hand, in a case where the tensile strength TS acquired in step S1 is less than 980 MPa, the break discernment standard calculation unit 2 reads Equations (4) to (6) and each correction term of Equation (8) where each of coefficients a1, a2, b1, b2, c1, c2, d1, d2, e1, e2, f1 and f2 is determined, from the storage unit 4, and calculates the break limit moment Mf, by substituting the tensile strength TS, the effective width W, the thickness t, the nugget diameter D, the total elongation el, and the element size e of the base material portion into each equation, among the parameters acquired in step S1.

(59) Next, the break discernment unit 3 determines whether or not the bending moment applied to the spot-welded portion in the collision deformation analysis of the element model of the hat type member exceeds the break limit moment Mf obtained in step S2, and outputs the determination result as the break prediction result of the spot-welded portion (step S3: third process).

(60) Specifically, in a case where the bending moment M1 applied to the spot-welded portion in the collision deformation analysis of the element model of the hat type member and the break limit moment Mf satisfies Equation (9), the break discernment unit 3 outputs a result indicating that there is a break, as the break prediction result.
M1/Mf1(9)

(61) The bending moment M1 applied to the spot-welded portion in the collision deformation analysis of an element model of a hat type member can be obtained from the result of the collision deformation analysis process of a hat type member based on LS-DYNA.

(62) In a case where the break prediction result output from the break discernment unit 3 indicates that there is a break, the arithmetic processing unit of the break prediction device deletes the spot-welded portion included in the element model of a hat type member, in the collision deformation analysis process of a hat type member based on the LS-DYNA, so as to inform the user that a break has occurred in the spot-welded portion.

(63) On the other hand, in a case where the break prediction result output from the break discernment unit 3 indicates that there is no break, the arithmetic processing unit of the break prediction device leaves the spot-welded portion included in the element model of the hat type member, in the collision deformation analysis process of a hat type member based on the LS-DYNA, so as to inform the user that there is no break in the spot-welded portion.

(64) FIG. 3A and FIG. 3B show a result of investigating the relationship between the load (maximum load) immediately before the start of load reduction, which is determined as a break in the present embodiment, and the experimental breaking load (maximum load) obtained by experiments by using the tensile test model of an L-shaped joint of which the base material portion is made of a shell element and the spot-welded portion is made of a solid element, in the tensile test of an L-shaped joint as an example. FIG. 3A shows a case where the tensile strength of an L-shaped joint which is an object to be analyzed is equal to or greater than 980 MPa, and FIG. 3B shows a case where the tensile strength of an L-shaped joint which is an object to be analyzed is less than 980 MPa.

(65) A present invention example and a comparative example in FIG. 3A are common in analysis performed under the condition that the tensile strength TS is 1057 MPa, the Young's modulus E is 205800 MPa, the effective width W is 40 mm, the thickness t is 1.6 mm, the nugget diameter D is 6.3 mm, the arm length L is 10 mm, and the total elongation el is 0.15.

(66) In each of the present invention example and the comparative example in FIG. 3A, the tensile test analysis is performed using an analysis model (model A) prepared so that the element size e (average value) of the base material portion is 3.0 mm, an analysis model (model B) prepared so that the element size e (average value) of the base material portion is 4.7 mm and an analysis model (model C) prepared so that the element size e (average value) of the base material portion is 5.8 mm.

(67) Whereas the predicted breaking load is calculated based on the break limit moment calculated using the functions (see Equations (1) to (3)) including the element size e of the base material portion as a variable in the present invention example in FIG. 3A, the predicted breaking load is calculated based on the break limit moment calculated using the function not including the element size e of the base material portion as a variable in the comparative example. Specifically, in the comparative example, the break limit moment is not corrected by the element size e of the base material portion, by inputting a fixed value of 5 mm to the correction term f(e) in any of model A, model B, and model C.

(68) A present invention example and a comparative example in FIG. 3B are common in the analysis performed under the condition that the tensile strength TS is 467 MPa, the effective width W is 50 mm, the thickness t is 1.6 mm, the nugget diameter D is 5.0 mm, and the total elongation el is 0.36.

(69) In each of the present invention example and the comparative example in FIG. 3B, the tensile test analysis is performed using an analysis model (model A) prepared so that the element size e (average value) of the base material portion is 3.0 mm, an analysis model (model B) prepared so that the element size e (average value) of the base material portion is 4.7 mm and an analysis model (model C) prepared so that the element size e (average value) of the base material portion is 5.8 mm.

(70) Whereas the predicted breaking load is calculated based on the break limit moment calculated using the functions (see Equations (4) to (6)) including the element size e of the base material portion as a variable in the present invention example in FIG. 3B, the predicted breaking load is calculated based on the break limit moment calculated using the function not including the element size e of the base material portion as a variable in the comparative example. Specifically, in the comparative example, the break limit moment is not corrected by the element size e of the base material portion, by inputting a fixed value of 5 mm to the correction term f(e) in any of model A, model B, and model C.

(71) As shown in FIG. 3A and FIG. 3B, in the example of the present invention, while the deviation between the predicted breaking load and the experimental breaking load (5.6 kN or 3.9 kN) is also small in any condition that the element size e of the base material portion of the tensile test model of an L-shaped joint is different, there is a case where the deviation between the predicted breaking load and the experimental breaking load is large depending on the element size e of the base material portion of the tensile test model of an L-shaped joint in the comparative example.

(72) In other words, according to the present invention example, the analysis results in FIG. 3A and FIG. 3B show that high break prediction accuracy can be stably obtained without depending on the element size e of the base material portion, when predicting a break of a spot-welded portion of an L-shaped joint which is an object to be analyzed (in particular, a break caused by the application of a moment to the spot-welded portion) by using a finite element method.

(73) As described above, according to the present embodiment, it is possible to stably obtain high break prediction accuracy without depending on the element size of the base material portion, when predicting a break of a joint portion (for example, a spot-welded portion) of an object to be analyzed (for example, a hat type member) including a pair of members joined to each other (in particular, a break caused by the application of a moment to the joint portion) by using a finite element method.

(74) As described above, the function of each component of the break prediction device according to the present embodiment (the parameter acquisition unit 1, the break discernment standard calculation unit 2, and the break discernment unit 3 in FIG. 1), and the break prediction method (first to third steps) according to the present embodiment can be realized by the arithmetic processing unit operating according to the program stored in the nonvolatile storage device of the computer. The program and the computer readable recording medium having the program recorded therein are included in the present embodiment.

(75) Specifically, the program is recorded in a recording medium such as a CD-ROM, or provided to a computer through various transmission media. As a recording medium having the program recorded therein, in addition to a CD-ROM, a flexible disk, a hard disk, a magnetic tape, a magneto-optical disk, a nonvolatile memory card, or the like can be used. On the other hand, as the transmission medium of a program, it is possible to use a communication medium in a computer network system for propagating and providing program information as a carrier wave. Here, the computer network is a local area network (LAN), a wide area network (WAN) such as the Internet, a wireless communication network or the like, and the communication medium is a wired line such as an optical fiber or a wireless line.

(76) Further, the program included in the present embodiment is not limited to the one in which the function of the present embodiment is realized by the computer executing the supplied program. For example, in a case where the function of the present embodiment is realized in cooperation with an operating system (OS) or other application software running on the computer, such program is included in the present embodiment. Further, in a case where all or part of the process of the supplied program is performed by a function expansion board or a function expansion unit of the computer to realize the function of the present embodiment, such a program is included in the present embodiment.

(77) For example, FIG. 4 is a schematic diagram showing an internal structure of a personal user terminal device. In FIG. 4, reference numeral 1200 denotes a personal computer (abbreviated as a PC) provided with a CPU 1201. The PC 1200 executes device control software stored in the ROM 1202 or a hard disk (abbreviated as HD) 1211 or supplied from a flexible disk drive (abbreviated as FD) 1212. The PC 1200 collectively controls each device connected to the system bus 1204.

(78) By the program stored in the CPU 1201, the ROM 1202 or the HD 1211 of the PC 1200, the procedures of the steps S1 to S3 in FIG. 2 of the present embodiment are realized.

(79) Reference numeral 1203 denotes a random access memory (RAM), which functions as a main memory, a work area, and the like of the CPU 1201. Reference numeral 1205 denotes a keyboard controller (abbreviated as KBC), and controls an instruction input from a keyboard (abbreviated as KB) 1209 or a device not shown.

(80) Reference numeral 1206 denotes a CRT controller (abbreviated as CRTC), which controls the display of a CRT display (abbreviated as CRT) 1210. Reference numeral 1207 denotes a disk controller (abbreviated as DKC). The DKC 1207 controls access to the HD 1211 and the FD 1212 that store a boot program, a plurality of applications, an editing file, a user file, and a network management program or the like. Here, the boot program is an activation program for starting the execution (operation) of the hardware and software of the PC 1200.

(81) Reference numeral 1208 denotes a network interface card (abbreviated as NIC), which bidirectionally communicates data with a network printer, another network device, or another PC through the LAN 1220.

(82) Instead of using the PC 1200 which is a personal user terminal device, a predetermined calculator or the like specialized for a break prediction device may be used.

(83) It should be noted that the present invention is not limited to the above embodiment, and the following modifications are conceivable.

(84) (1) In the above embodiment, a description has been made on a method for predicting a break of a joint portion of an object to be analyzed including a pair of members joined to each other, in particular, a break caused by the application of a moment to the joint portion by using a finite element method.

(85) However, it is generally known that not only moment but also shear force and axial force are added to a joint portion, in collision deformation analysis of an object to be analyzed. Therefore, in the break prediction method in the related art, a break caused by the moment, a break caused by the shear force, and a break caused by the axial force are predicted using separate break prediction equations, respectively.

(86) The break prediction method (program) according to the above embodiment can be used as a method (subroutine program) for predicting a break caused by a moment, out of the above three break modes, but in a case where the axial force applied to the joint portion is the compressed axial force, a break does not occur at the time of actual collision, the moment may occur even under the compressed axial force, and there is a possibility of erroneously predicting the occurrence of a moment break.

(87) Therefore, in the third step of the break prediction method according to the present embodiment, in the collision deformation analysis of an object to be analyzed, in a case where the axial force applied to the joint portion is a compressed axial force, it is preferable that a result indicating that there is no break is forcibly output as the break prediction result.

(88) In other words, in step S3, in a case where the axial force applied to a joint portion in the collision deformation analysis of an object to be analyzed is a compressed axial force, it is preferable that the break discernment unit 3 of the break prediction device according to the break prediction device forcibly outputs a result indicating that there is no break as the break prediction result.

(89) This makes it possible to prevent erroneous prediction that the moment that should not occur at the time of actual collision occurs in numerical analysis and a moment break occurs, in a case where the axial force applied to the joint portion is a compressed axial force.

(90) (2) In the above embodiment, the object to be analyzed in which a pair of members is spot-welded is exemplified, and the case of predicting the moment break of the spot-welded portion is exemplified. The joint portion in the present invention is not limited to the spot-welded portion, and the present invention can also be applied to, for example, a case of predicting the moment break occurring in the joint portion of the object to be analyzed including a pair of members joined by other welding methods such as spot welding or line welding.

(91) (3) In the above embodiment, the case where the object to be analyzed is a steel plate is exemplified, but the material of the object to be analyzed in the present invention is not limited to a steel plate, and the present invention can also be applied to a case of predicting the moment break of the object to be analyzed of which material is iron, aluminum, titanium, stainless steel, a composite material (metal-resin material or dissimilar metal material), carbon fiber or the like.

(92) (4) In the above embodiment, the case is exemplified where the break limit moment is calculated by separately using two types of functions (a function used in a case where the tensile strength is 980 MPa or more and a function used in a case where the tensile strength is less than 980 MPa), but three or more functions may be used separately depending on the tensile strength.

(93) (5) The present invention can be applied not only to structural members of vehicles such as hat type members, but also as a method for predicting the moment break of structural members such as various vehicles including railway vehicles, general machinery, or ships.

(94) (6) In the above embodiment, a description has been made on a method for predicting a break of a joint portion of an object to be analyzed including a pair of members joined to each other, in particular, a break caused by the application of a moment to the joint portion by using a finite element method. However, there may be users who request only the method (program) for calculating the break discernment standard.

(95) In response to such a request from the user, a break discernment standard calculation method obtained by deleting the third step from the break prediction method according to the above embodiment may be provided.

(96) That is, the break discernment standard calculation method is a method which calculates a break discernment standard used when predicting a break of a joint portion of an object to be analyzed including a pair of members joined to each other by using a finite element method, and includes a first step of acquiring at least an element size of a base material portion, from among parameters set in an element model for the object to be analyzed; and a second step of calculating, as a break discernment standard, a break limit moment defined by a function including the element size of the base material portion as a variable.

BRIEF DESCRIPTION OF THE REFERENCE SYMBOLS

(97) 1 PARAMETER ACQUISITION UNIT 2 BREAK DISCERNMENT STANDARD CALCULATION UNIT 3 BREAK DISCERNMENT UNIT 4 STORAGE UNIT