Search results

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.

It is of great significance to study the influence of subgrade filling on permafrost temperature field in permafrost area for the smooth construction and safe operation of railway.

The paper builds up the model for the hydrothermal coupling calculation of permafrost using finite element software COMSOL to study how permafrost temperature field changes in the short term after subgrade filling, on which basis it proposes the method of calculation for the concave distortion of freezing front in the subgrade-covered area.

The results show that the freezing front below the subgrade center sinks due to the thermal effect of subgrade filling, which will trigger hydrothermal erosion in case of sufficient moisture inflows, leading to the thawing settlement or the cracking of the subgrade, etc. The heat output of soil will be hindered the most in case of July filling, in which case the sinking and the distortion of the freezing front is found to be the most severe, which the recovery of the permafrost temperature field, the slowest, constituting the most unfavorable working condition. The concave distortion of the freezing front in the subgrade area increases with the increase in temperature difference between the filler and ground surface, the subgrade height, the subgrade width and the volumetric thermal capacity of filler, while decreases with the increase of the thermal conductivity of filler. Therefore, the filler chose for engineering project shall be of small volumetric thermal capacity, low initial temperature and high thermal conductivity whenever possible.

The concave distortion of the freezing front under different working conditions at different times after filling can be calculated using the method proposed.

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.

We consider, as discretization in space, the nonconforming mesh developed in SUSHI (Scheme Using Stabilization and Hybrid Interfaces) developed in Eymard et al. (2010) for a semi-linear heat equation. The time discretization is performed using a uniform mesh. We are concerned with a nonlinear scheme that has been studied in Bradji (2016) in the context of the general framework GDM (Gradient Discretization Method) (Droniou et al., 2018) which includes SUSHI. We provide sufficient conditions on the size of the spatial mesh and the time step which allow to prove a W1,∞(L2)-error estimate. This error estimate can be viewed as an improvement for the W1,2(L2)-error estimate proved in Bradji (2016). The W1,∞(L2)-error estimate we want to prove in this note was stated without proof in Bradji (2016, Remark 7.2, Page 1302). Its proof is based on a comparison with an appropriately chosen auxiliary finite volume scheme along with the derivation of some new estimates on its solution.

The main objective of the present research is to investigate the benefits of using geogrid reinforcement in minimizing the rate of deterioration of ballasted rail track geometry resting on soft clay and to explore the effect of load amplitude, load frequency, presence of geogrid layer in ballast layer and ballast layer thickness on the behavior of track system. These variables are studied both experimentally and numerically. This paper examines the effect of geogrid reinforced ballast laying on a layer of clayey soil as a subgrade layer, where a half full scale railway tests are conducted as well as a theoretical analysis is performed.

The experimental tests work consists of laboratory model tests to investigate the reduction in the compressibility and stress distribution induced in soft clay under a ballast railway reinforced by geogrid reinforcement subjected to dynamic load. Experimental model based on an approximate half scale for general rail track engineering practice is adopted in this study which is used in Iraqi railways. The investigated parameters are load amplitude, load frequency and presence of geogrid reinforcement layer. A half full-scale railway was constructed for carrying out the tests, which consists of two rails 800 mm in length with three wooden sleepers (900 mm × 90 mm × 90 mm). The ballast was overlying 500 mm thick clay layer. The tests were carried out with and without geogrid reinforcement, the tests were carried out in a well tied steel box of 1.5 m length × 1 m width × 1 m height. A series of laboratory tests were conducted to investigate the response of the ballast and the clay layers where the ballast was reinforced by a geogrid. Settlement in ballast and clay, was measured in reinforced and unreinforced ballast cases. In addition to the laboratory tests, the application of numerical analysis was made by using the finite element program PLAXIS 3D 2013.

It was concluded that the settlement increased with increasing the simulated train load amplitude, there is a sharp increase in settlement up to the cycle 500 and after that, there is a gradual increase to level out between, 2,500 and 4,500 cycles depending on the load frequency. There is a little increase in the induced settlement when the load amplitude increased from 0.5 to 1 ton, but it is higher when the load amplitude increased to 2 ton, the increase in settlement depends on the geogrid existence and the other studied parameters. Both experimental and numerical results showed the same behavior. The effect of load frequency on the settlement ratio is almost constant after 500 cycles. In general, for reinforced cases, the effect of load frequency on the settlement ratio is very small ranging between 0.5 and 2% compared with the unreinforced case.

Increasing the ballast layer thickness from 20 cm to 30 cm leads to decrease the settlement by about 50%. This ascertains the efficiency of ballast in spreading the waves induced by the track.

The paper proposes an efficient and insightful approach for solving neutral delay differential equations (NDDE) with high-frequency inputs. This paper aims to overcome the need to use a very small time step when high frequencies are present. High-frequency signals abound in communication circuits when modulated signals are involved.

The method involves an asymptotic expansion of the solution and each term in the expansion can be determined either from NDDE without oscillatory inputs or recursive equations. Such an approach leads to an efficient algorithm with a performance that improves as the input frequency increases.

