Engineering Computations最新文献

Purpose

The constitutive models determine the mechanical response to the defined loading based on model parameters. In this paper, the inverse problem is researched, i.e. the identification of the model parameters based on the loading and responses of the material. The conventional methods for determining the parameters of constitutive models often demand significant computational time or extensive model knowledge for manual calibration. The aim of this paper is to introduce an alternative method, based on artificial neural networks, for determining the parameters of a viscoplastic model.

Design/methodology/approach

An artificial neural network was proposed to determine nine material parameters of a viscoplastic model using data from three half-life hysteresis loops. The proposed network was used to determine the material parameters from uniaxial low-cycle fatigue experimental data of an aluminium alloy obtained at elevated temperatures and three different mechanical strain rates.

Findings

A reasonable correlation between experimental and numerical data was achieved using the determined material parameters.

Originality/value

This paper fulfils a need to research alternative methods of identifying material parameters.

Purpose

The purpose of this research is to establish a robust numerical framework for the calibration of macroscopic constitutive parameters, based on the analysis of polycrystalline RVEs with computational homogenisation.

Design/methodology/approach

This framework is composed of four building-blocks: (1) the multi-scale model, consisting of polycrystalline RVEs, where the grains are modelled with anisotropic crystal plasticity, and computational homogenisation to link the scales, (2) a set of loading cases to generate the reference responses, (3) the von Mises elasto-plastic model to be calibrated, and (4) the optimisation algorithms to solve the inverse identification problem. Several optimisation algorithms are assessed through a reference identification problem. Thereafter, different calibration strategies are tested. The accuracy of the calibrated models is evaluated by comparing their results against an FE2 model and experimental data.

Findings

In the initial tests, the LIPO optimiser performs the best. Good results accuracy is obtained with the calibrated constitutive models. The computing time needed by the FE2 simulations is 5 orders of magnitude larger, compared to the standard macroscopic simulations, demonstrating how this framework is suitable to obtain efficient micro-mechanics-informed constitutive models.

Originality/value

This contribution proposes a numerical framework, based on FE2 and macro-scale single element simulations, where the calibration of constitutive laws is informed by multi-scale analysis. The most efficient combination of optimisation algorithm and definition of the objective function is studied, and the robustness of the proposed approach is demonstrated by validation with both numerical and experimental data.

Purpose

This paper aims to construct positivity-preserving finite volume schemes for the three-dimensional convection–diffusion equation that are applicable to arbitrary polyhedral grids.

Design/methodology/approach

The cell vertices are used to define the auxiliary unknowns, and the primary unknowns are defined at cell centers. The diffusion flux is discretized by the classical nonlinear two-point flux approximation. To ensure the fully discrete scheme has positivity-preserving property, an improved discretization method for the convection flux was presented. Besides, a new positivity-preserving vertex interpolation method is derived from the linear reconstruction in the discretization of convection flux. Moreover, the Picard iteration method may have slow convergence in solving the nonlinear system. Thus, the Anderson acceleration of Picard iteration method is used to solve the nonlinear system. A condition number monitor of matrix is employed in the Anderson acceleration method to achieve better robustness.

Findings

The new scheme is applicable to arbitrary polyhedral grids and has a second-order accuracy. The results of numerical experiments also confirm the positivity-preserving of the discretization scheme.

Originality/value

1. This article presents a new positivity-preserving finite volume scheme for the 3D convection–diffusion equation. 2. The new discretization scheme of convection flux is constructed. 3. A new second-order interpolation algorithm is given to eliminate the auxiliary unknowns in flux expressions. 4. An improved Anderson acceleration method is applied to accelerate the convergence of Picard iterations. 5. This scheme can solve the convection–diffusion equation on the distorted meshes with second-order accuracy.

Purpose

This study focuses on the CFD numerical simulation and analysis of the vortex stacking problem at the top of the impeller of a high-speed fuel pump, mainly using LCS and entropy production theory to visualize the vortex at the top of the impeller as well as quantitatively analyzing the energy loss caused by the vortex at the top of the impeller. By combining the two methods, the two are well verified with each other that the stacking problem of the vortex at the top of the impeller and the location of the energy loss caused by the vortex are consistent with the vortex location. Such a method can reveal the problem of vortex buildup at the top of the lobe well, and provide a novel guidance idea for improving the performance of high-speed fuel pumps.

Design/methodology/approach

