Search results

Accurately evaluating fluid flow behaviors and determining permeability for deforming porous media is time-consuming and remains challenging. This paper aims to propose a novel machine-learning method for the rapid estimation of permeability of porous media at different deformation stages constrained by hydro-mechanical coupling analysis.

A convolutional neural network (CNN) is proposed in this paper, which is guided by the results of finite element coupling analysis of equilibrium equation for mechanical deformation and Boltzmann equation for fluid dynamics during the hydro-mechanical coupling process [denoted as Finite element lattice Boltzmann model (FELBM) in this paper]. The FELBM ensures the Lattice Boltzmann analysis of coupled fluid flow with an unstructured mesh, which varies with the corresponding nodal displacement resulting from mechanical deformation. It provides reliable label data for permeability estimation at different stages using CNN.

The proposed CNN can rapidly and accurately estimate the permeability of deformable porous media, significantly reducing processing time. The application studies demonstrate high accuracy in predicting the permeability of deformable porous media for both the test and validation sets. The corresponding correlation coefficients (R²) is 0.93 for the validation set, and the R² for the test set A and test set B are 0.93 and 0.94, respectively.

This study proposes an innovative approach with the CNN to rapidly estimate permeability in porous media under dynamic deformations, guided by FELBM coupling analysis. The fast and accurate performance of CNN underscores its promising potential for future applications.

Realistic composite vehicles with 2, 3, 5 and 9 axles, consisting of a truck with one or two trailers, are addressed in this paper by computational models for vehicle–bridge interaction analysis.

The vehicle–bridge interaction (VBI) models are formed by sets of 2-D rigid blocks interconnected by mass, damping and stiffness elements to simulate their suspension system. The passage of the vehicles is performed at different speeds. Several rolling surface profiles are admitted, considering the maintenance grade of the pavement. The spectral density functions are generated from an experimental database to form the longitudinal surface irregularity profiles. A computational code written in Phyton based on the finite element method was developed considering the Euler–Bernoulli beam model.

Several models of composite heavy vehicles are presented as manufactured and currently travel on major roads. Dynamic amplification factors are presented for each type of composite vehicle.

The VBI models for compound heavy vehicles are 2-D.

This work contributes to improving the safety and lifetime of the bridges, as well as the stability and comfort of the vehicles when passing over a bridge.

The structural response of the bridge is affected by the type and size of the compound vehicles, their speed and the conservative grade of the pavement. Moreover, one axle produces vibrations that can be superposed by the vibrations of the other axles. This effect can generate not usual dynamic responses.

This study aims to reduce the redundant weight of the anti-roll torsion bar brought by the traditional empirical design and improving its strength and stiffness.

Based on the finite element approach coupled with the improved beluga whale optimization (IBWO) algorithm, a collaborative optimization method is suggested to optimize the design of the anti-roll torsion bar structure and weight. The dimensions and material properties of the torsion bar were defined as random variables, and the torsion bar's mass and strength were investigated using finite elements. Then, chaotic mapping and differential evolution (DE) operators are introduced to improve the beluga whale optimization (BWO) algorithm and run case studies.

The findings demonstrate that the IBWO has superior solution set distribution uniformity, convergence speed, solution correctness and stability than the BWO. The IBWO algorithm is used to optimize the anti-roll torsion bar design. The error between the optimization and finite element simulation results was less than 1%. The weight of the optimized anti-roll torsion bar was lessened by 4%, the maximum stress was reduced by 35% and the stiffness was increased by 1.9%.

The study provides a methodological reference for the simulation optimization process of the lateral anti-roll torsion bar.

A novel method for modelling permanent magnets is investigated based on numerical approximations with rational functions. This study aims to introduce the AAA algorithm and other recently developed, cutting-edge mathematical tools, which provide outstandingly fast and accurate numerical computation of potentials and vector fields.

First, the AAA algorithm is briefly introduced along with its main variants and other advanced mathematical tools involved in the modelling. Then, the analysis of a circular Halbach array with a one-pole pair is carried out by means of the AAA-least squares method, focusing on vector potential and flux density in the bore and validating results by means of classic finite element software. Finally, the investigation is completed by a finite difference analysis.

