• Current Issue
  • Archive
  • Most Downloaded
    Select AllDeselectExport
    Display Method:
    2026,24(2):1-15, DOI: 10.6052/1672-6553-2025-089
    Abstract:
    The complex dynamics of neuronal firing in the brain provide a crucial basis for high order cognitive functions and the emergence of pathological states. Astrocytes play a key role in regulating neural activity on a timescale of seconds, influencing nervous system function through multiple mechanisms such as modulating neurotransmitter levels, ion concentrations, and energy metabolism, as well as responding to exogenous disturbances (e.g., temperature fluctuations and noise). Based on recent biological experimental studies on astrocytes, this paper reviews the latest progress in dynamic modeling of astrocyte regulated neural firing and elaborates on the contributions of astrocytes to synaptic information transmission and synaptic plasticity at the neural network level. This study offers theoretical support for an in depth understanding of astrocytes’roles in cognitive functions such as memory and attention. Additionally, the paper discusses the influence of astrocytes on the dynamic behavior of abnormal neural firing in neurological diseases (e.g., epilepsy), thereby providing potential clinical value for the prevention and treatment of related neurological diseases. Finally, considering current developments in artificial intelligence (AI) technology, the paper outlines future research directions for more comprehensively uncovering the regulatory role of astrocytes in nervous system dynamics through the integration of experimental data and dynamic modeling.
    2026,24(2):16-25, DOI: 10.6052/1672-6553-2025-087
    Abstract:
    Time scales are defined as any non empty closed subset of the real number field, which unifies the treatment of continuous and discrete systems. In this paper, the Herglotz type Vacco dynamics of nonholonomic systems are extended to the time scales, and its Noether symmetry and conservation law are investigated. Firstly, based on the Herglotz variational principle on time scales, the Herglotz type Vacco dynamics equations on time scales are established. Secondly, according to the invariance of Hamilton Herglotz action on time scales under infinitesimal transformations, the Noether symmetry of Herglotz type Vacco dynamics of nonholonomic systems on time scales is defined, and the corresponding Noether identities are presented. Finally, the Noether’s theorem of Herglotz type Vacco dynamics for nonholonomic systems on time scales is proven, and the corresponding conserved quantities are provided. At the conclusion of the paper, two examples are presented to demonstrate the results of theoretical analysis.
    2026,24(2):26-36, DOI: 10.6052/1672-6553-2025-086
    Abstract:
    A numerical model of a Nextel/Kevlar filled protective structure under hypervelocity impact is developed to investigate its structural response and optimize its configuration. Four key design parameters, namely the thicknesses of the front plate, Nextel and Kevlar layers, and rear plate, are chosen as variables. Latin hypercube sampling is employed to generate sample points, and high fidelity simulations are carried out using LS DYNA to obtain the projectile’s kinetic energy dissipation. A Kriging surrogate model is built to approximate the simulation response with reduced computational cost. Based on the surrogate, a multi objective optimization using the NSGA II algorithm is performed, aiming to minimize the areal density and maximize energy dissipation. The optimized structure achieves a 2.56% reduction in areal density and a 1.91% increase in energy loss.
    2026,24(2):37-44, DOI: 10.6052/1672-6553-2025-076
    Abstract:
    Due to the severity and complexity of the working environment, turbine blades often experience fatigue failure and fracture damage under large vibration amplitude. In order to reduce the risk of fatigue failure of turbine blades, the underplatform dampers are typically employed to suppress the blade vibration. However, most experimental studies use weights and wires to simulate the application of normal load, which cannot quantitatively and precisely apply normal load. To address this problem, this paper proposes a precise loading method for underplatform dampers. Moreover, experimental studies on the vibration reduction effects of underplatform dampers are conducted. The study investigates the influence of parameters such as normal load, excitation amplitudes, and friction area ratios on the vibration reduction characteristics of dampers. Finally, a MATLAB Python ABAQUS joint simulation method is used to analyze the vibration effect of the underplatform damper, which validates the effectiveness of the proposed experimental method.
    2026,24(2):45-53, DOI: 10.6052/1672-6553-2025-093
    Abstract:
    Modal localization, owing to its high sensitivity to structural perturbations, has demonstrated unique advantages in the field of microelectromechanical systems (MEMS) sensors. However, traditional modal localization resonators mostly rely on frequency swept excitation, which can only stimulate a single mode and suffer from poor amplitude stability in open loop operation as well as instability issues in dual closed loop driving. These limitations hinder the realization of real time measurement and fast response. To address this challenge, this paper investigates a coupled double beam resonator and proposes a modal localization sensing method based on broadband noise excitation. The proposed method drives multiple modes simultaneously using broadband noise and extracts the variation characteristics of the modal energy distribution through power spectral density analysis, enabling efficient and sensitive perturbation detection. A coupled dynamic model of a dual beam MEMS resonator under broadband noise excitation is established, and numerical simulations together with experimental studies are conducted to comparatively analyze the modal responses under frequency swept and noise driven excitations. The results show that the modal localization effect under broadband noise driven conditions is highly consistent with that under harmonic driving, thereby verifying the feasibility and effectiveness of the proposed method.
    2026,24(2):54-66, DOI: 10.6052/1672-6553-2025-091
    Abstract:
    This paper investigates the nonlinear dynamic behavior of PVC gel cylindrical shell under multiple electromechanical parameters. First, the nonlinear vibration equations of the cylindrical shell structure are theoretically derived based on the Gent hyperelastic material model, and the response and stability of the system under static and dynamic voltages are subsequently discussed. Studies on static voltage reveal that the cylindrical shell exhibits a critical voltage threshold because of the positive feedback effect of the electric field. Exceeding this threshold leads to instability and damage. And the critical voltage is greatly influenced by thickness and boundary conditions. When dynamic sinusoidal voltage is applied, the system shows complex nonlinear vibration characteristics. Through analyses of the time domain responses, phase trajectories,Poincaré sections, bifurcation characteristics and Lyapunov exponents, the presence of periodic vibrations and bifurcation phenomena is confirmed. Numerical simulations show that multi frequency resonance occurs when the excitation frequency changes. Such resonance induces amplitude jumps, which leads to structural damage. By combining nonlinear verification at the same frequency with phase diagram analysis, it is shown that the vibration state has regular consistency under specific parameters. This finding confirms that the vibration response can be controlled by adjusting voltage parameters.
    2026,24(2):67-73, DOI: 10.6052/1672-6553-2025-090
    Abstract:
    This paper studies the dynamic response of wind turbine systems (WTS) under non Gaussian stochastic excitation. Firstly, considering the limitations of traditional Gaussian noise in representing the actual wind speed and system uncertainty, the α stable Lévy noise with heavy tail and pulse characteristics is introduced to establish a more practical WTS stochastic dynamical model. Secondly, based on the theory of stochastic differential, the fractional Fokker Planck Kolmogorov (FPK) equation corresponding to WTS under the excitation of α stable Lévy noise is derived, which precisely describes the evolution law of the transient probability density function (PDF) for the system state. Finally, to effectively solve the fractional partial differential equation, a physics informed neural networks (PINNs) framework is proposed, which takes the physical control equation as the constrained embedding loss function, and can directly learn the space time continuous PDF solution without grid discretization. Numerical experiments show that the PINNs solution is highly consistent with the Monte Carlo simulation results, which verifies the accuracy of this method in solving fractional FPK equations. Meanwhile, PINNs shows much higher computational efficiency than traditional Monte Carlo methods.
    2026,24(2):74-83, DOI: 10.6052/1672-6553-2025-088
    Abstract:
    Wind and rain loads exert a significant influence on the safety of the transmission tower line system. In this paper, the dynamic characteristics of the tower line coupling system in a 110 kV transmission line under wind and rain loads are analyzed. Firstly, the finite element model of a 110 kV transmission tower line system is established in Ansys software, and the modal analysis is carried out. Combined with the theory of Davenport wind speed spectrum and rain load, the wind rain coupling load corresponding to the loading node is generated in Matlab. The dynamic response of the tower line system under this coupling load is studied, and the results of its dynamic response under different wind direction angles and with or without rain load are explored. The results show that the grounding wire increases the stiffness of the tower line system as a whole, but because of its strong ‘galloping effect’in the horizontal direction, the coupling effect in the horizontal direction and the along line direction is different. The wind direction angle of 90 ° is the most unfavorable wind direction angle of the tower line system, and the dynamic response of the tower line system reaches the maximum value. In the case of rainfall of 20 mm/h, the dynamic response range of the tower line system will be increased by about 5 % ~10 % considering the rain load. The effect of rain load on a single tower is almost negligible, which mainly increases the dynamic response of the tower line system by aggravating the galloping effect’of the ground wire.
    2026,24(2):84-94, DOI: 10.6052/1672-6553-2025-078
    Abstract:
    In response to the challenge of efficient prediction of the dynamic behavior of the dovetail connected bladed disk systems involving boundary nonlinearity, rotating effects, complex loads, etc., the fixed interface modal synthesis method is applied to the reduced order modeling of a rotating dovetail connected bladed disk system. By introducing thin layer solid elements on both groove and tenon contact surfaces to capture interference behavior, and considering the rotation induced stiffening and softening effects as well as the dovetail joint induced local load action, the reduced system level model including the disk substructure, dovetail joint zone, and blade substructure under aerodynamic excitation is then established. The influence of the modal truncation numbers of the blade and disk on the first three natural frequencies of the reduced system is discussed, and the effects of rotating speed and friction coefficient on the modal characteristics and vibration responses of both the full and reduced models are compared with each other. The results show that: (1) within the studied parameter range, the maximum deviation of the reduced model in predicting the first three natural frequencies compared to the full model does not exceed 0.6%, and the maximum deviation in predicting the critical speed does not exceed 0.1%; (2) the nonlinearity of the dovetail connection makes the response spectrum of the system to exhibit multiples of the excitation frequency, and a smaller friction coefficient induces slip between the tenon and the groove thus leading to a quasi linear component in the vibration response.
    Select AllDeselectExport
    Display Method:
    2014,12(3):269-273, DOI: 10.6052/1672-6553-2014-043
    [Abstract] (3842) [HTML] (0) [PDF 336.30 K] (22046)
    Abstract:
    Aiming at the problem that the longitudinal modes of structural system of rocket need to be identified from its integral modes in engineering, a method that automatically identifies the longitudinal modes of structural system of rocket was proposed according to the theory of modal effective mass. Taking the vibration characteristics of system with lumped mass as a computing example, applying the finite element software, the beam model of system with lumped mass was established, and the longitudinal modes of the system were automatically identified based on the method. Compared with the system modal information calculated by the method of modal analysis, this automatic identification method not only can accurately identify the longitudinal modes of vibrating system, but also has automatic and high efficiency identification feature. It provides a theoretical basis for the dynamic model of POGO vibrating system in liquid rockets and other model of engineering systems to be accurately and promptly established.
    2014,12(3):193-200, DOI: 10.6052/1672-6553-2014-046
    [Abstract] (3645) [HTML] (0) [PDF 748.46 K] (19552)
    Abstract:
    The spacecraft's ability to adapt to the harsh dynamics environment is critical for the whole space mission. Vibration test control technology is the key part of the dynamic environment test. The current progress, fundamental principles and key techniques development level of the spacecraft and vibration control algorithms overseas were analyzed. The basic ideas, effective ways and suggestions were given to domestic following research.
    2017,15(5):385-405, DOI: 10.6052/1672-6553-2017-039
    [Abstract] (2383) [HTML] (0) [PDF 1.91 M] (19403)
    Abstract:
    In this review article, the growth and related academic communications in the dynamics of multibody systems are firstly surveyed. Then, the recent advances in the numerical algorithms for solving the dynamic equations of flexible multibody systems, the contact/impact dynamics of flexible multibody systems and the deployment of flexible space structures are systematically reviewed, together with several open problems of concern. Finally, some suggestions are made for the prospective researches on the dynamics of flexible multibody systems.
    2014,12(1):18-23, DOI: 10.6052/1672-6553-2013-068
    [Abstract] (3450) [HTML] (0) [PDF 1.13 M] (19000)
    Abstract:
    The electromechanical coupling model of cantilevered piezoelectric harvester was developed by considering the nonlinearities of piezoelectric material, based on Hamilton theory, Rayleigh-Ritz method, Euler-Bernoulli beam theory and constant electrical field across the piezoelectric element. The response characteristics of the system were investigated numerically, and the influences of piezoelectric material nonlinear coefficient on the system response were analyzed. By exploring the nonlinear characteristics of the piezoelectric vibrator near the resonant frequency, the nature of the multi-solutions and jump phenomena in the resonance region was revealed. The results were verified experimentally. which provides a theoretical basis for the study of nonlinear mechanism of piezoelectric power generation system.
    2014,12(3):243-247, DOI: 10.6052/1672-6553-2014-054
    [Abstract] (2874) [HTML] (0) [PDF 1.07 M] (18590)
    Abstract:
    Random vibration test is very important to the aerospace equipment. The windage of ASD outside the frequency bandwidth was analyzed. The reason, premonition, affection and effective way were given. And the commonly used random vibration test and vibration test metrology standard for the ASD outside the frequency bandwidth were analyzed.
    2016,14(2):97-108, DOI: 10.6052/1672-6553-2015-009
    [Abstract] (2896) [HTML] (0) [PDF 1.87 M] (18427)
    Abstract:
    Based on the current research status of multi agents system control theory and technology, the paper makes a detailed overview for unmanned ground systems. From two aspects of behaviour and task cooperative control for multi agents, the relevant theory and application problem is discussed. Moreover, some existed open problems are presented and a possible future development is proposed. For unmanned ground systems, cooperative control will be of great importance in promoting social and military benefits and maximizing the executive function of ground mission.
    2015,13(5):361-366, DOI: 10.6052/1672-6553-2014-064
    [Abstract] (1586) [HTML] (0) [PDF 826.19 K] (18249)
    Abstract:
    Dynamic and control systems often contain uncertain parameters that may result in uncertain predictions. In the interest of quantifying the effects of parameter uncertainties on response variability, this paper develops a stochastic response surface based method for the sensitivity analysis of uncertain parameters. Stochastic response surfaces were firstly constructed to describe the explicit relationships between uncertain parameters and responses. Then partial derivations were performed on the mathematical expressions of stochastic response surfaces in order to obtain sensitivity indices that simultaneously embody the effects of parameter means and standard deviations. Lastly, the developed method has been verified against a numerical cantilever beam containing uncertain geometric and material parameters. The sensitivity analysis results were compared with those given by the analysis of variance method.
    2014,12(3):225-229, DOI: 10.6052/1672-6553-2014-051
    [Abstract] (2822) [HTML] (0) [PDF 1.38 M] (18174)
    Abstract:
    Based on the staggered solution procedure of ANSYS and CFX software, the fluid structure coupling response of projectile during tail slapping has been researched. Structural response was simulated by using FEM and flow field was simulated by using inhomogeneous model and SST turbulence model. Finally, the influences of fluid structure coupling effect have been analyzed and the change law of body stress has been given.
    2013,11(1):12-19, DOI: 10.6052/1672-6553-2013-003
    [Abstract] (2848) [HTML] (0) [PDF 530.95 K] (18047)
    Abstract:
    The dynamics of the whole aero-engine system has always been the important part that cannot be neglected in the research and design of the engine, as the framework of the engine, the vibration of the casing directly reflects the level of the whole aero engine vibration. In this paper, an analysis was made on the research of the problems and fault classification of aero-engine casing dynamics, and an overview was made on the research of present situation, development trend, problems and solutions of the domestic and foreign of casing dynamics, which expanded the present situation of the Inclusiveness problems of aero-engine casing dynamics. Finally, some proposals were put forward for the development of the casing dynamics suitable for our country aero-engine technology level.
    2008,6(4):301-306, DOI:
    [Abstract] (1585) [HTML] (0) [PDF 0.00 Byte] (17715)
    Abstract:
    The problem of P and PI feedback control to a time delay system was investigated,with the emphasis on the determination of the feedback gains that ensured the asymptotical stability of the delayed system. By means of Lambert W function,the feedback gain of P control can be expressed explicitly,so that the optimal feedback gain can be easily obtained. For the system under a PI control,the stable region of the feedback gains was determined on the basis of stability switches and D subdivision,and the optimal feedback gains that enabled the system to admit maximal stable margin were figured out numerically by using Lambert W function. From the viewpoint of computation,the present method is much simpler than the available methods.
    2012,10(1):21-26, DOI:
    [Abstract] (1143) [HTML] (0) [PDF 567.09 K] (17682)
    Abstract:
    A discrete finite dimensional dynamical model was built to describe the space large overall motion of tethered satellite system with an infinite dimensional viscoelastic tether in a long time. The tethered satellite system is a complex nolinear dynamic system. Considering the tether’s viscoelasticity, distributed mass and space form, the established improved bead model can meticulously describe the tether’s vertical and horizontal vibration. According to tether’s characteristic of tensile and not compressive, the slack tether unit model was set up to accurately reflect real stress of tether. The determination of the number of degrees of freedom of the system was studied. Based on numerical integral calculation, the dynamic response was obtained via numerical simulation of the deployment, retrievement and retainment process of tethered satellite system in a long time. The result is convergent. The simulation proves the important role of the stable equilibrium position in the dynamics of tethered space system.
    2014,12(2):183-187, DOI: 10.6052/1672-6553-2014-025
    [Abstract] (1663) [HTML] (0) [PDF 811.52 K] (17412)
    Abstract:
    PMSM due to little harmonic, high precision torque, commonly is used in the servo system and the high performance speed control system. In this paper,the physical model of PMSM is simplified and the mathematical model of the motor is established in order to facilitate research.Thispaper uses id=0 control manner which is the simplest manner is vector control methods,motor electromagnetic torque equation is established based on rotor field oriented vector control.The system model, speed and current control block are build and simulated with MATLAB/Simulink.Simulation results shows that the waveform is consistent with thoretical analysis,the model has fast response and small overshoot.The system runs stably with good dynamic and static characteristics. The simulation makes full use of modularization design. All the parameters and their influence on the system can be changed and observed.It also can easily validate the control strategies and select the most suitable one. So this kind of simulation is good for system design and adjusting and validating.
    2014,12(3):205-209, DOI: 10.6052/1672-6553-2014-059
    [Abstract] (2431) [HTML] (0) [PDF 467.15 K] (14921)
    Abstract:
    This paper investigated the passive nonlinear vibration control method used for energy absorbing in structures of spacecrafts. The structure and the dynamic model of the nonlinear energy sink which could adapt to the space environment were proposed. As nonlinear spring could not be acquired easily in reality, we proposed a new design for the NES based on employing an asymmetric NES force which was generated by two pairs of aligned permanent magnets. Then, the dynamic model for a cantilever beam structure coupled with nonlinear energy sink had been built theoretically. In addition, the passive vibration suppression effect of the nonlinear energy sink on the cantilever beam structure under transient excitation had been analyzed through Galerkin method and numerical analysis method. The results showed that the NES acquired up to 92% dissipation of the system energy imposed by shock excitation, hence the NES could adapt to the space environment and improve the reliability of space system.
    2014,12(3):283-288, DOI: 10.6052/1672-6553-2014-061
    [Abstract] (3507) [HTML] (0) [PDF 479.50 K] (14254)
    Abstract:
    A method was presented to analyze the nonlinear flutter. Based on this method, the flutter characteristics of the high aspect wing were illustrated. The numerical results show that the flutter speed is decreased when the first horizontal bending mode involved. Secondly, this study discussed how the main direction of the composite influenced the character of the nonlinear vibration and flutter, and established the method of the flutter clipping to the high aspect wing. And the result shows that the stiffness of structure can be changed by changing the main direction of the composite. It mainly changes the horizontal bending mode, makes the main direction tend to the trailing edge, and then makes the section line move to the leading edge. Further analyzing the nonlinear flutter reveals that it is the changing of the horizontal bending mode that causes the flutter speed change obviously. And by the section line of this mode moves ahead, the flutter speed will become larger. In the study, two examples were illustrated to validate its truthiness.
    2013,11(4):357-362, DOI: 10.6052/1672-6553-2013-041
    [Abstract] (1744) [HTML] (0) [PDF 350.50 K] (12288)
    Abstract:
    At present, the main task of designing a mounting system of automotive engine powertrain is to select appropriate stiffness, position and angle of mounting components so that free-vibration modal frequency of the mounting system can avert from the exciting-force frequency at the idle speed of the engine and the natural frequency of vibration of the vehicle body and that the decoupling degree of each mode shape is increased as far as possible, so as to improve the vibration-isolation effect of the mounting system. The design of a mounting system based on strict decoupling at predetermined frequencies is to make the modal frequencies of the designed mounting system completely equal to the frequencies predetermined in accordance with the frequency planning of automotive design, and to enable strict decoupling of each mode shape of each mode, i.e., the decoupling degree of vibration energy in every direction equals to 1. Based on a free-vibration equation for a mounting system, this paper presents an equation system for designing a mounting system with strict decoupling at predetermined frequencies, provides a solving method for this equation system by using the theory of generalized inverse matrix or method of constructing function, so as to provide an optimal design method more efficient and simpler than the current modal optimization method of mounting system. Relevant example has validated the correctness of equations and solving method of the strict-decoupling design at predetermined frequencies.
    2016,14(3):247-252, DOI: 10.6052/1672-6553-2015-052
    [Abstract] (2232) [HTML] (0) [PDF 2.18 M] (9536)
    Abstract:
    Tires are the only carrier of the contact between the vehicle and road surface. Their mechanical property is important for analysis and control on the vehicle dynamic response. At present, the tire simulation mostly focuses on the steady state model. But it can not describe the dynamic characteristic of the tire accurately. Therefore, it plays a significant role to add the dynamic tire model in the vehicle dynamics simulation. The tire friction model in the multi body dynamical software ADAMS is static, where the friction is regarded as a static value. However, in actual, the friction between the tire and road surface is dynamic, and it should be a dynamic function of the relative velocity and displacement. To this end, in this paper, the dynamic tire LuGre model using the Matlab/Simulink software is constructed. Through connecting the interface with Adams/Car, co simulation between the vehicle model and the simulink tire model is carried out in order to achieve the dynamic contact between tire and road and improve the accuracy of vehicle system analysis.
    2014,12(1):36-43, DOI: 10.6052/1672-6553-2013-110
    [Abstract] (2539) [HTML] (0) [PDF 2.22 M] (7984)
    Abstract:
    The coupling nonlinear dynamic model of dual rotor system was established by using finite element method, and then the critical speed of revolution and mode shape were calculated by using the software MATLAB. In addition, the unbalance responses of dual rotor system were studied, and the vibration performances in different speeds of dual rotor casing systems were obtained. The research provides a theoretical basis for the design of the dual rotors system in engineering.
    2014,12(3):253-258, DOI: 10.6052/1672-6553-2014-056
    [Abstract] (2416) [HTML] (0) [PDF 1.00 M] (6389)
    Abstract:
    Taking account of the structural stiffness and the low order vibration frequencies, two schemes of multi-objective topology optimization were proposed to obtain the best aircraft lifting-surface structural design. Based on penalized density theory, the scheme one (named as constrain method) is to convert the multi-objective optimization to single-objective optimization by considering the minimum structural mass as the objective with constraints of reference points displacements and the low order vibration frequencies. The scheme two (named as the combination of constrain method and criterion function method) settles the multi-objective optimization by defining combined compliance index (CCI) as the objective, with the constraints of volume fraction and the low order vibration frequencies. The CCI is the function of structural compliance and low order vibration frequencies. Numerical results demonstrate the proposed schemes not only realize reducing the structural mass but also raise the first and second order frequencies.
    2014,12(3):201-204, DOI: 10.6052/1672-6553-2014-048
    [Abstract] (2517) [HTML] (0) [PDF 296.79 K] (6089)
    Abstract:
    Non-linear factors cannot be avoided in the design of aircraft structures. In this paper, a two-degree-of-freedom airfoil and cubic stiffness nonlinearities in pitching degree-of-freedom operating in supersonic flight speed regimes has been analyzed. The averaging method and the theory of flutter were used to analyze the nonlinear dynamic system of the dualistic airfoil in the supersonic flow. Then the correctness of the theoretical calculation was verified by numerical calculation, and the analysis result was given.
    2014,12(1):62-66, DOI: 10.6052/1672-6553-2013-097
    [Abstract] (1973) [HTML] (0) [PDF 832.87 K] (5518)
    Abstract:
    A spline finite point method was presented to study the natural frequency of arch. The displacement mode shape function of the arch free vibration was simulated with a linear combination of cubic B spline. The free vibration frequency equation of arch structures was derived according to Hamilton principle, in which the effect of the dead load was considered. Meanwhile, the effect of the dead load on the natural frequency of arch structures was analyzed. The results show that the natural frequency of arch is reduced. The effect of influence depends on the stiffness of the arch itself. When the arch stiffness is certain, the bigger the rise span ration and the radius to thickness ration, the higher the effect of the dead load on the natural frequency of arch structures.

External Links Links

WeChat

Mobile website