Search

This is the first volume of the concurring series of Reports in Applied Mechanics, which is based on the outcome of the advanced project course TMPM07 in Applied Mechanics at Link¨oping University during the autumn of 2022. The course lay-up is based on several industrial related projects within the field of Solid Mechanics, concerning fatigue, topology optimisation, structural dimensioning, contacts etc, and Fluid Mechanics, concerning fluid dynamics, flow, aerodynamics, heat transfer etc. The students tackle industry relevant projects in close collaboration with industry from near and neighbouring regions and work in project groups to solve the given tasks within the time limit of the course. Close collaboration with the industry is necessary to define planning, update and feedback for further evaluation at the industry. Three projects were performed during the course of 2022, two within Solid Mechanics and one in Fluid Mechanics. The projects were all performed in tight collaboration with industry partners, and had a close application to real industrial problems. A good opportunity for the students to show-off all their gained knowledge and apply in the best possible way to make innovative solutions in the respective projects. Something they all managed to do with success!, Collection of student reports.

In this paper, the high-temperature constitutive behaviour of an additively manufactured ductile nickel-based superalloy is investigated and modelled, with application to thermomechanical fatigue, low-cycle fatigue and creep conditions at temperatures up to 800∘800∘C. Thermomechanical fatigue tests have been performed on smooth specimens in both in-phase and out-of-phase conditions at a temperature range of 100−800∘100−800∘C, and creep tests at 625∘625∘C, 700∘700∘C, 750∘750∘C and 800∘800∘C. Additionally, low-cycle fatigue tests at different strain ranges and load ratios have been performed at 700∘700, Funding Agencies|Linkoping University and Siemens Energy AB

In this study, a developed co-simulation method, which couples 1D-fluid and 3D-structural models, has been utilised to simulate wear in a hydraulic percussion unit. The effect of wear is generally detrimental on performance and lifetime for such units, but can also cause catastrophic failure and breakdown, requiring a total overhaul and replacement of core components. One experiment of standard straight impact was performed to investigate the tolerance against seizure. The percussion unit was operated at successively increasing operating pressures, and the level of wear was registered at each step, until seizure occurred. The co-simulation model was used to replicate the running conditions from the experiment to simulate the structural response to be used as input for the wear routine to calculate the wear depth. The wear pattern from the simulations corresponds well to the wear pattern from the experiment. Further, the effect of a misaligned impact on wear development was also studied, as this is a loading situation that typically occurs for hydraulic percussion units. The study demonstrates that the simulation method used has a potential for simulating wear and predicting seizure in hydraulic percussion units., Funding: Epiroc Tools & Attachments Division

The crack initiation life of a ductile additively manufactured nickel-based superalloy is studied and modeled for low-cycle fatigue and thermomechanical fatigue conditions up to 600 degrees C. Isothermal experiments were performed on smooth specimens at temperatures up to 600 degrees C with different applied strain ranges. Additionally, thermomechanical fatigue experiments at 100-450 degrees C and 100-600 degrees C were performed on smooth specimens under in-phase and out-of-phase conditions. A life prediction model accounting for the anisotropy was developed, where the temperature cycle is accounted with a Delta T$$ \Delta T $$-functionality, generating good agreements with the experiments. The model was also validated on notched specimens undergoing thermomechanical fatigue conditions at 100-500 degrees C using simplified notch correction methods., Funding Agencies|Linkoping University; Siemens Energy AB

The Ottosen-Stenstrom-Ristinmaa (OSR) incremental fatigue damage model, based on an moving endurance surface centered around a backstress, is adapted for high-cycle fatigue at stress-raisers in AA7050-T7 specimens. Fatigue experiments are carried out for circumferentially-notched, axisymmetric specimens subjected to constant-amplitude (CA) load. The OSR model parameters are fitted to CA fatigue data, showing fair agreement for one set of model parameters across different stress ratios, stress concentration factors, uniaxial stress and biaxial stress. To demonstrate predictive capability, the fatigue life is integrated for an aircraft load spectrum (TWIST), and compared with experimental fatigue life data for holeplate specimens in the literature., Funding Agencies|National Aeronautical Research Program 7 of Vinnova [2019-02778]

The crack growth behaviour of the alloy CM 247 LC is investigated for out-of-phase TMF and isothermal tests at the same temperature as the minimum temperature in the TMF tests. The results suggest that it is possible to characterise crack growth behaviour if experimental corrections for crack closure are accounted for. The repli -cation of these experimental tests using a numerical FE-solver results in similar crack growth behaviour, sug-gesting that the main mechanism in place is plasticity-induced crack closure. A pragmatic analytical model to characterise crack closure including hold time effects is proposed. The comparison of the response from this model with the experimental and numerical results suggests that the proposed analytical model is capable to approximate crack closure effects for cases where substantial creep deformation is to be expected., Funding Agencies|Link oping University; Siemens Energy

Crack closure effects are known to have a large impact on crack growth behaviour. In this work, tests were performed on Inconel 792 specimens under TMF loading conditions at 100–850 °C with extended hold times at tensile stress. The effective stress-intensity range was estimated experimentally using a compliance-based method leading to the conclusion that crack closure appears to have a primary impact on the crack growth behaviour for this material under the conditions studied. The crack closure behaviour for the tests was successfully modelled using numerical simulations, including creep., Ytterligare forskningsfinansiär: Siemens Industrial Turbomachinery AB through “Turbines for Future Energy Systems” (Turbiner för framtidens energisystem), Grant No. 44100-1

In this study, a previously developed co-simulation method has been expanded to also simulate the dynamic behaviour of sealing gap regions in hydraulic percussion units. This approach is based on a 1D system model representing the fluid components and a 3D finite element model representing the structural parts of a hydraulic hammer. The sealing gap is a fundamental feature of a percussion unit, where the reciprocating motion of the piston is generated by the valve mechanism of the sealing gap. When the gap is closed it will prevent fluid flow between regions of different pressure levels. However, a small leakage flow through the gap will always occur which size depends on the clearance and the position of the piston. The method proposed here will take the structural motion and deformation into consideration when calculating the leakage flow. The deformed state of the structure is approximated by a cylindrical surface, in a least square manner, and communicated through the co-simulation interface to the fluid simulation module, and then used when calculating the leakage flow. This method aims at a more accurate simulation of the leakage flow that will not only yield a more realistic description of the mechanism on the local level, but also a more accurate estimation of global parameters such as overall performance and efficiency. The results indicate that the simulated leakage flow will decrease when dynamic gaps are used in comparison to static gaps, which is a consequence of the deformed structure that will generate smaller clearances. The leakage flow for the dynamic gaps will even be lower than for the static perfectly concentric case, mainly due to the reduction of clearances. The results also indicate that the dynamic eccentricity does not have a major influence on the leakage flow. The outcome from this study highlights the potentials of the described co-simulation approach for analysing the dynamics of the sealing gaps in a hydraulic percussion unit (i.e. gap

The current paper presents work on identification and evaluation of a range of factors influencing accuracy and comparability of data generated by three laboratories carrying out stress-controlled thermo-mechanical fatigue crack growth tests. It addresses crack length measurements, heating methods and temperature measurement techniques. It also provides guidance for pre-cracking and use of different specimen geometries as well as Digital Image Correlation imaging for crack monitoring. The majority of the tests have been carried out on a coarse grain polycrystalline nickel-base superalloy using two phase angles, Out-of-Phase and In-Phase cycles with a triangular waveform and a temperature range of 400-750 degrees C., Funding Agencies|European UnionEuropean Union (EU) [686600]

Single-crystal nickel-base superalloys frequently experience two distinct fatigue crack growth modes. It has been observed that, under certain conditions, cracks transition from a path perpendicular to the loading direction to a crystallographic slip plane. As crystallographic cracking is associated with an increased fatigue crack growth rate, it is important to be able to predict when this transition occurs. In this work three different criteria for crystallographic cracking based on resolved anisotropic stress intensity factors are evaluated in a three-dimensional finite element context. The criteria were calibrated and evaluated using isothermal fatigue experiments on two different specimen geometries. It is suggested by the results, that a threshold value of a resolved shear stress intensity factor can act as a conservative criterion indicating cracking mode transition. Further, a trend hinting towards a loading frequency dependency could be observed., Funding Agencies|Linkoping University; Siemens Industrial Turbomachinery AB

DevTMF

Crack propagation rates from isothermal and in-phase high temperature thermomechanical fatigue crack propagation tests for Inconel 792 resulted in a rather large scatter band when plotted against the linear-elastic stress intensity function Delta K. Previous observations indicated that the residual strains generated during the first few cycles in the uncracked specimen can explain, to a large extent, the differences between the tests. An effective Delta K-eff was calculated for each test based on non-linear FE simulations. When the crack propagation rates for each test were plotted against the calculated Delta K-eff, the curves collapsed to a smaller scatter band., Funding Agencies|Swedish Energy AgencySwedish Energy Agency; Siemens Industrial Turbomachinery AB through "Turbines for Future Energy Systems" (Turbiner for framtidens energisystem) [44100-1]

The polycrystalline nickel-base superalloy RR1000 is used as turbine rotor material in Rolls-Royce aero engines and has to withstand a wide variety of load and temperature changes during operation. In order to maximize the potential of the material and to improve component design, it is of great interest to understand, and subsequently be able to accurately model the crack propagation caused by thermo-mechanical fatigue conditions. In this work, experimental data is analysed and used to inform unified modelling approaches in order to predict the crack propagation behaviour of RR1000 under a variety of stress-controlled thermo-mechanical fatigue conditions., Funding Agencies|European Unions Horizon 2020 research and innovation programme and Joint Undertaking Clean Sky 2 [686600]

The computational efficiency in analysing cyclically loaded structures is a highly prioritised issue for the gas turbine industry, as a cycle-by-cycle simulation of e.g. a turbine disc is far too time consuming. Hence, in this paper, the efficiency of two different procedures to handle computational expansive load cases, a numerical extrapolation and a parameter modification procedure, are evaluated and compared to a cycle-by-cycle simulation. For this, a local implementation approach was adopted, where a user-defined material subroutine is used for the cycle jumping procedures with good results. This in contrast to a global approach where the finite element simulation is restarted and mapping of the solution is performed at each cycle jump. From the comparison, it can be observed that the discrete parameter modification procedure is by margin the fastest one, but the accuracy depends on the material parameter optimisation routine. The extrapolation procedure can incorporate stability and/or termination criteria., Funding Agencies|Cleansky [686600]

Cracks in single-crystal nickel-base superalloys have been observed to switch cracking mode from Mode I to crystallographic cracking. The crack propagation rate is usually higher on the crystallographic planes compared to Mode I, which is important to account for in crack growth life predictions. In this paper, a method to evaluate the crystallographic fatigue crack growth rate, based on a previously developed crystallographic crack driving force parameter, is presented. The crystallographic crack growth rate was determined by evaluating heat tints on the fracture surfaces of the test specimens from the experiments. Complicated crack geometries including multiple crystallographic crack fronts were modelled in a three dimensional finite element context, The data points of the crystallographic fatigue crack growth rate collapse on a narrow scatter band for the crystallographic cracks indicating a correlation with the previously developed crystallographic crack driving force., Funding Agencies|Linkoping University; Siemens Industrial Turbomachinery AB

In this study, a previously developed co-simulation approach has been adopted to simulate the responses of an existing hydraulic hammer product. This approach is based on a 1D system model representing the fluid components and a 3D finite element model representing the structural parts of the hydraulic hammer. The simulation model was validated against four experiments with different running conditions. The corresponding set-ups were analysed using the co-simulation method in order to evaluate the overall responses. A parameter study was also performed involving the working pressure and the restrictor diameter, with the objective to validate that a parameter change in the simulation model will affect the input and output power in the same direction as in the experiments. The experimental responses used in the validation were time history data of fluid pressure, component position and acceleration, and structural stresses. The experiments result in high frequency and high amplitude excitations of the hydraulic hammer and thus require a model with a high resolution of the model dynamics. The conclusion of the validation is that the simulation model is able to replicate the experimental responses with high accuracy including the high frequency dynamics. The favourable outcome of the validation makes the described co-simulation approach promising as an efficient tool for a wide range of other applications where short time duration mechanisms need to be studied., Funding Agencies|Epiroc Hydraulic Attachment Tools Division

Fatigue crack propagation tests under both isothermal and non-isothermal thermomechanical fatigue conditions have been performed on wrought Haynes 230, a ductile combustor material. A number of specimens were thermally aged by pre-straining and subsequent furnace exposure for 3000 h at 600 degrees C. The tests were performed both under load and strain control, between room temperature and 600 degrees C. The thermally aged notched specimens show a decrease in the crack initiation life, similar to results previously reported for smooth test specimens at room temperature. For the crack growth rates, the effects of thermal ageing were less pronounced than for crack initiation. Further, the tests have been simulated using the finite element method to calculate the crack driving force, where the plasticity induced crack closure is handled with a full history description. A temperature dependent linear kinematic hardening plasticity law has been adopted for describing the material behaviour between room temperature and 600 degrees C. A post-processing tool was used in which the plasticity induced crack opening level was calculated, followed by a calculation of the effective Delta J range for each crack length. The adopted procedure yields good correlation between the different tests, under both isothermal and non-isothermal conditions., Funding Agencies|Siemens Industrial Turbomachinery AB, Finspang, Sweden

The inherent anisotropy of single-crystal nickel-base superalloys brings many difficulties in terms of modelling, evaluation and prediction of fatigue crack growth. Two models to predict on which crystallographic plane cracking will occur is presented. The models are based on anisotropic stress intensity factors resolved on crystallographic slip planes calculated in a three-dimensional finite-element context. The developed models have been compared to experiments on two different test specimen geometries. The results show that a correct prediction of the crystallographic cracking plane can be achieved. This knowledge is of great interest for the industry and academia to better understand and predict crack growth in single-crystal materials., Funding Agencies|Swedish Energy Agency; Siemens Industrial Turbomachinery AB through the Research Consortium of Materials Technology for Thermal Energy Processes [KME-702]

Single-crystal nickel-base superalloys are often used in the hot sections of gas turbines due to their good mechanical properties at high temperatures such as enhanced creep resistance. However, the anisotropic material properties of these materials bring many difficulties in terms of modelling and crack growth prediction. Cracks tend to switch cracking mode from Mode I cracking to crystallographic cracking. Crystallographic crack growth is often associated with a decrease in crack propagation life compared to Mode I cracking and this must be taken into account for reliable component lifing. In this paper a method to evaluate the crystallographic crack propagation rate related to a crystallographic crack driving force parameter is presented. The crystallographic crack growth rate is determined by an evaluation of heat tints on the fracture surface of a specimen subjected to fatigue loading. The complicated crack geometry including two crystallographic crack fronts is modelled in a three dimensional finite element context. The crack driving force parameter is determined by calculating anisotropic stress intensity factors along the two crystallographic crack fronts by finite-element simulations and post-processing the data in a fracture mechanics tool that resolves the stress intensity factors on the crystallographic slip planes in the slip directions. The evaluated crack propagation rate shows a good correlation for both considered crystallographic cracks fronts., Funding Agencies|Swedish Energy Agency; Siemens industrial Turbomachinery AB through Research Consortium of Materials Technology for Thermal Energy Processes [KME-702]