Contact detection between interacting blocks is of great importance to discontinuity-based numerical methods, such as DDA, DEM, and NMM. A rigorous contact theory is a prerequisite to describing the interactions of mu...Contact detection between interacting blocks is of great importance to discontinuity-based numerical methods, such as DDA, DEM, and NMM. A rigorous contact theory is a prerequisite to describing the interactions of multiple blocks. Currently, the penalty method, in which mathematical springs with high stiffness values are employed, is always used to calculate the contact forces. High stiffness values may cause numerical oscillations and limit the time step. Furthermore, their values are difficult to identify. The intention of this study is to present a two-scale contact model for the calculation of forces between colliding blocks. In this new model, a calculation step taken from the moment of contact will be divided into two time stages: the free motion time stage and the contact time stage. Actually, these two time stages correspond to two real physical processes. Based on this, we present a new numerical model that is intended to be more precise and useful in calculating the contact forces without mathematical springs. The propagation of the elastic wave during collision is of a characteristic length, which determines the volume of material involved in the contact force calculation. In conventional contact models, this range is always regarded as the length of one element, which may lead to an inaccurate calculation of contact forces. In fact, the real scale of this range is smaller than the length of a single element, and subdivided elements, which are refined according to the characteristic length and are presented in the new contact model.展开更多
Rock,concrete and other geo-materials,due to the presence of microstructural inhomogeneity,their fracture processes and damage characteristics are associated with the distribution of micro-cracks contained in the mate...Rock,concrete and other geo-materials,due to the presence of microstructural inhomogeneity,their fracture processes and damage characteristics are associated with the distribution of micro-cracks contained in the materials.In this study,by introducing a cohesive zone model based on fracture mechanics into the framework of deformable discrete element method,a continuous-discontinuous coupling analysis approach for simulating the fracture of quasi-brittle materials is proposed.The cohesive interface elements are inserted into certain engineering or research region.It is assumed that damage and fracture occur only in the interface elements,while bulk material is modeled to be elastic.The Mohr-Coulomb criterion with tension cut-off is adopted as the damage initiation criterion,and a scalar damage variable representing damage in the material is used to describe the rate at which the material stiffness is degraded.Cracks are simulated explicitly by the failure of the interface elements.Numerical simulations are performed in order to validate the suggested method.Partial applications are also listed.The results show that this method provides a simple but effective tool for the simulation of crack initiation and propagation,and it can reflect the whole process of quasi-brittle materials from small deformation to large deformation and failure.展开更多
In this paper,an approach of square coordinate transformation is proposed to approximate the spectral abscissa for continuous-time switched linear systems.By applying elementary transformations iteratively,a series of...In this paper,an approach of square coordinate transformation is proposed to approximate the spectral abscissa for continuous-time switched linear systems.By applying elementary transformations iteratively,a series of minimums of least μ1 matrix set measures are obtained,which are utilized to approximate the spectral abscissa of the switched system.The approach is developed into tractable numerical algorithms that provide upper bound estimates of the spectral abscissa.Numerical simulations show the effectiveness of the proposed method.展开更多
基金supported by the National Basic Research Program of China("973"Project)(Grant No.2015CB250903)the CAS Strategic Priority Research Program(B)(Grant No.XDB10030303)+1 种基金the National Key Technology Research and Development Program of the Ministry of Science and Technology of China(Grant No.2012BAK10B01)the Youth Science Fund of the National Natural Science Foundation of China(Grant No.11302230)
文摘Contact detection between interacting blocks is of great importance to discontinuity-based numerical methods, such as DDA, DEM, and NMM. A rigorous contact theory is a prerequisite to describing the interactions of multiple blocks. Currently, the penalty method, in which mathematical springs with high stiffness values are employed, is always used to calculate the contact forces. High stiffness values may cause numerical oscillations and limit the time step. Furthermore, their values are difficult to identify. The intention of this study is to present a two-scale contact model for the calculation of forces between colliding blocks. In this new model, a calculation step taken from the moment of contact will be divided into two time stages: the free motion time stage and the contact time stage. Actually, these two time stages correspond to two real physical processes. Based on this, we present a new numerical model that is intended to be more precise and useful in calculating the contact forces without mathematical springs. The propagation of the elastic wave during collision is of a characteristic length, which determines the volume of material involved in the contact force calculation. In conventional contact models, this range is always regarded as the length of one element, which may lead to an inaccurate calculation of contact forces. In fact, the real scale of this range is smaller than the length of a single element, and subdivided elements, which are refined according to the characteristic length and are presented in the new contact model.
基金supported by the National Basic Research Program of China(973 Program)(Grant No.2013CB035901)the National Natural Science Foundation of China(Grant No.51379161)+2 种基金the Program for New Century Excellent Talents in University(Grant No.NCET-10-0657)the Specialized Research Fund for the Doctoral Program of Higher Education of China(Grant No.20120141110008)the Fundamental Research Funds for the Central Universities(Grant No.2012206020207)
文摘Rock,concrete and other geo-materials,due to the presence of microstructural inhomogeneity,their fracture processes and damage characteristics are associated with the distribution of micro-cracks contained in the materials.In this study,by introducing a cohesive zone model based on fracture mechanics into the framework of deformable discrete element method,a continuous-discontinuous coupling analysis approach for simulating the fracture of quasi-brittle materials is proposed.The cohesive interface elements are inserted into certain engineering or research region.It is assumed that damage and fracture occur only in the interface elements,while bulk material is modeled to be elastic.The Mohr-Coulomb criterion with tension cut-off is adopted as the damage initiation criterion,and a scalar damage variable representing damage in the material is used to describe the rate at which the material stiffness is degraded.Cracks are simulated explicitly by the failure of the interface elements.Numerical simulations are performed in order to validate the suggested method.Partial applications are also listed.The results show that this method provides a simple but effective tool for the simulation of crack initiation and propagation,and it can reflect the whole process of quasi-brittle materials from small deformation to large deformation and failure.
基金supported by the National Key Basic Research Program(973 Plan)under Grant No.2014CB845302the National Natural Science Foundation of China under Grant No.61273121
文摘In this paper,an approach of square coordinate transformation is proposed to approximate the spectral abscissa for continuous-time switched linear systems.By applying elementary transformations iteratively,a series of minimums of least μ1 matrix set measures are obtained,which are utilized to approximate the spectral abscissa of the switched system.The approach is developed into tractable numerical algorithms that provide upper bound estimates of the spectral abscissa.Numerical simulations show the effectiveness of the proposed method.