In light of the above studies, the fatigue damage parameters for fatigue life prediction are determined by the results of the cyclic rain flow counting. The damage parameters of the SW model and modified FS model are as follows:
3.3.1. Parameters of the SW model
It is indicated by the tension–torsion multiaxial fatigue test that the amplitude of
on the maximum shear plane is very small under proportional loading, whereas the amplitude of
obviously increases with the phase angle under nonproportional loading [
36]. Furthermore, Shang and Wang point out that fatigue crack growth is a decohesion process along the crack tip shear band from the point of view of the micrometer, and the normal strain on the crack plane accelerates this behavior. This means that
and
on the critical plane are two important fatigue damage-controlled parameters, and the size of the normal strain excursion between adjacent turning points is one of the parameters affecting fatigue crack growth. Based on the above ideas, Shang and Wang suggest that an equivalent strain to the von Mises criterion is used as a multiaxial fatigue damage parameter:
where
is the equivalent strain,
is the normal strain excursion between adjacent turning points, and
is the maximum shear strain range on the critical plane.
The multiaxial fatigue damage formula relating to the Manson–Coffin equation is expressed as Equation (10), that is, the SW model:
In the case of proportional loading, Equation (10) reduces the equivalent strain approach form. In the case of uniaxial loading, Equation (10) reduces the Manson–Coffin equation to a uniaxial form. Thus, Shang and Wang proposed that Equation (10) could be used as a unified fatigue damage criterion under either proportional loading (including uniaxial and multiaxial loading) or nonproportional loading.
3.3.2. Parameters of the Modified FS model
In terms of multiaxial fatigue study, researchers have made great progress in the area of damage accumulation and life prediction based on the strain components or the combined strain with stress components on the critical plane. The other popular approach is that proposed by Fatemi and Socie. They introduced a maximum normal stress
on the maximum shear strain plane to reflect the effect of additional cyclic hardening under nonproportional loading [
32]. The FS damage parameter is expressed as
where
is an empirical constant,
is the yield strength, and
,
are elastic Poisson ratio and plastic Poisson ratio, respectively. The other symbols are as previously stated. It is worth noting here that the stress-correlated factor,
, found in Equation (11), is not constant. The maximum normal stress
is determined by the modified Ramberg–Osgood relation [
31]:
where
is the maximum normal strain,
and
are as previously stated, and
is the nonproportional additional hardening coefficient, which depends on the material.
is determined by [
37]
where
,
are the equivalent stress amplitude under 90° out-of-phase loading and the equivalent stress amplitude under proportional loading, respectively.
is the nonproportionality of the strain path determined by [
31]
where
is the swept area of the
polar coordinate space and
is the circle area with a radius of the maximum shear strain during one cycle.
It is believed that Equation (12) can account for the nonproportionality of both material additional hardening and the strain path. However, it is worth noting that not all metallic materials exhibit nonproportional additional hardening characteristics under nonproportional loading. The value of
is zero for materials without a nonproportional additional hardening effect, and the value of the modified coefficient
in Equation (12) is then zero, even under nonproportional loading. This means that Equation (12) is similar to the Ramberg–Osgood relation, i.e., Equation (2). That is to say, Equation (12) cannot accurately reflect the effect of additional hardening when applied to fatigue life prediction for a material without a nonproportional additional hardening effect. To solve this problem, a new modified coefficient
is introduced in Equation (11) to reflect the effect of additional hardening by our research team. In addition, it is suggested that
can be replaced by
. Therefore, our previously proposed modified FS model is expressed as
where
is the maximum normal stress on the critical plane determined by the Ramberg–Osgood relation. The other symbols are as previously stated except
, which is defined as the nonproportional factor, and is given by [
35]
It can be observed from Equation (16) that the value of is zero for materials without a nonproportional additional hardening effect, whereas the effect of nonproportionality of the strain path is not ignored. That is to say, can account for the effect of the nonproportionality of both the strain path and material additional hardening. In other words, the proposed modified FS model for life prediction appears to be relatively accurate in adopting the modified coefficient , and may be used for different materials.
The nonproportionality of the strain path cannot be adequately accounted for if determined by Equation (14). This is because the value of
may be negative or zero when
is less than or equal to
, resulting in limitations of application. To solve this problem, we carry out a second amendment for the nonproportional factor
.
is determined by the introduction of the moment of inertia [
34]. As shown in
Figure 15,
is the random strain path simplified by its enveloping line, i.e., the equivalent convex path
, and the maximum chord length of
is the diameter of the maximum circle path
. It is worth noting that the midpoint of the maximum chord length is used as the origin of the
coordinate system, the x-axis represents the in-phase proportional loading path, and the y-axis represents the antiphase proportional loading path.
is calculated using the following formula:
where
,
, and
are respectively given by
where
is the moment of inertia of
to the x-axis in Equation (18), which indicates the deviation of the point from the in-phase proportional loading path. In other words, the moment of inertia can account for the nonproportional contribution of the strain path. Any point in area
surrounded by
contributes to nonproportional additional hardening, and
represents the distance from any point in area
to the x-axis. Similar to
,
is the moment of inertia of
to the x-axis.
is the circle area with diameter of the maximum chord length.
is used to express the distance from any point in area
to the x-axis in Equation (19). Because an equivalent convex path may appear even under different random strain path loading, the contribution of the random strain path is taken into account by the parameter
.
is the area surrounded by the random strain path
.
and
are the perimeter and the maximum strain range of
, respectively.
The strain path of the dangerous position at the weld zone of NWB_2 is shown in
Figure 16. The dotted line represents the original coordinate position and the new coordinate position is determined by the midpoint of the maximum chord length of the strain path marked with the blue line. Its equivalent convex path and maximum circular path are marked with the red and black lines, respectively.
It can be seen from
Figure 16 that the strain path changes constantly and has obvious randomness, while the determination of
is convenient for engineering applications when the convex hull in the MATLAB program is adopted, and its physical significance is clear.
In light of the above study, the nonproportional parameters in the modified FS model are shown in
Table 6.