Next Article in Journal
A Review Study on Total Ankle Replacement
Next Article in Special Issue
Field Inspection of High-Density Polyethylene (HDPE) Storage Tanks Using Infrared Thermography and Ultrasonic Methods
Previous Article in Journal
Networked Control System Based on PSO-RBF Neural Network Time-Delay Prediction Model
Previous Article in Special Issue
Acoustic Emission and Deep Learning for the Classification of the Mechanical Behavior of AlSi10Mg AM-SLM Specimens
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Model of Self-Magnetic Flux Leakage and Its Application in Estimating the Depth Direction of a Fatigue Crack

by
Jinyi Lee
1,2,3,4,*,
Dabin Wang
2 and
I Dewa Made Oka Dharmawan
2,3,*
1
Department of Electronics Engineering, Chosun University, Gwangju 61452, Republic of Korea
2
Department of Control, Instrumentation, and Robot Engineering, Graduated School of Chosun University, Gwangju 61452, Republic of Korea
3
Interdisciplinary Program in IT-BIO Convergence System, Chosun University, Gwangju 61452, Republic of Korea
4
IT-Based Real-Time NDT Center, Chosun University, Gwangju 61452, Republic of Korea
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2023, 13(1), 533; https://doi.org/10.3390/app13010533
Submission received: 10 November 2022 / Revised: 22 December 2022 / Accepted: 29 December 2022 / Published: 30 December 2022
(This article belongs to the Special Issue Advances in Nondestructive Testing and Evaluation)

Abstract

:

Featured Application

Compact proposal of theoretical model of self-magnetic flux leakage based on magnetic metal memory for estimating depth direction of a fatigue crack.

Abstract

In this study, theoretical models were proposed to explain the changes in self-magnetic flux density (SMFD) due to fatigue cracks in the presence and absence of external magnetic fields. Three theoretical models were proposed: rotation domain model (RDM), concentration domain model (CDM), and vertical domain model (VDM), considering the deformation and non-deformation possibilities. To prove the theoretical model, fatigue cracks with different depth angles were fabricated through fatigue testing and EDM processing on the CT specimens. In addition, tunnel magnetoresistance (TMR) sensors were used to evaluate the 3-axis distribution of SMFD. Comparing the simulation and experimental results, similar tendencies of the occurrence and depth angle of fatigue cracks and their effect on the distribution of SMFD were observed. According to the RDM, the distribution of SMFD occurs in the direction of the crack length (y-direction), while the CDM explains that the SMFD does not occur if the fatigue crack is in a direction perpendicular to the surface. In addition, the VDM shows that SMFDs occur in a direction perpendicular to the crack length (x-direction) and the specimen surface (z-direction). Interestingly, these trends agree with the experimental results, which confirms the validity of the theoretical model and thus can be used to estimate the depth direction of a fatigue crack.

1. Introduction

Large-scale building infrastructures, such as power generation facilities, railroads, bridges, and high-rise buildings, require extensive budgets and time investments. Therefore, social and economic benefits from investments can only be ensured when these infrastructures are operated with safety measures while maintaining the target performance throughout the design life. On the other hand, the damage-tolerant design states that a structure can be continuously operated even after the design life has elapsed if the target performance, safety, and economic feasibility are maintained. In this case, non-destructive (NDT) testing plays an important role to support defect detection. In damage tolerance engineering [1,2,3,4,5], determining tolerance size and appropriate maintenance and replacement are crucial factors. However, the stress concentration coefficients depend on the type of defects, such as voids, corrosion, fatigue cracks, delamination, and shelling, and consequently, the allowable tolerance sizes are different [6,7,8]. Therefore, it is essential to quantitatively evaluate the size and shape of the defects detected by NDT.
One of the challenges of NDT is detecting fatigue cracks that consist of a very thin crack layer. The initiation of fatigue crack on the metal and alloy is a convergence phenomenon that is influenced by many parameter interactions, physically and chemically. Widely known as damage accumulation, crack nucleation and growth occur when the local stress in the atomic displacement exceeds the yield strength of the material [9,10]. The accumulation of this displacement leads to continuous faults, such as deformation, dislocation motions, and evolution, in each atomic layer and initiates micro voids which transform into micro defects with high tendencies to grow with the applied loads. The fatigue crack can be initiated on the material surface or at an internal layer [11]. If fatigue cracks develop in a direction close to the normal to the surface, the structure will fracture rapidly. Furthermore, if fatigue cracks develop in a direction inclined parallel to the surface, the probability of delaminating or shelling increases. Therefore, it is important to evaluate the angle between the fatigue crack and the surface of the structure because this phenomenon can significantly reduce the life span of components and increase the probability of unexpected in-service failures [12]. Many approaches to NDT have been used to detect fatigue defects, one of which is surface defect detection.
In the case of a surface defect in NDT, large mechanical structures are detected using several methods such as visual, penetration, and magnetic particle testing, which can effectively determine the presence and length of a defect. However, it is difficult to implement these methods to evaluate the depth of the material defect [13]. A study on magnetic metal memory (MMM) has been reported in leakage magnetic flux testing which is able to detect an internal defect by measuring the self-magnetic flux leakage (SMFL) [14,15,16,17]. Many related studies have reported the development of this technology for solid ferromagnetic materials and assemblies of ferromagnetic materials, through experimental data models [18,19,20], the probability evaluation model [21], and the magnetic theoretical model. It was found that the SMFD signal around the defect changed depending on the defect types, defect position, and cross-sectional area [17].
Several studies have presented a theoretical equation for the leakage flux generated around a crack using a dipole model [22,23,24,25,26,27,28,29,30,31]. A formula was presented for the magnetization of a specimen by an external magnetic field, assuming that the magnetic dipole moment per unit area occurs at the crack surface. The formula has been used as a model that simulates the qualitative (and partially quantitative) change in MFL based on the presence of defects [22,23,24,25,26,27] and sensor designs [28,29], and is used for quantitative evaluation [22,23,24,25,26,27,29,30,31]. However, it is difficult to apply it to very narrow cracks, such as fatigue cracks. In addition, the change in the magnetic flux density (MFD) with the change in stress remains unclear, where conventional studies on MMM states that the magnetic flux density changes when stress changes or fatigue cracks occur, even when no external magnetic field is applied [32,33,34,35,36]. Therefore, it is uncoherent to use a magnetic dipole that assumes that an external magnetic field is applied which make the specimen in a magnetized state. Therefore, there is a need for a theoretical model to explain the SMFL phenomenon that plays the same role as the magnetic dipole model introduced in the MFLT and SMFL.
Furthermore, to formulate the theoretical model of MMM, several phenomena are considered. Wilson et al. [32] reported that, for applied stress measurements in the stress region of a steel sample, the magnetic field component parallel to the applied stress has a considerably greater correlation than the normal component. This is because of the magnetic field distribution in the stress region, which is related to the magnetic field distribution generated in the sample after the stress is released, and the pattern was confirmed experimentally. On the other hand, Berkache et al. [33] compared the stress analysis using the finite element method (FEM), measured the 3-axial magnetic flux density distribution, and examined the correlation. The stress distribution, which uses the parameter of the 3-axial magnetic flux density, was set as a qualitative relationship in the third polynomial formula. As a result, it was reported that it is possible to examine stress using the MMM. However, a quantitative evaluation of stress is accompanied by a large error. Kim et al. [34] derived a formula in which a change in stress accompanying mechanical deformation causes a change in SMFD. In addition, it was experimentally verified that a change in the SMFD distribution occurred owing to the change in the stress in the weld zone. A system for measuring the residual stress of a wide welding line with a high spatial resolution, using 64 GMR sensor elements arranged at intervals of 0.6 mm, was proposed. In addition to the review of stress performance and its effect on the MFD, the residual magnetic fields have also been considered. Roskosz et al. [35] investigated the residual magnetic field as a diagnostic signal for the reliable evaluation of the material. Accordingly, the possibility of defining the plastic deformation region of a component, determining the approximate distribution of stresses and strains, and defining the remaining lifespan of the component was presented. Finally, Dubov et al. [36] investigated the applicability of the metal magnetic memory (MMM) method for the stress–strain state (SSS) and NDT evaluation of gas and petroleum pipelines. Using the MMM method, it was suggested that the stress concentration zone (SCZ), where the development process of fatigue damage is intensively strengthened, can be detected and used to remove this zone. Unfortunately, despite the several studies on MMM presented above, it is challenging to develop a theoretical model to explain fatigue cracking and evaluate the depth direction of fatigue cracks by measuring the self-magnetization, and deciding whether the edge of the ferromagnetic metal structure is magnetized or the magnetic flux density shifted near fatigue crack despite the absence of external magnetic field.
Dealing with this gap, the change in magnetic domain caused by SSS and fatigue cracking with a narrow width was reviewed and modeled in this study. A valid domain model was selected to evaluate the depth direction inside fatigue cracks. The depth direction of the fatigue cracks was evaluated by measuring the MFD distribution based on the self-magnetization of the subject. Carbon steel specimens (SA 106) were subjected to fatigue crack tests using a compact tension (CT) specimen, and a small volume of 8 × 16 × 6.4 mm (width × length × height) containing the crack was cut with different cutting angles of 30, 60, and 90°. The MFD distribution on the surface of the specimen was measured using a 3-axis tunnel magneto-resistance (TMR) sensor. An evaluation algorithm for the depth direction was derived using the 3-axial SMFD distribution.

