







There are many approaches to the durability calculation that are used in engineering practice. At the same time the existing accident studies show that the leading position is still hold by fatigue failures. This means that there is still no universal approach to fatigue problem solution, and the existing approaches have their limitations. In addition, there is lack of information about the comparison between the precision of the obtained results using different approaches. In this paper different fatigue life calculation methods, like nominal stress, hot spot stress, notch stress and fracture mechanics are used to calculate the durability of T-type welded joint. The obtained results are compared with the fatigue test ones and the approaches, which give the closest results, are found.
Time varying working loads are typical for metal constructions of chassis frames, material handling machines, ship hulls etc. According to accident studies for offshore structures [1], that took place in the North Sea, for period from 1972 to 1992, all reasons have been split into several groups according to their significance:
In spite of the existence of different guides and approaches that have being used for fatigue design the significant part of failures caused by fatigue reveals the imperfection of using analysis methods. That is why the development of a new methodology is the pressing issue.
Modern fatigue design approaches are based on stress information about designing joint received from the finite element analysis of a structure. This gives the possibility of using the local stress in the probable area of the fatigue crack appearance instead of using nominal stress in the joint and broadens horizons for further enhancements.
Metal fatigue phenomena has been attracting a lot of researchers‘ interest for a long time and with the welding invention this interest even increased. The main problem was that all of researches solved particular problems (i.e. the effect of mean stress on the durability etc.) but there was no general practical approach with thorough step by step recommendations for the practicing engineers how to perform the analysis. The situation is changed during last decade when International Institute of Welding [2, 3, 4, 5], British Standard [6, 7], DNV [8, 9] have represented researches that are summarized in particular guides for the fatigue analysis with detailed description of practical utilization of the approaches, starting from mesh description and finishing with recommendations about what type of S-N curve to use.
With the aforementioned guides in the place the question of the analysis result validation has appeared. Thus, many researches have their goal to compare the fatigue experiment and analysis results [10, 11, 12, 13]. The main problem in our opinion is that in those researches only one method of the analysis is compared with the test results. But at the same time in engineering practice at least four of them are frequently used:
In this paper the comparison between main analytical approaches and test results for the fatigue life assessment has been done. This comparison could help to the practicing engineer to decide which approach to the durability analysis is more accurate for designing of similar joints.
For the analysis the T-type welded joint (Fig. 1) is chosen. Despite the fact that this type of connection is typical for a chassis frame, it is not covered in the researches. All the existing analysis, done for the T weld connection [10, 12, 13], have their welded gusset plate serving for stress concentration purpose only, when in the T-weld connection that is studied, the force and moment are transmitted to the main plate (crossbeam) through the gusset plate (longeron).
In the following chapters the durability of the joint is obtained using testing and different analysis approaches. The results are discussed in chapter “Discussion of the obtained results”.
The article objective is to define the approaches that give the closest result of fatigue life assessment to ones taken from fatigue test for T-type welded joint of a chassis frame [14].

Specimens have been tested using symmetric stress cycle (R= -1). The crossbeam was fixed using 4 holes of 10 mm in diameter and the 2 forces were applied using the 2 holes of 14 mm in diameter in longeron.
The fact of the crossbeam vertical deformation amplitude increasing beyond 30 % has been used as a collapse criterion to stop the fatigue tests.
The six joints have been tested on 6 different stress levels (Table 1). The fatigue curve of Weibull type has been used
![]()
Using linear interpolation the following parameters of Eq. 1 have been found: mw = -2.489; Cw=3.3319.
Table 1. Fatigue test results for T-weld joint crossbeam to longeron

Based on Eq. 1 the fatigue life for stress amplitude with 50% failure probability is 425 100 cycles.

Fatigue life with failure probability of 2.3 % has been calculated using next equation:
![]()
where d – standard deviation amount below mean value, zP=2.3%= zP=97.7%=2 (quantile for failure probability of 2.3%)
lg σN – standard deviation of lgN, 0.178, p. 20 [2] for the specimen amount n<10.
Fatigue life with failure probability of 97.7 % has been calculated using next equation:
![]()

Fig. 3. Test machine
Traditionally beam theory for nominal stress calculation is used for S-N curve. But that stress is not representative for current joint because the fracture happens not in the crossbeam outer layers but in the area of welding seam transition to the longeron (Areas 1 and 2, Fig. 1).

