An Efficient Algorithm for Calculating Aircraft RCS Based on the Geometrical Characteristics

Taking into account the influences of scatterer geometrical shapes on induced currents, an algorithm, termed the sparse-matrix method （SMM）, is proposed to calculate radar cross section （RCS） of aircraft configuration. Based on the geometrical characteristics and the method of moment （MOM）, the SMM points out that the strong current coupling zone could be predefined according to the shape of scatterers. Two geometrical parameters, the surface curvature and the electrical space between the field position and source position, are deducted to distinguish the dominant current coupling. Then the strong current coupling is computed to construct an impedance matrix having sparse nature, which is solved to compute RCS. The efficiency and feasibility of the SMM are demonstrated by computing electromagnetic scattering of some kinds of shapes such as a cone-sphere with a gap, a bi-arc column and a stealth aircraft configuration. The numerical results show that： （1） the accuracy of SMM is satisfied, as compared with MOM, and the computational time it spends is only about 8% of the MOM; （2） with the electrical space considered, making another allowance for the surface curvature can reduce the computation time by 9.5%.

A Density Functional Theory Study of the Geometric and Electronic Structure of MgF_2 (110) Surface

Abstract＂ Ab initio density functional theory （DFT） was employed to study geometric and electronic structure of MgF2 （110） surface. Three different clean surface models have been considered. The results show that the surface terminated with one-layer F has the smallest relaxation and the lowest surface energy, which indicates this model is the most energetically favorable structure of MgF2（110） surface. Furthermore, the electronic properties are also discussed from the point of density of states and charge density. Analysis of electronic structure shows that the band gap of the surface is significantly narrowed with respect to the bulk. The electrons of the surface exhibit strong locality and larger effective mass.

Influence of joint spacing and rock characteristics on the toppling stability of cut rock slope through a simplified limit equilibrium method

Toppling failure of rock mass/soil slope is an important geological and environmental problem.Clarifying its failure mechanism under different conditions has great significance in engineering.The toppling failure of a cutting slope occurred in a hydropower station in Kyushu,Japan illustrates that the joint characteristic played a significant role in the occurrence of rock slope tipping failure.Thus,in order to consider the mechanical properties of jointed rock mass and the influence of geometric conditions,a simplified analytical approach based on the limit equilibrium method for modeling the flexural toppling of cut rock slopes is proposed to consider the influence of the mechanical properties and geometry condition of jointed rock mass.The theoretical solution is compared with the numerical solution taking Kyushu Hydropower Station in Japan as one case,and it is found that the theoretical solution obtained by the simplified analysis method is consistent with the numerical analytical solution,thus verifying the accuracy of the simplified method.Meanwhile,the Goodman-Bray approach conventionally used in engineering practice is improved according to the analytical results.The results show that the allowable slope angle may be obtained by the improved Goodman-Bray approach considering the joint spacing,the joint frictional angle and the tensile strength of rock mass together.

THE NONLINEAR STABILITY OF PLANE PARALLEL SHEAR FLOWS WITH RESPECT TO TILTED PERTURBATIONS

The nonlinear stability of plane parallel shear flows with respect to tilted perturbations is studied by energy methods.Tilted perturbation refers to the fact that perturbations form an angleθ∈(0,π/2)with the direction of the basic flows.By defining an energy functional,it is proven that plane parallel shear flows are unconditionally nonlinearly exponentially stable for tilted streamwise perturbation when the Reynolds number is below a certain critical value and the boundary conditions are either rigid or stress-free.In the case of stress-free boundaries,by taking advantage of the poloidal-toroidal decomposition of a solenoidal field to define energy functionals,it can be even shown that plane parallel shear flows are unconditionally nonlinearly exponentially stable for all Reynolds numbers,where the tilted perturbation can be either spanwise or streamwise.

Stability and Time-Step Constraints of Implicit-Explicit Runge-Kutta Methods for the Linearized Korteweg-de Vries Equation

This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.

TOPOLOGY SYNTHESIS OF GEOMETRICALLY NONLINEAR COMPLIANT MECHANISMS USING MESHLESS METHODS

This paper presents a new method for topology optimization of geometrical nonlinear compliant mechanisms using the element-free Galerkin method （EFGM）. The EFGM is employed as an alternative scheme to numerically solve the state equations by fully taking advantage of its capability in dealing with large displacement problems. In the meshless method, the imposition of essential boundary conditions is also addressed. The popularly studied solid isotropic material with the penalization （SIMP） scheme is used to represent the nonlinear dependence between material properties and regularized discrete densities. The output displacement is regarded as the objective function and the adjoint method is applied to finding the sensitivity of the design functions. As a result, the optimization of compliant mechanisms is mathematically established as a nonlinear programming problem, to which the method of moving asymptotes （MMA） belonging to the sequential convex programming can be applied. The availability of the present method is finally demonstrated with several widely investigated numerical examples.