2. Magnetic Domain Models for Stress and Fatigue Cracks

Definition and statements of magnetic dipole moment ( m ), magnetization ( M ), magnetic flux density ( B ), vector magnetic potential ( A ), and magnetization surface current density ( J ms ) are referred to as reference [37] in this paper. A magnetic dipole moment m , is assumed to be generated by electrons rotating and orbiting around the nucleus of a ferromagnetic metal. Regions with magnetic moments in the same direction are called magnetic domains and the boundaries of these regions are called domain walls. Magnetization M is defined by Equation (1). Here, v is the micro-volume in the 3-dimensional model, and can be replaced by the micro-area s in the 2-dimensional model.
M = lim v 0 k = 1 v m k   v ,   ( A / m )
Figure 1a,b shows the rotation domain model (RDM) in which the magnetization is rotated. Figure 1c,d show a concentration domain model (CDM), which is arranged in the direction in which the magnetization is concentrated at the center. Figure 1e,f show the vertical domain model (VDM), where the magnetization is arranged in a direction perpendicular to the surface. RDM, CDM, and VDM are separate concepts from the domain theory of magnetism [38]. Figure 1a,c), 1e show a two-dimensional magnetic domain model with magnetizations M 1 , M 2 , M 3 , and M 4 in each magnetic domain with areas S1, S2, S3, and S4 in a state where no external magnetic field H is applied ( H = 0 ). In contrast, the equivalent magnetization M eqv can be expressed as weighted average from the total surface area obtained by multiplying the magnetization by the corresponding area of the magnetic domain S k divided by total area of magnetic domain, as shown in Equation (2), when H = 0 , M eqv = 0 . In Equation (2), n is the number of magnetic domains. In this study, we specified n = 4 or n = 8 correspond to the number of models arrayed to the system.
M eqv = 1 S k = 1 n S k M k ,   ( A / m )
S = k = 1 n S k ,   ( m 2 )
When an external magnetic field H is applied ( H 0 ), as shown in Figure 1d,e, the magnetic domain wall moves, and the area of the magnetic domain (S2 or S1) with a magnetic moment in the same direction as H increases. In addition, the area of magnetization (S4 or S3) with a magnetic moment opposite to H becomes smaller. Consequently, in the state where H 0 , M eqv 0 . However, in the VDM of Figure 1f, even when an external magnetic field H is applied ( H 0 ), M eqv = 0 , so in the two-dimensional magnetic domain model, the changes of RDM and CDM are justifiable. The equivalent magnetization M eqv can be expressed as the product of magnetic susceptibility χ m , which is an intrinsic property of a material, and the external magnetic field H , which is expressed in Equation (4).
M eqv = χ m H ,   ( A / m )
The magnetic susceptibility is not a constant entity, but corresponds to a parameter for an external magnetic field, which causes magnetic hysteresis. In other words, even after the external magnetic field is removed, the residual magnetization phenomenon of M eqv 0 appears. In this state, for M eqv = 0 , an external magnetic field in the reverse direction, i.e., a coercive force, must be applied.
The MFD B outside the ferromagnetic metal can be expressed as the superposition of the external magnetic field H and the equivalent magnetization, as in Equation (5). Here, μ 0 , μ r , and μ are the permeability in free space, relative permeability, and permeability, respectively. According to Equation (5), when H = 0 , M eqv = 0 .
B = μ 0 H + μ 0 M eqv = μ 0 1 + χ m H = μ 0 μ r H = μ H ,   ( Wb / m 2 )
μ 0 = 4 π × 10 7 ,   ( H / m )
Although it is not a magnetic hysteresis phenomenon, as shown in Figure 2, when a mechanical load is applied and deformation occurs, a phenomenon in which M eqv 0 may occur at H = 0 . Figure 2a–c show conceptual diagrams applying RDM, CDM, and VDM, respectively, with n = 8 to a cantilever beam. Figure 2d–f show the deformation that occurred when a load F was applied to the tip of the cantilever beam. The area Sk of each magnetic domain was deformed into Sk′ by the load F . The strain ε α β in each area can be expressed using Equation (7). In addition, the stress σ α β in each area can be expressed using Equation (8). In this study, both the tensile and shear stresses are expressed as σ α β , (Equation (8)), where E is the elastic modulus and α   or   β refers to the x or y direction
S k = ε α β , k S k ,   ( m 2 )
σ α β , k = E   ε α β , k ,   ( N / m 2 )
In the state of Figure 2d,e, when H = 0 , the equivalent magnetization M eqv and the magnetic flux density B can be expressed by Equations (9) and (10), respectively. S is the sum of the total area. In addition, Equation (10) can be expressed as Equations (11) and (12) by substituting in Equations (7) and (8). This formula indicates that a change in strain or stress causes a change in magnetic flux density.
M eqv = 1 S k = 1 n S k M k ,   ( A / m )
B = μ 0 M eqv = μ 0 S k = 1 n S k M k ,   ( Wb / m 2 )
B = μ 0 S k = 1 n ε α β , k S k M k ,   ( Wb / m 2 )
B = μ 0 SE k = 1 n σ α β , k S k M k ,   ( Wb / m 2 )
Figure 3 shows an example of a 3-dimensional representation of the RDM model in Figure 2a by extending it in the direction perpendicular to the xy plane, that is, the z-direction. This gap can be expressed by separating regions S1–S4 and S5–S8 in Figure 2. Hereafter, the gap is collectively referred to as a crack. Figure 3b shows a magnetic domain model with cracks in the direction perpendicular to the surface of the specimen, and Figure 3c shows a magnetic domain model with slanted cracks while maintaining hexahedronal shape. Figure 3d shows a case in which the top surface of the magnetic domain is deformed to be horizontal to the xy plane and has an inclined crack. In each figure, a n denotes the surface unit vector. In addition, the unit vectors in the x, y, and z directions are expressed as a x , a y , and a z . The magnetization surface current density (MSCD), J ms , represented by the dotted arrow, is related to the magnetization ( M k ), represented by the solid arrow and Equation (13).
In addition, the relationship between the equivalent MSCD ( J eqv , mn ) and the equivalent magnetization ( M eqv , mn ) can be expressed by Equation (14). In the case where there is no external magnetic field ( H = 0 ) and when no deformation occurs ( ε α β , k = 0 ), M eqv , 14 = 0 and M eqv , 58 = 0 . This corresponds to the case shown in Figure 3a, where the equivalent magnetization M eqv , 35 by two adjacent magnetizations M 3 and M 5 is 0 . Moreover, in the case of Figure 3b–d, assuming that there are continuous magnetic domains in the ±x and ±y directions, M eqv , 14 = 0 , M eqv , 58 = 0 . However, because the equivalent magnetization expressed in Equation (2) is a model limited to adjacent magnetic domains bordering the magnetic domain wall, the equivalent magnetization M eqv , 35 by M 3 and M 5 is not 0 . At this time, MSCDs J 3 and J 5 , are located at the interface of the crack.
J ms , k = M k × a n , k ,   ( A / m )
J eqv , mn = M eqv , mn × a eqv , mn ,   ( A / m )
Figure 3. Simplified 3-D magnetic domains: (a) without gap, (b) with gap, (c) with gradient gap without deformation, (d) with gradient gap with deformation.
Figure 3. Simplified 3-D magnetic domains: (a) without gap, (b) with gap, (c) with gradient gap without deformation, (d) with gradient gap with deformation.
Applsci 13 00533 g003
Figure 4, Figure 5 and Figure 6 shows the simplified magnetization surface current densities on the RDM, CDM, and VDM. Figure 4 shows a schematic of J k , edge and M k , edge in the xz plane by RDM at the crack tip. When the width of the crack is w, Equation (13) can be expressed as Equation (15): Here, the intensity of M k , edge is assumed to be M 0 . w and θ represent the crack width and inclination angle, respectively. P x , h z represents a specific position at the top of the magnetic domain model. D is the depth of the crack, and h is the lift-off. Equation (15) is a general expression expressing all of Figure 3b–d.
J RDM , k , edge = M 0 sin θ a x + cos θ a z ,   ( A / m )
Figure 5 shows a schematic of J k , edge and M k , edge in the xz plane by CDM at the crack tip. According to the CDM domain model in the state in which the domain is not deformed, as shown in Figure 3c, Equation (13) can be expressed as Equation (16). This is generally represented as in Figure 5b,c. According to the CDM in the state in which the domain is deformed, as shown in Figure 5d, Equation (13) can be expressed as Equation (17).
J CDM , k , edge = 0 ,   ( A / m )
J CDM , k , edge = M 0 sin θ a y ,   ( A / m )
Figure 5. Simplified magnetization surface current densities on the concentration domain model: (a) without gap, (b) with gap, (c) with gradient gap without deformation, (d) with gradient gap with deformation.
Figure 5. Simplified magnetization surface current densities on the concentration domain model: (a) without gap, (b) with gap, (c) with gradient gap without deformation, (d) with gradient gap with deformation.
Applsci 13 00533 g005aApplsci 13 00533 g005b
Figure 6 is a schematic of J k , edge and M k , edge in the xz-plane by VDM at the crack tip. According to the VDM in the state in which the magnetic domain is not deformed, as shown in Figure 3c, Equation (13) can be expressed as Equation (18). This is generally represented by Figure 6b,c. According to the VDM in the state in which the magnetic domain is deformed, as shown in Figure 6d, Equation (13) can be expressed as Equation (19).
J VDM , k , edge = M 0 a y ,   ( A / m )
J VDM , k , edge = M 0 cos θ a y ,   ( A / m )
Figure 6. Simplified magnetization surface current densities on the vertical domain model: (a) without gap, (b) with gap, (c) with gradient gap without deformation, (d) with gradient gap with deformation.
Figure 6. Simplified magnetization surface current densities on the vertical domain model: (a) without gap, (b) with gap, (c) with gradient gap without deformation, (d) with gradient gap with deformation.
Applsci 13 00533 g006
The vector magnetic potential (VMP) ( A ) and MFD ( B ) are expressed by Equations (20) and (21), respectively. In RDM, if the magnetic domain model is not subjected to an external load, A can be expressed by Equations (22)–(24), where R 3 and R 5 represent the position vectors from the crack wall to a specific position, P(x,h).
d A = μ 0 4 π C J R dl ,   ( Wb / m )
d B = × d A ,   ( Wb / m 2 )
d A RDM x , h = μ 0 M 0 4 π C 1 R 3 + 1 R 5 sin θ a x + cos θ a z dy
R 3 = h z 2 + w 2 + x 2
R 5 = h z 2 + w 2 x 2
The unit area per VMPs, A RDM , a , A RDM , b , and A RDM , c for each case of RDM shown in Figure 4a, b and c can be expressed by Equations (25)–(27), respectively. In Equation (27), w = 0 and θ = 0° can be expressed by Equations (25) and (26), respectively. Therefore, Equation (27) is a general formula that includes the cases with and without cracks. In addition, as shown in Figure 3d, inclined cracks caused by the undeformed and deformed RDMs are expressed in the same manner as in Equations (27) and (28), respectively. Therefore, MFD ( B ) in the RDM can be interpreted using Equations (29) and (30).
d A RDM , a x , z = 0
d A RDM , b x , h = μ 0 M 0 4 π 1 R 3 + 1 R 5 a z
d A RDM , c x , h = μ 0 M 0 4 π 1 R 3 + 1 R 5 sin θ a x + cos θ a z
d A RDM , d x , h = μ 0 M 0 4 π 1 R 3 + 1 R 5 sin θ a x + cos θ a z
d B RDM x , h = μ 0 M 0 4 π w 2 + x R 3 3 w 2 x R 5 3 cos θ + h z R 3 3 + h z R 5 3 sin θ a y
B RDM x , h = 0 D d B RDM dz
In CDM, A is interpreted using Equations (31)–(34). The B CDM , c by inclined cracks in the non-deformed CDM, is expressed by Equation (35). The B CDM , d for cracks analyzed using the deformed CDM is interpreted as Equation (36) and (37).
d A CDM , a x , z = 0
d A CDM , b x , z = 0
d A CDM , c x , z = 0
d A CDM , d x , z = μ 0 M 0 4 π 1 R 3 + 1 R 5 sin θ a y
d B CDM , c x , z = 0
d B CDM , d x , z = μ 0 M 0 4 π h z R 3 3 + h z R 5 3 sin θ a x + w 2 + x R 3 3 w 2 x R 5 3 sin θ a z
B CDM , d x , h = 0 D d B CDM , d dz
Equations (38)–(41) show A in VDM. The MFDs ( B VDM , c , B VDM , d ) generated around the slanted crack caused by the non-deformed and deformed VDM are interpreted by Equations (42) and (43), and Equations (44) and (45), respectively.
d A VDM , a x , z = 0
d A VDM , b x , z = μ 0 M 0 4 π 1 R 3 + 1 R 5 a y
d A VDM , c x , z = μ 0 M 0 4 π 1 R 3 + 1 R 5 a y
d A VDM , d x , z = μ 0 M 0 4 π 1 R 3 + 1 R 5 cos θ a y
d B VDM , c x , z = μ 0 M 0 4 π h z R 3 3 + h z R 5 3 a x + w 2 + x R 3 3 w 2 x R 5 3 a z
B VDM , c x , h = 0 D d B VDM , c dz
d B VDM , d x , z = μ 0 M 0 4 π h z R 3 3 + h z R 5 3 cos θ a x + w 2 + x R 3 3 w 2 x R 5 3 cos θ a z
B VDM , d x , h = 0 D d B VDM , d dz
In the RDM by Equation (29), when cracks occur, only the MFD of the y-component exists regardless of the deformation of the magnetic domain. In the undeformed CDM model, Equation (35), the MFD does not change, even if cracks occur. However, in the deformed CDM model according to Equation (36), when an inclined crack occurs, the MFDs in the x- and z-components change. In the undeformed VDM model using Equation (42), the MFD of the x- and z-components changes regardless of the inclination direction of the crack. In addition, in the deformed VDM model in the inclined crack according to Equation (44), the MFDs in the x- and z-components change. Moreover, the intensity of the MFD increases as the inclination angle θ of the crack approaches 0°.
Figure 7 shows a graph expressing the components in each direction expressed by Equation (29) according to the depth direction of the crack. To understand the qualitative trend, it was assumed that M 0 = 4 π / μ 0 , h = 0.2 mm, w = 0.01 mm, and D = 5 mm. The y-component of B RDM , c according to Figure 7 has a positive sign on the left and a negative sign on the right, near the crack center. As the angle of the crack increases, the negative maximum value on the right-side increases. In B CDM , d , according to Figure 8, only the x- and z-components appear, and when θ = 0 ° , the magnetic flux density is 0. In addition, in B CDM , c by Equation (35), the magnetic flux density does not change. In B VDM , c according to Figure 9, only the x- and z-components appear, and the magnitude of the absolute value increases with the angle. Because the length in the oblique direction increases for the same depth, a larger value appears as the inclination angle increases. However, the x-component indicates the maximum value in the vicinity of the crack and shows a symmetrical distribution. Further, the z-component has a negative maximum value on the left, centered on the maximum value of the x-component, indicating a distribution opposite to the y-component of B RDM , c (Figure 7). In B VDM , d shown in Figure 10, only the x- and z-components appear, and the x- and z-components that decrease in inverse proportion to the size of the angle can be observed. Furthermore, these results will form the theoretical model for estimating the depth direction of a fatigue crack. This trend is compared to the experimental results mentioned in the next section.

