Structural Crack Identification Based on the Variational Mode Decomposition北大核心CSCD

In order to enrich the bridge damage detection method and further improve the accuracy of bridge damage identification, a detection method for simply supported beams with cracks under dynamic loads was proposed not based on the complete finite element model. Under the premise of not blocking traffic, the method only needs to analyze and deal with the acceleration responses of the simply supported beam span, which reduces the mounting, dismounting and maintenance of sensors in practical engineering. At the same time, based on the model, an analytical formula of the acceleration at the midspan of the simply supported cracked beam was derived. Based on the theoretical derivation, the instantaneous energy and the mean energy difference were constructed through the variational mode decomposition and the Hilbert transform, and these 2 crack identification indexes were used to effectively identify small cracks with a crack depth ratio of only 5%. Then the influences of different wheel loads, environmental noises and damage degrees on detection results were studied. The results show that: ① the instantaneous frequency has a better recognition effect for crack positions; ② the mean energy difference is sensitive to crack depth ratio δ and the wheel load magnitude; ③ this method has strong noise robustness. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

AN EFFECTIVE CONTINUOUS ALGORITHM FOR APPROXIMATE SOLUTIONS OF LARGE SCALE MAX-CUT PROBLEMS

An effective continuous algorithm is proposed to find approximate solutions of NP-hardmax-cut problems.The algorithm relaxes the max-cut problem into a continuous nonlinearprogramming problem by replacing n discrete constraints in the original problem with onesingle continuous constraint.A feasible direction method is designed to solve the resultingnonlinear programming problem.The method employs only the gradient evaluations ofthe objective function,and no any matrix calculations and no line searches are required.This greatly reduces the calculation cost of the method,and is suitable for the solutionof large size max-cut problems.The convergence properties of the proposed method toKKT points of the nonlinear programming are analyzed.If the solution obtained by theproposed method is a global solution of the nonlinear programming problem,the solutionwill provide an upper bound on the max-cut value.Then an approximate solution to themax-cut problem is generated from the solution of the nonlinear programming and providesa lower bound on the max-cut value.Numerical experiments and comparisons on somemax-cut test problems(small and large size)show that the proposed algorithm is efficientto get the exact solutions for all small test problems and well satisfied solutions for mostof the large size test problems with less calculation costs.     

Identification of Pipeline Inner Wall Geometry Based on the POD⁃RBF Method北大核心CSCD

Based on the proper orthogonal decomposition⁃radial basis function (POD⁃RBF), a geometric identification method for pipeline inner wall was proposed to solve the internal corrosion detection problem of natural gas and oil pipelines. In view of the static magnetic field, the simplified finite element model for the pipelines was established, and the variable⁃geometry sample library was constructed, to realize the response prediction of arbitrary geometry by the POD⁃RBF. The proposed method achieves reduced⁃order analysis and avoids repeated solution of the stiffness matrix due to the geometrical change during the identification process. Hence, it can significantly improve the computation efficiency. Finally, the grey wolf optimization (GWO) algorithm was used to optimize the objective function and avoid the calculation of the sensitivity in the process of geometry change. The numerical examples show that, the proposed method has high efficiency and accuracy in the geometric identification of the pipeline inner wall, with good stability even under introduced noises. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Forced Vibration Analysis of Euler-Bernoulli Double-Beam Systems by Means of Green’s Functions北大核心CSCD

Double-curved-beam (DCB) systems are usually seen in many engineering fields. Compared to straight double-beam systems, DCB systems are more efficient in noise and vibration control problems. To obtain closed-form solutions of steady-state forced vibrations of DCB systems, the classical Euler-Bernoulli curved beam (ECB) model was employed to model vibration equations for the DCB systems. Green’s functions and the Laplace transform methods were used to get the closed-form solutions to the vibration equations for the DCB systems. These solutions apply to arbitrary boundary conditions. Numerical tests were conducted to verify the present solutions with related results from previous literatures. Effects of some important geometric and physical parameters on vibration responses and the interaction between the elastic layer stiffness and the DCB system, were discussed. The results show that, the DCB system will degenerate to a straight double-beam system when the 2 radii approach infinity, moreover, the DCB system can be simplified as one comprising a straight beam and a curved beam. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

A Self-Adaptive Alternating Direction Multiplier Method for Frictionless Elastic Contact Problems北大核心CSCD

