Search results

This study implements the fourth-order phase field method (PFM) for modeling fracture in brittle materials. The weak form of the fourth-order PFM requires C¹ basis functions for the crack evolution scalar field in a finite element framework. To address this, non-Sibsonian type shape functions that are nonpolynomial types based on distance measures, are used in the context of natural neighbor shape functions. The capability and efficiency of this method are studied for modeling cracks.

The weak form of the fourth-order PFM is derived from two governing equations for finite element modeling. C⁰ non-Sibsonian shape functions are derived using distance measures on a generalized quad element. Then these shape functions are degree elevated with Bernstein-Bezier (BB) patch to get higher-order continuity (C¹) in the shape function. The quad element is divided into several background triangular elements to apply the Gauss-quadrature rule for numerical integration. Both fourth-order and second-order PFMs are implemented in a finite element framework. The efficiency of the interpolation function is studied in terms of convergence and accuracy for capturing crack topology in the fourth-order PFM.

It is observed that fourth-order PFM has higher accuracy and convergence than second-order PFM using non-Sibsonian type interpolants. The former predicts higher failure loads and failure displacements compared to the second-order model due to the addition of higher-order terms in the energy equation. The fracture pattern is realistic when only the tensile part of the strain energy is taken for fracture evolution. The fracture pattern is also observed in the compressive region when both tensile and compressive energy for crack evolution are taken into account, which is unrealistic. Length scale has a certain specific effect on the failure load of the specimen.

Fourth-order PFM is implemented using C¹ non-Sibsonian type of shape functions. The derivation and implementation are carried out for both the second-order and fourth-order PFM. The length scale effect on both models is shown. The better accuracy and convergence rate of the fourth-order PFM over second-order PFM are studied using the current approach. The critical difference between the isotropic phase field and the hybrid phase field approach is also presented to showcase the importance of strain energy decomposition in PFM.

Fused Filament Fabrication (FFF) is an extrusion-based manufacturing process using fused thermoplastics. Despite its low cost, the FFF is not extensively used in high-value industrial sectors mainly due to parts' anisotropy (related to the deposition strategy) and residual stresses (caused by successive heating cycles). Thus, this study aims to investigate the process improvement and the optimization of the printed parts.

In this work, a meshless technique – the Radial Point Interpolation Method (RPIM) – is used to numerically simulate the viscoplastic extrusion process – the initial phase of the FFF. Unlike the FEM, in meshless methods, there is no pre-established relationship between the nodes so the nodal mesh will not face mesh distortions and the discretization can easily be modified by adding or removing nodes from the initial nodal mesh. The accuracy of the obtained results highlights the importance of using meshless techniques in this field.

Meshless methods show particular relevance in this topic since the nodes can be distributed to match the layer-by-layer growing condition of the printing process.

Using the flow formulation combined with the heat transfer formulation presented here for the first time within an in-house RPIM code, an algorithm is proposed, implemented and validated for benchmark examples.

Bayesian cubature (BC) has emerged to be one of most competitive approach for estimating the multi-dimensional integral especially when the integrand is expensive to evaluate, and alternative acquisition functions, such as the Posterior Variance Contribution (PVC) function, have been developed for adaptive experiment design of the integration points. However, those sequential design strategies also prevent BC from being implemented in a parallel scheme. Therefore, this paper aims at developing a parallelized adaptive BC method to further improve the computational efficiency.

By theoretically examining the multimodal behavior of the PVC function, it is concluded that the multiple local maxima all have important contribution to the integration accuracy as can be selected as design points, providing a practical way for parallelization of the adaptive BC. Inspired by the above finding, four multimodal optimization algorithms, including one newly developed in this work, are then introduced for finding multiple local maxima of the PVC function in one run, and further for parallel implementation of the adaptive BC.

The superiority of the parallel schemes and the performance of the four multimodal optimization algorithms are then demonstrated and compared with the k-means clustering method by using two numerical benchmarks and two engineering examples.