3. Materials and Methods

3.1. Specimen with Fatigue Crack

Fatigue crack specimens were prepared to investigate the validity of the proposed magnetic domain model and the change in MFD distribution according to the crack inclination angle. Figure 11 shows a schematic of the extraction of fatigue specimens from compact test (CT) specimens with fatigue cracks. A fatigue crack with a depth of 3–7 mm was created by applying a repeated load of 350–3500 N, 630,000–1,080,000 times to a CT specimen made of SS400 carbon steel with the dimensions 32 × 30.72 × 12.8 mm. As shown in Figure 12, small specimens of 6.4 × 16 × 8 mm were machined at different cutting angles by electric discharge machining (EDM). Consequently, each specimen contained fatigue cracks inclined at 0, 30, and 60°. Figure 13 shows the micrographs of the surface and side of each specimen. Each specimen was polished using 200, 400, 800, 1200, and 2000 grades SiC pads. Finally, fine polishing was performed for 20 s using polycrystalline diamond paste with a lattice size of 1 μm and a polishing cloth. After micro-polishing was completed, an etching solution prepared by diluting 3 g of nitric acid in 65 g of ethanol was applied to the observation surface for 5 s, followed by ultrasonic cleaning with distilled water for 15 s and washing with ethanol for 10 s. After repeating this chemical etching thrice, it was photographed with an optical microscope, and the 3-axial MFD was measured at a lift-off of 0.2 mm.

