Based on the discrete element software PFC 2D , the biaxial test of clay is simulated. Contact bond model simulates the physical and mechanical behaviors of the clay, and is used to research the influences between the seven microscopic parameters, such as particle radius, particle radius ratio, stiffness ratio, friction coefficient, contact modulus, contact normal bond strength, contact shear bond strength and the macroscopic parameters of the clay, such as Young's modulus, Poisson's ratio, cohesion and internal friction angle, respectively. In order to make calibration process of clay efficiently, screen out the two most significant microscopic parameters for each response by design of experiment-the Plackett-Burman design method, and obtain the linear regression equation of macroscopic parameters. After that, estimate the nonlinear relationship between the two microscopic parameters and macroscopic parameters by the Response Surface Methods-Central Composite Design method, and obtain an optimum set of microscopic parameters. KEYWORDS: clay; microscopic parameter; the Plackett-Burman design method; Central Composite Design method
INTRODUCITONThe discrete element method (DEM) is used to simulate the mechanical behavior of rock by integrating the system ideal into a set of structural elements such as springs, beams, or separate particles bonded together through the contact point [1] . While the continuum model represents damages indirectly by the empirical relation, the discrete element method is used to represent the failure directly by the individual structural elements or the bond failure. DEM has been successfully applied to the modeling of soil behavior, especially the soil failure and nonlinear behavior. The particle flow code is mainly used for the study of particles or system analysis which can be simplified as particles, there are two kinds of two-dimensional (PFC
2D) and three-dimensional (PFC 3D ), which belongs to discrete element category.Different with the traditional finite element numerical software, discrete element numerical software PFC 2D can't assign value to macroscopic parameters directly (such as Young's modulus, Poisson's ratio, cohesion and internal friction angle). A series of experiments (such as uniaxial, biaxial test or Brazil splitting test, etc.) are carried out on the model which has explicit microscopic parameters, and these macroscopic parameters are obtained, this process is usually called the calibration process. So far the only way to find the