AAA methods for field analysis prove to be strikingly fast and accurate. Results are in excellent agreement with those provided by the finite element model, and the very good agreement with those from finite differences suggests future improvements. They are also easy programming; the MATLAB code is less than 200 lines. This indicates they can provide an effective tool for rapid analysis.

AAA methods in magnetostatics are novel, but their extension to analogous physical problems seems straightforward. Being a meshless method, it is unlikely that local non-linearities can be considered. An aspect of particular interest, left for future research, is the capability of handling inhomogeneous domains, i.e. solving general interface problems.

The authors use cutting-edge mathematical tools for the modelling of complex physical objects in magnetostatics.

This study investigates the impact of the fire decay phase on structural damage using the sectional analysis method. The primary objective of this work is to forecast the non-dimensional capacity parameters for the axial and flexural load-carrying capacity of reinforced concrete (RC) sections for heating and the subsequent post-heating phase (decay phase) of the fire.

The sectional analysis method is used to determine the moment and axial capacities. The findings of sectional analysis and heat transfer for the heating stage are initially validated, and the analysis subsequently proceeds to determine the load capacity during the fire’s heating and decay phases by appropriately incorporating non-dimensional sectional and material parameters. The numerical analysis includes four fire curves with different cooling rates and steel percentages.

The study’s findings indicate that the rate at which the cooling process occurs after undergoing heating substantially impacts the axial and flexural capacity. The maximum degradation in axial and flexural capacity occurred in the range of 15–20% for cooling rates of 3 °C/min and 5 °C/min as compared to the capacity obtained at 120 min of heating for all steel percentages. As the fire cooling rate reduced to 1 °C/min, the highest deterioration in axial and flexural capacity reached 48–50% and 42–46%, respectively, in the post-heating stage.

The established non-dimensional parameters for axial and flexural capacity are limited to the analysed section in the study owing to the thermal profile, however, this can be modified depending on the section geometry and fire scenario.

The study primarily focusses on the degradation of axial and flexural capacity at various time intervals during the entire fire exposure, including heating and cooling. The findings obtained showed that following the completion of the fire’s heating phase, the structural capacity continued to decrease over the subsequent post-heating period. It is recommended that structural members' fire resistance designs encompass both the heating and cooling phases of a fire. Since the capacity degradation varies with fire duration, the conventional method is inadequate to design the load capacity for appropriate fire safety. Therefore, it is essential to adopt a performance-based approach while designing structural elements' capacity for the desired fire resistance rating. The proposed technique of using non-dimensional parameters will effectively support predicting the load capacity for required fire resistance.

The fire-resistant requirements for reinforced concrete structures are generally established based on standard fire exposure conditions, which account for the fire growth phase. However, it is important to note that concrete structures can experience internal damage over time during the decay phase of fires, which can be quantitatively determined using the proposed non-dimensional parameter approach.

The “chiralisation” of Euclidean polygonal tessellations is a novel, recent method which has been used to design new auxetic metamaterials with complex topologies and improved geometric versatility over traditional chiral honeycombs. This paper aims to design and manufacture chiral honeycombs representative of four distinct classes of 2D Euclidean tessellations with hexagonal rotational symmetry using fused-deposition additive manufacturing and experimentally analysed the mechanical properties and failure modes of these metamaterials.

Finite Element simulations were also used to study the high-strain compressive performance of these systems under both periodic boundary conditions and realistic, finite conditions. Experimental uniaxial compressive loading tests were applied to additively manufactured prototypes and digital image correlation was used to measure the Poisson’s ratio and analyse the deformation behaviour of these systems.

The results obtained demonstrate that these systems have the ability to exhibit a wide range of Poisson’s ratios (positive, quasi-zero and negative values) and stiffnesses as well as unusual failure modes characterised by a sequential layer-by-layer collapse of specific, non-adjacent ligaments. These findings provide useful insights on the mechanical properties and deformation behaviours of this new class of metamaterials and indicate that these chiral honeycombs could potentially possess anomalous characteristics which are not commonly found in traditional chiral metamaterials based on regular monohedral tilings.