3.2. Measurement of Magnetic Flux Density with 3-Axis TMR Sensor Scanning

Figure 14a shows the experimental setup for measuring the MFD distribution on the upper part of each specimen using a 3-axis TMR sensor. A microscope and a 3-axis TMR sensor were installed on an optical board. Figure 14b shows an enlarged photograph of the red square in Figure 14a. The specimen was fixed with a key to reproduce the absolute positions of the optical images and TMR scanning before and after polishing. In addition, the specimen fixed by the key measured the MFD distribution with a spatial resolution of 4 μm and acquired an optical microscope image using a motorized XY stage. The MFD distribution was measured using a 3-axis TMR sensor (TMR2305, Multidimension Co., Paris, France). The sensitivity of each sensor was 25 mV/V/Oe, and the saturation magnetic field was ±10 Oe. As a result of amplification by 46 dB (200 times) by the signal processing circuit, 1 V is the output for an MFD of 0.2 G. The distance between the surface of the specimen and the sensor surface and the lift-off was maintained at approximately 0.2 mm.

4. Experimental Results

Training, testing, and analysis processes were performed using MATLAB R2021a software on a Dell ALIENWARE PC with an Intel® Core™ i9-9900HQ 3.6 GHz of CPU and 32 GB of RAM. Several positive results were achieved in this study and a theoretical model of magnetic metal memory was proposed to estimate the depth direction of a fatigue crack (see Figure 7, Figure 8 and Figure 9). To provide a comparison, an experiment was conducted using a 3-axis TMR sensor on CT scan specimens with fatigue cracks on three variations of depth angle: 0 o , 30 o , and 60 o . In addition, to maintain the validity of the measurement, the experiment was conducted on a wide demagnetized aluminum plate to reduce the probability of an external magnetic field disturbing the recorded dataset.
Figure 15 shows the raw output signal of the 3-axis TMR sensor in the x- ( V raw , X ), y- ( V raw , Y ), and z- ( V raw , Z ) directions. To map the distribution of the magnetic fields, the contour lines are added to all of the variations. In addition, the coordination of the specimen boundaries and the centerline line of the fatigue cracks are represented by the brown box and thick pink line on the vertical axis. The upper part of the rectangle is the surface side of the CT specimen, and the lower part is the center surface which is cut by EDM processes. To provide clear data for supporting the analysis, signal processing was conducted. The idea is to delete the background distribution of SMFD data by subtracting the raw data ( V raw x , y ) from the background distribution ( V back x , y ). The V back x , y was gathered by measuring the distribution of SMFD in the air by removing the specimen from the setup. The signal processing processes can be described mathematically as follow:
V 1 x , y = V raw x , y V back x , y
V raw x , y = V raw , X x , y a x + V raw , Y x , y a y + V raw , Z x , y a Z
V back x , y = V back , X x , y a x + V back , Y x , y a y + V back , Z x , y a Z
where, Figure 16 shows the distribution of Equation (46) and corresponds to the background signal in the x- ( V back , x ), y- ( V back , Y ), and z- ( V back , Z ) directions. Interestingly, the overall trends in each case in Figure 15 and Figure 16 are very similar, indicating that the external magnetic field has negligible effect on the experimental values. As shown in Figure 16a,d,g, the X-component, V 1 , X , has a maximum value near the center of the crack, regardless of the depth angle of the crack, and has a symmetrical shape. As shown in Figure 16c,f,i, the Z-component, V 1 , Z , represents the distribution of the different signs on the left and right sides around the crack. This phenomenon can be explained by the VDM in Figure 9 and Figure 10. Moreover, as shown in Figure 16b,e,h, the Y-component, V 1 , Y , does not differ significantly along the center line of the crack, indicating that the VDM model is valid.
To investigate the quantitative trend of MFD based on the angle direction to the depth of cracks, Figure 17, Figure 18 and Figure 19 show the MFD distribution along a line passing through the center of the specimen (y = 5.1 mm) in Figure 15 and Figure 16. In each graph, V 2 x , y is the value expressed by Equation (49). V 1 ,   avr x , y is the average value expressed by Equation (50). This variable was introduced to observe the correlation between the position of x and the crack inclination angle when V 2 x , y = 0 . In addition, M, N, and r are the number of data points in the x-axis direction (M = 5000), number of data points in the y-axis direction (N = 2500), and spatial resolution (r = 0.004 mm), respectively. The thick solid line near the center of the longitudinal axis (x = 9.3–9.5 mm) is the central position of the crack observed under the microscope. The solid vertical lines at x = 1.4 mm and x = 15.9 mm indicate the positions of both corners of the specimen.
V 2 x , y = V 1 x , y V 1 ,   avr x , y
V 1 , avr x , y = 1 MN i = 1 M j = 1 N V 1 , X ir , jr a x + i = 1 M j = 1 N V 1 , Y ir , jr a y + i = 1 M j = 1 N V 1 , Z ir , jr a Z
Comparing V raw , x x , 5.1 in Figure 17a and V 1 , x x , 5.1 in Figure 17b, the background signal of about −0.9 V (approximately −0.18 G) is applied uniformly. As the overall depth angle increased, the peak value also increased. This means that B CDM , c x , z in Figure 9 exhibits a more accurate trend than B CDM , d x , z in Figure 10, i.e., for the same depth, the area of magnetization widens according to the inclination angle, which means that a larger MFD distribution is obtained as a result.
In contrast, the position of the peak was deflected to the right (direction with an acute angle) according to the depth angle. Therefore, if the central position of the crack can be specified by penetrant inspection, the crack direction to depth based on the correlation with the peak position can be estimated. However, it is difficult to estimate the crack direction if only the MFD distribution is provided without knowing the exact location of the crack. Furthermore, according to V 2 , x x , 5.1 in Figure 17c, the left–right symmetry collapses with respect to the peak according to the depth direction. To understand the left-right symmetry of V 2 , x , an angle estimation function defined by Equation (51), x 2 , x , was proposed. x 2 , x is defined as the distance between x L V 2 , x ; 0 , the left peak position, and x R V 2 , x ; 0 , the right peak position close to the x position of the peak, x C V 2 , x ; peak . The relationship between x 2 , x and the crack angle is discussed in Figure 20, along with another angle estimation function, Σ V 2 , z .
x 2 , x = x C V 2 , x ; peak x L V 2 , x ; 0 x C V 2 , x ; peak x R V 2 , x ; 0
Figure 17. MFD distribution in the X direction on a line passing through the center of the specimen (y = 5.1 mm): (a) raw data; with (b) background data subtracted; and (c) average value subtracted.
Figure 17. MFD distribution in the X direction on a line passing through the center of the specimen (y = 5.1 mm): (a) raw data; with (b) background data subtracted; and (c) average value subtracted.
Applsci 13 00533 g017
Figure 18. MFD distribution in the Y direction on a line passing through the center of the specimen (y = 5.1 mm): (a) raw data; with (b) background data subtracted and (c) average value subtracted.
Figure 18. MFD distribution in the Y direction on a line passing through the center of the specimen (y = 5.1 mm): (a) raw data; with (b) background data subtracted and (c) average value subtracted.
Applsci 13 00533 g018
Figure 19. MFD distribution in the Z direction on a line passing through the center of the specimen (y = 5.1 mm): (a) raw data; with (b) background data subtracted and (c) subtracted by average value.
Figure 19. MFD distribution in the Z direction on a line passing through the center of the specimen (y = 5.1 mm): (a) raw data; with (b) background data subtracted and (c) subtracted by average value.
Applsci 13 00533 g019
Comparing V raw , y x , 5.1 in Figure 18a and V 1 , y x , 5.1 in Figure 18b, the background signal of approximately −0.5 V (at approximately −0.1 G) is applied uniformly. An MFD with different signs to the left and right around the crack occurred, and the overall absolute value increased as the crack direction to depth increased. This can be explained by the trend of the y-component of B RDM x , z in Figure 7.
Comparing V raw , Z x , 5.1 in Figure 19a and V 1 , z x , 5.1 in Figure 19b, a background signal of about 0.05 V (corresponding to about 0.01 G) is applied uniformly. MFDs with different signs are observed on the left and right sides of the cracks. At this time, the left and right signs are opposite to those shown in Figure 18. This is consistent with the theoretical consideration of RDM shown in Figure 9. The overall absolute values of V 1 , z x , 5.1 increased as the overall inclination angle increased. This can be explained by the tendency of the z-component of B RDM x , z , as shown in Figure 9. From the above results, the domain model for evaluating fatigue cracks is validated as a model in which the RDM where the y-component occurs and the VDM where the x- and z-components overlap. In addition, a model in which magnetic domains are not deformed when cracks occur is appropriate because the change in the MFD increases based on the depth direction of fatigue cracks at the same depth. When stress is generated, it can be expressed as Equation (52). Here, C1, C2, and C3 are constants that represent the respective ratios when the stress change, RDM, and VDM overlap.
B = C 1 μ 0 SE k = 1 n σ α β , k S k M k + C 2 B RDM , c + C 3 B VDM , c ,   ( Wb / m 2 )
However, according to V 2 , z x , 5.1 ) in Figure 19c, the position and area of the positive peak on the right and the negative peak on the left with respect to the zero point vary according to the inclination angle. To understand the left-right symmetry of V 2 , z , an angle estimation function defined by Equation (53)–(55), Σ V 2 , z , is proposed.
Σ V 2 , z is defined as the difference between the integral up to the negative peak ( Σ V L V 2 , z ; negative   peak ) and the integral up to the positive peak ( V 2 , z ; positive   peak ) from the x position on V 2 , z = 0 .
Σ V 2 , z = Σ V R Σ V L
Σ V R = x V 2 , z = 0 x V 2 , z = positive   peak V 2 , z dV
Σ V L = x V 2 , z = negative   peak x V 2 , z = 0 V 2 , z dV
Figure 20 shows the relationship between x 2 , x and Σ V 2 , z expressed by Equations (51) and (53)–(55), and the crack angle. In the case of x 2 , x , the lowest value appears when θ = 30°, making it difficult to evaluate the crack angle quantitatively. However, according to Σ V 2 , z reflecting the peak position ( x C ) in the distribution of V 2 , x , the depth direction of fatigue crack can be quantitatively evaluated using empirical Formula (57).
Σ V L = 0.1778 θ 2 + 2.781 θ + 209
θ = 8.671 × 10 5 Σ V L 2 + 0.1806 Σ V L 33.96
Figure 20. Relationship between the angle of fatigue crack and angle estimation functions.
Figure 20. Relationship between the angle of fatigue crack and angle estimation functions.
Applsci 13 00533 g020

