Ship To Shore cranes, Rubber Tyred Cranes, Rail Mounted Cranes

Ship board cranes, offshore cranes acc DNV, Lloyds, RMRS, ABS

palf
palf 2

Fatigue Life Assessment Approaches Comparison Based on Typical Welded Joint of Chassis Frame

Abstract

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.

Introduction

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:

  • fatigue 25%;
  • structure collision with a ship 24%;
  • dropping objects 9%;
  • corrosion6%

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:

  • nominal stress approach;
  • hot spot stress approach;
  • notch stress approach;
  • fracture mechanics approach.

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”.

Fatigue test 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].

Crossbeam to longeron T-type welded joint from 93571 ODAZ trailer chassis frame (1 – crossbeam; 2 – longeron)

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

The fatigue curve of Weibull type

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.

Nominal stress in crossbeam vs the number of stress cycles (S-N curve) obtained from fatigue tests

Fatigue life with failure probability of 2.3 % has been calculated using next equation:

Fatigue life with failure probability

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:

Fatigue life with failure probability95%

test machine

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).

Crossbeam stress calculation using finite elements of beam and shell types

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.

Nominal Stress approach

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).

Nominal Stress approach joint classification

Hot spot stress approach

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 finite elements. In this case welding seam is to be create in such ways:
    • model without welding seams;
    • using oblique shell elements to model welding seams;
    • using shell element with increased thickness for welding seams modelling;
  • Solid modelling with volume finite elements. Idealized welding seam shape is 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:

the hot spot stress

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.

 Finite element model

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).

Stress concentrator area meshing. Concentrators are circled by red line

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

Sy stress graphical plots for the boundary conditions

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:

The general equation

where

Δσhs= σhs_maxhs_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

“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).

Welding seam modelling with oblique shell elements

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.

Example of welding seam modelling with oblique shell elements

Fig. 10. Example of welding seam modelling with oblique shell elements

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

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.

Solid model of the crossbeam-longeron welding connection

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

Durability corresponding to different failure probabilities

Notch Stress approach

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:

equation

In addition to the weld toe modelling radius (Fig. 13) the approach specifies the welding seam geometry creation method, finite element size etc.

Welding seam modelling requirements

Fig. 13. Welding seam modelling requirements

Crossbeam-longeron welding connection model for Notch Stress Analysis

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.

Crossbeam-longeron welding connection sub-model

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

First principal stress graphical plot for subassembly

Table 5. Principal stress variation during cycle and durability assessment

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.

Fracture Mechanics based approach

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:

Fracture Mechanics based approach

According to [7] either of two types of the crack growth relationship (Fig.17) could be used.

Crack growth relationship (taken from [7])

Using Eq. 7 the crack length – stress cycle relationship could be obtained:

 stress cycle relationship

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:

  • Define the SIF variation as the approximation  . To do this the models of the joint with different crack lengths are created and for each crack length the SIF is calculated (calculation results are shown in Table 6 and Table 7)
  • Substitute the obtained approximation into the integral equation Eq. 8 and integrate

The initial limit

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.

Crack modelling in the welding seam vicinity. The finite elements with shifted nodes have been used

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

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:

First model for equivalent SIF calculation

Second model for equivalent SIF calculation:

Second model for equivalent SIF calculation:

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

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

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.

Total durability

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.

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

Table 7. Crack growth modelling results

BENDING MODELS

crack length relation

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

Second Model for equivalent SIF calculation with two stage crack growth relationship

Fatigue life assessment results for crossbeam to longeron welding connection using different methods

Discussion of the obtained results

  • It has been found that for the case of Hot Spot stress approach analysis without weld seam modelling the local orientation of 1st principal stress near the gusset plate to main plate connection ends is not perpendicular to the welding seam and that is the reason for using stress component perpendicular to the seam. At the same time for the cases where the welding seam is modelled (both shell and solid models) the 1st principal stress is perpendicular to the welding seam. Thus, the local stress strain state in models without modelled seams does not reflect the reality and the fatigue analysis based on local stress in these areas is not correct.
  • The thickness correction for 4 mm plate, applied with „Hot Spot“ stress approach, when the higher FAT class is used, gives significant overestimation of the joint durability.
  • For the case when weld seam is NOT modelled the lower stress is , but for models with welding seam . As the result the stress range for the models without welding seam is approximately twice bigger than for model with seam modelled.

Conclusion