To the best of the authors’ knowledge, the authors have analysed for the first time the high strain behaviour and failure modes of chiral metamaterials based on Euclidean multi-polygonal tessellations.

This study establishes a method for predicting ground vibrations caused by railway tunnels in unsaturated soils with spatial variability.

First, an improved 2.5D finite-element-method-perfect-matching-layer (FEM-PML) model is proposed. The Galerkin method is used to derive the finite element expression in the ub-pl-pg format for unsaturated soil. Unlike the ub-v-w format, which has nine degrees of freedom per node, the ub-pl-pg format has only five degrees of freedom per node; this significantly enhances the calculation efficiency. The stretching function of the PML is adopted to handle the unlimited boundary domain. Additionally, the 2.5D FEM-PML model couples the tunnel, vehicle and track structures. Next, the spatial variability of the soil parameters is simulated by random fields using the Monte Carlo method. By incorporating random fields of soil parameters into the 2.5D FEM-PML model, the effect of soil spatial variability on ground vibrations is demonstrated using a case study.

The spatial variability of the soil parameters primarily affected the vibration acceleration amplitude but had a minor effect on its spatial distribution and attenuation over time. In addition, ground vibration acceleration was more affected by the spatial variability of the soil bulk modulus of compressibility than by that of saturation.

Using the 2.5D FEM-PML model in the ub-pl-pg format of unsaturated soil enhances the computational efficiency. On this basis, with the random fields established by Monte Carlo simulation, the model can calculate the reliability of soil dynamics, which was rarely considered by previous models.

Thermal protection of a flange is critical for preventing tower icing and collapse of wind turbines (WTs) in extremely cold weather. This study aims to develop a novel thermal protection system for the WTs flanges using an electrical heat-tracing element.

A three-dimensional model and the Poly-Hexacore mesh structure are used, and the fluid-solid coupling method was validated and then deployed to analyze the heat transfer and convection process. Intra-volumetric heat sources are applied to represent the heat generated by the heating element, and the dynamic boundary conditions are considered. The steady temperature and temperature uniformity of the flange are the assessment criteria for the thermal protection performance of the heating element.

Enlarging the heating area and increasing the heating power improved the flange's temperature and temperature uniformity. A heating power of 4.9 kW was suitable for engineering applications with the lowest temperature nonuniformity. Compared with continuous heating, the increased temperature nonuniformity was buffered, and the electrical power consumption was reduced by half using pulse heating. Pulse heating time intervals of 1, 3 and 4 h were determined for the spring, autumn and winter, respectively.

The originality of this study is to propose a novel electrical heat-tracing thermal protection system for the WTs flanges. The effect of different arrangements, heating powers and heating strategies was studied, by which the theoretical basis is provided for a stable and long-term utilization of the WT flange.

The use of waste polyethylene terephthalate (PET) plastics derived from shredded bottles in concrete is not formalized yet, especially in reinforced members such as beams and columns. The disposal of plastic wastes in concrete is a viable alternative to manage those wastes while minimizing the environmental impacts associated to recycling, carbon dioxide emissions and energy consumption.

This paper evaluates the suitability of 2D deterministic and stochastic finite element (FE) modeling to predict the shear strength behavior of reinforced concrete (RC) beams without stirrups. Different concrete mixtures prepared with 1.5%–4.5% PET additions, by volume, are investigated.

Test results showed that the deterministic and stochastic FE approaches are accurate to assess the maximum load of RC beams at failure and corresponding midspan deflection. However, the crack patterns observed experimentally during the different stages of loading can only be reproduced using the stochastic FE approach. This later method accounts for the concrete heterogeneity due to PET additions, allowing a statistical simulation of the effect of mechanical properties (i.e. compressive strength, tensile strength and Young’s modulus) on the output FE parameters.

Data presented in this paper can be of interest to civil and structural engineers, aiming to predict the failure mechanisms of RC beams containing plastic wastes, while minimizing the experimental time and resources needed to estimate the variability effect of concrete properties on the performance of such structures.

The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.

An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.

The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.

The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.