Multimodal behavior of acquisition function for BC is comprehensively investigated. All the local maxima of the acquisition function contribute to adaptive BC accuracy. Parallelization of adaptive BC is realized with four multimodal optimization methods.

The purpose of this study is to model the theory of the low-velocity impact (LVI) process on sandwich beams consisting of flexible cores and face sheets reinforced with functionally graded carbon nanotubes (CNTs).

A series of parameters derived from molecular dynamics are used to consider the size scale in the mixture rule for the combination of CNTs and resin. A procedure involving the use of the first-order shear deformation theory of the beam is used to provide the displacement field of the sandwich beam. The energy method and subsequently the generalized Lagrange method are used to derive the motion equations. Due to the use of Hertz’s nonlinear theory to calculate the contact force, the equations of motion are nonlinear. Validation of the problem is carried out by comparing natural frequencies with other papers.

The influence of a series of parameters such as CNTs distributions pattern in the face sheets, the influence of the CNTs volume fraction and the influence of the core thickness to the face sheets thickness ratio in the issue of LVI on sandwich beams with clamped-clamped boundary conditions is investigated. The result shows that the type of CNTs pattern in the face sheet and the CNTs volume fraction have a very important effect on the answer to the problem, which is caused by the change in the value of the Young’s modulus of the beam at the contact surface. Changes in the core thickness to the face sheets thickness ratio has little effect on the impact response.

Considering the important application of sandwich structures in vehicles, aviation and ships, in this research, sandwich beams consisting of flexible core and CNTs-reinforced face sheets are investigated under LVI.

This study aims to assess the robustness and accuracy of the face-centred finite volume (FCFV) method for the simulation of compressible laminar flows in different regimes, using numerical benchmarks.

The work presents a detailed comparison with reference solutions published in the literature –when available– and numerical results computed using a commercial cell-centred finite volume software.

The FCFV scheme provides first-order accurate approximations of the viscous stress tensor and the heat flux, insensitively to cell distortion or stretching. The strategy demonstrates its efficiency in inviscid and viscous flows, for a wide range of Mach numbers, also in the incompressible limit. In purely inviscid flows, non-oscillatory approximations are obtained in the presence of shock waves. In the incompressible limit, accurate solutions are computed without pressure correction algorithms. The method shows its superior performance for viscous high Mach number flows, achieving physically admissible solutions without carbuncle effect and predictions of quantities of interest with errors below 5%.

The FCFV method accurately evaluates, for a wide range of compressible laminar flows, quantities of engineering interest, such as drag, lift and heat transfer coefficients, on unstructured meshes featuring distorted and highly stretched cells, with an aspect ratio up to ten thousand. The method is suitable to simulate industrial flows on complex geometries, relaxing the requirements on mesh quality introduced by existing finite volume solvers and alleviating the need for time-consuming manual procedures for mesh generation to be performed by specialised technicians.

Ensuring the safety of structures is important. However, when a structure possesses both an implicit performance function and an extremely small failure probability, traditional methods struggle to conduct a reliability analysis. Therefore, this paper proposes a reliability analysis method aimed at enhancing the efficiency of rare event analysis, using the widely recognized Relevant Vector Machine (RVM).

Drawing from the principles of importance sampling (IS), this paper employs Harris Hawks Optimization (HHO) to ascertain the optimal design point. This approach not only guarantees precision but also facilitates the RVM in approximating the limit state surface. When the U learning function, designed for Kriging, is applied to RVM, it results in sample clustering in the design of experiment (DoE). Therefore, this paper proposes a FU learning function, which is more suitable for RVM.

Three numerical examples and two engineering problem demonstrate the effectiveness of the proposed method.

By employing the HHO algorithm, this paper innovatively applies RVM in IS reliability analysis, proposing a novel method termed RVM-HIS. The RVM-HIS demonstrates exceptional computational efficiency, making it eminently suitable for rare events reliability analysis with implicit performance function. Moreover, the computational efficiency of RVM-HIS has been significantly enhanced through the improvement of the U learning function.

