High-frequency oscillation(HFO)of gridconnected wind power generation systems(WPGS)is one of the most critical issues in recent years that threaten the safe access of WPGS to the grid.Ensuring the WPGS can damp HFO is...High-frequency oscillation(HFO)of gridconnected wind power generation systems(WPGS)is one of the most critical issues in recent years that threaten the safe access of WPGS to the grid.Ensuring the WPGS can damp HFO is becoming more and more vital for the development of wind power.The HFO phenomenon of wind turbines under different scenarios usually has different mechanisms.Hence,engineers need to acquire the working mechanisms of the different HFO damping technologies and select the appropriate one to ensure the effective implementation of oscillation damping in practical engineering.This paper introduces the general assumptions of WPGS when analyzing HFO,systematically summarizes the reasons for the occurrence of HFO in different scenarios,deeply analyses the key points and difficulties of HFO damping under different scenarios,and then compares the technical performances of various types of HFO suppression methods to provide adequate references for engineers in the application of technology.Finally,this paper discusses possible future research difficulties in the problem of HFO,as well as the possible future trends in the demand for HFO damping.展开更多
The purpose of this paper is to study the oscillation of second-order half-linear neutral differential equations with advanced argument of the form(r(t)((y(t)+p(t)y(τ(t)))')^(α))'+q(t)yα(σ(t))=0,t≥t_(0),w...The purpose of this paper is to study the oscillation of second-order half-linear neutral differential equations with advanced argument of the form(r(t)((y(t)+p(t)y(τ(t)))')^(α))'+q(t)yα(σ(t))=0,t≥t_(0),when∫^(∞)r^(−1/α)(s)ds<∞.We obtain sufficient conditions for the oscillation of the studied equations by the inequality principle and the Riccati transformation.An example is provided to illustrate the results.展开更多
In this paper, we consider the following second order retarded differential equations x″(t)+cx′(t)=qx(t-σ)-lx(t-δ) (1) x″(t)+p(t)x(t-τ)=0 (2) We give some sufficient conditions for the oscillation of all solutio...In this paper, we consider the following second order retarded differential equations x″(t)+cx′(t)=qx(t-σ)-lx(t-δ) (1) x″(t)+p(t)x(t-τ)=0 (2) We give some sufficient conditions for the oscillation of all solutions of Eq. (1) in the case where q, ι, σ, δ are positive numbers and c is a real number. And also, we study the asymptotic behavior of the nonoscillatory solutions. If necessary, we give some examples to illustrate our results. At last, we study Eq. (2) with some conditions on p(t).展开更多
A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO mode...A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO model are obtained and the asymptotic behaviour of solution for corresponding problem is considered.展开更多
An analytical technique, namely the homotopy analysis method (HAM), is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small phys...An analytical technique, namely the homotopy analysis method (HAM), is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small physical parameters at all, and thus valid for both weakly and strongly nonlinear problems. In addition, HAM is different from all other analytic techniques in providing a simple way to adjust and control convergence region of the series solution by means of an auxiliary parameter h. In the present paper, a periodic analytic approximations for nonlinear oscillations with parametric excitation are obtained by using HAM, and the results are validated by numerical simulations.展开更多
The EI Nino-southern oscillation (ENSO) is an interannual phenomenon involved in tropical Pacific ocean- atmosphere interactions. In this paper, we develop an asymptotic method of solving the nonlinear equation usin...The EI Nino-southern oscillation (ENSO) is an interannual phenomenon involved in tropical Pacific ocean- atmosphere interactions. In this paper, we develop an asymptotic method of solving the nonlinear equation using the ENSO model. Based on a class of the oscillator of the ENSO model, a approximate solution of the corresponding problem is studied employing the perturbation method.展开更多
In this paper,a nonlinear time transformation method is presented for the analysis of strong nonlinear oscillation systems.This method can be used to study the limit cycle behavior of the autonomous systems and to ana...In this paper,a nonlinear time transformation method is presented for the analysis of strong nonlinear oscillation systems.This method can be used to study the limit cycle behavior of the autonomous systems and to analyze the forced vibration of a strong nonlinear system.展开更多
In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
The utilization of thin plate systems based on acoustic vibration holds significant importance in micro-nano manipulation and the exploration of nonlinear science. This paper focuses on the analysis of an actual thin ...The utilization of thin plate systems based on acoustic vibration holds significant importance in micro-nano manipulation and the exploration of nonlinear science. This paper focuses on the analysis of an actual thin plate system driven by acoustic wave signals. By combining the mechanical analysis of thin plate microelements with the Bubnov–Galerkin integral method, the governing equation for the forced vibration of a square thin plate is derived. Notably,the reaction force of the thin plate vibration system is defined as f=α|w|, resembling Hooke’s law. The energy function and energy level curve of the system are also analyzed. Subsequently, the amplitude–frequency response function of the thin plate oscillator is solved using the harmonic balance method. Through numerical simulations, the amplitude–frequency curves are analyzed for different vibration modes under the influence of various parameters. Furthermore, the paper demonstrates the occurrence of conservative chaotic motions in the thin plate oscillator using theoretical and numerical methods. Dynamics maps illustrating the system’s states are presented to reveal the evolution laws of the system. By exploring the effects of force fields and system energy, the underlying mechanism of chaos is interpreted. Additionally, the phenomenon of chaos in the oscillator can be controlled through the method of velocity and displacement states feedback, which holds significance for engineering applications.展开更多
High voltage fracturing technology was widely used in the field of reservoir reconstruction due to its advantages of being clean, pollution-free, and high-efficiency. However, high-frequency circuit oscillation occurs...High voltage fracturing technology was widely used in the field of reservoir reconstruction due to its advantages of being clean, pollution-free, and high-efficiency. However, high-frequency circuit oscillation occurs during the underwater high voltage pulse discharge process, which brings security risks to the stability of the pulse fracturing system. In order to solve this problem, an underwater pulse power discharge system was established, the circuit oscillation generation conditions were analyzed and the circuit oscillation suppression method was proposed. Firstly, the system structure was introduced and the charging model of the energy storage capacitor was established by the state space average method. Next, the electrode high-voltage breakdown model was established through COMSOL software, the electrode breakdown process was analyzed according to the electron density distribution image, and the plasma channel impedance was estimated based on the conductivity simulation results. Then the underwater pulse power discharge process and the circuit oscillation generation condition were analyzed, and the circuit oscillation suppression strategy of using the thyristor to replace the gas spark switch was proposed. Finally, laboratory experiments were carried out to verify the precision of the theoretical model and the suppression effect of circuit oscillation. The experimental results show that the voltage variation of the energy storage capacitor, the impedance change of the pulse power discharge process, and the equivalent circuit in each discharge stage were consistent with the theoretical model. The proposed oscillation suppression strategy cannot only prevent the damage caused by circuit oscillation but also reduce the damping oscillation time by77.1%, which can greatly improve the stability of the system. This research has potential application value in the field of underwater pulse power discharge for reservoir reconstruction.展开更多
Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices cal...Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices call for an efficient and reliable method to analyze the dynamic response of DE. This remains to be a challenge since the resultant vibration equation of DE, for example, the vibration of a DE balloon considered here is highly nonlinear with higher-order power terms and time-dependent coefficients. Previous efforts toward this goal use largely the numerical integration method with the simple harmonic balance method as a supplement. The numerical integration and the simple harmonic balance method are inefficient for large parametric analysis or with difficulty in improving the solution accuracy. To overcome the weakness of these two methods,we describe formulations of the incremental harmonic balance(IHB) method for periodic forced solutions of such a unique system. Combined with an arc-length continuation technique, the proposed strategy can capture the whole solution branches, both stable and unstable, automatically with any desired accuracy.展开更多
The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equation...The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equations were expanded with respect to a snmll parameter to get the zeroth- and first-order approximate equations. By using the differenl2al quadrature method with only a few grid points, the high-accurate numerical results were obtained.展开更多
The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially i...The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.展开更多
We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, ...We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.展开更多
In the measurement of G with the angular acceleration method,the improved correlation method developed by Wu et al.(Wu W H,Tian Y,Luo J,Shao C G,Xu J H and Wang DH 2016 Rev.Sci.Instrum.87 094501) is used to accurate...In the measurement of G with the angular acceleration method,the improved correlation method developed by Wu et al.(Wu W H,Tian Y,Luo J,Shao C G,Xu J H and Wang DH 2016 Rev.Sci.Instrum.87 094501) is used to accurately estimate the amplitudes of the prominent harmonic components of the gravitational background signal with time-varying frequency.Except the quadratic slow drift,the angular frequency of the gravitational background signal also includes a cosine oscillation coming from the useful angular acceleration signal,which leads to a deviation from the estimated amplitude.We calculate the correction of the cosine oscillation to the amplitude estimation.The result shows that the corrections of the cosine oscillation to the amplitudes of the fundamental frequency and second harmonic components obtained by the improved correlation method are within respective errors.展开更多
For differential equations with piecewise constant arguments of advanced type, numerical stability and oscillations of Runge-Kutta methods are investigated. The necessary and sufficient conditions under which the nume...For differential equations with piecewise constant arguments of advanced type, numerical stability and oscillations of Runge-Kutta methods are investigated. The necessary and sufficient conditions under which the numerical stability region contains the analytic stability region are given. The conditions of oscillations for the Runge-Kutta methods are obtained also. We prove that the Runge-Kutta methods preserve the oscillations of the analytic solution. Moreover, the relationship between stability and oscillations is discussed. Several numerical examples which confirm the results of our analysis are presented.展开更多
This paper studies a delayed air-sea coupled oscillator describing the physical mechanism of El Nino Southern Oscillation. The approximate expansions of the delayed differential equation's solution are obtained succe...This paper studies a delayed air-sea coupled oscillator describing the physical mechanism of El Nino Southern Oscillation. The approximate expansions of the delayed differential equation's solution are obtained successfully by the modified variational iteration method. The numerical results illustrate the effectiveness and correctness of the method by comparing with the exact solution of the reduced model.展开更多
In this paper, a six-order nonlinear dynamic model with three degrees of freedom is presented for the study of the 'fishtail' motion of a Single Point Mooring System. The effect of parameter variations on the ...In this paper, a six-order nonlinear dynamic model with three degrees of freedom is presented for the study of the 'fishtail' motion of a Single Point Mooring System. The effect of parameter variations on the equilibrium state of the system is analyzed. In order to study the stability of the equilibrium state, the mooring-line length l is chosen as a bifurcation parameter, so that all eigenvalues of the Jacobian matrix under different parameters can be worked out, and then the Hopf-bifurcation point can be found. Finally, the Hopf-bifurcation periodic solution of the system is computed.展开更多
In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is deriv...In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.展开更多
The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elasti...The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elastic case. However, the correspon- dence principle becomes invalid when the materials exhibit ageing. To deal with this problem, a second-order two-scale (SOTS) computational method in the time domain is presented to predict the ageing linear viscoelastic performance of composite materials with a periodic structure. First, in the time domain, the SOTS formulation for calcu- lating the effective relaxation modulus and displacement approximate solutions of the ageing viscoelastic problem is formally derived. Error estimates of the displacement ap- proximate solutions for SOTS method are then given. Numerical results obtained by the SOTS method are shown and compared with those by the finite element method in a very fine mesh. Both the analytical and numerical results show that the SOTS computational method is feasible and efficient to predict the ageing linear viscoelastic performance of composite materials with a periodic structure.展开更多
基金supported in part by the Fundamental Research Funds for the Central Universities under Grant 2682023CX019National Natural Science Foundation of China under Grant U23B6007 and Grant 52307141Sichuan Science and Technology Program under Grant 2024NSFSC0115。
文摘High-frequency oscillation(HFO)of gridconnected wind power generation systems(WPGS)is one of the most critical issues in recent years that threaten the safe access of WPGS to the grid.Ensuring the WPGS can damp HFO is becoming more and more vital for the development of wind power.The HFO phenomenon of wind turbines under different scenarios usually has different mechanisms.Hence,engineers need to acquire the working mechanisms of the different HFO damping technologies and select the appropriate one to ensure the effective implementation of oscillation damping in practical engineering.This paper introduces the general assumptions of WPGS when analyzing HFO,systematically summarizes the reasons for the occurrence of HFO in different scenarios,deeply analyses the key points and difficulties of HFO damping under different scenarios,and then compares the technical performances of various types of HFO suppression methods to provide adequate references for engineers in the application of technology.Finally,this paper discusses possible future research difficulties in the problem of HFO,as well as the possible future trends in the demand for HFO damping.
基金This research is supported by the Shandong Provincial Natural Science Foundation of China(ZR2017MA043).
文摘The purpose of this paper is to study the oscillation of second-order half-linear neutral differential equations with advanced argument of the form(r(t)((y(t)+p(t)y(τ(t)))')^(α))'+q(t)yα(σ(t))=0,t≥t_(0),when∫^(∞)r^(−1/α)(s)ds<∞.We obtain sufficient conditions for the oscillation of the studied equations by the inequality principle and the Riccati transformation.An example is provided to illustrate the results.
文摘In this paper, we consider the following second order retarded differential equations x″(t)+cx′(t)=qx(t-σ)-lx(t-δ) (1) x″(t)+p(t)x(t-τ)=0 (2) We give some sufficient conditions for the oscillation of all solutions of Eq. (1) in the case where q, ι, σ, δ are positive numbers and c is a real number. And also, we study the asymptotic behavior of the nonoscillatory solutions. If necessary, we give some examples to illustrate our results. At last, we study Eq. (2) with some conditions on p(t).
基金Project supported by the National Natural Science Foundation of China (Grant Nos 90111011 and 10471039), the National Key Basic Research Special Foundation of China (Grant Nos 2003CB415101-03 and 2004CB418304), the Key Basic Research Foundation of the Chinese Academy of Sciences (Grant No KZCX3-SW-221) and in part by E-Institutes of Shanghai Municipal Education Commission (Grant No N.E03004).
文摘A class of coupled system for the E1 Nifio-Southern Oscillation (ENSO) mechanism is studied. Using the method of variational iteration for perturbation theory, the asymptotic expansions of the solution for ENSO model are obtained and the asymptotic behaviour of solution for corresponding problem is considered.
文摘An analytical technique, namely the homotopy analysis method (HAM), is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small physical parameters at all, and thus valid for both weakly and strongly nonlinear problems. In addition, HAM is different from all other analytic techniques in providing a simple way to adjust and control convergence region of the series solution by means of an auxiliary parameter h. In the present paper, a periodic analytic approximations for nonlinear oscillations with parametric excitation are obtained by using HAM, and the results are validated by numerical simulations.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11071205 and 11101349), the “Strate- gic Priority Research Program-Climate Change: Carbon Budget and Relevant Issues” of the Chinese Academy of Sciences, China (Grant No. XDA01020304), the Natural Science Foundation from the Education Bureau of Anhui Province, China (Grant No. KJ2011A135), and the Natural Science Foundation of Jiangsu Province, China (Grant No. BK2011042).
文摘The EI Nino-southern oscillation (ENSO) is an interannual phenomenon involved in tropical Pacific ocean- atmosphere interactions. In this paper, we develop an asymptotic method of solving the nonlinear equation using the ENSO model. Based on a class of the oscillator of the ENSO model, a approximate solution of the corresponding problem is studied employing the perturbation method.
基金The project partly supported by the Foundation of Zhongshan University Advanced Research Center
文摘In this paper,a nonlinear time transformation method is presented for the analysis of strong nonlinear oscillation systems.This method can be used to study the limit cycle behavior of the autonomous systems and to analyze the forced vibration of a strong nonlinear system.
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 61973172, 62003177, 62103204, 62003175, and 61973175)the Joint Fund of the Ministry of Education for Equipment Pre-research (Grant No. 8091B022133)General Terminal IC Interdisciplinary Science Center of Nankai University。
文摘The utilization of thin plate systems based on acoustic vibration holds significant importance in micro-nano manipulation and the exploration of nonlinear science. This paper focuses on the analysis of an actual thin plate system driven by acoustic wave signals. By combining the mechanical analysis of thin plate microelements with the Bubnov–Galerkin integral method, the governing equation for the forced vibration of a square thin plate is derived. Notably,the reaction force of the thin plate vibration system is defined as f=α|w|, resembling Hooke’s law. The energy function and energy level curve of the system are also analyzed. Subsequently, the amplitude–frequency response function of the thin plate oscillator is solved using the harmonic balance method. Through numerical simulations, the amplitude–frequency curves are analyzed for different vibration modes under the influence of various parameters. Furthermore, the paper demonstrates the occurrence of conservative chaotic motions in the thin plate oscillator using theoretical and numerical methods. Dynamics maps illustrating the system’s states are presented to reveal the evolution laws of the system. By exploring the effects of force fields and system energy, the underlying mechanism of chaos is interpreted. Additionally, the phenomenon of chaos in the oscillator can be controlled through the method of velocity and displacement states feedback, which holds significance for engineering applications.
基金financially supported by the National Science and Technology Major Project(No.2016ZX05034004)。
文摘High voltage fracturing technology was widely used in the field of reservoir reconstruction due to its advantages of being clean, pollution-free, and high-efficiency. However, high-frequency circuit oscillation occurs during the underwater high voltage pulse discharge process, which brings security risks to the stability of the pulse fracturing system. In order to solve this problem, an underwater pulse power discharge system was established, the circuit oscillation generation conditions were analyzed and the circuit oscillation suppression method was proposed. Firstly, the system structure was introduced and the charging model of the energy storage capacitor was established by the state space average method. Next, the electrode high-voltage breakdown model was established through COMSOL software, the electrode breakdown process was analyzed according to the electron density distribution image, and the plasma channel impedance was estimated based on the conductivity simulation results. Then the underwater pulse power discharge process and the circuit oscillation generation condition were analyzed, and the circuit oscillation suppression strategy of using the thyristor to replace the gas spark switch was proposed. Finally, laboratory experiments were carried out to verify the precision of the theoretical model and the suppression effect of circuit oscillation. The experimental results show that the voltage variation of the energy storage capacitor, the impedance change of the pulse power discharge process, and the equivalent circuit in each discharge stage were consistent with the theoretical model. The proposed oscillation suppression strategy cannot only prevent the damage caused by circuit oscillation but also reduce the damping oscillation time by77.1%, which can greatly improve the stability of the system. This research has potential application value in the field of underwater pulse power discharge for reservoir reconstruction.
基金the National Natural Science Foundation of China(Nos.11702215 and11972277)the Natural Science Basic Research Plan in Shaanxi Province of China(Nos.2017JQ5062 and 2018JQ1029)。
文摘Dielectric elastomer(DE) is suitable in soft transducers for broad applications,among which many are subjected to dynamic loadings, either mechanical or electrical or both. The tuning behaviors of these DE devices call for an efficient and reliable method to analyze the dynamic response of DE. This remains to be a challenge since the resultant vibration equation of DE, for example, the vibration of a DE balloon considered here is highly nonlinear with higher-order power terms and time-dependent coefficients. Previous efforts toward this goal use largely the numerical integration method with the simple harmonic balance method as a supplement. The numerical integration and the simple harmonic balance method are inefficient for large parametric analysis or with difficulty in improving the solution accuracy. To overcome the weakness of these two methods,we describe formulations of the incremental harmonic balance(IHB) method for periodic forced solutions of such a unique system. Combined with an arc-length continuation technique, the proposed strategy can capture the whole solution branches, both stable and unstable, automatically with any desired accuracy.
文摘The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equations were expanded with respect to a snmll parameter to get the zeroth- and first-order approximate equations. By using the differenl2al quadrature method with only a few grid points, the high-accurate numerical results were obtained.
基金the National Supercomputer Center in Tianjin for their patient assistance in providing the compilation environment.We thank the editor,Huajian Yao,for handling the manuscript and Mingming Li and another anonymous reviewer for their constructive comments.The research leading to these results has received funding from National Natural Science Foundation of China projects(Grant Nos.92355302 and 42121005)Taishan Scholar projects(Grant No.tspd20210305)others(Grant Nos.XDB0710000,L2324203,XK2023DXC001,LSKJ202204400,and ZR2021ZD09).
文摘The thermal evolution of the Earth’s interior and its dynamic effects are the focus of Earth sciences.However,the commonly adopted grid-based temperature solver is usually prone to numerical oscillations,especially in the presence of sharp thermal gradients,such as when modeling subducting slabs and rising plumes.This phenomenon prohibits the correct representation of thermal evolution and may cause incorrect implications of geodynamic processes.After examining several approaches for removing these numerical oscillations,we show that the Lagrangian method provides an ideal way to solve this problem.In this study,we propose a particle-in-cell method as a strategy for improving the solution to the energy equation and demonstrate its effectiveness in both one-dimensional and three-dimensional thermal problems,as well as in a global spherical simulation with data assimilation.We have implemented this method in the open-source finite-element code CitcomS,which features a spherical coordinate system,distributed memory parallel computing,and data assimilation algorithms.
基金Project supported by the National Natural Science Foundation of China(Grant No.11471262)the National Basic Research Program of China(Grant No.2012CB025904)the State Key Laboratory of Science and Engineering Computing and the Center for High Performance Computing of Northwestern Polytechnical University,China
文摘We study the hyperbolic–parabolic equations with rapidly oscillating coefficients. The formal second-order two-scale asymptotic expansion solutions are constructed by the multiscale asymptotic analysis. In addition, we theoretically explain the importance of the second-order two-scale solution by the error analysis in the pointwise sense. The associated explicit convergence rates are also obtained. Then a second-order two-scale numerical method based on the Newmark scheme is presented to solve the equations. Finally, some numerical examples are used to verify the effectiveness and efficiency of the multiscale numerical algorithm we proposed.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11575160,11175160,11275075,and 11511130011)
文摘In the measurement of G with the angular acceleration method,the improved correlation method developed by Wu et al.(Wu W H,Tian Y,Luo J,Shao C G,Xu J H and Wang DH 2016 Rev.Sci.Instrum.87 094501) is used to accurately estimate the amplitudes of the prominent harmonic components of the gravitational background signal with time-varying frequency.Except the quadratic slow drift,the angular frequency of the gravitational background signal also includes a cosine oscillation coming from the useful angular acceleration signal,which leads to a deviation from the estimated amplitude.We calculate the correction of the cosine oscillation to the amplitude estimation.The result shows that the corrections of the cosine oscillation to the amplitudes of the fundamental frequency and second harmonic components obtained by the improved correlation method are within respective errors.
文摘For differential equations with piecewise constant arguments of advanced type, numerical stability and oscillations of Runge-Kutta methods are investigated. The necessary and sufficient conditions under which the numerical stability region contains the analytic stability region are given. The conditions of oscillations for the Runge-Kutta methods are obtained also. We prove that the Runge-Kutta methods preserve the oscillations of the analytic solution. Moreover, the relationship between stability and oscillations is discussed. Several numerical examples which confirm the results of our analysis are presented.
基金supported by the National Natural Science Foundation of China (Grant Nos.41105063 and 61070041)
文摘This paper studies a delayed air-sea coupled oscillator describing the physical mechanism of El Nino Southern Oscillation. The approximate expansions of the delayed differential equation's solution are obtained successfully by the modified variational iteration method. The numerical results illustrate the effectiveness and correctness of the method by comparing with the exact solution of the reduced model.
文摘In this paper, a six-order nonlinear dynamic model with three degrees of freedom is presented for the study of the 'fishtail' motion of a Single Point Mooring System. The effect of parameter variations on the equilibrium state of the system is analyzed. In order to study the stability of the equilibrium state, the mooring-line length l is chosen as a bifurcation parameter, so that all eigenvalues of the Jacobian matrix under different parameters can be worked out, and then the Hopf-bifurcation point can be found. Finally, the Hopf-bifurcation periodic solution of the system is computed.
基金Project supported by the National Natural Science Foundation of China(Nos.11971303 and 11871330)。
文摘In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.
基金Project supported by the National Natural Science Foundation of China(No.11471262)
文摘The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elastic case. However, the correspon- dence principle becomes invalid when the materials exhibit ageing. To deal with this problem, a second-order two-scale (SOTS) computational method in the time domain is presented to predict the ageing linear viscoelastic performance of composite materials with a periodic structure. First, in the time domain, the SOTS formulation for calcu- lating the effective relaxation modulus and displacement approximate solutions of the ageing viscoelastic problem is formally derived. Error estimates of the displacement ap- proximate solutions for SOTS method are then given. Numerical results obtained by the SOTS method are shown and compared with those by the finite element method in a very fine mesh. Both the analytical and numerical results show that the SOTS computational method is feasible and efficient to predict the ageing linear viscoelastic performance of composite materials with a periodic structure.