Based on CFD numerical simulation and analysis, this study mainly uses LCS and entropy production theory to visualize the top vortex of the impeller. Through the combination of the two methods, the accumulation problem of the top vortex of the impeller and the location of the energy loss caused by the vortex can be well revealed.

Findings

(1) The CFD numerical simulation analysis of the high-speed fuel pump is carried out, and the test is conducted to verify the numerical simulation results. The inlet and outlet pressure difference? P is used as the validation index, and the error analysis shows that the error between numerical simulation and test results is within 10%, which meets our requirements. Therefore, we carry out the next analysis with the help of CFD numerical simulation. By analyzing the full working condition simulation, its inlet and outlet differential pressure? P and efficiency? Are evaluated. It is found that its differential pressure decreases with the flow rate and its efficiency reaches its maximum at Qv = 9.87 L/s with a maximum efficiency of 78.32%. (2) We used the LCS in the analysis of vortices at the top of the impeller blades of a high-speed fuel pump. One of the metrics used to describe the LCS in fluid dynamics is the FTLE. The high FTLE region represents the region with the highest and fastest particle trajectory stretching velocity in the fluid flow. We performed a cross-sectional analysis of the FTLE field on the different height surfaces of the impeller on 25% Plane, 50% Plane, and 75% Plane, respectively. And a quarter turn of the rotor rotation was analyzed as a cycle divided into 8 moments. It is found that on 25% Plane, the vortex at the top of the lobe is not obvious, but there are high FTLE values on the shroud surface. On 50% Plane, the lobe top vortex is relatively obvious and the number of vortices is three. The vortex pattern remains stable with the rotating motion of the rotor. At 75% Plane, the lobe top vortex is more visible and its number of vortices increases to about 5 and the vortex morphology is relatively stable. The FTLE ri

Purpose

With the development of fracturing technology, the research of multi-well hydrofracturing becomes the key issue. Frac-hits in multi-well hydrofracturing has an important effect on fracture propagation and final production of fractured well; in the process of hydrofracturing, there are many implement parameters that can affect frac-hits, and previous studies in this area have not systematically targeted the influence of a single parameter on multi-well hydrofracturing. Therefore, it is of great significance to study the occurrence rule and influence of frac-hits for optimizing the design of fracturing wells.

Design/methodology/approach

Based on the proposed numerical models, the effects of different fracturing implement parameters (perforation cluster spacing, well spacing and injection rate) on frac-hits are compared in numerical cases. Through the analysis of fracture network, stress field and microseismic, the effects of different fracturing implement parameters on frac-hits and connections are compared.

Findings

The simulation results show that the effect of perforation cluster spacing and well spacing on frac-hits is greater than that of injection rate. Smaller well spacing makes it easier for fractures between adjacent wells to interact with each other, which increases the risk of frac-hits and reduces the risk of fracture connections. Smaller perforation cluster spacing results in larger individual fracture lengths and greater deflection angles, which makes the possibility of frac-hits and connections greater. The lower the injection rate, the lower the probability of frac-hits.

Originality/value

In this study, the influence of different fracturing implement parameters on frac-hits and connections in multi-well hydrofracturing is studied, and the mechanism of frac-hits and connections is analyzed through fracture network, stress field and microseismic analysis. Different simulation results are compared to optimize fracturing well parameter design and provide reference for engineering application.

Purpose

Firepower conflicts usually decay the firepower plan's enforceability, thus incurring high survival risks. Previous studies have shown little attention to avoiding firepower conflicts during the weapon target assignment process. This research proposes a new constrained optimization model named Firepower Conflict Free WTA (FCFWTA) and designs a Survival Evolution (SE) strategy for Artificial Fish Swarm Algorithm (AFSA) to solve the complex constrained WTA problem. In this way, commanders can get more reliable firepower assignment decision support.

Design/methodology/approach

A new constrained optimization model named Firepower Conflict Free WTA (FCFWTA) is constructed. FCFWTA unifies firepower decision variables for different kinds of weapons and takes the firing time point as a clue for firepower conflict checking. The objective function of FCFWTA is the weighted sum of the minimum threat value rest rate (RRTV), maximum hit efficiency (HE) and minimum latest interception time percentage (PLT). Since previous algorithms do not check and resolve intermediate results during optimization, an adapted strategy named Survival Evolution is designed. It enables making full use of the limited firepower without adjusting the coordination scenario in execution.

Findings