The strength and stiffness of steel deteriorate rapidly at elevated temperatures. Thus, the characteristics of steel structures exposed to fire have been concerned in recent years. Most studies on the fire response of steel structures were conducted at uniformly distributed temperatures. This study aims to evaluate the buckling capacity of steel H-beams subjected to different loading conditions under non-uniform heating.

A numerical investigation was conducted employing finite element analysis software, ABAQUS. A comparison between the numerical analysis results and the experimental data from previous studies was conducted to verify the beam model. Simply supported beams were loaded with several loading conditions including one end moment, end equal moments, uniformly distributed load and concentrated load at midspan. The effects of initial imperfections were considered. The buckling capacities of steel beams under fire using the existing fire design code and the previous study were also generated and compared.

The results showed that the length-to-height ratio and loading conditions have a great effect on the buckling resistance of steel beams under fire. The capacity of steel beams under non-uniform temperature distribution using the existing fire design code and the previous study can give unconservative values or too conservative values depending on loading conditions. The maximum differences of unconservative and conservative values are −44.5 and 129.2% for beams subjected to end equal moments and one end moment, respectively.

This study provides the buckling characteristics of steel beams under non-uniform temperature considering the influences of initial imperfections, length-to-height ratios, and loading conditions. This study will be beneficial for structural engineers in properly evaluating structures under non-uniform heating conditions.

The purpose of this paper is to design a simple and accurate a-posteriori Lagrangian-based error estimator is developed for the class of backward differentiation formula (BDF) algorithms with variable time step size, and the adaptive time-stepping in BDF algorithms is demonstrated for efficient time-dependent simulations in fluid flow and heat transfer.

The Lagrange interpolation polynomial is used to predict the time derivative, and then the accurate primary result is obtained by the Gauss integral, which is applied to evaluate the local error. Not only the generalized formula of the proposed error estimator is presented but also the specific expression for the widely applied BDF1/2/3 is illustrated. Two essential executable MATLAB functions to implement the proposed error estimator are appended for practical applications. Then, the adaptive time-stepping is demonstrated based on the newly proposed error estimator for BDF algorithms.

The validation tests show that the newly proposed error estimator is accurate such that the effectivity index is always close to unity for both linear and nonlinear problems, and it avoids under/overestimation of the exact local error. The applications for fluid dynamics and coupled fluid flow and heat transfer problems depict the advantage of adaptive time-stepping based on the proposed error estimator for time-dependent simulations.

In contrast to existing error estimators for BDF algorithms, the present work is more accurate for the local error estimation, and it can be readily extended to practical applications in engineering with a few changes to existing codes, contributing to efficient time-dependent simulations in fluid flow and heat transfer.

The purpose of this study is to solve two unique but difficult partial differential equations: the foam drainage equation and the nonlinear time-fractional fisher’s equation. Through our methods, we aim to provide accurate solutions and gain a deeper understanding of the intricate behaviors exhibited by these systems.

In this study, we use a dual technique that combines the Aboodh residual power series method and the Aboodh transform iteration method, both of which are combined with the Caputo operator.

We develop exact and efficient solutions by merging these unique methodologies. Our results, presented through illustrative figures and data, demonstrate the efficacy and versatility of the Aboodh methods in tackling such complex mathematical models.

Owing to their fractional derivatives and nonlinear behavior, these equations are crucial in modeling complex processes and confront analytical complications in various scientific and engineering contexts.

The purpose of the paper is the simulation of nonuniform transmission lines.

The method involves a Magnus expansion and a numerical Laplace transform. The method involves a judicious arrangement of the governing equations so as to enable efficient simulation.

The results confirm an effective and efficient numerical solver for inclusion of nonuniform transmission lines in circuit simulation.

The work combines a Magnus expansion and numerical Laplace transform algorithm in a novel manner and applies the resultant algorithm for the effective and efficient simulation of nonuniform transmission lines.