Having analysed obtained results for crossbeam – longeron welding connection and compared them with the fatigue test following conclusion has been done:

  1. Fatigue life assessment based on Nominal stress approach could be utilized only if the geometry and boundary conditions (type of joint fixation and applying loads) of the analysing joint comply with the one from the existing schemes of the Codes, for which data has been originally obtained by fatigue testing. The biggest problem is that the Codes do not cover all possible types of boundary conditions. For example, in the case of crossbeam-longeron joint analysis this method could not be used because the appropriate loading scheme could not be found in the Standard.
  2. The closest to the fatigue test results are given by the fracture mechanics approach based on equivalent Stress Intensity Factor calculated acc. Eq. 11 in combination with:
    1. one stage crack growth relationship (difference with test is 3.75%; the result is NOT conservative as the calculated durability is more than the test results);
    2. two stage crack growth relationship (difference with test is -16.75 %; the result is conservative as the calculated durability is less than the test results).
  3. The worst correlation with the test shows the “Hot Spot” stress-based approach without the seam modelling.
  4. “Notch Stress” analysis result is close to the one obtained using “Hot Spot” stress analysis.
  5. Regarding to the “Notch Stress” approach its main merit is that only this method among described above could predict the durability for the cases where the crack initiates from the weld root. Thus, sometimes it is only one option for analysis.

Acknowledgement

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.

References  

  • Review of Repairs to Offshore Structures and Pipelines, Publication 94/102 Marine Technology Directorate, UK, 1994.
  • Fricke, Guideline for the fatigue design of welded joints and components XIII-1539-96/XV-845-96, 2010.
  • Niemi, W. Fricke and S.J. Maddox, Fatigue Analysis of Welded Components. Designer guide to the structural hot-spot stress approach, IIW-1430-00, 2007.
  • Fricke, IIW Guideline for the Assessment of Weld Root Fatigue IIW-Doc. XIII-2380r3-11/XV-1383r3-11, 2012
  • Guideline for the Fatigue Assessment by Notch Stress Analysis for Welded Structures IIW-Doc. XIII-2240r2-08/XV-1289r2-08 Wolfgang Fricke, 2010
  • BS 7608-1993 Fatigue design and assessment of steel structures, 1993.
  • BS 7910 Guide on methods for assessing the acceptability of flaws in metallic structures, 1999.
  • RP-203: Fatigue Design of Offshore Steel Structures, DNVGL-RP-0005:2014-06
  • Fatigue assessment of ship structures DNVGL-CG-0129, 2015
  • Baik, K. Yamada, T. Ishikawa, Fatigue Strength of Fillet Welded Joint subjected to Plate Bending, Steel Structures, 8 (2008) 163-169
  • Lassen, N. Recho, Fatigue life analysis of welded structures, ISTE, 2006
  • Comparison of hot spot stress evaluation methods for welded structures J-M. Lee, J-K. Seo, M-H. Kim, S-B. Shin, M-S. Han, June-Soo Park, and M. Mahendran, Inter J Nav Archit Oc Engng, (2010) 2:200~210
  • Lotsberg, Fatigue design of plated structures using finite element analysis, Ships and Offshore Structures, (2006), 45-54
  • General classification creation: Doctor of technical science thesis: 05.02.02/ Konoplyov A. V; Odessa National Polytechnical University. – O., 2013. – 39p. [inUkrainian] (Експериментально-розрахункові методи визначення границі витривалості деталей машин. Створення їх єдиної класифікації [Текст] : автореф. дис. … д-ра техн. наук : 05.02.02 / Конопльов Анатолій Васильович ; Одес. нац. політехн. ун-т. – О., 2013. – 39 с. : рис.)
  • Draper, Modern Metal Fatigue Analysis, Birchwood Park, Warrington : EMAS Pub, 2008.

Measured strain signal as a Boundary conditions for Finite Element submodel with application to crane structure residual life assessment

Abstract

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.

Introduction

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:

  • Crane’s line diagram, where all structure elements are represented with their center lines or point masses; in addition, the elements’ sections geometrical parameters are also given (width, height, web/flange thickness etc.).
  • Crane elements’ weight and the Centers of Gravities coordinates.
  • Material properties of the crane elements’ grades.

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

  • Find out the crane critical elements by combining the results from the crane structural surveying and the fatigue damage check method based on coercive force measurements.
  • Creating FEM models of the structure critical elements using geometrical measurements.
  • Measure the stress fluctuation during crane operation using strain gauges.
  • Process the measured signal and convert it to the FE model nodal forces.
  • Cut the material from the most critical element and test it to find the current, degraded properties of the grade.
  • Calculate the residual life of the critical elements based on crack growth modelling, using the material parameters obtained from the test and the “measured” boundary conditions.

More detailed description of each step is given farther on.

Finding crane structure critical elements

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.

The Grab Ship Unloader general arrangement

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):

  • Boom, the section between Boom/Girder hinge and Boom/forestay hinge.
  • Girder, between Land side and Sea side legs, just above the hopper where the crane unloads its grab.
  • Sea side TGSB middle part, where the pulley support is fixed to it.

Creating FEM models of the structure critical elements

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.

TGSB submodel with crack

Measured stress for FEM submodels’ Boundary Conditions

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.

Box type cross section beam with applied forces and moments

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:           