A self⁃adaptive alternating direction multiplier method was designed for frictionless elastic contact problems. An augmented Lagrange function was introduced for the variational formulation of the problem with an auxiliary variable, to deduce a minimization problem and an equivalent saddle⁃point problem. Then the alter⁃ nating direction multiplier method was used to solve the problem. To enhance the performance of the algo⁃ rithm, a self⁃adaptive rule based on the iterative function on the boundary was proposed to automatically select the proper penalty parameter. The advantage of this algorithm is that, each iteration only needs to solve a linear variational problem and explicitly calculate the auxiliary variable and the Lagrange multiplier. The convergence of the algorithm was analyzed theoretically. The numerical results illustrate the feasibility and effectiveness of the proposed method. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Construction and analysis of the quadratic finite volume methods on tetrahedral meshes

We construct and analyze a family of quadratic finite volume method(FVM) schemes over tetrahedral meshes.In order to prove the stability and the error estimate,we propose the minimum V-angle condition on tetrahedral meshes,and the surface and volume orthogonal conditions on dual meshes.Through the technique of element analysis,the local stability is equivalent to a positive definiteness of a 9 X 9 element matrix,which is difficult to analyze directly or even numerically.With the help of the surf...     

Finite Integral Transform Solutions for Free Vibrations of Rectangular Thin Plates With Mixed Boundary Constraints北大核心CSCD

Analytical solutions, with unique research value, can serve as benchmarks for empirical formulas and numerical methods, a tool for rapid parameter analysis and optimization, and a theoretical basis for experimental designs. Conventional analytical methods, e.g., the Lévy solution method, are only applicable to mechanical problems of plates and shells with opposite simply-supported edges, which, however, may fail to obtain analytical solutions for the issues with complex boundary constraints. In recent years, the finite integral transform method for plate and shell problems was developed to deal with non-Lévy-type plates and shells, but it is still infeasible to solve the mixed boundary constrains-induced complex boundary value problems of higher-order partial differential equations. Herein, for the first time, the finite integral transform method was combined with the sub-domain decomposition technique to solve the free vibrations of rectangular thin plates with mixed boundary constraints. The rectangular plate was first divided into 2 sub-domains according to the mixed boundary constraints, and the 2 sub-domains were solved analytically with the finite integral transform method. Finally, the continuity conditions were introduced to obtain the analytical solution of the original problem. Based on the side spot-welded cantilever plates commonly used in engineering, the free vibration problem of a rectangular thin plate with 1 edge subjected to clamped-simply supported constraints and the other 3 edges free, was analyzed. The obtained natural frequencies and mode shapes are in good agreement with those from the finite element method as well as the solutions in literature, thus verifying the accuracy of the proposed method. The solution procedure of the finite integral transform method can be implemented based on the governing equations without any assumption of the solution form. Therefore, this strict analytical method is widely applicable to complex boundary value problems of higher-order partial differential equations for such mechanical problems of plates and shells. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

A REDUCED FE FORMULATION BASED ON POD METHOD FOR HYPERBOLIC EQUATIONS

A proper orthogonal decomposition(POD) method was successfully used in the reduced-order modeling of complex systems.In this paper,we extend the applications of POD method,namely,apply POD method to a classical finite element(FE) formulation for second-order hyperbolic equations with real practical applied background,establish a reduced FE formulation with lower dimensions and high enough accuracy,and provide the error estimates between the reduced FE solutions and the classical FE solutions and the implementation of algorithm for solving reduced FE formulation so as to provide scientific theoretic basis for service applications.Some numerical examples illustrate the fact that the results of numerical computation are consistent with theoretical conclusions.Moreover,it is shown that the reduced FE formulation based on POD method is feasible and efficient for solving FE formulation for second-order hyperbolic equations.     

关于一类二阶矩阵方程的敏度分析

In this paper, the sensitivity of the solution for a class of quadratic matrix equation which arises in the analysis of structural systems and vibration problems is discussed. With Brouwer fixed piont theory, the perturbation of the quadratic matrix equation is analyzed and two computational perturbation bounds are derived. Then a Rice condition number of some kind of solutions is given using the analytic expansion method. Two examples are presented in the last part.     

ELASTIC MEMBRANE EQUATION WITH MEMORY TERM AND NONLINEAR BOUNDARY DAMPING:GLOBAL EXISTENCE,DECAY AND BLOWUP OF THE SOLUTION