Slope Stability Evaluations by Limit Equilibrium and Finite Element Methods Applied to a Railway in the Moroccan Rif

Since 1930, the analysis of slope stability is done according to the limit equilibrium approach. Several methods were developed of which certain remain applicable because of their simplicity. However, major disadvantages of these methods are (1) they do not take into account the soil behavior and (2) the complex cases cannot be studied with precision. The use of the finite elements in calculations of stability has to overcome the weakness of the traditional methods. An analysis of stability was applied to a slope, of complex geometry, composed of alternating sandstone and marls using finite elements and limit equilibrium methods. The calculation of the safety factors did not note any significant difference between the two approaches. Various calculations carried out illustrate perfectly benefits that can be gained from modeling the behavior by the finite elements method. In the finite elements analysis, the shape of deformations localization in the slope is nearly circular and confirms the shape of the failure line which constitutes the basic assumption of the analytical methods. The integration of the constitutive laws of soils and the use of field’s results tests in finite elements models predict the failure mode, to better approach the real behavior of slope soil formations and to optimize its reinforcement.

THE METHOD OF MULTIPLE SCALES APPLIED TO THE NONLINEAR STABILITY PROBLEM OF A TRUNCATED SHALLOW SPHERICAL SHELL OF VARIABLE THICKNESS WITH THE LARGE GEOMETRICAL PARAMETER

Using the modified method of multiple scales, the nonlinear stability of a truncated shallow spherical shell of variable thickness with a nondeformable rigid body at the center under compound loads is investigated. When the geometrical parameter k is larger, the uniformly valid asymptotic solutions of this problem are obtained and the remainder terms are estimated.

H-stability of the Runge-Kutta methods with general variable stepsize for system of pantograph equations with two delay terms

This paper deals with H-stability of the Runge-Kutta methods with a general variable stepsize for the system of pantograph equations with two delay terms. It is shown that the Runge-Kutta methods with a regular matrix A are H-stable if and only if the modulus of the stability function at infinity is less than 1.

ASYMPTOTICAL STABILITY OF NEUTRAL REACTION-DIFFUSION EQUATIONS WITH PCAS AND THEIR FINITE ELEMENT METHODS

This paper focuses on the analytical and numerical asymptotical stability of neutral reaction-diffusion equations with piecewise continuous arguments.First,for the analytical solutions of the equations,we derive their expressions and asymptotical stability criteria.Second,for the semi-discrete and one-parameter fully-discrete finite element methods solving the above equations,we work out the sufficient conditions for assuring that the finite element solutions are asymptotically stable.Finally,with a typical example with numerical experiments,we illustrate the applicability of the obtained theoretical results.

THE STABILITY OF θ-METHODS FOR PANTOGRAPH DELAY DIFFERENTIAL EQUATIONS

This paper deals with the numerical solution of initial value problems for pantograph differential equations with variable delays. We investigate the stability of one leg θ-methods in the numerical solution of these problems. Sufficient conditions for the asymptotic stability of θ-methods are given by Fourier analysis and Ergodic theory.

Evaluation of the effect of geometrical parameters on stope probability of failure in the open stoping method using numerical modeling

Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes which are subjected to stress-induced failure.The probability of failure(POF)depends on many factors,of which the geometry of an open stope is especially important.In this study,a methodology is proposed to assess the effect of stope geometrical parameters on the POF,using numerical modelling.Different ranges for each input parameter are defined according to previous surveys on open stope geometry in a number of Canadian underground mines.A Monte-Carlo simulation technique is combined with the finite difference code FLAC3D,to generate model realizations containing stopes with different geometrical features.The probability of failure(POF)for different categories of stope geometry,is calculated by considering two modes of failure;relaxation-related gravity driven(tensile)failure and rock mass brittle failure.The individual and interactive effects of stope geometrical parameters on the POF,are analyzed using a general multi-level factorial design.Finally,mathematical optimization techniques are employed to estimate the most stable stope conditions,by determining the optimal ranges for each stope’s geometrical parameter.

Comprehensive analysis of slope stability and determination of stable slopes in the Chador-Malu iron ore mine using numerical and limit equilibrium methods

