HR: 0800h
AN: NG41B-0432    [Abstracts]
TI: 3D Modeling of fracture of brittle rock using Lattice Solid Model with rotational mechanism
AU: * Wang, Y
EM: wangyc@quakes.uq.edu.au
AF: QUAKES, Earth Systems Science Computational Centre,, ESSCC,The University of Queensland, St. Lucia., Brisbane, Qld 4072 Australia
AU: * Wang, Y
EM: wangyc@quakes.uq.edu.au
AF: Australian Computational Earth Systems Simulator, ACcESS,The University of Queensland, St. Lucia., Brisbane, Qld 4072 Australia
AU: Abe, S
EM: steffen@quakes.uq.edu.au
AF: QUAKES, Earth Systems Science Computational Centre,, ESSCC,The University of Queensland, St. Lucia., Brisbane, Qld 4072 Australia
AU: Abe, S
EM: steffen@quakes.uq.edu.au
AF: Australian Computational Earth Systems Simulator, ACcESS,The University of Queensland, St. Lucia., Brisbane, Qld 4072 Australia
AU: Latham, S
EM: slatham@esscc.uq.edu.au
AF: QUAKES, Earth Systems Science Computational Centre,, ESSCC,The University of Queensland, St. Lucia., Brisbane, Qld 4072 Australia
AU: Latham, S
EM: slatham@esscc.uq.edu.au
AF: Australian Computational Earth Systems Simulator, ACcESS,The University of Queensland, St. Lucia., Brisbane, Qld 4072 Australia
AU: Mora, p
EM: morap@esscc.uq.edu.au
AF: QUAKES, Earth Systems Science Computational Centre,, ESSCC,The University of Queensland, St. Lucia., Brisbane, Qld 4072 Australia
AU: Mora, p
EM: morap@esscc.uq.edu.au
AF: Australian Computational Earth Systems Simulator, ACcESS,The University of Queensland, St. Lucia., Brisbane, Qld 4072 Australia
AB: The Lattice Solid Model (LSMearth) is a particle based model similar to the Discrete Element Model (DEM). The current LSMerath includes only radial interaction between two linked particles, involving only translational motions of particles. In this study, we extend the LSM by introducing full rigidity interactions between particles and full degrees of freedom for a single particle. In the new model, for each particle we introduce six degrees of freedom: 3 for translational motion, and 3 for orientation. Six kinds of relative motions are permitted, and six interactions are transferred, i.e., radial, two shearing forces, twisting and two bending torques. Particle motion is decomposed into translational motion of the center of mass and rotation about the center. The former is solved using conventional Molecular Dynamics algorithms. The latter is integrated using Fincham's leap-frog algorithm using quaternion representation of orientations. The relative rotation between two particles is decomposed into two sequence-independent rotations. Using such decomposition, all interactions due to the relative translational and rotational motions between interactive rigid bodies can be uniquely determined. We carried out several tests on 3-D rock failure under uni-axial compression and frictional instability between two blocks. Compared with the simulations without the single particle rotational mechanism, the new simulation results match more closely with experimental results
DE: 5104 Fracture and flow
DE: 3210 Modeling
SC: Nonlinear Geophysics [NG]
MN: 2004 AGU Fall Meeting