5. Conclusions

In this study, theoretical models of magnetic metal memory (MMM) have been proposed. The idea was to utilize these models to estimate the depth direction of a fatigue crack. Three models were used: the rotation domain model (RDM), concentration domain model (CDM), and vertical domain model (VDM). In addition, case studies on the deformation and non-deformation conditions when cracks occur were also studied theoretically. In general, the results found that using the magnetic domain model, the self-magnetic flux density (SMFD) changes owing to the stress–strain state (SSS) and fatigue cracking, even when no external magnetic field is applied to the ferromagnetic specimen. This indicates that MMM can be used to detect internal defects by measuring the self-magnetic flux leakage (SMFL) on the specimen. This result was observed on both the theoretical and experimental approaches, which supports the model validity.
In the theoretical model, if the magnetic domain is not deformed when the fatigue cracking occurs, RDM explains that the SMFL will appear on the longitudinal component of the crack (y-direction). Meanwhile, the SMFL of the longitudinal component (y-direction) of the crack occurs according to the deformed RDM, accompanied by fatigue cracking. On the other hand, CDM explains that if the magnetic domain is not deformed when fatigue cracks occur, the SMFL components perpendicular to the crack length (x-and z-directions) will exist, while the longitudinal component of the crack (y-direction) does not occur. Compared to the CDM, VDM pursues a similar tendency on each SMFL component. However, VDM found that the larger the depth direction of crack angle, the larger the SMFD intensity of each component.
Furthermore, to reinforce the results from theoretical models, the magnetic flux density was measured using a 3-axis TMR sensor in a state where no external magnetic field was applied to the specimen with fatigue cracks. It was found that although no external magnetic field was applied, the SMFD of the x-, y-, and z-components changed. Furthermore, in the identical tendency with the VDM, as the inclination angle of the crack increased, the magnitude of the magnetic flux density of each component increased. Therefore, it is justifiable to explain the change in the self-magnetic flux density due to the occurrence of fatigue cracks using a magnetic domain model in which the undeformed RDM and VDM are superimposed. Therefore, by using the integral value of the z-component reflecting the peak position of the x-component, the depth angle of the cracks can be quantitatively evaluated.
Finally, it can be concluded that both the theoretical models and experimental results have similar trends. This supports the model validity, and therefore, the model can be used to estimate the depth direction of a fatigue crack.

