A Boundary Element Method Code in 2-D using Matlab for Fracture and Fault Modeling
Personnel: David D. Pollard, Gaurav Chopra, Kurt Sternlof
Sponsor: Stanford Rock Fracture Project
Description: BEM2D is a computer program written by Gaurav Chopra in MATLAB® which models the quasi-static linear elastic deformation of an infinite medium in two dimensions. The code is designed to model the behavior of fractures and faults in nature. Each fracture/fault is referred to as a "structure" in the code and these structures are composed of many different linear "elements". The normal and shear stiffnesses of elements or the displacement discontinuities across elements may be specified as boundary conditions (inputs given by the user). The medium surrounding the elements is postulated to be homogenous, isotropic and linear with respect to elastic properties while the elements can possess non-linear stress displacement relations. The user may define any geometry of the structures (resolved as straight segments) and there may be any number of structures which contain any number of elements. The external loading is specified as a uniform stress state. The structures are allowed to propagate according to a criterion, provided by the user. The stress state at any point in the elastic whole space results from the superposed contributions of every element, as well as the remotely applied stress state. The simulation results from the code may be compared to field observations and changes in the propagation criterion or changes in the element stiffness made to find better agreement.