Computer Science

The local and overall responses of nonlinear composites are classically investigated by the Finite Element Method. We propose an alternate method based on Fourier series which avoids meshing and which makes direct use of microstructure images. It is based on the exact expression of the Green function of a linear elastic and homogeneous comparison material. First, the case of elastic nonhomogeneous constituents is considered and an iterative procedure is proposed to solve the Lippman-Schwinger equation which naturally arises in the problem. Then, the method is extended to non-linear constituents by a step-by-step integration in time. The accuracy of the method is assessed by varying the spatial resolution of the microstructures. The flexibility of the method allows it to serve for a large variety of microstructures. (C) 1998 Elsevier Science S.A.

While crack nucleation and propagation in the brittle or quasi-brittle regime can be predicted via variational or material-force-based phase field fracture models, these models often assume that the underlying elastic response of the material is non-polar and yet a length scale parameter must be introduced to enable the sharp cracks represented by a regularized implicit function. However, many materials with internal microstructures that contain surface tension, micro-cracks, micro-fracture, inclusion, cavity or those of particulate nature often exhibit size-dependent behaviors in both the path-independent and path-dependent regimes. This paper is intended to introduce a unified treatment that captures the size effect of the materials in both elastic and damaged states. By introducing a cohesive micropolar phase field fracture theory, along with the computational model and validation exercises, we explore the interacting size-dependent elastic deformation and fracture mechanisms exhibits in materials of complex microstructures. To achieve this goal, we introduce the distinctive degradation functions of the force-stress-strain and couple-stress-micro-rotation energy-conjugated pairs for a given regularization profile such that the macroscopic size-dependent responses of the micropolar continua is insensitive to the length scale parameter of the regularized interface. Then, we apply the variational principle to derive governing equations from the micropolar stored energy and dissipative functionals. Numerical examples are introduced to demonstrate the proper way to identify material parameters and the capacity of the new formulation to simulate complex crack patterns in the quasi-static regime.

In this paper, a new computational framework based on the topology derivative concept is presented for evaluating stochastic topological sensitivities of complex systems. The proposed framework, designed for dealing with high dimensional random inputs, dovetails a polynomial dimensional decomposition (PDD) of multivariate stochastic response functions and deterministic topology derivatives. On one hand, it provides analytical expressions to calculate topology sensitivities of the first three stochastic moments which are often required in robust topology optimization (RTO). On another hand, it offers embedded Monte Carlo Simulation (MCS) and finite difference formulations to estimate topology sensitivities of failure probability for reliability-based topology optimization (RBTO). For both cases, the quantification of uncertainties and their topology sensitivities are determined concurrently from a single stochastic analysis. Moreover, an original example of two random variables is developed for the first time to obtain analytical solutions for topology sensitivity of moments and failure probability. Another 53-dimension example is constructed for analytical solutions of topology sensitivity of moments and semi-analytical solutions of topology sensitivity of failure probabilities in order to verify the accuracy and efficiency of the proposed method for high-dimensional scenarios. Those examples are new and make it possible for researchers to benchmark stochastic topology sensitivities of existing or new algorithms. In addition, it is unveiled that under certain conditions the proposed method achieves better accuracies for stochastic topology sensitivities than for the stochastic quantities themselves.

Based on previous work for the static problem, in this paper we first derive one form of dynamic finite-strain shell equations for incompressible hyperelastic materials that involve three shell constitutive relations. In order to single out the bending effect as well as to reduce the number of shell constitutive relations, a further refinement is performed, which leads to a refined dynamic finite-strain shell theory with only two shell constitutive relations (deducible from the given three-dimensional (3D) strain energy function) and some new insights are also deduced. By using the weak formulation of the shell equations and the variation of the 3D Lagrange functional, boundary conditions and the two-dimensional (2D) shell virtual work principle are derived. As a benchmark problem, we consider the extension and inflation of an arterial segment. The good agreement between the asymptotic solution based on the shell equations and that from the 3D exact one gives verification of the former. The refined shell theory is also applied to study the plane-strain vibrations of a pressurized artery, and the effects of the axial pre-stretch, pressure and fibre angle on the vibration frequencies are investigated in detail.