Author Contributions

Conceptualization, J.L.; methodology, J.L., D.W. and I.D.M.O.D.; software, J.L. and I.D.M.O.D.; validation, J.L. and I.D.M.O.D.; formal analysis, J.L.; investigation, J.L., D.W. and I.D.M.O.D.; resources, J.L. and D.W.; data curation, J.L., D.W. and I.D.M.O.D.; writing—original draft preparation, J.L.; writing—review and editing, J.L. and I.D.M.O.D.; visualization, J.L., D.W. and I.D.M.O.D.; supervision, J.L. and I.D.M.O.D.; visualization project administration, J.L.; funding acquisition, J.L. This paper was prepared with the contributions of all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. NRF-2019R1A2C2006064).

Data Availability Statement

The code and dataset will be made available on request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

MFDMagnetic flux density
SMFDSelf-magnetic flux density
SMFLSelf-magnetic flux leakage
FEMFinite element method
RDMRotation domain model
CDMConcentration domain model
VDMVertical domain model
CTCompact tension
MMMMagnetic metal memory
TMRTunnel magneto-resistance
NDTNon-destructive testing
SSSStress–strain state
SCZStress concentration zone
MSCDMagnetization surface current density
VMPVector magnetic potential
EDMElectric discharge machining

References

  1. Sierakowski, R.L.; Newaz, G.M. Damage Tolerance in Advanced Composites; Routledge: London, UK, 2018; ISBN 9781315140735. [Google Scholar]
  2. Gosselin, S.R.; Simonen, F.A.; Heasler, P.G.; Doctor, S.R. Fatigue Crack Flaw Tolerance in Nuclear Power Plant Piping—A Basis for Improvements to ASME Code Section XI Appendix L (NUREG/CR-6934); Pacific Northwest National Laboratory: Richland, WA, USA, 2007; p. 99352. [Google Scholar]
  3. Zerbst, U.; Lundén, R.; Edel, K.-O.; Smith, R.A. Introduction to the Damage Tolerance Behaviour of Railway Rails—A Review. Eng. Fract. Mech. 2009, 76, 2563–2601. [Google Scholar] [CrossRef] [Green Version]
  4. Valiente, A.; Iordachescu, M.; Aznar, A.; de Abreu, M. Damage Tolerance of the Prestressing Steel in a Bridge Stay-Cable after Thirty Years of Service in a Graded, Extremely Adverse Environment. Constr. Build. Mater. 2021, 285, 122945. [Google Scholar] [CrossRef]
  5. De Abreu, M.; Iordachescu, M.; Valiente, A. Effects of Hydrogen Assisted Stress Corrosion on Damage Tolerance of a High-Strength Duplex Stainless Steel Wire for Prestressing Concrete. Constr. Build. Mater. 2014, 66, 38–44. [Google Scholar] [CrossRef] [Green Version]
  6. Pao, P.S.; Holtz, R.L.; Jones, H.N.; Feng, C.R. Effect of Environment on Fatigue Crack Growth in Ultrafine Grain Al–Mg. Int. J. Fatigue 2009, 31, 1678–1683. [Google Scholar] [CrossRef]
  7. Pao, P.S.; Holtz, R.L. Effect of Crack Tip Stress Concentration Factor on Fracture Resistance in Vacuum Environment; Naval Research Laboratory: Washington, DC, USA, 2015. [Google Scholar]
  8. Jones, R. Fatigue Crack Growth and Damage Tolerance. Fatigue Fract. Eng. Mater. Struct. 2014, 37, 463–483. [Google Scholar] [CrossRef]
  9. Stachurski, Z.H.; Wang, G.; Tan, X. Fracture Mechanics of Metallic Glasses. In An Introduction to Metallic Glasses and Amorphous Metals; Elsevier: Amsterdam, The Netherlands, 2021; pp. 313–335. [Google Scholar] [CrossRef]
  10. Roeder, R.K. Mechanical Characterization of Biomaterials. In Characterization of Biomaterials; Newnes: Southampton, UK, 2013; pp. 49–104. [Google Scholar] [CrossRef]
  11. Starke, E.A.; Staley, J.T. Application of Modern Aluminium Alloys to Aircraft. In Fundamentals of Aluminium Metallurgy: Production, Processing and Applications; Woodhead Publishing: Sawston, UK, 2011; pp. 747–783. [Google Scholar] [CrossRef]
  12. Kong, Y.; Bennett, C.J.; Hyde, C.J. A Review of Non-Destructive Testing Techniques for the in-Situ Investigation of Fretting Fatigue Cracks. Mater. Des. 2020, 196, 109093. [Google Scholar] [CrossRef]
  13. Dharmawan, I.D.M.O.; Lee, J.; Sim, S. Defect Shape Classification Using Transfer Learning in Deep Convolutional Neural Network on Magneto-Optical Nondestructive Inspection. Appl. Sci. 2022, 12, 7613. [Google Scholar] [CrossRef]
  14. Wang, H.; Dong, L.; Wang, H.; Ma, G.; Xu, B.; Zhao, Y. Effect of Tensile Stress on Metal Magnetic Memory Signals during On-Line Measurement in Ferromagnetic Steel. NDT E Int. 2021, 117, 102378. [Google Scholar] [CrossRef]
  15. Roskosz, M. Metal Magnetic Memory Testing of Welded Joints of Ferritic and Austenitic Steels. NDT E Int. 2011, 44, 305–310. [Google Scholar] [CrossRef]
  16. Dubov, A.A. A Study of Metal Properties Using the Method of Magnetic Memory. Met. Sci. Heat Treat. 1997, 39, 401–405. [Google Scholar] [CrossRef]
  17. Meng, Q.; Pan, P.; Yang, X.; Song, J.; Wang, J. Self-Magnetic Flux Leakage-Based Detection and Quantification for High-Strength Steel Wires of Bridge Cables Considering Corrosion-Fatigue Coupling Effect. J. Magn. Magn. Mater. 2022, 561, 169641. [Google Scholar] [CrossRef]
  18. Karthik, M.M.; Terzioglu, T.; Hurlebaus, S.; Hueste, M.B.; Weischedel, H.; Stamm, R. Magnetic Flux Leakage Technique to Detect Loss in Metallic Area in External Post-Tensioning Systems. Eng. Struct. 2019, 201, 109765. [Google Scholar] [CrossRef]
  19. Xia, R.; Zhang, H.; Zhou, J.; Liao, L.; Zhang, Z.; Yang, F. Probability Evaluation Method of Cable Corrosion Degree Based on Self-Magnetic Flux Leakage. J. Magn. Magn. Mater. 2021, 522, 167544. [Google Scholar] [CrossRef]
  20. Ni, Y.; Zhang, Q.; Xin, R. Magnetic Flux Detection and Identification of Bridge Cable Metal Area Loss Damage. Measurement 2021, 167, 108443. [Google Scholar] [CrossRef]
  21. Yao, K.; Wang, Z.D.; Deng, B.; Shen, K. Experimental Research on Metal Magnetic Memory Method. Exp. Mech. 2012, 52, 305–314. [Google Scholar] [CrossRef]
  22. Rao, B.P.C. Magnetic Flux Leakage Testing: Basics. J. Non Destr. Test. Eval. 2012, 11, 7–17. [Google Scholar]
  23. Dutta, S.M.; Ghorbel, F.H.; Stanley, R.K. Dipole Modeling of Magnetic Flux Leakage. IEEE Trans. Magn. 2009, 45, 1959–1965. [Google Scholar] [CrossRef]
  24. Lee, J.; Lyu, S.; Nam, Y. An Algorithm for the Characterization of Surface Crack by Use of Dipole Model and Magneto-Optical Non-Destructive Inspection System. KSME Int. J. 2000, 14, 1072–1080. [Google Scholar] [CrossRef]
  25. Lee, J.Y.; Jun, J.W.; Hwang, J.S.; Lee, S.H. Development of Numerical Analysis Software for the NDE by Using Dipole Model. Key Eng. Mater. 2007, 353, 2383–2386. [Google Scholar] [CrossRef]
  26. Le, M.; Lee, J.; Shoji, T. A Simulation of Magneto-Optical Eddy Current Imaging. NDT E Int. 2011, 44, 783–788. [Google Scholar] [CrossRef]
  27. Le, M.; Kim, J.; Vu, H.; Shin, K.; Lee, J. Fast Simulation of Alternating Magnetic Fields around a Crack in a Pipeline System Using a Dipole Model Method. Int. J. Numer. Model. Electron. Netw. Devices Fields 2014, 27, 580–589. [Google Scholar] [CrossRef]
  28. Lee, J.; Choe, E.; Pham, C.-T.; Le, M. Measurement of Thinned Water-Cooled Wall in a Circulating Fluidized Bed Boiler Using Ultrasonic and Magnetic Methods. Appl. Sci. 2021, 11, 2498. [Google Scholar] [CrossRef]
  29. Shi, Y.; Zhang, C.; Li, R.; Cai, M.; Jia, G. Theory and Application of Magnetic Flux Leakage Pipeline Detection. Sensors 2015, 15, 31036–31055. [Google Scholar] [CrossRef] [Green Version]
  30. Le, M.; Lee, J.; Jun, J.; Kim, J. Estimation of Sizes of Cracks on Pipes in Nuclear Power Plants Using Dipole Moment and Finite Element Methods. NDT E Int. 2013, 58, 56–63. [Google Scholar] [CrossRef]
  31. Minkov, D.; Lee, J.; Shoji, T. Study of Crack Inversions Utilizing Dipole Model of a Crack and Hall Element Measurements. J. Magn. Magn. Mater. 2000, 217, 207–215. [Google Scholar] [CrossRef]
  32. Barrans, S.; Tian, G.; Wilson, J.W. Residual Magnetic Field Sensing for Stress Measurement and Defect Detection. In Proceedings of the 9th European Conference on NDT (ECNDT 2006), Berlin, Germany, 25–29 September 2006. [Google Scholar]
  33. Berkache, A.; Lee, J.; Wang, D.; Sim, S. Distribution of Magnetic Flux Density under Stress and Its Application in Nondestructive Testing. Appl. Sci. 2022, 12, 7612. [Google Scholar] [CrossRef]
  34. Kim, J.; Le, M.; Park, J.; Seo, H.; Jung, G.; Lee, J. Measurement of Residual Stress Using Linearly Integrated GMR Sensor Arrays. J. Mech. Sci. Technol. 2018, 32, 623–630. [Google Scholar] [CrossRef]
  35. Roskosz, M.; Rusin, A.; Kotowicz, J. The Metal Magnetic Memory Method in the Diagnostics of Power Machinery Componen. J. Achiev. Mater. Manuf. Eng. 2010, 43, 362–370. [Google Scholar]
  36. Dubov, A.; Kolokolnikov, S. Assessment of the Material State of Oil and Gas Pipelines Based on the Metal Magnetic Memory Method. Weld. World 2012, 56, 11–19. [Google Scholar] [CrossRef]
  37. Cheng, D.K. Fundamentals of Engineering Electromagnetics, 1st ed.; Pearson Education: London, UK, 2014. [Google Scholar]
  38. DiGiuseppe, M.; Howes, C.T.; Speijer, J.; Stewart, C.; van Bemmel, H.; Vucic, R.; Wraight, V. Physics 11 Study Guide, 1st ed.; Nelson: Toronto, ON, Canada, 2012; Volume 1, p. 476. [Google Scholar]
