The circular membrane solar array has attracted extensive attention due to its high storage ratio and strong power supply capability. In order to adjust the tension of large film structures, a tension adjusting device composed of ropes and springs is usually introduced, and its mechanical characteristics are highly nonlinear, with the effects rarely studied yet. Aimed at the tension adjustment, a mechanism model was proposed. The nonlinear dynamics equation for the 2DOF system was established with the Lagrangian energy method. With an engineering prototype as the example, the responses of the tension mechanism with unsymmetrical ribs under resonance excitation were studied. The results show that, the change of the excitation amplitude has an important influence on the characteristics of the beat response of the system. Consequently, the responses of the system exhibit chaotic, almost periodic and multifold periodic phenomena. The research makes an important reference to the parameter design of tension mechanisms.

The radial basis function partition of unity (RBF-PU) method  was applied to obtain the numerical solution of 2D nonlocal diffusion and peridynamic problems. The main idea is to partition the original domain into several patches, use the RBF approximation on each local domain, and then give weighting to obtain the global approximation of the unknown function. The radial basis function method based on the strong form of the equation has many advantages, such as avoiding an additional layer of integral calculation, no need to deal with intersections of neighborhoods with the mesh, and easiness of implementation. The numerical results show that, this method can solve nonlocal diffusion equations and peridynamic equations accurately and efficiently.

The grazing-induced non-smooth dynamical behaviors of single-degree-of-freedom cantilever beam systems with bilateral elastic constraints were studied. Firstly, based on the dynamical equations for the cantilever beam under elastic impacts and the definition of grazing points, the existence condition for the bilateral grazing periodic motion was analyzed. Secondly, the zero-velocity Poincaré section was selected to derive the high-order discontinuous mapping with parameters near bilateral grazing orbits. Then a new composite piecewise normal form mapping was established through combination of the smooth flow mapping and the high-order discontinuous mapping. Finally, the validity of the high-order mapping was verified through comparison of the bifurcation diagram of the low-order mapping with that of the high-order mapping, and the grazing dynamics of the cantilever beam under elastic impacts were further revealed through numerical simulation.

The superconducting thin film is a kind of multilayer structure prepared by chemical coating. As a conductive functional structure material with excellent performance, its structural integrity is directly related to the current-carrying capacity. During the preparation of superconducting thin films, it is hard to avoid the interface cracks between the superconducting layer and the metal substrate. In this case, along with the current-carrying operation, the strength of the interface crack in an external magnetic field makes a key problem. Therefore, based on the theory of flux through the thin film and the linear elastic fracture, an analytical model was established for the strength of the interface crack between the superconducting film and the substrate. The effects of the viscous flux flow on the stress field and the energy release rate at the crack tip were obtained. The results show that, the higher the flux flow velocity is, the greater the stress and the energy release rate at the crack tip of the interface will be, which will lead to crack propagation along the interface. The work is helpful for the analysis of interface cracks mentioned above.

The anti-plane problem of asymmetric collinear interface cracks emanating from a circular hole in 1D hexagonal quasicrystal bi-materials was studied. With the Stroh formula and the complex function method, the complex potential functions under the coupling action of the phonon field and the phason field were obtained. The analytical expressions of the stress intensity factor (SIF) and the energy release rate (ERR) at the crack tip were given. The effects of the circular hole radius and the crack length on the SIF, and the effects of the coupling coefficient, the phonon field stress and the phason field stress on the ERR, were discussed. The results show that, the SIF tends to be stable with the increase of the right crack length for a constant circular hole radius. For a certain phason field stress value, the ERR reaches the minimum value, which indicates that a specific phason field stress can inhibit the crack growth.