We describe a computational framework for simulating suspensions of rigid particles in Newtonian Stokes flow. One central building block is a collision-resolution algorithm that overcomes the numerical constraints arising from particle collisions. This algorithm extends the well-known complementarity method for non-smooth multi-body dynamics to resolve collisions in dense rigid body suspensions. This approach formulates the collision resolution problem as a linear complementarity problem with geometric `non-overlapping' constraints imposed at each timestep. It is then reformulated as a constrained quadratic programming problem and the Barzilai-Borwein projected gradient descent method is applied for its solution. This framework is designed to be applicable for any convex particle shape, e.g., spheres and spherocylinders, and applicable to any Stokes mobility solver, including the Rotne-Prager-Yamakawa approximation, Stokesian Dynamics, and PDE solvers (e.g., boundary integral and immersed boundary methods). In particular, this method imposes Newton's Third Law and records the entire contact network. Further, we describe a fast, parallel, and spectrally-accurate boundary integral method tailored for spherical particles, capable of resolving lubrication effects. We show weak and strong parallel scalings up to 8× 10 4 particles with approximately 4× 10 7 degrees of freedom on 1792 cores. We demonstrate the versatility of this framework with several examples, including sedimentation of particle clusters, and active matter systems composed of ensembles of particles driven to rotate.

In this paper, based on the idea of self-adjusting steepness based schemes[5], a two-dimensional calculation method of steepness parameter is proposed, and thus a two-dimensional self-adjusting steepness based limiter is constructed. With the application of such limiter to the over-intersection based remapping framework, a low dissipation remapping method has been proposed that can be applied to the existing ALE method.

In this paper, we present a simple artificial damping method to enhance the robustness of total Lagrangian smoothed particle hydrodynamics (TL-SPH). Specifically, an artificial damping stress based on the Kelvin-Voigt type damper with a scaling factor imitating a von Neumann-Richtmyer type artificial viscosity is introduced in the constitutive equation to alleviate the spurious oscillation in the vicinity of the sharp spatial gradients. After validating the robustness and accuracy of the present method with a set of benchmark tests with very challenging cases, we demonstrate its potentials in the field of bio-mechanics by simulating the deformation of complex stent structures.

This paper presents a simplified implementation of the arc-length method for computing the equilibrium paths of nonlinear structural mechanics problems using the finite element method. In the proposed technique, the predictor is computed by extrapolating the solutions from two previously converged load steps. The extrapolation is a linear combination of the previous solutions; therefore, it is simple and inexpensive. Additionally, the proposed extrapolated predictor also serves as a means for identifying the forward movement along the equilibrium path without the need for any sophisticated techniques commonly employed for explicit tracking. The ability of the proposed technique to successfully compute complex equilibrium paths in static structural mechanics problems is demonstrated using seven numerical examples involving truss, beam-column and shell models. The computed numerical results are in excellent agreement with the reference solutions. The present approach does not require prohibitively small increments for its success.

We demonstrate the ability of a stabilized finite element method, inspired by the weighted Nitsche approach, to alleviate spurious traction oscillations at interlaminar interfaces in multi-ply multi-directional composite laminates. In contrast with the standard (penalty-like) method, the stabilized method allows the use of arbitrarily large values of cohesive stiffness and obviates the need for engineering approaches to estimate minimum cohesive stiffness necessary for accurate delamination analysis. This is achieved by defining a weighted interface traction in the stabilized method, which allows a gradual transition from penalty-like method for soft elastic contact to Nitsche-like method for rigid contact. We conducted several simulation studies involving constant strain patch tests and benchmark delamination tests under mode-I, mode-II and mixed-mode loadings. Our results show clear evidence of traction oscillations with the standard method with structured and perturbed finite element meshes, and that the stabilized method alleviates these oscillations, thus illustrating its robustness.

Tensile instability, often observed in smoothed particle hydrodynamics (SPH), is a numerical artifact that manifests itself by unphysical clustering or separation of particles. The instability originates in estimating the derivatives of the smoothing functions which, when interact with material constitution may result in negative stiffness in the discretized system. In the present study, a stable formulation of SPH is developed where the kernel function is continuously adapted at every material point depending on its state of stress. Bspline basis function with a variable intermediate knot is used as the kernel function. The shape of the kernel function is then modified by changing the intermediate knot position such that the condition associated with instability does not arise. While implementing the algorithm the simplicity and computational efficiency of SPH are not compromised. One-dimensional dispersion analysis is performed to understand the effect adaptive kernel on the stability. Finally, the efficacy of the algorithm is demonstrated through some benchmark elastic dynamics problems.