One of the critical aspects in mine design is slope stability analysis and the determination of stable slopes. In the Chador- Malu iron ore mine, one of the most important iron ore mines in central Iran, it was considered vital to perform a comprehensive slope stability analysis. At first, we divided the existing rock hosting pit into six zones and a geotechnical map was prepared. Then, the value of MRMR （Mining Rock Mass Rating） was determined for each zone. Owing to the fact that the Chador-Malu iron ore mine is located in a highly tectonic area and the rock mass completely crushed, the Hoek-Brown failure criterion was found suitable to estimate geo-mechanical parameters. After that, the value of cohesion （c） and friction angle （tp） were calculated for different geotechnical zones and relative graphs and equations were derived as a function of slope height. The stability analyses using numerical and limit equilibrium methods showed that some instability problems might occur by increasing the slope height. Therefore, stable slopes for each geotechnical zone and prepared sections were calculated and presented as a function of slope height.

THE STABILITY OF LINEAR MULTISTEP METHODS FOR SYSTEMS OF DELAY DIFFERENTIAL EQUATIONS

This paper deals with the numerical solution of initial value problems for systems of differential equations with a delay argument. The numerical stability of a linear multistep method is investigated by analysing the solution of the lest equation y’(t)=Ay(t) + By(1-t),where A,B denote constant complex N×N-matrices,and t】0.We investigate carefully the characterization of the stability region.

Geometric interpretation of several classical iterative methods for linear system of equations and diverse relaxation parameter of the SOR method

Two kinds of iterative methods are designed to solve the linear system of equations, we obtain a new interpretation in terms of a geometric concept. Therefore, we have a better insight into the essence of the iterative methods and provide a reference for further study and design. Finally, a new iterative method is designed named as the diverse relaxation parameter of the SOR method which, in particular, demonstrates the geometric characteristics. Many examples prove that the method is quite effective.

A FIXED POINT APPROACH TO THE STABILITY OF A GENERALIZED APOLLONIUS TYPE QUADRATIC FUNCTIONAL EQUATION

Using the fixed point method, this article proves the Hyers-Ulam-Rassias stability of a generalized Apollonius type quadratic functional equation in Banach spaces. The conditions of these results are demonstrated by the quadratic functional equation of Apollonius type.

THE NUMERICAL STABILITY OF THE BLOCK θ-METHODS FOR DELAY DIFFERENTIAL EQUATIONS

This paper focuses on the numerical stability of the block θ methods adapted to differential equations with a delay argument. For the block θ methods, an interpolation procedure is introduced which leads to the numerical processes that satisfy an important asymptotic stability condition related to the class of test problems y′(t)=ay(t)+by(t-τ) with a,b∈C, Re(a)<-|b| and τ>0. We prove that the block θ method is GP stable if and only if the method is A stable for ordinary differential equations. Furthermore, it is proved that the P and GP stability are equivalent for the block θ method.

Numerical Analysis of Balanced Methods for the Impulsive Stochastic Differential Equations

Positive results are proved here about the ability of balanced methods to reproduce the mean square stability of the impulsive stochastic differential equations. It is shown that the balanced methods with strong convergence can preserve the mean square stability with the sufficiently small stepsize. Weak variants and their mean square stability are also considered. Several numerical experiments are given for illustration and show that the fully implicit methods are superior to those of the explicit methods in terms of mean-square stabilities for relatively large stepsizes especially.

Stability Loss of Rotating Elastoplastic Discs of the Specific Form

A method of calculating a possible stability loss by a rotating circular annular disc of variable thickness is suggested within the theory of perfect plasticity with the help of small parameter method. A characteristic equation for a critical radius of a plastic zone is obtained as a first approximation. The formula for the critical angular velocity, determining the stability loss of the disc according to the self-balanced form, is derived. The method using which we can take into account the disc’s geometry and loading parameters is also specified. The efficiency of the proposed method is shown in Section 5 while considering an illustrative example. The values of critical angular velocity of rotating are found numerically for different parameters of the disc.

Transverse Vibration and Stability Analysis of Circular Plate Subjected to Follower Force and Thermal Load

Transverse vibration and stability analysis of circular plate subjected to follower force and thermal load are analyzed.Based on the thin plate theory in involving the variable temperature,the differential equation of transverse vibration for the axisymmetric circular plate subjected to follower force and thermal load is established.Then,the differential equation of vibration and corresponding boundary conditions are discretized by the differential quadrature method.Meanwhile,the generalized eigenvalue under three different boundary conditions are calculated.In this case,the change curve of the first order dimensionless complex frequency of the circular plate subjected to the follower force in the different conditions with the variable temperature coefficient and temperature load is analyzed.The stability and corresponding critical loads of the circular plate subjected to follower force and thermal load with simply supported edge,clamped edge and free edge are discussed.The results provide theoretical basis for improving the dynamic stability of the circular plate.