In this paper we consider the Elastic membrane equation with memory term and nonlinear boundary damping.Under some appropriate assumptions on the relaxation function h and with certain initial data,the global existence of solutions and a general decay for the energy are established using the multiplier technique.Also,we show that a nonlinear source of polynomial type is able to force solutions to blow up in finite time even in presence of a nonlinear damping.     

The Panel Method for Rigid Section Group Added Mass Coefficients and Its Application to Fuel Assemblies北大核心CSCD

A numerical method based on the panel method was developed to calculate the added fluid mass coefficients of the rigid section group with arbitrarily complex shapes, and was successfully applied to the PWR fuel assemblies. The variation law of added mass coefficients with position deviations was analyzed in the seismic test of 1×5 fuel assemblies. The results show that, this method is suitable for the calculation of the added mass coefficients of rigid section groups with complex continuous boundaries. Compared with the gap between assemblies, the gap between baffles and assemblies has a dominant influence on the added mass coefficient. Regardless of the position deviation, the sum of the added mass coefficients of all assemblies and baffles in the assumed motion direction is approximately equal to –1, and that in the perpendicular direction is approximately equal to 0. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

EXISTENCE OF PERIODIC SOLUTIONS TO A SYSTEM WITH FUNCTIONAL RESPONSE ON TIME SCALES

This paper investigates the existence of periodic solutions of a three-species food-chain diffusive system with Beddington-DeAngelis functional responses and time delays in a two-patch environment on time scales. By using a continuation theorem based on coincidence degree theory, we obtain sufficient criteria for the existence of periodic solutions for the system. Moreover, when the time scale T is chosen as R or Z, the existence of the periodic solutions of the corresponding continuous and discrete models follows. Therefore, the method is unified to provide the existence of the desired solutions for continuous differential equations and discrete difference equations.     

Effects of Cone Angles on Nonlinear Vibration Responses of Functionally Graded Shells北大核心CSCD

The nonlinear vibration responses of functionally graded materials (FGMs) shells with different cone angles under external loads were studied. Firstly, the Voigt model was employed to describe the physical properties along the thickness direction of FGMs conical shells. Then, the motion equations were derived based on the 1st-order shear deformation theory, the von Kármán geometric nonlinearity and Hamilton’s principle. Next, the Galerkin method was applied to discretize the motion equations and the governing equations were simplified into a 1DOF nonlinear vibration differential equation under Volmir’s assumption. Finally, the nonlinear motion equations were solved with the harmonic balance method and the Runge-Kutta method, and the amplitude frequency response characteristic curves of the FGMs conical shells were obtained. The effects of different material distribution functions and different ceramic volume fraction exponents on the amplitude frequency response curves of conical shells were discussed. The bifurcation diagrams of conical shells with different cone angles, as well as time process diagrams and phase diagrams for different excitation amplitudes, were described. The motion characteristics were characterized by Poincaré maps. The results show that, the FGMs conical shells present the nonlinear characteristics of hardening springs. The chaotic motions of the FGMs conical shells are restrained and not prone to motion instability with the increase of the cone angle. The FGMs conical shell present a process from the periodic motion to the multi-periodic motion and then to chaos with the increase of the excitation amplitude. © 2022 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

AN ITERATIVE HYBRIDIZED MIXED FINITE ELEMENT METHOD FOR ELLIPTIC INTERFACE PROBLEMS WITH STRONGLY DISCONTINUOUS COEFFICIENTS

An iterative algorithm is proposed and analyzed based on a hybridized mized finite element method for numerically solving two-phase generalized Stefan interface problems with strongly discontinuous solutions,conormal derivatives,and coefficients.This algorithm iteratively solves small problems for each single phase with good accuracy and exchange information at the interface to advance the iteration until convergence ,following the idea of Schwarz Alternating Methods,Error estimates are derived to show that this algorithm always converges provided that relaxation parameters are suitably chosen,Numeric exper-iments with matching and non-matching grids at the interface from different phases are performed to show the accuracy of the method for capturing discontinuities in the solutions and coefficients.In contrast to standard numerical methods,the accuracy of our method does not seem to deteriorate as the coefficient discontinuity increases.     

A Quasi-Interpolation Satisfying Quadratic Polynomial Reproduction with Radial Basis Functions