Figure 1. Magnetic domain models; (a) rotational model without and (b) with external magnetic field, (c) concentration model without and (d) with external magnetic field, (e) vertical model without and (f) with external magnetic field.
Figure 1. Magnetic domain models; (a) rotational model without and (b) with external magnetic field, (c) concentration model without and (d) with external magnetic field, (e) vertical model without and (f) with external magnetic field.
Applsci 13 00533 g001
Figure 2. Distributions of magnetization; (a) on the rotational domain model (RDM) without and (b) with external force, (c) on the concentration domain model (CDM) without and (d) with external force, and (e) on the vertical domain model (VDM) without and (f) with external force.
Figure 2. Distributions of magnetization; (a) on the rotational domain model (RDM) without and (b) with external force, (c) on the concentration domain model (CDM) without and (d) with external force, and (e) on the vertical domain model (VDM) without and (f) with external force.
Applsci 13 00533 g002
Figure 4. Simplified magnetization surface current densities on the rotational domain model: (a) without gap, (b) with gap, (c) with gradient gap without deformation, (d) with gradient gap with deformation.
Figure 4. Simplified magnetization surface current densities on the rotational domain model: (a) without gap, (b) with gap, (c) with gradient gap without deformation, (d) with gradient gap with deformation.
Applsci 13 00533 g004
Figure 7. Distribution of magnetic flux densities with rotation domain model, B RDM , c and B RDM , d ; (a) x-, (b) y-, (c) z-directional component.
Figure 7. Distribution of magnetic flux densities with rotation domain model, B RDM , c and B RDM , d ; (a) x-, (b) y-, (c) z-directional component.
Applsci 13 00533 g007
Figure 8. Distribution of magnetic flux densities with concentration domain model, B CDM , d ; (a) x-, (b) y-, (c) z-directional component.
Figure 8. Distribution of magnetic flux densities with concentration domain model, B CDM , d ; (a) x-, (b) y-, (c) z-directional component.
Applsci 13 00533 g008
Figure 9. Distribution of magnetic flux densities with undeformed vertical domain model; (a) x-, (b) y-, (c) z-directional component.
Figure 9. Distribution of magnetic flux densities with undeformed vertical domain model; (a) x-, (b) y-, (c) z-directional component.
Applsci 13 00533 g009
Figure 10. Distribution of magnetic flux densities with deformed vertical domain model, B VDM , d ; (a) x-, (b) y-, (c) z-directional component.
Figure 10. Distribution of magnetic flux densities with deformed vertical domain model, B VDM , d ; (a) x-, (b) y-, (c) z-directional component.
Applsci 13 00533 g010
Figure 11. Specimen configuration. (a) fatigue CT experiment; (b) EDM cutting; (c) gradient fatigue crack.
Figure 11. Specimen configuration. (a) fatigue CT experiment; (b) EDM cutting; (c) gradient fatigue crack.
Applsci 13 00533 g011
Figure 12. Cutting and mounting of the specimens. (a) θ = 0 ° ; (b) θ = 30 ° ; (c) θ = 60 ° .
Figure 12. Cutting and mounting of the specimens. (a) θ = 0 ° ; (b) θ = 30 ° ; (c) θ = 60 ° .
Applsci 13 00533 g012
Figure 13. Microscopic images of each fatigue cracks. (ac) top-view; (df) side-view; (a,d) θ = 0 ° ; (b,e) θ = 30 ° ; (c,f) θ = 60 ° .
Figure 13. Microscopic images of each fatigue cracks. (ac) top-view; (df) side-view; (a,d) θ = 0 ° ; (b,e) θ = 30 ° ; (c,f) θ = 60 ° .
Applsci 13 00533 g013
Figure 14. Experimental setup. (a) experimental set-up on the optical board; (b) 3-axis TMR sensor scanning on a specimen.
Figure 14. Experimental setup. (a) experimental set-up on the optical board; (b) 3-axis TMR sensor scanning on a specimen.
Applsci 13 00533 g014
Figure 15. Distribution of 3-axial TMR signals. (a) X-, (b) Y-, (c) Z-directional component on the 0° specimen; (d) X-, (e) Y-, (f) Z-directional component on the 30° specimen; and (g) X-, (h) Y-, (i) Z-directional component on the 60° specimen.
Figure 15. Distribution of 3-axial TMR signals. (a) X-, (b) Y-, (c) Z-directional component on the 0° specimen; (d) X-, (e) Y-, (f) Z-directional component on the 30° specimen; and (g) X-, (h) Y-, (i) Z-directional component on the 60° specimen.
Applsci 13 00533 g015aApplsci 13 00533 g015b
Figure 16. Distribution of 3-axial TMR signals subtracted background data. (a) X-, (b) Y-, (c) Z-directional component on the 0° specimen; (d) X-, (e) Y-, (f) Z-directional component on the 30° specimen; and (g) X-, (h) Y-, (i) Z-directional component on the 60° specimen.
Figure 16. Distribution of 3-axial TMR signals subtracted background data. (a) X-, (b) Y-, (c) Z-directional component on the 0° specimen; (d) X-, (e) Y-, (f) Z-directional component on the 30° specimen; and (g) X-, (h) Y-, (i) Z-directional component on the 60° specimen.
Applsci 13 00533 g016
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lee, J.; Wang, D.; Dharmawan, I.D.M.O. Theoretical Model of Self-Magnetic Flux Leakage and Its Application in Estimating the Depth Direction of a Fatigue Crack. Appl. Sci. 2023, 13, 533. https://doi.org/10.3390/app13010533

AMA Style

Lee J, Wang D, Dharmawan IDMO. Theoretical Model of Self-Magnetic Flux Leakage and Its Application in Estimating the Depth Direction of a Fatigue Crack. Applied Sciences. 2023; 13(1):533. https://doi.org/10.3390/app13010533

Chicago/Turabian Style

Lee, Jinyi, Dabin Wang, and I Dewa Made Oka Dharmawan. 2023. "Theoretical Model of Self-Magnetic Flux Leakage and Its Application in Estimating the Depth Direction of a Fatigue Crack" Applied Sciences 13, no. 1: 533. https://doi.org/10.3390/app13010533

APA Style

Lee, J., Wang, D., & Dharmawan, I. D. M. O. (2023). Theoretical Model of Self-Magnetic Flux Leakage and Its Application in Estimating the Depth Direction of a Fatigue Crack. Applied Sciences, 13(1), 533. https://doi.org/10.3390/app13010533

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop