Behaviour Of Materials Research Articles

Comparison of elements and state-variable transfer methods for quasi-incompressible material behaviour in the particle finite element method

AbstractThe Particle Finite Element Method (PFEM) is attractive for the simulation of large deformation problems, e.g. in free-surface fluid flows, fluid–structure interaction and in solid mechanics for geotechnical engineering and production processes. During cutting, forming or melting of metal, quasi-incompressible material behaviour is often considered. To circumvent the associated volumetric locking in finite element simulations, different approaches have been proposed in the literature and a stabilised low-order mixed formulation (P1P1) is state-of-the-art. The present paper compares the established mixed formulation with a higher order pure displacement element (TRI6) under 2d plane strain conditions. The TRI6 element requires specialized handling, involving the deletion and re-addition of edge-mid-nodes during triangulation remeshing. The robustness of both element formulations is analysed along with different state-variable transfer schemes, which are not yet widely discussed in the literature. The influence of the stabilisation factor in the P1P1 element formulation is investigated, and an equation linking this factor to the Poisson ratio for hyperelastic materials is proposed.

Effects of Rubber Core on the Mechanical Behaviour of the Carbon-Aramid Composite Materials Subjected to Low-Velocity Impact Loading Considering Water Absorption.

The large-scale use of composite materials reinforced with carbon-aramid hybrid fabric in various outdoor applications, which ensures increased mechanical resistance including in impact loadings, led to the need to investigate the effects of aggressive environmental factors (moisture absorption, temperature, thermal cycles, ultra-violet rays) on the variation of their mechanical properties. Since the literature is still lacking in research on this topic, this article aims to compare the low-velocity impact behaviour of two carbon-aramid hybrid composite materials (with and without rubber core) and to investigate the effects of water absorption on impact properties. The main objectives of this research were as follows: (i) the investigation of the mechanical behavior in tests for two impact energies of 25 J and 50 J; (ii) comparison of the results obtained in terms of the force, displacement, velocity, and energy related to the time; (iii) analysis of the water absorption data; (iii) low-velocity impact testing of wet specimens after saturation; (iv) comparison between the impact behaviour of the wet specimens with that of the dried ones. One of the main findings was that for the wet specimens without rubber core, absorbed impact energy was 16% less than the one recorded for dried specimens at an impact energy of 50 J. The failure modes of the dried specimens without rubber core are breakage for both carbon and aramid fibres, matrix cracks, and delamination at matrix-fibre interfaces. The degradation for the wet specimens with rubber core is much more pronounced because the decrease in the absorbed impact energy was 53.26% after 10,513 h of immersion in water and all the layers were broken.

Strength and fabric anisotropy of granular materials under true triaxial configurations using DEM

This paper investigates the strength and fabric anisotropy of granular materials under true triaxial configurations using the discrete element method. A clump logic based on the multi-sphere approach was utilized to replicate realistic particle shapes. The evolutions of deviatoric stress and volumetric strains were found to be independent of mean effective stress, and intermediate principal stress parameter (b). From a macroscopic viewpoint, all samples exhibited initial strain hardening followed by softening behaviour, stabilizing at large deviatoric strains due to the dense state of assemblies. The peak state friction angles showed dependency on the b-value, aligning closely with experimental findings. Microscopic parameters such as mechanical coordination number and sliding fraction were found to be approximately independent of b-values. Additionally, non-coaxiality was observed for various b-values except for specific shear modes, aligning with previous findings using spherical particles. These results highlight the capability of clumped particles to simulate the true mechanical behaviour of granular materials and contribute to our understanding of their complex response under different loading conditions.

Frictional strength of siliciclastic sediment mixtures in fault stability assessment

Frictional strength of fault zones is a key parameter for evaluation of fault stability and reactivation. We measure friction using the direct shear test (DST) for (i) sand-clay mixes mimicking fault gouges, and (ii) strength of fault zone interfaces. The sand-clay mixing ratio is linked to the established Shale Gouge Ratio (SGR) used in fault seal analysis, suggesting an approach for linking the measured frictional properties to the established subsurface fault characterization methods and risk assessment. We use powdered caprock material from the Draupne Formation mixed with sand to prepare the fault gouge and discuss application of results for fault zones on the Horda Platform, Norwegian North Sea.Shearing of fault gouge show a systematic decrease in residual friction coefficient with increasing clay content from 0.6 (φ = 31°) in pure sand (SGR 0%) to 0.4 (φ = 22°) for clay rich mixtures (SGR 100%). Interface testing mimicking fault gouge on sand surface is less systematic and shows residual friction coefficients ranging from 0.53 to 0.6 (φ = 28–31°) for SGR 0–50%. Detailed interpretation of the shear testing results shows changes in drainage properties, volumetric changes during shearing and shear responses to normal stresses indicating threshold values for sand versus clay dominated material behaviour. However, the results are non-conclusive on the question if a linear variation of friction with clay content or a threshold between sand dominated versus clay dominated friction provides the best approach for linking friction to fault clay content.

Tracking the movement of quartz sand particles with neural networks

Tracking the movement of particles is of paramount significance for studying the impact of particle breakage on the macroscopic mechanical behaviour of granular materials and remains a persistent challenge to date. This paper presents a novel particle tracking method that integrates Pointon and PointNetLK networks, enabling an accurate tracking of both intact and broken Leighton Buzzard sand (LBS) particles. Initially, morphological information of LBS particles was extracted from X-ray micro-tomography (CT) data collected from a miniature triaxial test. Various image processing techniques were applied to the raw CT images to achieve a realistic three-dimensional (3D) reconstruction. Subsequently, particle point cloud data was processed through sampling, Gaussian noise injection, and grouping for training and testing the Poynton and PointNetLK networks. Next, the correspondences among particles across different scans were established by PointConv, and the transformation matrix between two mutually matched particles was predicted using PointNetLK. Finally, an examination of the changes in the spatial distribution and morphological parameters of both tracked and untracked particles throughout the shearing process was conducted and followed by an analysis of particle kinematics.

Non-linear mechanical behaviour of thermoplastic elastomeric materials and its vulcanizate under tension/tension fatigue deformation by fourier transform rheological studies

Fourier transform (FT) rheology opens up novel frontiers in understanding non-linear mechanical behaviours of polymeric materials under sustained long-term dynamic load. The prediction of the exact point of appearance of a crack in a sample under dynamic mechanical strains of large amplitude and the fatigue analysis has been made possible to such degrees of precision, previously not possible through conventional rheological and mechanical analysis. In this work, the fatigue behaviour of a thermoplastic elastomeric material (TPE) and a thermoplastic vulcanizate (TPV) from thermoplastic polyurethane (TPU) and epichlorohydrin− ethylene oxide−allyl glycidyl ether (GECO) rubber is investigated by Fourier transform studies under oscillatory strain-controlled tensional test to understand the linear and non-linear mechanical behaviour. Fatigue analysis is performed through the fingerprint harmonics of the material’s stress response, i.e., the second and third harmonics, ( I 2/1 and I 3/1) to establish the stress nonlinearity, asymmetry, and the formation of macrocracks. Furthermore, these higher harmonics are fitted with the Neo-Hooke and Mooney-Rivlin model to fundamentally understand the mode of fatigue failure, and a close agreement between theoretical fitting and experimental outcomes is established. Strain-life curves were utilized for the thorough investigation of the fatigue behaviour of the specimens and more than 9-fold enhancement in the strain life of the TPV is observed over TPU at a strain amplitude of ∼1.05 indicating an increasing ductility. Finally, the necessity of FT rheology was further emphasised as it is observed that I 2/1 and I 3/1 harmonics are more sensitive towards fatigue response as compared to the storage modulus and the complex modulus. Henceforth, the FT rheology analysis has enabled the current study to efficiently establish the ductility of each individual specimens and for prediction of the dynamic service lifetime of the developed materials.

Accounting for friction in the mechanical testing of athletics tracks

This investigation deals with the problem of identifying the mechanical behaviour of rubbers from compression tests, performed on specimens having unfavorable geometry. A typical situation is that of flat specimens obtained from high-friction sports surfaces. To this purpose, experimental tests were conducted, aimed at measuring friction under various conditions and evaluating its effect on the compressive behavior of different rubber samples. The experimental results have been interpreted in view of an existing analytical model proposed by Gent and coworkers. The method was shown to be valid within a relatively broad range of conditions (in terms of materials, lubrication and aspect ratio). Its application allowed the creation of virtual “frictionless” curves, by rescaling experimental data for the stiffening factor predicted by Gent model. These curves represent more closely the intrinsic material behaviour, removing the large frictional contribution present in the experimental tests, and can be used as a more reliable input for numerical simulations.

Post-fire structural material behaviour of high- and ultra-high-strength steels – Experimental study and aids for assessment for further use and potential reuse after fires

A thorough understanding and realistic modelling of the essential post-fire material behaviour is crucial for a safe assessment of the residual structural performance of steel components and the reliable further usage of steel structures after fires. This paper presents a comprehensive experimental study of the constitutive material behaviour after heating-cooling cycles of two quenched and tempered high-strength steels (S690QL, S960QL) as well as two thermomechanically rolled steels (S960M, S1100M). The results are compared to both the behaviour of grade S355J2+N mild steels and a comprehensive database of existing post-fire tests. Except for high- and ultra-high-strength steels, which are exposed to very high temperatures above the A1-phase-transformation-temperature during fires, the assumption of reversibility of the mechanical properties for structural designs of steel structures has been confirmed. At high and very high steel temperatures, the microstructure changes and the strength properties after fire are significantly lower than those of virgin high-strength steel specimens. The study thus provides the basis for reliable structural designs and future normative guidance of the post-fire structural steel designs.

Supramolecular J- and H-Aggregation of Naphthalimide-Conjugated Dipeptide in Mixed-Solvent Systems and its Application in Cell Imaging.

In this work, a core-substituted NMI-conjugated dipeptide (4MNLV) was extensively studied in mixed solvent systems to explore the polarity effect on the self-assembly pattern and their photophysical property. 4MNLV adopted J- or H- type aggregation pattern depending upon the polarity index of the solvent system chosen. The self-assembly process was achieved through the anti-solvent effect. UV-vis study suggested that if the stock solution of 4MNLV was diluted with a relatively more polar solvent (compared to the stock solvent), then the system acquired J- type of aggregation pattern by showing a red-shift in their absorption maxima (λmax). Conversely, when the stock was diluted by a relatively less polar solvent, H-type of aggregation was observed, where blue shift of λmax was noticed. The emission spectra and the lifetime of the self-assembled materials were also influenced by the chosen solvent system. The chirotopic behaviour of these self-assembled materials was studied through CD spectroscopy. Morphological study indicated the formation of helical nanofibrillar structures. The bright green fluorescence of these highly biocompatible naphthalimide-peptide conjugate was used for cell imaging application, indicating its futuristic scope.

Dynamic impact mechanical properties of red sandstone based on digital image correlation method

ABSTRACT In the process of underground blasting excavation, the safety and stability of rock mass are significantly influenced by dynamic load. To study the full-field strain and displacement behaviour of rock materials subjected to dynamic compressive force, the impact tests were done on red sandstone specimens using an electromagnetically driven Hopkinson pressure bar. At the same time, the data were collected by a high-speed camera, and non-contact measurements of the red sandstone were carried out with the digital image correlation method. The resulting images were screened and analysed, based on which the stress-strain curve, full-field strain cloud diagram and displacement cloud diagram were obtained. Results show that the dynamic peak stress increases with the increase of impact velocity. During the impact process, substantial internal uneven deformation occurs in the red sandstone close to the transmission bar’s end face. Under the impact load, the cracks initiate from the upper and lower sides near the transmission bar’s end face, followed by formation of cracks in the middle part. These cracks expand throughout the rock to form through-cracks, eventually resulting in complete failure. For rock materials, the evolution of their strain field during the failure process can reflect the initiation and propagation of their internal cracks. Therefore, the use of digital image correlation method proves effective for studying the deformation and failure mechanisms of rock. The conclusions obtained in this study provide a valuable reference for the macro-meso deformation and failure mechanisms of geotechnical materials.

X-IGA Used for Orthotropic Material Crack Growth.

In this paper, we propose a new approach for numerically simulating the growth of cracks in unidirectional composite materials, termed extended isogeometric analysis, evaluating the maximum stress intensity factor and T-stress. To validate our approach, we used a small anisotropic plate with two edge cracks, beginning with formulating the governing equations based on the energy integral method, Stroh's Formula, and the Elastic Law describing the behaviour of anisotropic materials, while considering boundary conditions and initial states. A MATLAB code was developed to solve these equations numerically and to post-process the tensile stress and the stress intensity factor (SIF) in the first mode. The results for the SIF closely match those obtained using the extended finite element method (X-FEM), with a discrepancy of only 0.0021 Pa·m0.5. This finding underscores the credibility of our approach. The extended finite element method has demonstrated robustness in predicting crack propagation in composite materials in recent years, leading to its adoption by several widely used software packages in various industries.

Modelling the effects of void swelling and spatial heterogeneity on the fracture of metallic materials by crystal plasticity

To theoretically analyse how void swelling and spatial heterogeneity impact the fracture behaviour of metallic materials, a coupled porous crystal plasticity framework has been developed. The proposed framework is then applied to analyse the fracture behaviour of metallic materials under different swelling conditions and types of spatial heterogeneity. It is revealed that at low porosity, the spatial heterogeneity of voids will have limited influence on the fracture behaviour of the material. However, at high porosity, owing to void swelling, a concentrated distribution of voids near grain boundaries can lead to a strong deformation localization, which can be seen as an indicator of possible quasi-brittle fracture. It is further discovered that the inner pressure of the voids tends to reduce the yield stress of the material and promote localized plastic deformation. As a result, such pressure will enhance the impact of concentrated distribution of voids, which can lead to further embrittlement of the material. The present study reveals the key role of void swelling and spatial heterogeneity in regulating the fracture behaviour of metallic materials, a step towards a better understanding over the failure mechanisms of metallic materials under extreme conditions.

Computer vision for enhanced quantification of FEA of ballistic impact

In finite element (FE) modelling, the visual assessment of images of contour plots is widely used for qualitative comparison and analysis. However, this approach has limitations in conducting in-depth and comprehensive analysis, which is rigorously made based on objective, massive, and quantitative data. To tackle the problems, we propose a methodology based on computer vision with k-means clustering algorithms, providing an effective quantitative strategy to improve the efficiency and accuracy of the analysis. To demonstrate the application of this methodology, three levels of FE models on impact, viz. macro (homogenous UD laminate model), meso (yarn-level fabric model), and micro (fibre-level yarn model), are developed, where the images with both simple and complex morphologies and backgrounds are all be considered. The contour plots of stress and transverse deformation are clustered into k categories based on the RGB values of the pixels in the images. Notably, the quantification of one image using the proposed methodology is completed in seconds, and the efficiency is enhanced by approximately 40 times compared with that using the manual method, without the accuracy sacrificed. This methodology enables the researchers to efficiently and precisely understand the ballistic behaviours of the textile materials from the perspectives of fibres, yarn, fabrics, etc., and has great potential to apply in other research fields, where computational simulation is widely used.

Seismic Assessment of an Old Historical Unreinforced Masonry Power House Building

Abstract Design of unreinforced masonry buildings is a major concern in seismically active zones in Nepal. The ChandraJyoti Jalvidhyut Power House, a historical landmark in Nepal, represents the nation’s pioneering effort in hydropower generation and signifies its rich cultural heritage. This study comprehensively assesses the structural assessment and seismic performance deficiencies of the ChandraJyoti Jalvidhyut Power House. Material tests on bricks with the individual properties and behaviour of construction materials were done. Utilizing a finite element model in DIANA software V 10.5, the structure’s response under pushover loading conditions is simulated, providing insights into its structural behaviour. Additionally, fragility curves are developed to correlate earthquake intensity with the probability of damage, offering essential information for preservation strategies and seismic risk mitigation. The compressive strength test of bricks, with a value of 9.6 MPa, falls short of the 10 MPa minimum requirement specified by NBC 203:2015, indicating insufficient load-bearing capacity. The building’s pushover analysis shows that the longitudinal direction experiences higher lateral forces (base shear of 945.87 kN) and greater drift ratios exceeding 0.4%, which increases the likelihood of cracks compared to the transverse direction (base shear of 826.77 kN). The analysis also highlights significant nonlinear behavior in both directions, with maximum roof displacements of 35.18 mm (transverse) and 34.4 mm (longitudinal), reflecting the impact of structural irregularities.

A Discrete Element Method Calibration Approach for Large Ore-flow Applications

The major constraint for large discrete element method (DEM) simulations is the high demand for computation power. To keep the model in a feasible range of simulation time and to conduct the numerical calculations efficiently, the implemented material can be simplified and upscaled. After this, a calibration of the material is inevitable and sets the link between the simulation model and reality. Therefore, an alternative calibration approach is established in this contribution to mimic the flow behaviour of the original material with multi-spherical DEM particles. The calibration approach utilises flow zone dimensions, which were observed in block-caving operations. Additionally, a passive draw-point is used to calibrate the angle of repose. Applying the calibration approach for Kiruna iron ore showed that, with the help of the calibration approach, the material in the simulation could successfully reproduce the targeted flow zone and angle of repose.

Study on the Yield Surface of a Fibre-Reinforced Stabilised Soil – Effect of the Number of Loading Cycles

Cyclic loading may induce changes in the geomechanical behaviour of materials that should be characterised. This work studies the impact of the number of loading cycles on the mechanical behaviour of a fibre-reinforced stabilised soil focusing on its behaviour before failure (yield surface). To this end, an experimental testing program based on triaxial tests was performed on samples not subjected to a cycling loading stage, as well as on samples previously subjected to a cycling loading stage varying the number of loading cycles from 1,000 to 100,000. The results were studied in terms of the accumulated permanent axial strain and the yield surface of the composite material. It was observed that increasing the number of loading cycles led to a rise in the accumulated permanent axial strain and in the undrained resilient modulus. The results also showed an expansion of the yield surface during the first 1,000 loading cycles (the yield occurs later due to the partial mobilization of the tensile strength of the fibres during the cyclic stage) but its shape is maintained. The results also showed a progressive reduction in the yield loci with the increase in the number of loading cycles, reflecting the greater degradation of the solid matrix induced by the accumulated permanent axial strains.

High temperature deformation behaviour and recrystallization mechanism of TiC/Ti6Al4V gradient material fabricated by laser directed energy deposition

Metal-ceramic gradient materials are of considerable interest due to their unique capability to integrate the ductility of metals with the mechanical strength and high-temperature resistance of ceramics, thereby reducing the risk of cracking induced by thermal stress. In this article, TiC/Ti6Al4V gradient materials were prepared through the laser directed energy deposition (LDED) process. Subsequently, high-temperature compression tests were performed at various temperatures to investigate the deformation behaviour and recrystallization mechanism of the gradient materials. The Zener-Holomon constitutive equation was then used to calculate the gradient material at 750–900 °C and 0.001–0.1 s−1 to obtain a deformation activation energy of 362.36 KJ·mol−1. The EBSD was used to analyze the microstructural evolution of the gradient material. The results indicate a significant reduction in the average grain size of the specimens after high-temperature compression, and a more random crystal orientation due to the occurrence of dynamic recrystallisation (DRX). The primary softening mechanism of TiC/Ti6Al4V gradient materials is discontinuous dynamic recrystallisation, with TiC primarily affecting the nucleation of discontinuous dynamic recrystallisation (DDRX) and the subsequent grain growth. The nucleation becomes easier with increasing TiC content. However, the TiC hinders the migration of grain boundaries, resulting in a smaller DDRX grain size. After high-temperature compression, the main slip system of the crystal is converted from Pri to Bas. Pyr is difficult to activate due to its high critical resolved shear stress.

Flow Analysis of Centrifugal Compressor using ANSYS

Abstract: The aerodynamic and structural designs of the impeller are vital to the functioning of the entire compressor stage in centrifugal compressors. In contrast to previous years, there has been a greater demand for industrial centrifugal compressors to operate with greater efficiency and operational range. Although the efficiency of low-pressure-ratio centrifugal compressors has nearly achieved its maximum, intermediate and high-pressure ratio machines still lag behind the state-of-the-art in terms of efficiency. Maintaining maximum efficiency while extending the working range is a challenging task in centrifugal compressor design. Furthermore, a single design must be usable globally due to product globalisation, even in the face of technological disparities and variances in electricity frequency of up to 10 percent in some nations. Often operating at different rotational speeds, centrifugal compressors add another layer of complexity to the impeller structural design. A full-3D impeller design that is optimised for both structural integrity and aerodynamic performance is shown in this study .For air conditioning systems of the future, centrifugal compressor performance must be improved. This research focusses on systematic numerical simulation, a computationally efficient method of optimising centrifugal compressor performance. ANSYS simulations are used in conjunction with a thorough literature review on the aerodynamic behaviour of materials used in centrifugal compressors. For the purpose of performing flow simulations under subsonic circumstances with suitable boundary conditions and grid independence checks, geometric modelling was completed using CATIAV5R22 and imported into ANSYS. Tests were conducted using various aluminium materials and mass flow rates, and enhanced efficiency findings were confirmed by comparison with experimental data. Materials Processing and Characterisations is in charge of peer review.

A software-agnostic benchmark for DEM simulation of cohesive and non-cohesive materials

Although DEM (Discrete Element Method) was introduced >40 years ago, there are still various challenges in applying it with a certain level of confidence. To develop a realistic material model, calibration, validation, particle shape representation and scaling are well-known challenges that have been studied by various researchers over the past two decades. In addition, once a realistic DEM material model is developed and published in line with open science principles, it is unknown to what extent the calibrated values can be used across software packages. Although a benchmark between 8 open source codes was recently published, it did not cover cohesive materials, rolling friction models, realistic calibrated material behaviour, and commercial software that is widely used in industries.The aim of this study is to investigate the portability of input parameters between different codes and to identify if software independent calibration is possible. We propose a framework to assist DEM users in industry and academia to reach a software independent DEM simulation when required.To compare the results between software packages, the framework considers the following aspects: 1) identical implementation of contact models, 2) single contact modelling, 3) bulk level simulation, and 4) identical post-processing. This framework is demonstrated for two contact models (Hertz-Mindlin and Edinburgh Elasto Plastic Adhesive) with rolling friction model and two commercial software packages (EDEM and PFC) but is applicable for any other (open-source) software. Moreover, two realistic calibrated material models are included to cover a range of material characteristics, from free-flowing incompressible materials to cohesive compressible materials. This study shows that individual particle level simulations give identical results, bulk level simulations show differences for both contact models. In general, users should use parameter values from other software packages with caution, especially where critical or sensitive applications are modelled. This paper also highlights the use of novel computation techniques, such as the GPU engine, to achieve practical computation times when modelling industrial applications.

Numerical Study on the Effect of Geometrical Changes on the Deformation Behaviour of Recycled Aluminium Alloy AA6061 undergoing High Velocity Impact using Hyperworks-Radioss

Numerical analysis is important to predict material behaviour without require an experimental implementation. This manuscript focuses on a numerical prediction establishment of a direct recycled aluminium alloys AA6061 undergoing Taylor Cylinder Impact test using Johnson-Cook model in HyperWorks Radioss. The numerical setup was first validated against the experimental data at the velocity range of 179 to 212 m/s. Good agreement between simulation and experimental data was obtained within the range that exhibits a mushrooming shape fracture mode. A parametric study was then conducted to study the deformation behaviour of the selected recycled aluminum alloys within the validated range at various geometrical settings. The analysis was made by focusing on the post-impact configuration of the projectile at different impact velocities in terms of residual length, deformed diameter, and the final length-to-diameter ratio. It was found that a broader projectile experienced a less significant reduction in its final length (Lf/Lo goes from 0.87 to 0.9 for projectile diameter 9mm to 34mm) and a smaller increase in the deformed diameter compared to a thinner projectile (Df/Do goes from 1.18 to 1.12 for projectile diameter 9mm to 34mm). It was found that a thinner projectile experienced more diameter expansion than length reduction post impact. In addition, a longer projectile experienced more residual length reduction (Lf/Lo goes from 0.92 to 0.87 for projectile length 2mm to 14mm) and more radial deformation compared to the one with a smaller initial length (Df/Do goes from 1.06 to 1.18 for projectile length 2mm to 14mm). All projectiles showed more significant changes on the deformed diameter compared to the changes in residual length post-impact. The results helped the understanding of a critical aspect of the deformation behaviour of recycled aluminum alloy AA6061 more effectively compared to experimental work implementation.