Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction ...Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints is presented. In particular, a wheeled multi-body system is considered. Here, the state transition of stick-slip between wheel and ground is transformed into a nonlinear complementarity problem （NCP）. An iterative algorithm for solving the NCP is then presented using an event-driven method. Dynamical equations of the multi-body system with holonomic and nonholonomic constraints are given using Routh equations and a con- straint stabilization method. Finally, an example is used to test the proposed numerical method. The results show some dynamical behaviors of the wheeled multi-body system and its constraint stabilization effects.展开更多
Through detailed three-dimensional （3D） finite element （FE） calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are o...Through detailed three-dimensional （3D） finite element （FE） calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are obtained near the crack front with aspect ratios （a/c） of 0.2, 0.4, 0.5, 0.6, 0.8 and 1.0. Tz decreases from an approximate value of Poisson ratio v at the crack tip to zero with increasing normalized radial distances （r/a） in the normal plane of the crack front line, and increases gradually when the elliptical parameter angle φ changes from 0° to 90°at the same r/a. With a/c rising to 1.0, Tz is getting nearly independent of φ and is only related to r/a. Based on the present FE calculations for Tz, empirical formulas for Tz are obtained to describe the 3D distribution of Tz for embedded center-elliptical cracks using the least squares method in the range of 0.2 ≤ a/c ≤ 1.0. These Tz results together with the corresponding stress intensity factor K are well suitable for the analysis of the 3D embedded centerelliptical crack from field, and a two-parameter K-Tz principle is proposed.展开更多
This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zer...This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.展开更多
An explicit Bargmann symmetry constraint is computed and its associated binary nonlinearization of Lax pairs is carried out for the super NLS-MKdV hierarchy. Under the obtained symmetry constraint, the n-th flow of th...An explicit Bargmann symmetry constraint is computed and its associated binary nonlinearization of Lax pairs is carried out for the super NLS-MKdV hierarchy. Under the obtained symmetry constraint, the n-th flow of the super NLS-MKdV hierarchy is decomposed into two super finite-dimensional integrable Hamiltonian systems, defined over the super-symmetry manifold R4N|2N with the corresponding dynamical variables x and tn. The integrals of motion required for Liouville integrability are explicitly given.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.11372018 and 11572018)
文摘Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints is presented. In particular, a wheeled multi-body system is considered. Here, the state transition of stick-slip between wheel and ground is transformed into a nonlinear complementarity problem （NCP）. An iterative algorithm for solving the NCP is then presented using an event-driven method. Dynamical equations of the multi-body system with holonomic and nonholonomic constraints are given using Routh equations and a con- straint stabilization method. Finally, an example is used to test the proposed numerical method. The results show some dynamical behaviors of the wheeled multi-body system and its constraint stabilization effects.
基金The project supported by the National Natural Science Foundation of China （50275073）
文摘Through detailed three-dimensional （3D） finite element （FE） calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are obtained near the crack front with aspect ratios （a/c） of 0.2, 0.4, 0.5, 0.6, 0.8 and 1.0. Tz decreases from an approximate value of Poisson ratio v at the crack tip to zero with increasing normalized radial distances （r/a） in the normal plane of the crack front line, and increases gradually when the elliptical parameter angle φ changes from 0° to 90°at the same r/a. With a/c rising to 1.0, Tz is getting nearly independent of φ and is only related to r/a. Based on the present FE calculations for Tz, empirical formulas for Tz are obtained to describe the 3D distribution of Tz for embedded center-elliptical cracks using the least squares method in the range of 0.2 ≤ a/c ≤ 1.0. These Tz results together with the corresponding stress intensity factor K are well suitable for the analysis of the 3D embedded centerelliptical crack from field, and a two-parameter K-Tz principle is proposed.
基金supported by the Development of airborne gravity gradiometer(No.2017YFC0601601)open subject of Key Laboratory of Petroleum Resources Research,Institute of Geology and Geophysics,Chinese Academy of Sciences(No.KLOR2018-8)
文摘This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.
文摘An explicit Bargmann symmetry constraint is computed and its associated binary nonlinearization of Lax pairs is carried out for the super NLS-MKdV hierarchy. Under the obtained symmetry constraint, the n-th flow of the super NLS-MKdV hierarchy is decomposed into two super finite-dimensional integrable Hamiltonian systems, defined over the super-symmetry manifold R4N|2N with the corresponding dynamical variables x and tn. The integrals of motion required for Liouville integrability are explicitly given.