Normal stress distribution from My bending moment:

Normal stress distribution from My bending moment

Shear stress distribution from Fz shear force:

 Shear stress plot for Fz

Normal stress distribution from Mz bending moment:

Normal stress plot for Mz
Normal stress plot for Mz

Shear stress distribution from Fy shear force:

Shear stress plot for Fy
Shear stress plot for Fy

Stress distribution over the cross section
Stress distribution over the cross section.

Stress superposition for all forces and moments is shown on Fig. 13.

Strain gauge minimum amount and positions

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:

  • the normal stress is represented by polynomial with the maximum order of 1;
  • the shear stress is represented by polynomial with the maximum order of 2.

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.

Minimum amount of the strain gauges for box section type beam
Minimum amount of the strain gauges for box section type beam

Application

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.

Strain gauges’ positions on Boom, Girder and TGSB
Strain gauges’ positions on Boom, Girder and TGSB

Strain gauge signal snippet for TGSB critical element during crane operation
Strain gauge signal snippet for TGSB critical element during crane operation

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.

Beam for strain gauge initial signal calibration
Beam for strain gauge initial signal calibration

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.

Measured signal conversion volt-strain-stress
Measured signal conversion volt-strain-stress

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.

The stress (normal and shear) distribution reconstruction along the plate edge
The stress (normal and shear) distribution reconstruction along the plate edge

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.

Nodal forces calculation based on the reconstructed stress
Nodal forces calculation based on the reconstructed stress

Boom critical element with applied boundary conditions
Boom critical element with applied boundary conditions

Material tests

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

  • Uniaxial tension in order to find the possible changes of yield stress value due to the accumulated fatigue damage.

 Uniaxial test machine and specimen
Uniaxial test machine and specimen

  • Impact toughness test to find the possible embrittlement of the grade due to the accumulated fatigue damage

Impact toughness test machine and specimen
Impact toughness test machine and specimen

  • Fractography for all mentioned above tested specimens

Fracture features of the specimen destruction tested for impact strength (with different resolution)
Fracture features of the specimen destruction tested for impact strength (with different resolution)

It is found that

  • Based on yield stress and impact toughness values there is no significant degradation of plasticity of the grade during the crane operational life.
  • At the same time, fractography found a large number of deep laminations in the steel. The laminations appear due to weakened cohesion between grade matrix elements, nonmetallic inclusions etc. Those laminations significantly decrease the robustness of the crane element and could be a reason of the crane collapse due to uncontrolled fatigue crack propagation.

Critical elements’ residual life assessment based on crack growth FEM analysis

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

14
14

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

Advantages of the approach

  • The proposed approach allows to perform the FEM analysis for the crane elements in the cases when the structural drawings are not available.
  • The measured stress contains all the information about forces that appear during the real operational process (i.e. skew load on the trolley wheels, forces from the grab swaying on ropes, dynamic forces during closing the grab etc.). This allows to calculate the residual life more precise for the considered crane. At the same time for FEA analysis of the whole FE model those forces are unknown and are taken according to the existing standards where certain exaggerations are used to represent the crane worst operational conditions that in real life happen seldom and for the particular crane under consideration could not happen at all.

Disadvantages of the approach

  • The analysis could be done only for the current operational parameters of the crane. For the case when it is needed to calculate the residual life for the case of increased SWL or mechanism accelerations, the FEA of the whole crane and as the result the full package of structural drawings is needed.

References

  • Starikov, Beljatinskij, Perentkovskis, & Klimenko. The use of magnetic coercivity method to diagnose crane metalware. TRANSPORT, 26 (3), (2011), 255-262.
  • A. Starikov, Nikiforov Yu. A. (2012). Residual life assessment for the metallic structures of hoisting machines. Strengs of Materials, 44(1), (2012), 108-113.
  • Starykov, M., Methodology of residual life assessment in the metallic structures of lifting machines. TRANSPORT, 28 (3), (2013), 238-243.
  • Practical evaluation of fatigue and stress state, and residual life of metal by non-destructive method for measuring magnetic characteristic “The coercive force” – A case study / R. Solomakha, R. D. Pittala, G. Bezlyudko, B. V. Baskaran // Asia Pacific Conference on Non-Destructive Testing (14th APCNDT), Mumbai, India, November 18–22, 2013. – P. 1–9.
  • Information on http://www.cae.tntech.edu/~chriswilson/FEA/ANSYS/ANSYS_LEFM02.pdf
  • Information on https://en.wikipedia.org/wiki/Mohr%27s_circle
  • Draper, Modern Metal Fatigue Analysis, Birchwood Park, Warrington : EMAS Pub, 2008.
  • Standard Practice for Cycle counting in Fatigue Analysis. ASTM E 1049-85 (Reapproved 1997). (n.d.).
  • BS 7910 Guide on methods for assessing the acceptability of flaws in metallic structures, 1999.