An example shall indicate the salient features of the method. Its improved performance shall be shown when the input frequency increases. The example is chosen as it is similar to that in literature concerned with partial element equivalent circuit (PEEC) circuits (Bellen et al., 1999). Its structure shall also be shown to enable insights into the behaviour of the system governed by the differential equation.

The method is novel in its application to NDDE as arises in engineering applications such as those involving PEEC circuits. In addition, the focus of the method is on a technique suitable for high-frequency signals.

This paper studies the lateral stability regulation of intelligent electric vehicle (EV) based on model predictive control (MPC) algorithm.

Firstly, the bicycle model is adopted in the system modelling process. To improve the accuracy, the lateral stiffness of front and rear tire is estimated using the real-time yaw rate acceleration and lateral acceleration of the vehicle based on the vehicle dynamics. Then the constraint of input and output in the model predictive controller is designed. Soft constraints on the lateral speed of the vehicle are designed to guarantee the solved persistent feasibility and enforce the vehicle’s sideslip angle within a safety range.

The simulation results show that the proposed lateral stability controller based on the MPC algorithm can improve the handling and stability performance of the vehicle under complex working conditions.

The MPC schema and the objective function are established. The integrated active front steering/direct yaw moments control strategy is simultaneously adopted in the model. The vehicle’s sideslip angle is chosen as the constraint and is controlled in stable range. The online estimation of tire stiffness is performed. The vehicle’s lateral acceleration and the yaw rate acceleration are modelled into the two-degree-of-freedom equation to solve the tire cornering stiffness in real time. This can ensure the accuracy of model.

The purpose of this paper is to investigate problems in performing stable lane changes and to find a solution to reduce energy consumption of autonomous electric vehicles.

An optimization algorithm, model predictive control (MPC) and Karush–Kuhn–Tucker (KKT) conditions are adopted to resolve the problems of obtaining optimal lane time, tracking dynamic reference and energy-efficient allocation. In this paper, the dynamic constraints of vehicles during lane change are first established based on the longitudinal and lateral force coupling characteristics and the nominal reference trajectory. Then, by optimizing the lane change time, the yaw rate and lateral acceleration that connect with the lane change time are limed. Furthermore, to assure the dynamic properties of autonomous vehicles, the real system inputs under the restraints are obtained by using the MPC method. Based on the gained inputs and the efficient map of brushless direct-current in-wheel motors (BLDC IWMs), the nonlinear cost function which combines vehicle dynamic and energy consumption is given and the KKT-based method is adopted.

The effectiveness of the proposed control system is verified by numerical simulations. Consequently, the proposed control system can successfully achieve stable trajectory planning, which means that the yaw rate and longitudinal and lateral acceleration of vehicle are within stability boundaries, which accomplishes accurate tracking control and decreases obvious energy consumption.

This paper proposes a solution to simultaneously satisfy stable lane change maneuvering and reduction of energy consumption for autonomous electric vehicles. Different from previous path planning researches in which only the geometric constraints are involved, this paper considers vehicle dynamics, and stability boundaries are established in path planning to ensure the feasibility of the generated reference path.

Performing accurate numerical simulations of electrical drives, the precise knowledge of the local magnetic material properties is of utmost importance. Due to the various manufacturing steps, e.g. heat treatment or cutting techniques, the magnetic material properties can strongly vary locally, and the assumption of homogenized global material parameters is no longer feasible. This paper aims to present the general methodology and two different solution strategies for determining the local magnetic material properties using reference and simulation data.

The general methodology combines methods based on measurement, numerical simulation and solving an inverse problem. Therefore, a sensor-actuator system is used to characterize electrical steel sheets locally. Based on the measurement data and results from the finite element simulation, the inverse problem is solved with two different solution strategies. The first one is a quasi Newton method (QNM) using Broyden's update formula to approximate the Jacobian and the second is an adjoint method. For comparison of both methods regarding convergence and efficiency, an artificial example with a linear material model is considered.

The QNM and the adjoint method show similar convergence behavior for two different cutting-edge effects. Furthermore, considering a priori information improved the convergence rate. However, no impact on the stability and the remaining error is observed.

The presented methodology enables a fast and simple determination of the local magnetic material properties of electrical steel sheets without the need for a large number of samples or special preparation procedures.

The authors characterize the conditions under which a country may eventually split and when it splits within an infinite horizon multi-stage differential game.

In contrast to the existing literature, the authors do not assume that after splitting, players will adopt Markovian strategies. Instead, the authors assume that while the splitting country plays Markovian, the remaining coalition remains committed to the collective control of pollution and plays open-loop.

Within a full linear-quadratic model, the authors characterize the optimal strategies. The authors later compare with the outcomes of the case where the splitting country and the remaining coalition play both Markovian. The authors highlight several interesting results in terms of the implications for long-term pollution levels and the duration of coalitions under heterogenous strategies as compared to Markovian behavior.

In this paper, the authors have illustrated the richness of the simplications of enlarging the set of strategies in terms of the emergence of coalitions, their duration and the implied welfare levels per player. Varying only three parameters (the technological gap, pollution damage and coalition payoff share distribution across players), the authors have been able to generate, among other findings, quite different rankings of welfare per player depending on whether the remaining coalitions after split play Markovian or stay precommited to the pre-splitting period decisions.