The vibration suppression performance of a damping composite structure depends on the material layout and the damping material properties. A topology optimization method was proposed for damping material microstructures with varied volume constraints, to obtain the damping material microstructure with desired properties under the smallest material consumption. Based on the homogenization method, a 3D finite element model for the damping material was established, and the effective elastic matrix of the damping material was formulated. The Hashin-Shtrikman bounds theory was used inversely to estimate the volume fraction bound of the damping material corresponding to the desired effective modulus, and a movement criterion for volume constraint bounds of damping materials was constructed. Then the optimization problem of achieving the desired properties of damping materials with microstructures was converted to another problem of maximizing the desired modulus under volume constraints, and a topology optimization model for the damping material microstructure was established. The optimality criteria method was employed to update the design variables, and the optimized topology configurations of damping material microstructures were obtained. The feasibility and effectiveness of the proposed method were verified with several numerical examples, and the influences of the initial configurations, the mesh density and Young’s modulus on the microstructure configurations of the damping material were also discussed.

The peridynamic differential operator (PDDO) theory was introduced to solve the 2D transient heat conduction problems. The heat conduction equation and the boundary condition were reformulated from the local differential form to the non-local integral form. Then the Lagrangian multiplier method and the variational analysis were used, and a non-local model for 2D transient heat conduction problems was established. Through error and convergence analysis, the accuracy of this model was verified in comparison with the results of other numerical methods. The model was further applied to solve the 2D transient heat conduction problems of plates with irregular boundaries and micro-defects (cracks and holes). The results show high accuracy, wide applicability and good convergence of this method, which provides new insights into the 2D transient heat conduction problems.

Aimed at the phenomenon of competing for water resources between young vegetation and adult vegetation in arid and semi-arid areas, a vegetation-soil water dynamic model with intraspecfic competition delay was established. The conditions for the existence of an unique vegetation survival equilibrium and the local stability of the vegetation extinction equilibrium were analyzed. The generating conditions for Hopf bifurcating periodic solutions of non-spatial and spatial systems were given, respectively. The periodic oscillation pattern appearing in the vegetation evolution with time was numerically simulated. Through the parameter sensitivity analysis, the rainfall and the vegetation growth rate were found to have significant influences on the generation, the amplitude and the period of this pattern, while the effects of evaporation was found to be the least significant. The results indicate that, the rainfall and the vegetation properties have profound impacts on the evolution and development of vegetation in arid and semi-arid areas. The introduction of spatial diffusion inhibits the occurrence of this pattern, but doesn’t affect the amplitude and the period. The work explains the phenomenon of vegetation periodic oscillation widely observed in nature, and provides theoretical supports for the sustainable development of the vegetation system.

A class of multi-objective fractional semi-infinite optimization problems with uncertain data were investigated. Firstly, a robust optimization model corresponding to the uncertain multi-objective optimization problem was introduced. Then the optimization model was converted to a multi-objective optimization problem with the Dinkelbach method. In turn, by means of the scalarization method, the corresponding scalarization optimization problem was built, and the relationship between robust solutions to the multi-objective optimization problem and its corresponding scalarization optimization problem was described. Finally, through a robust-type sub-differential constraint qualification, the robust optimality condition for approximate quasi-efficient solutions to the multi-objective fractional optimization problem was established.

Based on the deterministic network infectious disease model, a stochastic network infectious disease model under the influence of white noise was established, and the existence and uniqueness of the global solution to the model were proved. With the theory of stochastic differential equations, sufficient conditions for stochastic extinction and persistence of infectious diseases were obtained. The results show that, white noise has a great impact on the transmission dynamics of network infectious diseases. White noise can effectively suppress the spread of infectious diseases, and large white noise can even make the original persistent infectious diseases become extinct. Finally, the theoretical results were verified through numerical simulations.

With the central difference scheme to discretize the Riemann-Liouville time fractional derivatives and by means of the finite point method to establish discrete algebraic equation systems, a meshless finite point method was proposed for the numerical analysis of the fractional Cable equation. The error estimation of the method was derived and discussed in detail. Numerical examples verify the efficiency and convergence of the method and confirm the theoretical results.