Granular sands are characterized and modeled here by explicitly exploiting the discretecontinuum duality of granular matter. Grain-scale kinematics, obtained by shearing a sample under triaxial compression, are coupled with a recently proposed multiscale computational framework to model the behavior of the material without resorting to phenomenological evolution (hardening) laws. By doing this,...

Abstract. Since the multi-scale finite element analysis (FEA) requires large computation time, development of the parallel computing technique for the multi-scale analysis is inevitable. A parallel elastic/crystalline viscoplastic FEA code based on a crystallographic homogenization method has been developed using PC cluster. The homogenization scheme is introduced to compute macro-continuum pla...

This paper presents the formulation of numerical algorithms for the solution of the closest-point projection equations that appear in typical implementations of return mapping algorithms in elastoplasticity. The main motivation behind this work is to avoid the poor global convergence properties of a straight application of a Newton scheme in the solution of these equations, the socalled Newton–...

The theory and algorithm for the Material Point Method (MPM) are documented, with a detailed discussion on the treatments of boundary conditions and shock wave problems. A step-by-step solution scheme is written based on direct inspection of the two-dimensional MPM code currently used at the University of Missouri-Columbia (which is, in turn, a legacy of the University of New Mexico code). To t...

Abstract This paper deals with the structural modelling of fibre networks a focus on description populations initially crimped fibres. It presents systematic approach introducing appropriate strain measures for single fibres based deformation decomposition and by transferring knowledge from field elastoplasticity. On this basis, example, often used Biot-type $\lambda-\lambda_{\mathrm{w}}$ <mml:...

A new ordinary state-based peridynamic model (OSB-PD) in 2D consistent with J2 plasticity using a novel decomposition for force and extension states is presented. strategy testing the consistency of an OSB-PD formulation elastoplasticity introduced. In contrast other similar models, elasto-plastic objective works large rotations. Two rate-independent yield functions equivalent to associated flo...

Abstract Fictitious domain methods, such as the finite cell method, simplify discretization of a significantly. This is because mesh does not need to conform interest. However, generation simplified, broken cells with discontinuous integrands must be integrated using special quadrature schemes. The moment fitting very efficient scheme for integrating since number integration points generated mu...