Skip to main content
eScholarship
Open Access Publications from the University of California

Scholarly Works (5 results)

  • Thesis
  • Peer Reviewed
UC Santa Barbara Electronic Theses and Dissertations (2017)

Computer simulation of fracture in materials with nonlinear mechanical response can be computationally expensive. These simulations often require a large number of degrees of freedom, and the nonlinearity in the problem can pose difficulties when computing solutions. This work focuses on two material models. The first model consists of rigid bricks interacting through nonlinear cohesive springs. Fracture in the material occurs through the rupture of the cohesive springs. The second, more complicated, model consists of deformable elements interacting through nonlinear cohesive springs.

In the first model, we assume the bricks are under a quasi-static loading scenario. With this assumption, the problem can be solved using a global Monte Carlo minimization algorithm to minimize the energy of the system. The energy in the system comes from the deformation and rupture of the nonlinear cohesive springs. Since these simulations have a high computational cost, we have developed a GPU-based (Graphics Processing Unit) Monte Carlo minimization algorithm that offers a significant speedup compared to a conventional multithreaded CPU-based algorithm.

With the second model, we have dynamic simulations with explicit time discretization. In this case we compute the force, acceleration, velocity, and position explicitly. The force in the system comes from both the deformation of the elements as well as the deformation of the nonlinear cohesive springs. We have developed explicit, CPU-based methods and implicit-explict methods on both CPUs and GPUs. Our implicit-explict GPU-based method achieves substantial performance improvement compared to the explicit, CPU-based method.

We present our GPU-based implementation of AES (Advanced Encryption Standard), which is used in the Monte Carlo minimization algorithm to generate random numbers. Our implementation is substantially faster than CPU-based implementation of AES. It is also faster than previous GPU implementations of AES.

    Cover page: Algorithms and Software for High-Performance Fracture Simulation on GPU Architectures
    • Article
    • Peer Reviewed
    UC Santa Barbara Previously Published Works (2015)
    Stiff ceramic platelets (or bricks) that are aligned and bonded to a second ductile phase with low volume fraction (mortar) are a promising pathway to produce stiff, high-toughness composites. For certain ranges of constituent properties, including those of some synthetic analogs to nacre, one can demonstrate that the deformation is dominated by relative brick motions. This paper describes simulations of fracture that explicitly track the motions of individual rigid bricks in an idealized microstructure; cohesive tractions acting between the bricks introduce elastic, plastic and rupture behaviors. Results are presented for the stresses and damage near macroscopic cracks with different brick orientations relative to the loading orientation. The anisotropic macroscopic initiation toughness is computed for small-scale yielding conditions and is shown to be independent of specimen geometry and loading configuration. The results are shown to be in agreement with previously published experiments on synthetic nacre.
      • Article
      • Peer Reviewed
      UC Santa Barbara Previously Published Works (2018)
        • Article
        • Peer Reviewed
        UC Santa Barbara Previously Published Works (2016)
        Simulation of the nonlinear mechanical response of materials with explicit representation of microstructural features is extremely challenging. These models typically involve a very large number of degrees of freedom, and are prone to convergence difficulties when searching for roots to nonlinear equilibrium equations. We focus on an idealized material model that is motivated by the microstructure of synthetic nacre: individual ‘bricks’ (representing ceramic platelets) interact through nonlinear cohesive springs (representing a small volume fraction of polymer that bonds the platelets). The model simulates composite fracture through rupture of the cohesive springs. The problem is cast in terms of energy minimization and is essentially described by ‘nearest neighbor’ interactions. The principal focus of this paper is to illustrate the computational gains achievable by the strategic marriage of robust, global Monte Carlo minimization algorithms to the graphics processing unit architecture, and to describe how they were realized on the Nvidia GPU. Results comparing the computation times for graphics processing unit and central processing unit implementations demonstrate that a new adaptive version of the simulated annealing algorithm yields a speedup of approximately 5 times, whereas the graphics processing unit implementation yields a speed-up of about 16 times over conventional four-core central processing unit implementations. The resulting speed enhancement for adaptive graphics processing unit minimization of a factor of 80 enables a far broader range of simulations than has previously been possible. Simulations involving as many as 300,000 bricks can be performed in hours, as compared to weeks required by central processing unit implementation. Many aspects of this approach are translatable to other physical problems involving energy minimization in systems with large numbers of degrees of freedom.
          • Article
          • Peer Reviewed
          UC Santa Barbara Previously Published Works (2015)
          This paper examines the effect of non-uniform microstructures on the macroscopic fracture properties of idealized brick and mortar composites, which consist of rigid bricks bonded with elastic-plastic mortar that ruptures at finite strain. A simulation tool that harnesses the parallel processing power of graphics processing units (GPUs) was used to simulate fracture in virtual specimens, whose microstructures were generated by sampling a probability distribution of brick sizes. In the simulations, crack advance is a natural outcome of local ruptures in the cohesive zones bonding the bricks: the macroscopic initiation toughness for small-scale yielding is inferred by correlating the critical load needed to advance a pre-defined crack with an associated far-field energy release rate. Quantitative connections between the statistical parameters defining heterogeneous brick distributions and the statistics of initiation toughness are presented. The nature of crack tip damage and stresses ahead of the crack tip are illustrated as a function of brick size variability. The results offer quantitative insights that can be used to identify microstructural targets for process development, notably specific brick size distributions that still provide macroscopic toughening.
            Top