In this paper,a new quasi-interpolation with radial basis functions which satis- fies quadratic polynomial reproduction is constructed on the infinite set of equally spaced data.A new basis function is constructed by making convolution integral with a constructed spline and a given radial basis function.In particular,for twicely differ- entiable function the proposed method provides better approximation and also takes care of derivatives approximation.     

Exact Confidence Limits for the Parameter of an Exponential Distribution in the Accelerated Life Tests under Type-ⅠCensoring

Life data frequently arise in many reliability studies,such as accelerated life tests studies.This paper considers the part of life data where failure and censoring observations may exist.To develop statistical methods and theory for the analysis of these data,a new approach was proposed to obtain the exact lower and upper confidence limits for the mean life of the exponential distribution with Type-I censoring data.It is assumed that the acceleration factor is a random variable,and that the distribution of the acceleration factor is known from some empirical information or the meta analysis.A method for constructing the lower and upper confidence limits for the parameter based on an ordering relation among the sample space was proposed.Simulation studies and analyses of two examples suggest that the proposed method performed well.     

Elastic Analysis of Anisotropic Rotating Sandwich Circular Ring With a Functionally Graded Transition Region北大核心CSCD

The elastic analysis of anisotropic rotating sandwich ring with a functionally graded transition region was carried out. Like the shell sandwich structure in nature, the ring is composed of 3 well⁃bonded regions, of which the inner and outer regions are made of homogeneous anisotropic materials, and the intermediate transi⁃ tion region is made of a material with arbitrary⁃gradient properties along the radial direction. Based on the boundary conditions and the continuity conditions at the interface, the 2nd Fredholm integral equation for the radial stress was obtained with the integral equation method, then the stress and displacement fields of the sandwich ring structure were obtained through numerical solution. The distributions of the stress and displace⁃ ment fields in the sandwich ring structure were given. Different gradient changes encountered in engineering practice can be solved only through substitution of the corresponding function model. The effectiveness and ac⁃ curacy of the integral equation method were verified through comparison of the numerical solutions with the ex⁃ act ones for a special power function gradient variation form. The more general Voigt function model was adopt⁃ ed for the intermediate transition region, and the influences of the anisotropy degree, the gradient parameter, and the thickness on the stress and displacement fields were analyzed. The proposed Fredholm integral equation method provides a powerful tool for the optimal design of anisotropic functionally graded materials and sand⁃ wich ring structures. The numerical results make a theoretical guidance for the safety design of anisotropic functionally graded sandwich ring structures. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Convergence Rate to Stationary Solutions for Boltzmann Equation with External Force

For the Boltzmann equation with an external force in the form of the gradient of a potential function in space variable, the stability of its stationary solutions as local Maxwellians was studied by S. Ukai et al. (2005) through the energy method. Based on this stability analysis and some techniques on analyzing the convergence rates to stationary solutions for the compressible Navier-Stokes equations, in this paper, we study the convergence rate to the above stationary solutions for the Boltzmann equation which is a fundamental equation in statistical physics for non-equilibrium rarefied gas. By combining the dissipation from the viscosity and heat conductivity on the fluid components and the dissipation on the non-fluid component through the celebrated H-theorem, a convergence rate of the same order as the one for the compressible Navier-Stokes is obtained by constructing some energy functionals.     

A second order isoparametric finite element method for elliptic interface problems

A second order isoparametric finite element method (IPFEM) is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the immersed curved interface is discontinuous. Based on an initial Cartesian mesh, a mesh optimization strategy is presented by employing curved boundary elements at the interface, and an incomplete quadratic finite element space is constructed on the optimized mesh. It turns out that the number of curved boundary elements is far less than that of the straight one, and the total degree of freedom is almost the same as the uniform Cartesian mesh. Numerical examples with simple and complicated geometrical interfaces demonstrate the efficiency of the proposed method.     

OPTIMAL CONVERGENCE RATES OF LANDAU EQUATION WITH EXTERNAL FORCING IN THE WHOLE SPACE

In this paper, we combine the method of constructing the compensating function introduced by Kawashima and the standard energy method for the study on the Landau equation with external forcing. Both the global existence of solutions near the time asymptotic states which are local Maxwellians and the optimal convergence rates are obtained. The method used here has its own advantage for this kind of studies because it does not involve the spectrum analysis of the corresponding linearized operator.