Fig. 4. Crossbeam stress calculation using finite elements of beam and shell types
Using the shell finite elements gives realistic results. Maximum stress in crossbeam for the beam finite element (Fig 3. c) is 81.5 MPa, and for shell finite element (Fig 3. f) is 159 MPa.
Moreover, stress state of crossbeam in the area of welding seam is not more uniaxial one but complex i.e. all three principal stresses have non zero magnitudes.
The first step of nominal stress analysis [6] is to find among the variety of joint types with boundary conditions (showed in standard) the one that corresponds to the designing joint. But for currently calculating T-type welded connection the similar joint type does not exist. For the first look joint 5.3 (class F2, Fig. 5. a.), clause 2, Table 1, [6] could be taken, but its boundary conditions are different from analysing connection: unlike to the join from the standard the gusset plate (longeron) does not takes any load. That is why it can not be used further on. The joint on Fig 5. b can not be used for calculating either, because its boundary conditions differ from designing joint’s ones. It is also not clear stress in which element is taken for nominal (loading scheme is not shown).

This approach [3] allows calculating the joint fatigue life using its stress-strain state data obtained from the finite element analysis. The following joint modelling techniques are suggested to be used:
Modelling using shell elements
Model without welding seams
According to IIW Recommendations [3] welded element durability is to be calculated based on stress that acts in the weld toe. However, because of using linear elastic metal behaviour and the fact that the real weld profile is unknown on design stage, there is no possibility to use directly the stress read from welding toe. Instead, it has been proposed to use stress extrapolated value based on stress in the welding seam vicinity, so called Structural Stress.
For our case (model consists of 4 node linear shell finite elements with edge of 1.6 mm near the stress concentration point) the hot spot stress is given by:
![]()
where σ0.4·t – stress value at the distance of 0.4·t from the weld toe (the first extrapolation point);
σ1·t – stress value at the distance of 1·t from the weld toe (the second extrapolation point);
t – longeron thickness, 4 mm.
The finite element model of T-welded connection is shown in Fig. 6.

Fig. 6. Finite element model
The minimum thickness of the plate the approach is applicable for is 5 mm.
Area of the stress concentration has been meshed using two techniques (Fig. 7).

Results of finite element analysis are shown on Fig. 8; hot spot stress extrapolation calculation is put into Table 2.

In currently overlooking standard the fatigue life assessment is based on principal stress with the biggest range during cycle. However, if the angle between this stress direction and normal to the welding seam line is more than 60 degrees, the stress perpendicular to the welding seam must be used. In our case Sy is used.
Hot spot stress approach is much easier to use in comparison with the nominal stress approach because it is based only on two S-N curves to assess the fatigue life in a “hot spots”. They are known as FAT 90 and FAT 100. The numbers that come after letters „FAT“ indicate stress level in MPa that corresponds to fatigue life for 2·106 cycle durability. The general equation for these S-N curves is as follows:
![]()
where
Δσhs= σhs_max-σhs_min – stress range in the «hot spot», σhs_max – maximum hot spot stress of a cycle, σhs_min – minimum hot spot stress of a cycle;
m – index of power, 3.0;
С – coefficient, 2·1012;
N – life cycle
Table 2. “Hot spot” stress approximation and durability assessment

* Durability corresponding to different failure probabilities than other than 2.3% are calculated acc. Eq. 2 and Eq. 3.
Plane model with shell finite elements. Welding seam is modelled by oblique shell elements
The main concept of welding seam modelling is shown in Fig. 9 and meshed model – in Fig. 10 a).

Fig. 9. Welding seam modelling with oblique shell elements
For this case first principal stress is perpendicular to the welding seam. That is why it is used for the further analysis.

Fig. 10. Example of welding seam modelling with oblique shell elements
Table 3. “Hot spot” stress approximation and durability assessment

Solid model with volume finite elements
Solid model of the crossbeam-longeron welding connection is shown in Fig. 11. To reduce the computation time during model stress analysis only one half of the model has been created. 20 node Solid finite element with decreased integration and edge size of 4 mm is used.
The distances from the weld toe to the extrapolation points are the same (0.4·t to the first (nearest to weld) extrapolation point and 1·t to the second extrapolation point). Stress analyses result is shown in Fig. 12.