The proposed method offers significant advantages in two aspects. Firstly, it effectively enhances the optimization results of WTA in the absence of firepower conflicts. Evidence from Figure. 6 confirms that without the proposed method, there is a high likelihood of generating invalid outcomes. After implementing firepower conflict check and resolution, there is a substantial degradation in the objective function value. Secondly, the method excels at equitably distributing firepower among multiple targets while also enhancing the overall interception probability, irrespective of the varying complexities presented by different scenarios. This ability to maintain balance and efficiency is crucial for tackling defense-related issues.

Research limitations/implications

Specifically, SE is tailored for MWMT problem under time and space constraints. This approach diverges significantly from conventional MWMT research, which typically focuses solely on ammunition quantity or firing range. Consequently, the primary objective was to verify the efficacy of this method. Test results indicated that SE does not exhibit uniform performance across different algorithms; while it significantly enhances the efficacy with PSO and AFSA, its influence is considerably diminished when applied to GA. It might be attributed to the inherent randomness associated with crossover and mutation, which can increase the likelihood of firepower conflicts, coupled with SE's reorganization of the chromosome.

Originality/value

The work described was

Purpose

Parameter estimation in complex engineering structures typically necessitates repeated calculations using simulation models, leading to significant computational costs. This paper aims to introduce an adaptive multi-output Gaussian process (MOGP) surrogate model for parameter estimation in time-consuming models.

Design/methodology/approach

The MOGP surrogate model is established to replace the computationally expensive finite element method (FEM) analysis during the estimation process. We propose a novel adaptive sampling method for MOGP inspired by the traditional expected improvement (EI) method, aiming to reduce the number of required sample points for building the surrogate model. Two mathematical examples and an application in the back analysis of a concrete arch dam are tested to demonstrate the effectiveness of the proposed method.

Findings

The numerical results show that the proposed method requires a relatively small number of sample points to achieve accurate estimates. The proposed adaptive sampling method combined with the MOGP surrogate model shows an obvious advantage in parameter estimation problems involving expensive-to-evaluate models, particularly those with high-dimensional output.

Originality/value

A novel adaptive sampling method for establishing the MOGP surrogate model is proposed to accelerate the procedure of solving large-scale parameter estimation problems. This modified adaptive sampling method, based on the traditional EI method, is better suited for multi-output problems, making it highly valuable for numerous practical engineering applications.

Purpose

In the present study, the application of a recent damage plasticity model is presented for nonlinear dynamic analysis of the Koyna gravity dam. This is a single surface isotropic damage plasticity concrete model, which is based on the decomposition of stresses and was proposed in a previous study. The theoretical aspects of the model are initially reviewed, and a few preliminary verification examples are illustrated. Thereafter, the HHT-α (i.e. Hilber–Hughes–Taylor) algorithm is presented for nonlinear dynamic analysis of concrete gravity dams.

Design/methodology/approach

Based on the prepared tools, nonlinear behavior of the Koyna Dam is studied by applying the invoked damage plasticity model. For this purpose, three cases are considered for the present study. Case A, which is based on the linear model, is mainly used for comparative purposes. The other two cases (B and C) correspond to the nonlinear (i.e. damage plasticity) model. The basic data for these two cases are similar. However, the employed damping algorithms are different and correspond to constant and variable damping algorithms, respectively. This means that the damping matrix is either kept constant or updated for all iterations of different time increments through the course of analysis.

Findings

The time histories of horizontal displacement at the dam crest were initially compared for the three cases: the linear Case A, and two nonlinear Cases B and C. It was observed that nonlinear cases’ responses begin to deviate from the corresponding linear case after the time of about 4.3 s. However, the amount of change for Case C (i.e. variable damping) was much greater than for Case B (i.e. constant damping). This was manifested initially in the peaks of response. It was also noticed that the period of response changed slightly for Case B in comparison with the linear Case A, while this change was significant for Case C. The obtained tensile and compressive damages were subsequently compared for the two nonlinear cases. For constant damping Case B, it was noticed that tensile damage occurred in the D/S face kink and on the U/S face slightly at a lower elevation. Moreover, it had a scattered nature. However, in variable damping Case C, it was noticed that tensile damage was much more localized and acted similar to a discrete crack. Of course, both cases also show tensile damages at the dam’s heel. In regard to compressive damages, it is observed that low values are occurring for both nonlinear cases as expected.

Originality/value

The application of a recent single surface isotropic damage plasticity concrete model is presented for nonlinear dynamic analysis of the Koyna gravity dam. The nonlinear response of the dam is investigated for two different damping algorithms. Moreover, the influence of variable characteristic length is also investigated i