This paper investigates the magnetohydrodynamic （MHD） boundary layer flow of an incompressible upper-convected Maxwell （UCM） fluid over a porous stretching surface. Similarity transformations are used to reduce th...This paper investigates the magnetohydrodynamic （MHD） boundary layer flow of an incompressible upper-convected Maxwell （UCM） fluid over a porous stretching surface. Similarity transformations are used to reduce the governing partial differential equations into a kind of nonlinear ordinary differential equations. The nonlinear prob- lem is solved by using the successive Taylor series linearization method （STSLM）. The computations for velocity components are carried out for the emerging parameters. The numerical values of the skin friction coefficient are presented and analyzed for various parameters of interest in the problem.展开更多
Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this cha...Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.展开更多
In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the ...In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.展开更多
A high-speed train travelling from the open air into a narrow tunnel will cause the“sonic boom”at tunnel exit.When the maglev train’s speed reaches 600 km/h,the train-tunnel aerodynamic effect is intensified,so a n...A high-speed train travelling from the open air into a narrow tunnel will cause the“sonic boom”at tunnel exit.When the maglev train’s speed reaches 600 km/h,the train-tunnel aerodynamic effect is intensified,so a new mitigation method is urgently expected to be explored.This study proposed a novel asymptotic linear method(ALM)for micro pressure wave(MPW)mitigation to achieve a constant gradient of initial c ompression waves(ICWs),via a study with various open ratios on hoods.The properties of ICWs and MPWs under various open ratios of hoods were analyzed.The results show that as the open ratio increases,the MPW amplitude at the tunnel exit initially decreases before rising.At the open ratio of 2.28%,the slope of the ICW curve is linearly coincident with a supposed straight line in the ALM,which further reduces the MPW amplitude by 26.9%at 20 m and 20.0%at 50 m from the exit,as compared to the unvented hood.Therefore,the proposed method effectively mitigates MPW and quickly determines the upper limit of alleviation for the MPW amplitude at a fixed train-tunnel operation condition.All achievements provide a ne w potential measure for the adaptive design of tunnel hoods.展开更多
In this paper we present a filter-successive linearization method with trust region for solutions of nonlinear semidefinite programming. Such a method is based on the concept of filter for nonlinear programming introd...In this paper we present a filter-successive linearization method with trust region for solutions of nonlinear semidefinite programming. Such a method is based on the concept of filter for nonlinear programming introduced by Fletcher and Leyffer in 2002. We describe the new algorithm and prove its global convergence under weaker assumptions. Some numerical results are reported and show that the new method is potentially efficient.展开更多
Feedforward multi layer neural networks have very strong mapping capability that is based on the non linearity of the activation function, however, the non linearity of the activation function can cause the multiple ...Feedforward multi layer neural networks have very strong mapping capability that is based on the non linearity of the activation function, however, the non linearity of the activation function can cause the multiple local minima on the learning error surfaces, which affect the learning rate and solving optimal weights. This paper proposes a learning method linearizing non linearity of the activation function and discusses its merits and demerits theoretically.展开更多
In this paper, we develop a novel alternating linearization method for solving convex minimization whose objective function is the sum of two separable functions. The motivation of the paper is to extend the recent wo...In this paper, we develop a novel alternating linearization method for solving convex minimization whose objective function is the sum of two separable functions. The motivation of the paper is to extend the recent work Goldfarb et al.(2013) to cope with more generic convex minimization. For the proposed method,both the separable objective functions and the auxiliary penalty terms are linearized. Provided that the separable objective functions belong to C1,1(Rn), we prove the O(1/?) arithmetical complexity of the new method. Some preliminary numerical simulations involving image processing and compressive sensing are conducted.展开更多
An efficient approach is proposed for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations.The concentrated plastic hinges,described by the Bouc-Wen model,are as...An efficient approach is proposed for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations.The concentrated plastic hinges,described by the Bouc-Wen model,are assumed to occur at the two ends of a linear-elastic beam element.The auxiliary differential equations governing the plastic rotational displacements and their corresponding hysteretic displacements are replaced with linearized differential equations.Then,the two sets of equations of motion for the original nonlinear system can be reduced to an expanded-order equivalent linearized equation of motion for equivalent linear systems.To solve the equation of motion for equivalent linear systems,the nonstationary random vibration analysis is carried out based on the explicit time-domain method with high efficiency.Finally,the proposed treatment method for initial values of equivalent parameters is investigated in conjunction with parallel computing technology,which provides a new way of obtaining the equivalent linear systems at different time instants.Based on the explicit time-domain method,the key responses of interest of the converged equivalent linear system can be calculated through dimension reduction analysis with high efficiency.Numerical examples indicate that the proposed approach has high computational efficiency,and shows good applicability to weak nonlinear and medium-intensity nonlinear systems.展开更多
In this paper, two kinds of chaotic systems are controlled respectively with and without time-delay to eliminate their chaotic behaviors. First of all, according to the first-order approximation method and the stabili...In this paper, two kinds of chaotic systems are controlled respectively with and without time-delay to eliminate their chaotic behaviors. First of all, according to the first-order approximation method and the stabilization condition of the linear system, one linear feedback controller is structured to control the chaotic system without time-delay, its chaotic behavior is eliminated and stabilized to its equilibrium. After that, based on the first-order approximation method, the Lyapunov stability theorem, and the matrix inequality theory, the other linear feedback controller is structured to control the chaotic system with time-delay and make it stabilized at its equilibrium. Finally, two numerical examples are given to illustrate the correctness and effectiveness of the two linear feedback controllers.展开更多
To further understand the performance of the energy harvesters under the influence of the wind force and the random excitation,this paper investigates the stochastic response of the bio-inspired energy harvesters subj...To further understand the performance of the energy harvesters under the influence of the wind force and the random excitation,this paper investigates the stochastic response of the bio-inspired energy harvesters subjected to Gaussian white noise and galloping excitation,simulating the flapping pattern of a seagull and its interaction with wind force.The equivalent linearization method is utilized to convert the original nonlinear model into the Itôstochastic differential equation by minimizing the mean squared error.Then,the second-order steady-state moments about the displacement,velocity,and voltage are derived by combining the moment analysis theory.The theoretical results are simulated numerically to analyze the stochastic response performance under different noise intensities,wind speeds,stiffness coefficients,and electromechanical coupling coefficients,time domain analysis is also conducted to study the performance of the harvester with different parameters.The results reveal that the mean square displacement and voltage increase with increasing the noise intensity and wind speed,larger absolute values of stiffness coefficient correspond to smaller mean square displacement and voltage,and larger electromechanical coupling coefficients can enhance the mean square voltage.Finally,the influence of wind speed and electromechanical coupling coefficient on the stationary probability density function(SPDF)is investigated,revealing the existence of a bimodal distribution under varying environmental conditions.展开更多
The long-term Mittag-Leffler stability of solutions to multi-term timefractional diffusion equations with constant coefficients was rigorously established,which demonstrated that the algebraic decay rate of the soluti...The long-term Mittag-Leffler stability of solutions to multi-term timefractional diffusion equations with constant coefficients was rigorously established,which demonstrated that the algebraic decay rate of the solution,characterized by||u_(n)||L^(2)(Ω)=O(t^(−αs)) as t→∞,is determined by the minimum order α_(s) of the time-fractional derivatives.Building on this foundational result,this article pursues two primary objectives.First,we introduce a strongly A-stable fractional linear multistep method and derive the numerical stability region for the governing equation.Second,we rigorously prove the long-term decay rate of the numerical solution through a detailed singularity analysis of its generating function.Notably,the numerical decay rate||u_(n)||L^(2)(Ω)=O(t_(n)^(−α_(s)) as t_(n)→∞aligns precisely with the continuous case.Theoretical findings are further validated through comprehensive numerical simulations,underscoring the robustness of our proposed method.展开更多
The analysis of natural vibration characteristics has become one of important steps of the manufacture and dynamic design in the aerospace industry. This paper presents a new scenario called virtual cutting in the con...The analysis of natural vibration characteristics has become one of important steps of the manufacture and dynamic design in the aerospace industry. This paper presents a new scenario called virtual cutting in the context of the transfer matrix method of linear multibody systems closed- loop topology for computing the free vibration characteristics of elastically coupled flexible launch vehicle boosters. In this approach, the coupled system is idealized as a triple-beam system-like structure coupled by linear translational springs, where a non-uniform free-free Euler-Bemoulli beam is used. A large thrust-to-weight ratio leads to large axial accelera- tions that result in an axial inertia load distribution from nose to tail. Consequently, it causes the development of significant compressive forces along the length of the launch vehicle. Therefore, it is important to take into account this effect in the transverse vibration model. This scenario does not need the global dynamics equations of a system, and it has high computational efficiency and low memory requirements. The validity of the presented scenario is achieved through com- parison to other approaches published in the literature.展开更多
This paper deals with the stability of linear multistep methods for multidimensional differential systems with distributed delays. The delay-dependent stability of linear multistep methods with compound quadrature rul...This paper deals with the stability of linear multistep methods for multidimensional differential systems with distributed delays. The delay-dependent stability of linear multistep methods with compound quadrature rules is studied. Several new sufficient criteria of delay-dependent stability are obtained by means of the argument principle. An algorithm is provided to check delay-dependent stability. An example that illustrates the effectiveness of the derived theoretical results is given.展开更多
Significant wave height is an important criterion in designing coastal and offshore structures.Based on the orthogonality principle, the linear mean square estimation method is applied to calculate significant wave he...Significant wave height is an important criterion in designing coastal and offshore structures.Based on the orthogonality principle, the linear mean square estimation method is applied to calculate significant wave height in this paper.Twenty-eight-year time series of wave data collected from three ocean buoys near San Francisco along the California coast are analyzed.It is proved theoretically that the computation error will be reduced by using as many measured data as possible for the calculation of significant wave height.Measured significant wave height at one buoy location is compared with the calculated value based on the data from two other adjacent buoys.The results indicate that the linear mean square estimation method can be well applied to the calculation and prediction of significant wave height in coastal regions.展开更多
We consider the interior inverse scattering problem for recovering the shape of a penetrable partially coated cavity with external obstacles from the knowledge of measured scattered waves due to point sources.In the f...We consider the interior inverse scattering problem for recovering the shape of a penetrable partially coated cavity with external obstacles from the knowledge of measured scattered waves due to point sources.In the first part,we obtain the well-posedness of the direct scattering problem by the variational method.In the second part,we establish the mathematical basis of the linear sampling method to recover both the shape of the cavity,and the shape of the external obstacle,however the exterior transmission eigenvalue problem also plays a key role in the discussion of this paper.展开更多
In this paper, we consider the inverse scattering by chiral obstacle in electromagnetic fields, and prove that the linear sampling method is also effective to determine the support of a chiral obstacle from the noisy ...In this paper, we consider the inverse scattering by chiral obstacle in electromagnetic fields, and prove that the linear sampling method is also effective to determine the support of a chiral obstacle from the noisy far field data.展开更多
The stability analysis of linear multistep methods for the numerical solutions of the systems of generalized neutral delay differential equations is discussed. The stability behaviour of linear multistep methods was a...The stability analysis of linear multistep methods for the numerical solutions of the systems of generalized neutral delay differential equations is discussed. The stability behaviour of linear multistep methods was analysed for the solution of the generalized system of linear neutral test equations, After the establishment of a sufficient condition for asymptotic stability of the solutions of the generalized system, it is shown that a linear multistep method is NGP(G)-stable if and only if it is A-stable.展开更多
A class of preconditioned iterative methods, i.e., preconditioned generalized accelerated overrelaxation （GAOR） methods, is proposed to solve linear systems based on a class of weighted linear least squares problems...A class of preconditioned iterative methods, i.e., preconditioned generalized accelerated overrelaxation （GAOR） methods, is proposed to solve linear systems based on a class of weighted linear least squares problems. The convergence and comparison results are obtained. The comparison results show that the convergence rate of the preconditioned iterative methods is better than that of the original methods. Furthermore, the effectiveness of the proposed methods is shown in the numerical experiment.展开更多
We discuss estimates for the rate of convergence of the method of successive subspace corrections in terms of condition number estimate for the method of parallel subspace corrections.We provide upper bounds and in a ...We discuss estimates for the rate of convergence of the method of successive subspace corrections in terms of condition number estimate for the method of parallel subspace corrections.We provide upper bounds and in a special case,a lower bound for preconditioners defined via the method of successive subspace corrections.展开更多
C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolati...C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolation （NNI）, with interpolation realized to nodal function and nodal gradient values, so that the essential boundary conditions （EBCs） can be imposed directly in a Galerkin scheme for partial differential equations （PDEs）. In the present paper, C^1 NEM for strain gradient linear elasticity is constructed, and sev- eral typical examples which have analytical solutions are presented to illustrate the effectiveness of the constructed method. In its application to microstructures, the size effects of bending stiffness and stress concentration factor （SCF） are studied for microspeciem and microgripper, respectively. It is observed that the size effects become rather strong when the width of spring for microgripper, the radius of circular perforation and the long axis of elliptical perforation for microspeciem come close to the material characteristic length scales. For the U-shaped notch, the size effects decline obviously with increasing notch radius, and decline mildly with increasing length of notch.展开更多
文摘This paper investigates the magnetohydrodynamic （MHD） boundary layer flow of an incompressible upper-convected Maxwell （UCM） fluid over a porous stretching surface. Similarity transformations are used to reduce the governing partial differential equations into a kind of nonlinear ordinary differential equations. The nonlinear prob- lem is solved by using the successive Taylor series linearization method （STSLM）. The computations for velocity components are carried out for the emerging parameters. The numerical values of the skin friction coefficient are presented and analyzed for various parameters of interest in the problem.
基金funded by the National Key R&D Program of China(Grant No.2022YFC2903904)the National Natural Science Foundation of China(Grant Nos.51904057 and U1906208).
文摘Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.
基金Supported by the Beijing Municipal Science&Technology Commission(Z211100004421012),the Key Reaserch and Development Pro⁃gram of China(2022YFF0605902)。
文摘In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.
基金Project(24A0006)supported by the Key Project of Scientific Research Fund of Hunan Provincial Department of Education,ChinaProject(2024JJ5430)supported by the Natural Science Foundation of Hunan Province,ChinaProjects(2024JK2045,2023RC3061)supported by the Science and Technology Innovation Program of Hunan Province,China。
文摘A high-speed train travelling from the open air into a narrow tunnel will cause the“sonic boom”at tunnel exit.When the maglev train’s speed reaches 600 km/h,the train-tunnel aerodynamic effect is intensified,so a new mitigation method is urgently expected to be explored.This study proposed a novel asymptotic linear method(ALM)for micro pressure wave(MPW)mitigation to achieve a constant gradient of initial c ompression waves(ICWs),via a study with various open ratios on hoods.The properties of ICWs and MPWs under various open ratios of hoods were analyzed.The results show that as the open ratio increases,the MPW amplitude at the tunnel exit initially decreases before rising.At the open ratio of 2.28%,the slope of the ICW curve is linearly coincident with a supposed straight line in the ALM,which further reduces the MPW amplitude by 26.9%at 20 m and 20.0%at 50 m from the exit,as compared to the unvented hood.Therefore,the proposed method effectively mitigates MPW and quickly determines the upper limit of alleviation for the MPW amplitude at a fixed train-tunnel operation condition.All achievements provide a ne w potential measure for the adaptive design of tunnel hoods.
基金supported by National Natural Science Foundation of China (Grant No. 10871098)Science Foundation of Jiangsu Province (Grant No. BK2006214)
文摘In this paper we present a filter-successive linearization method with trust region for solutions of nonlinear semidefinite programming. Such a method is based on the concept of filter for nonlinear programming introduced by Fletcher and Leyffer in 2002. We describe the new algorithm and prove its global convergence under weaker assumptions. Some numerical results are reported and show that the new method is potentially efficient.
文摘Feedforward multi layer neural networks have very strong mapping capability that is based on the non linearity of the activation function, however, the non linearity of the activation function can cause the multiple local minima on the learning error surfaces, which affect the learning rate and solving optimal weights. This paper proposes a learning method linearizing non linearity of the activation function and discusses its merits and demerits theoretically.
基金supported by National Natural Science Foundation of China(Grant Nos.11301055 and 11401315)Natural Science Foundation of Jiangsu Province(Grant No.BK2009397)the Fundamental Research Funds for the Central Universities(Grant No.ZYGX2013J103)
文摘In this paper, we develop a novel alternating linearization method for solving convex minimization whose objective function is the sum of two separable functions. The motivation of the paper is to extend the recent work Goldfarb et al.(2013) to cope with more generic convex minimization. For the proposed method,both the separable objective functions and the auxiliary penalty terms are linearized. Provided that the separable objective functions belong to C1,1(Rn), we prove the O(1/?) arithmetical complexity of the new method. Some preliminary numerical simulations involving image processing and compressive sensing are conducted.
基金Fundamental Research Funds for the Central Universities under Grant No.2682022CX072the Research and Development Plan in Key Areas of Guangdong Province under Grant No.2020B0202010008。
文摘An efficient approach is proposed for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations.The concentrated plastic hinges,described by the Bouc-Wen model,are assumed to occur at the two ends of a linear-elastic beam element.The auxiliary differential equations governing the plastic rotational displacements and their corresponding hysteretic displacements are replaced with linearized differential equations.Then,the two sets of equations of motion for the original nonlinear system can be reduced to an expanded-order equivalent linearized equation of motion for equivalent linear systems.To solve the equation of motion for equivalent linear systems,the nonstationary random vibration analysis is carried out based on the explicit time-domain method with high efficiency.Finally,the proposed treatment method for initial values of equivalent parameters is investigated in conjunction with parallel computing technology,which provides a new way of obtaining the equivalent linear systems at different time instants.Based on the explicit time-domain method,the key responses of interest of the converged equivalent linear system can be calculated through dimension reduction analysis with high efficiency.Numerical examples indicate that the proposed approach has high computational efficiency,and shows good applicability to weak nonlinear and medium-intensity nonlinear systems.
基金Supported by the National Natural Science Foundation of China (61863022)the Natural Science Foundation of Gansu Province(20JR10RA329)Scientific Research and Innovation Fund Project of Gansu University of Chinese Medicine in 2019 (2019KCYB-10)。
文摘In this paper, two kinds of chaotic systems are controlled respectively with and without time-delay to eliminate their chaotic behaviors. First of all, according to the first-order approximation method and the stabilization condition of the linear system, one linear feedback controller is structured to control the chaotic system without time-delay, its chaotic behavior is eliminated and stabilized to its equilibrium. After that, based on the first-order approximation method, the Lyapunov stability theorem, and the matrix inequality theory, the other linear feedback controller is structured to control the chaotic system with time-delay and make it stabilized at its equilibrium. Finally, two numerical examples are given to illustrate the correctness and effectiveness of the two linear feedback controllers.
文摘To further understand the performance of the energy harvesters under the influence of the wind force and the random excitation,this paper investigates the stochastic response of the bio-inspired energy harvesters subjected to Gaussian white noise and galloping excitation,simulating the flapping pattern of a seagull and its interaction with wind force.The equivalent linearization method is utilized to convert the original nonlinear model into the Itôstochastic differential equation by minimizing the mean squared error.Then,the second-order steady-state moments about the displacement,velocity,and voltage are derived by combining the moment analysis theory.The theoretical results are simulated numerically to analyze the stochastic response performance under different noise intensities,wind speeds,stiffness coefficients,and electromechanical coupling coefficients,time domain analysis is also conducted to study the performance of the harvester with different parameters.The results reveal that the mean square displacement and voltage increase with increasing the noise intensity and wind speed,larger absolute values of stiffness coefficient correspond to smaller mean square displacement and voltage,and larger electromechanical coupling coefficients can enhance the mean square voltage.Finally,the influence of wind speed and electromechanical coupling coefficient on the stationary probability density function(SPDF)is investigated,revealing the existence of a bimodal distribution under varying environmental conditions.
文摘The long-term Mittag-Leffler stability of solutions to multi-term timefractional diffusion equations with constant coefficients was rigorously established,which demonstrated that the algebraic decay rate of the solution,characterized by||u_(n)||L^(2)(Ω)=O(t^(−αs)) as t→∞,is determined by the minimum order α_(s) of the time-fractional derivatives.Building on this foundational result,this article pursues two primary objectives.First,we introduce a strongly A-stable fractional linear multistep method and derive the numerical stability region for the governing equation.Second,we rigorously prove the long-term decay rate of the numerical solution through a detailed singularity analysis of its generating function.Notably,the numerical decay rate||u_(n)||L^(2)(Ω)=O(t_(n)^(−α_(s)) as t_(n)→∞aligns precisely with the continuous case.Theoretical findings are further validated through comprehensive numerical simulations,underscoring the robustness of our proposed method.
基金supported by the Research Fund for the Doctoral Program of Higher Education of China(Grants 20113219110025,20133219110037)the National Natural Science Foundation of China(Grants 11102089,61304137)the Program for New Century Excellent Talents in University(NCET-10-0075)
文摘The analysis of natural vibration characteristics has become one of important steps of the manufacture and dynamic design in the aerospace industry. This paper presents a new scenario called virtual cutting in the context of the transfer matrix method of linear multibody systems closed- loop topology for computing the free vibration characteristics of elastically coupled flexible launch vehicle boosters. In this approach, the coupled system is idealized as a triple-beam system-like structure coupled by linear translational springs, where a non-uniform free-free Euler-Bemoulli beam is used. A large thrust-to-weight ratio leads to large axial accelera- tions that result in an axial inertia load distribution from nose to tail. Consequently, it causes the development of significant compressive forces along the length of the launch vehicle. Therefore, it is important to take into account this effect in the transverse vibration model. This scenario does not need the global dynamics equations of a system, and it has high computational efficiency and low memory requirements. The validity of the presented scenario is achieved through com- parison to other approaches published in the literature.
基金Project supported by the National Natural Science Foundation of China(No.11471217)
文摘This paper deals with the stability of linear multistep methods for multidimensional differential systems with distributed delays. The delay-dependent stability of linear multistep methods with compound quadrature rules is studied. Several new sufficient criteria of delay-dependent stability are obtained by means of the argument principle. An algorithm is provided to check delay-dependent stability. An example that illustrates the effectiveness of the derived theoretical results is given.
基金support for this study was provided by the National Natural Science Foundation of China (No.40776006)Research Fund for the Doctoral Program of Higher Education of China (Grant No.20060423009)the Science and Technology Development Program of Shandong Province (Grant No.2008GGB01099)
文摘Significant wave height is an important criterion in designing coastal and offshore structures.Based on the orthogonality principle, the linear mean square estimation method is applied to calculate significant wave height in this paper.Twenty-eight-year time series of wave data collected from three ocean buoys near San Francisco along the California coast are analyzed.It is proved theoretically that the computation error will be reduced by using as many measured data as possible for the calculation of significant wave height.Measured significant wave height at one buoy location is compared with the calculated value based on the data from two other adjacent buoys.The results indicate that the linear mean square estimation method can be well applied to the calculation and prediction of significant wave height in coastal regions.
基金supported by the Natural Science Foundation of Xinjiang Uygur Autonomous Region of China(2019D01A05)supported by the NSFC(11571132)。
文摘We consider the interior inverse scattering problem for recovering the shape of a penetrable partially coated cavity with external obstacles from the knowledge of measured scattered waves due to point sources.In the first part,we obtain the well-posedness of the direct scattering problem by the variational method.In the second part,we establish the mathematical basis of the linear sampling method to recover both the shape of the cavity,and the shape of the external obstacle,however the exterior transmission eigenvalue problem also plays a key role in the discussion of this paper.
文摘In this paper, we consider the inverse scattering by chiral obstacle in electromagnetic fields, and prove that the linear sampling method is also effective to determine the support of a chiral obstacle from the noisy far field data.
文摘The stability analysis of linear multistep methods for the numerical solutions of the systems of generalized neutral delay differential equations is discussed. The stability behaviour of linear multistep methods was analysed for the solution of the generalized system of linear neutral test equations, After the establishment of a sufficient condition for asymptotic stability of the solutions of the generalized system, it is shown that a linear multistep method is NGP(G)-stable if and only if it is A-stable.
基金supported by the National Natural Science Foundation of China (No. 11071033)the Fundamental Research Funds for the Central Universities (No. 090405013)
文摘A class of preconditioned iterative methods, i.e., preconditioned generalized accelerated overrelaxation （GAOR） methods, is proposed to solve linear systems based on a class of weighted linear least squares problems. The convergence and comparison results are obtained. The comparison results show that the convergence rate of the preconditioned iterative methods is better than that of the original methods. Furthermore, the effectiveness of the proposed methods is shown in the numerical experiment.
文摘We discuss estimates for the rate of convergence of the method of successive subspace corrections in terms of condition number estimate for the method of parallel subspace corrections.We provide upper bounds and in a special case,a lower bound for preconditioners defined via the method of successive subspace corrections.
基金supported by the SDUST Spring Bud (2009AZZ021)Taian Science and Technology Development (20112001)
文摘C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolation （NNI）, with interpolation realized to nodal function and nodal gradient values, so that the essential boundary conditions （EBCs） can be imposed directly in a Galerkin scheme for partial differential equations （PDEs）. In the present paper, C^1 NEM for strain gradient linear elasticity is constructed, and sev- eral typical examples which have analytical solutions are presented to illustrate the effectiveness of the constructed method. In its application to microstructures, the size effects of bending stiffness and stress concentration factor （SCF） are studied for microspeciem and microgripper, respectively. It is observed that the size effects become rather strong when the width of spring for microgripper, the radius of circular perforation and the long axis of elliptical perforation for microspeciem come close to the material characteristic length scales. For the U-shaped notch, the size effects decline obviously with increasing notch radius, and decline mildly with increasing length of notch.