Table 4. “Hot spot” stress approximation and durability assessment

This approach [4, 5] demands solid model creation and volume finite element mesh using. For the plate thickness less than 5 mm the Notch Radius of 0.05mm instead of 1mm has to be used, Special attention must be paid to a weld seam modelling particularly in the area where welding seam material merges to the main metal (Fig. 14 b) because the stress in this area is used for the fatigue life estimation. Only one S-N curve uses for this analysis (FAT 630) which equation takes a form of:
![]()
In addition to the weld toe modelling radius (Fig. 13) the approach specifies the welding seam geometry creation method, finite element size etc.

Fig. 13. Welding seam modelling requirements

Fig. 14. Crossbeam-longeron welding connection model for Notch Stress Analysis
Due to the high level of detail needed for welding area modelling the scope of problem increases with the growth of the joint complexity. That is why calculation time could increase from i.e. 20 minutes to several days. In this case the Sub-modelling feature is very useful. It helps to create more dense mesh and retrieve more precise solution for the smaller part of a model.
For crossbeam-longeron joint welding seam area sub-model of a fatigue crack initiation is shown in Fig 15.

Fig. 15. Crossbeam-longeron welding connection sub-model

Table 5. Principal stress variation during cycle and durability assessment

* Durability corresponding to different failure probabilities than other than 2.3% are calculated acc. Eq. 2 and Eq. 3. According [5] standard deviation of the lgN=0.206.
The central idea of the approach [2, 3] consists in the using Paris equation for assessment of the joint fatigue stress cycles number till failure:

According to [7] either of two types of the crack growth relationship (Fig.17) could be used.
![Crack growth relationship (taken from [7])](https://i0.wp.com/crane-fem-consulting.com/wp-content/uploads/2019/06/Crack-growth-relationship-taken-from-7.jpg?resize=506%2C327&ssl=1)
Using Eq. 7 the crack length – stress cycle relationship could be obtained:

After solving integral equation Eq. 8 the stress cycle number could be defined (N=N2-N1) that is needed for crack growth from length 2a1 to 2a2.
As per fracture mechanics theory a crack starts to grow if SIF range exceeds some threshold value (), which is different for different grades. Only SIF ranges more than this threshold are considered in analysis.
The failure criterion for the fatigue testing of the crossbeam-longeron welding connection is the 30% of longeron deformation range increasing. This corresponds to the crack length of L=2a=35.5 mm.
The method of solving Eq. 8 is as follows:

The initial limit, a1 corresponds to SIF threshold value of the material (170 for R=-1, acc. (48 c), 8.2.3.6 [7]). Final limit, a2=17,75 mm comes from the failure criterion during test.
As the life of crack initiation for welded joints is a small part of the total life [15], we will neglect it. The minimum crack length is defined for each case based on threshold SIF.

After analysis it became clear, that SIFs for all three modes are nonzero. Next Eq. 10 and Eq. 11 have been used to calculate the effective SIF, corresponding to the complex loading, that takes into consideration SIFs for all three different modes. Linear elastic material model has been used.
Table 6. Crack growth modelling results

As all three SIF are not equal to 0 the equivalent SIF has to be used for further analysis.
First model for equivalent SIF calculation:

Second model for equivalent SIF calculation:

First Model for equivalent SIF calculation with one stage crack growth relationship

The SIF approximation is shown on the graph above as a trend line equation

where m – index of power, 3, clause 8.3.3.5, [7],
А – coefficient of proportionality, 5.21·10-13, clause 8.3.3.5 [7],
a1 for this case equals to 0.9 mm.
First Model for equivalent SIF calculation with two stage crack growth relationship
Total durability would consist of durability for two stages (stage A and stage B). For the Mean Curve (Table 10 [7]) the stage A/Stage B transition point is 196, which corresponds to a=1.15 mm.

where A1=4.8·10-18, m1=5.1, A2=5.86·10-13, m2=2.88,
For the Mean Curve + 2SD (Table 10 [7]) The stage A/Stage B transition point is 144, which is smaller than the threshold value and that why during the Stage A the crack will not propagate.

Table 7. Crack growth modelling results


Fig. 20. Approximation of SIF range vs. crack length relation (the polynomial approximation is shown above the trend line)
The SIF approximation is shown on the graph above as a trend line equation



Having analysed obtained results for crossbeam – longeron welding connection and compared them with the fatigue test following conclusion has been done:
I would like to express my gratitude to the following my colleagues from Liebherr Mining Equipment for their valuable comments to may paper: James Witfield PE, Dr. Vladimir Pokras, Michael Karge.
Special appreciation is to my teacher, Doctor of technical science, Prof. Konoplyov A. V., for his permission to use the fatigue test results he has carried out.
The general flowchart of a crane structure residual life assessment could be split into two basic stages. On the first stage crane structure critical elements are determined usually using the Finite Element Analysis (FEA) for the given operational parameters (i.e. safe working load, maximum outreach and back reach, speed of mechanisms etc.). On the second stage the fatigue life of the critical elements is calculated using the same Finite Element (FE) model.
The key factor for such approach is the availability of the crane structural drawings, which allows to create precise FE model. And here the main problem appears that to find the drawings for the old machines is sometimes either a big problem or even impossible (i.e. the manufacturer company does not exist anymore). Similar problem could be challenging even for new cranes. The approach described in this paper allows to overcome the drawing absence problem by using the following modifications of the aforementioned procedure: critical elements are found using Non Destructive Testing (NDT); FE models of the critical areas are created using local structure measurements; FE analysis is performed using strain gauge measurements for Boundary Conditions during the crane operation.
The main flow of the proposed approach is that the residual life could be calculated for the crane operational parameters, the measurements have been performed for. Despite the paper uses crane structure as an example, the approach is general and it could be transferred on different types of steel structures.
In the engineering practice the mostly used approach for stress analysis is Finite Element Analysis (FEA). In order to perform the FEA for the crane structure the following information about the crane must be presented:
This information could be taken from the structural drawings of the crane. Normally the crane owner may have the drawings, but in reality, it is quite hard to find the drawings especially for the machine that has been in operation for more than 20 years. Of course, one more way to find this information is to contact the crane manufacturer. However, it could lead at least to additional expenses to buy those drawings or sometimes the company does not even exist anymore and it is impossible to enquire the documents.
In this paper the approach that allows to overcome the structural drawing absence problem is proposed. As an example, the approach is used for 38-year-old Quay Grab Unloader residual life estimation as the FEA is a basic part of such type of project.
The idea of the approach is to find several critical elements of the crane, create the FE submodels for each of them based on geometrical measurements, performed on site. And instead of taking the boundary conditions (the forces that represent the action of the rest of the crane on the submodel) for these submodels from the whole crane Finite Element Analysis, get them from the processing of strain gauge signal measured during crane operation.
The residual life assessment approach for the considered Grab Unloader consists of the following steps
More detailed description of each step is given farther on.
The crane critical elements have been found using the following two methods separately: structural survey; crane structure accumulated fatigue damage measurement. The obtained results have been analyzed and combined.
Structural survey. During the structural surveying the fatigue cracks of significant length have been found in the sea side Trolley Girder Support Beam (TGSB), near its center, where pulleys are attached Fig. 1. Thus, the critical element based on structural survey is TGSB central part.

Crane structure accumulated fatigue damage measurement. The approach is based on measured coercive force parameter. The measured coercive force is in proportion with the fatigue damage accumulated in the element. For each grade exists the coercive force critical magnitude that shows that the further operation is dangerous and could lead to structure failure. For the simplicity the whole range of the coercive force is split into several intervals: Reliable operation, Controllable operation, Critical Operation [1, 2, 3, 4]
The method could be used for the residual life estimation for the tested grades only, where the critical value of the coercive force is found. For the Grab Unloader grade the study of changing coercive force vs number of stress cycles does not exist. That is why only comparison analysis that shows the critical elements could be done instead of finding the residual life directly for each element. The critical elements are the ones with the biggest coercive force magnitude.
Based on the measured coercive force parameter it appears that the critical elements are Boom and Girder.
Combining the surveying results the following critical elements have been taken for the further residual life assessment (see Fig. 2):

Based on the geometrical measurements performed on site (i.e. member cross section height, width, flange and web thicknesses, bulkheads’ and stiffeners’ geometry etc.) the FE models of the critical elements have been created. As the further analysis is planned to utilize the Linear Fracture Mechanics approach one or two cracks have been introduced acc. [5] to each submodel (there are several models for each critical element where the only difference is the crack length). The cracks’ positions are based on preliminary analysis of the element without crack and the cracks were introduced in the areas with the biggest stress range. Each crack is oriented to be normal to the first principal strain. The linear material behavior model is used.



Concept. Firstly, we have to find out the normal and shear stress distribution in a cross section which is used for the further boundary condition application. Then we could decide about the minimum amount of the strain gauges, their axis directions and positions along the section.
Stress distribution in cross section. For the normal and shear stress distribution the thin wall section beam theory is used. In crane structures the most frequently used types of cross sections are box and I beam types. As soon as all three critical elements of the Grab Unloader have box cross section type, all further analysis is done for this particular cross section type. This approach could be easily used for another cross section type (i.e. I beam). The typical box cross section with internal forces and moments is shown in Fig. 6.

The appropriate stress type (normal and/or shear) is calculated for each load component and the resultant stress distribution is found as a stress superposition. The following rules for the stress signs are used acc. [6]:
For normal stress: tension stress is assumed to be positive; compression stress is assumed to be negative.
For shear stress: positive shear stress acts on positive faces of the material element in the positive direction of an axis. Also, positive shear stress acts on negative faces of the material element in the negative direction of an axis. A positive face has its normal vector in the positive direction of an axis, and a negative face has its normal vector in the negative direction of an axis.
Normal stress from normal to the cross section force Fx:






Stress superposition for all forces and moments is shown on Fig. 13.
On this stage we could decide about the minimum amount of strain gauges we need to reconstruct the stress distribution based on measurements. Having analyzed the stress distribution it could be concluded that:
In other words, we need to have measurements from at least two point along each edge to reconstruct the normal stress distribution and at least three points to reconstruct shear stress distribution.

The stress has been measured during 26 crane operational cycles using the modern 64-channel data acquisition device, Fig. 15. The device has several position sensors that were located along the Boom and Girder to show the trolley position during the measurements.

The strain gauges have been attached from the inner side of the structure. The gauges, positioned along and transversely to the critical element axis, were used; the positions are shown in Fig. 16 by arrows that are the gauge center lines.


The snippet of measured signals from the strain gauges vs time are shown in Fig. 17.
The calibration coefficient to convert measured signal in volts to strain has been taken from laboratory tests (for each channel) using the specially designed beam (Fig. 18) with varying width to provide constant strain all over its length.

On the first step strain gauge signal in volts was measured for test weights of 0, 0.5 kg, 1 kg, 1.5 kg and 2 kg. Then the corresponding strain magnitudes were calculated for all applied test loads using Strength of Material approach. This allowed to plot point in strain vs voltage axes.
On the second step the calibration coefficient was calculated using linear interpolation for strain vs voltage signal plot.
Using calibration coefficient Initial signal from the strain gauge (electrical voltage) has been converted to strain (Fig. 19, a; Step 1), then to stress (Fig. 19, b; Step 2), filtered out [7] and the stress cycles were extracted using Range-Pair Counting Method algorithm [8], Fig. 19.

Based on the stress at the points where the strain gauges are placed the stress distribution diagram along each plate edge has been reconstructed for normal (Step 3) and shear (Step 4) stresses, Fig. 20.

Finally, knowing number of FE along each edge, plate thickness and FE edge length, the stress to the FE edge has been converted to nodal forces, Fig. 21. In order to automatize the aforementioned steps C++ based programs and ANSYS APDL language macros have been used.


The following material tests have been performed to find out the degraded properties of the Ship Unloader grade:



It is found that
The residual life has been calculated acc. [9] based on linear fracture analysis, which utilizes Paris low:

The cycles with SIF amplitude which is more than 0 have been considered in analysis. Also cycles with compression stress in the crack tip have been excluded.
Having analyzed all three submodels (Boom, Girder and TGSB) it was found that the element with the lowest value of the residual life is the TGSB with the critical area of the sheave support flange connection. The residual life of the unloader after repair, based on the fatigue analysis of its critical element, is approximately 9.7 years. It has been assumed that the crane will work with the same regime (same Safe Working Load, mechanisms accelerations etc.).
Summary
Disadvantages of the approach
References