Elastoviscoplastic constitutive frameworks for generalized continua

Summary.  A unifying thermomechanical constitutive framework for generalized continua including additional degrees of freedom or/and the second gradient of displacement is presented. Based on the analysis of the dissipation, state laws, flow rules and evolution equations are proposed for Cosserat, strain gradient and micromorphic continua. The case of the gradient of internal variable approach is also incorporated by regarding the nonlocal internal variable as an actual additional degree of freedom. The consistency of the continuum thermodynamical framework is ensured by the introduction of a viscoplastic pseudo–potential of dissipation, thus extending the classical class of so–called standard material models to generalized continua. Variants of the higher order and higher grade theories are also reported based on the explicit introduction of the plastic strain tensor as additional degree of freedom. Within this new class of models, called here gradient of strain models, one recognizes the fact that, in a second grade theory for instance, the plastic part of the strain gradient can be identified with the gradient of plastic strain. Simple examples dealing with bending and shearing of Cosserat or second grade media are given to illustrate two types of extensions of classical J ₂-plasticity : single-criterion and multi-mechanism generalized elastoplasticity. Finally, formulations at finite deformation of the proposed models are provided focusing on proper decompositions of Cosserat curvature, strain gradient and gradient of micromorphic deformation into elastic and viscoplastic parts. Received March 11, 2002 Published online: January 16, 2003     

Efficient generation of pareto‐optimal topologies for compliance optimization

In multi‐objective optimization, a design is defined to beit pareto‐optimal if no other design exists that is better with respect to one objective, and as good with respect to other objectives. In this paper, we first show that if a topology is pareto‐optimal, then it must satisfy certain properties associated with the topological sensitivity field, i.e. no further comparison is necessary. This, in turn, leads to a deterministic, i.e. non‐stochastic, method for efficiently generating pareto‐optimal topologies using the classic fixed‐point iteration scheme. The proposed method is illustrated, and compared against SIMP‐based methods, through numerical examples. In this paper, the proposed method of generating pareto‐optimal topologies is limited to bi‐objective optimization, namely compliance–volume and compliance–compliance. The future work will focus on extending the method to non‐compliance and higher dimensional pareto optimization. Copyright © 2011 John Wiley & Sons, Ltd.     

A gradient‐only line search method for the conjugate gradient method applied to constrained optimization problems with severe noise in the objective function

A new implementation of the conjugate gradient method is presented that economically overcomes the problem of severe numerical noise superimposed on an otherwise smooth underlying objective function of a constrained optimization problem. This is done by the use of a novel gradient‐only line search technique, which requires only two gradient vector evaluations per search direction and no explicit function evaluations. The use of this line search technique is not restricted to the conjugate gradient method but may be applied to any line search descent method. This method, in which the gradients may be computed by central finite differences with relatively large perturbations, allows for the effective smoothing out of any numerical noise present in the objective function. This new implementation of the conjugate gradient method, referred to as the ETOPC algorithm, is tested using a large number of well‐known test problems. For initial tests with no noise introduced in the objective functions, and with high accuracy requirements set, it is found that the proposed new conjugate gradient implementation is as robust and reliable as traditional first‐order penalty function methods. With the introduction of severe relative random noise in the objective function, the results are surprisingly good, with accuracies obtained that are more than sufficient compared to that required for engineering design optimization problems with similar noise. Copyright © 2004 John Wiley & Sons, Ltd.     

Stochastic inverse heat conduction using a spectral approach

An adjoint‐based functional optimization technique in conjunction with the spectral stochastic finite element method is proposed for the solution of an inverse heat conduction problem in the presence of uncertainties in material data, process conditions and measurement noise. The ill‐posed stochastic inverse problem is restated as a conditionally well‐posed L₂ optimization problem. The gradient of the objective function is obtained in a distributional sense by defining an appropriate stochastic adjoint field. The L₂ optimization problem is solved using a conjugate‐gradient approach. Accuracy and effectiveness of the proposed approach is appraised with the solution of several stochastic inverse heat conduction problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Phase‐field modeling of ductile fracture at finite strains: A robust variational‐based numerical implementation of a gradient‐extended theory by micromorphic regularization

This work provides a robust variational‐based numerical implementation of a phase field model of ductile fracture in elastic–plastic solids undergoing large strains. This covers a computationally efficient micromorphic regularization of the coupled gradient plasticity‐damage formulation. The phase field approach regularizes sharp crack surfaces within a pure continuum setting by a specific gradient damage modeling with geometric features rooted in fracture mechanics. It has proven immensely successful with regard to the analysis of complex crack topologies without the need for fracture‐specific computational structures such as finite element design of crack discontinuities or intricate crack‐tracking algorithms. The proposed gradient‐extended plasticity‐damage formulation includes two independent length scales that regularize both the plastic response as well as the crack discontinuities. This ensures that the damage zones of ductile fracture are inside of plastic zones or vice versa and guarantees on the computational side a mesh objectivity in post‐critical ranges. The proposed setting is rooted in a canonical variational principle. The coupling of gradient plasticity to gradient damage is realized by a constitutive work density function that includes the stored elastic energy and the dissipated work due to plasticity and fracture. The latter represents a coupled resistance to plasticity and damage, depending on the gradient‐extended internal variables that enter plastic yield functions and fracture threshold functions. With this viewpoint on the generalized internal variables at hand, the thermodynamic formulation is outlined for gradient‐extended dissipative solids with generalized internal variables that are passive in nature. It is specified for a conceptual model of von Mises‐type elasto‐plasticity at finite strains coupled with fracture. The canonical theory proposed is shown to be governed by a rate‐type minimization principle, which fully determines the coupled multi‐field evolution problem. This is exploited on the numerical side by a fully symmetric monolithic finite element implementation. An important aspect of this work is the regularization towards a micromorphic gradient plasticity‐damage setting by taking into account additional internal variable fields linked to the original ones by penalty terms. This enhances the robustness of the finite element implementation, in particular, on the side of gradient plasticity. The performance of the formulation is demonstrated by means of some representative examples. Copyright © 2016 John Wiley & Sons, Ltd.     

Dynamic response of intrinsic continua for use in biological and molecular modeling: Explicit finite element formulation

An intrinsic beam formulation has recently appeared (AIAA J. 2003; 41 (6):1131–1137), which presents the three‐dimensional equations of motion governing spatial and temporal changes in a one‐dimensional continua's curvature, strain, rectilinear velocity, and angular velocity. The formulation would suggest several computational advantages over more‐traditional displacement‐based continua approaches: low‐order interpolation functions can describe generally curved and twisted continua configurations; inter‐element displacements, slopes, strains, and curvatures can be matched; and finite rotational variables and their complexities are absent. Here, we present a completed intrinsic continua finite element development and critical analysis, which follows from an earlier preliminary treatment as applied to carbon nanotubes (Int. J. Solids Struct. 2007; 44 :874–894). Modeling of nodal displacements and rotations are included. Explicit time stepping, with desired high‐frequency damping, is accomplished using an implementation of the generalized‐α method. Zero‐energy modes inherent in the formulation are also identified and rectified. Finally, we document very good agreement between results predicted with the intrinsic continua finite element simulator and results generated using more‐traditional finite element simulations. Copyright © 2009 John Wiley & Sons, Ltd.     

A robust dual algorithm for topology design of structures in discrete variables

Dual optimization algorithms for the topology optimization of continuum structures in discrete variables are gaining popularity in recent times since, in topology design problems, the number of constraints is small in comparison to the number of design variables. Good topologies can be obtained for the minimum compliance design problem when the perimeter constraint is imposed in addition to the volume constraint. However, when the perimeter constraint is relaxed, the dual algorithm tends to give bad results, even with the use of higher‐order finite element models as we demonstrate in this work. Since, a priori, one does not know what a good value of the perimeter to be specified is, it is essential to have an algorithm which generates good topologies even in the absence of the perimeter constraint. We show how the dual algorithm can be made more robust so that it yields good designs consistently in the absence of the perimeter constraint. In particular, we show that the problem of checkerboarding which is frequently observed with the use of lower‐order finite elements is eliminated. Copyright © 2001 John Wiley & Sons, Ltd.     

A finite element displacement formulation for gradient elastoplasticity

We present a second gradient elastoplastic model for strain‐softening materials based entirely on a finite element displacement formulation. The stress increment is related to both the strain increment and its Laplacian. The displacement field is the only field needed to be discretized using a C¹ continuity element. The required higher‐order boundary conditions arise naturally from the displacement field. The model is developed to regularize the ill‐posedness caused by strain‐softening material behaviour. The gradient terms in the constitutive equations introduce an extra material parameter with dimensions of length allowing robust modelling of the post‐peak material behaviour leading to localization of deformation. Mesh insensitivity is demonstrated by modelling localization of deformation in biaxial tests. It is shown that both the thickness and inclination of the shear‐band zone are insensitive to the mesh directionality and refinement and agree with the expected theoretical and experimental values. Copyright © 2001 John Wiley & Sons, Ltd.     

Task assignment problems in distributed computing systems by simulated annealing

Abstract

The stochastic, heuristic search algorithm called simulated annealing is considered for the problems of static task assignment in distributed computing systems. The purposes of task assignment problems are to assign modules of programs over a set of interconnected processors in order to both maximize the utilization of processors and minimize interprocessor communication costs. This problem has been proven to be NP‐hard. Although simulated annealing has been applied to a broad class of combinatorial optimization problems, but it requires a long computation time in order to converge to the globally optimal solution. In this paper, we design a very efficient annealing schedule with good move generation strategies and use the concept of specific heat and the frozen condition to obtain near‐optimal solutions for task assignment problems with a significantly large reduction in the number of iterations.     

An X‐FEM approach for large sliding contact along discontinuities

The extended finite element method (X‐FEM) has been developed to minimize requirements on the mesh design in a problem with a displacement discontinuity. This advantage, however, still remains limited to the small deformation hypothesis when considering sliding discontinuities. The approach presented in this paper proposes to couple X‐FEM with a Lagrangian large sliding frictionless contact algorithm. A new hybrid X‐FEM contact element was developed with a contact search algorithm allowing for an update of contacting surfaces pairing. The stability of the contact formulation is ensured by an algorithm for fulfilling Ladyzhenskaya‐Babuska‐Brezzi (LBB) condition. Several 2D simple examples are presented in this paper in order to prove its efficiency and stability. Copyright © 2008 John Wiley & Sons, Ltd.     

Stochastic optimization using a sparse grid collocation scheme

In computational sciences, optimization problems are frequently encountered in solving inverse problems for computing system parameters based on data measurements at specific sensor locations, or to perform design of system parameters. This task becomes increasingly complicated in the presence of uncertainties in boundary conditions or material properties. The task of computing the optimal probability density function (PDF) of parameters based on measurements of physical fields of interest in the form of a PDF, is posed as a stochastic optimization problem. This stochastic optimization problem is solved by dividing it into two problems—an auxiliary optimization problem to construct stochastic space representations from the PDF of measurement data, and a stochastic optimization problem to compute the PDF of problem parameters. The auxiliary optimization problem is solved using a downhill simplex method, whilst a gradient based approach is employed for solving the stochastic optimization problem. The gradients required for stochastic optimization are defined, using appropriate stochastic sensitivity problems. A computationally efficient sparse grid collocation scheme is utilized to compute the solution of these stochastic sensitivity problems. The implementation discussed, requires minimum intrusion into existing deterministic solvers, and it is thus applicable to a variety of problems. Numerical examples involving stochastic inverse heat conduction problems, contamination source identification problems and large deformation robust design problems are discussed.     

Topology optimization of continuum structures with stress constraints and uncertainties in loading

Topology optimization using stress constraints and considering uncertainties is a serious challenge, since a reliability problem has to be solved for each stress constraint, for each element in the mesh. In this paper, an alternative way of solving this problem is used, where uncertainty quantification is performed through the first‐order perturbation approach, with proper validation by Monte Carlo simulation. Uncertainties are considered in the loading magnitude and direction. The minimum volume problem subjected to local stress constraints is formulated as a robust problem, where the stress constraints are written as a weighted average between their expected value and standard deviation. The augmented Lagrangian method is used for handling the large set of local stress constraints, whereas a gradient‐based algorithm is used for handling the bounding constraints. It is shown that even in the presence of small uncertainties in loading direction, different topologies are obtained when compared to a deterministic approach. The effect of correlation between uncertainties in loading magnitude and direction on optimal topologies is also studied, where the main observed result is loss of symmetry in optimal topologies.     

Active control of FGM plates subjected to a temperature gradient: Modelling via finite element method based on FSDT

An efficient finite element formulation based on a first‐order shear deformation theory (FSDT) is presented for the active control of functionally gradient material (FGM) plates with integrated piezoelectric sensor/actuator layers subjected to a thermal gradient; this is accomplished using both static and dynamic piezothermoelastic analyses. The formulation based on FSDT can be applied to a range of relatively thin‐to‐moderately thick plates. A constant displacement‐cum‐velocity feedback control algorithm coupling the direct and inverse piezoelectric effects is applied to provide active feedback control of the integrated FGM plate in a self‐monitoring and self‐controlling system. Numerical results for the control of bending and torsional deflections and/or vibrations are presented for a FGM plate comprising zirconia and aluminium. The effects of constituent volume fraction and the influence of feedback control gain on the static and dynamic responses of the FGM plates are examined in detail. Copyright © 2001 John Wiley & Sons, Ltd.     

A stabilized finite element formulation for monolithic thermo‐hydro‐mechanical simulations at finite strain

An adaptively stabilized monolithic finite element model is proposed to simulate the fully coupled thermo‐hydro‐mechanical behavior of porous media undergoing large deformation. We first formulate a finite‐deformation thermo‐hydro‐mechanics field theory for non‐isothermal porous media. Projection‐based stabilization procedure is derived to eliminate spurious pore pressure and temperature modes due to the lack of the two‐fold inf‐sup condition of the equal‐order finite element. To avoid volumetric locking due to the incompressibility of solid skeleton, we introduce a modified assumed deformation gradient in the formulation for non‐isothermal porous solids. Finally, numerical examples are given to demonstrate the versatility and efficiency of this thermo‐hydro‐mechanical model. Copyright © 2015 John Wiley & Sons, Ltd.     

A finite element method based on C0‐continuous assumed gradients

This article presents an alternative approach to assumed gradient methods in FEM applied to three‐dimensional elasticity. Starting from nodal integration (NI), a general C0‐continuous assumed interpolation of the deformation gradient is formulated. The assumed gradient is incorporated using the principle of Hu‐Washizu. By dual Lagrange multiplier spaces, the functional is reduced to the displacements as the only unknowns. An integration scheme is proposed where the integration points coincide with the support points of the interpolation. Requirements for regular finite element meshes are explained. Using this interpretation of NI, instabilities (appearance of spurious modes) can be explained. The article discusses and classifies available strategies to stabilize NI such as penalty methods, SCNI, α‐FEM. Related approaches, such as the smoothed finite element method, are presented and discussed. New stabilization techniques for NI are presented being entirely based on the choice of the assumed gradient interpolation, i.e. nodal‐bubble support, edge‐based support and support using tensor‐product interpolations. A strategy is presented on how the interpolation functions can be derived for various element types. Interpolation functions for the first‐order hexahedral element, the first‐order and the second‐order tetrahedral elements are given. Numerous examples illustrate the strengths and limitations of the new schemes. Copyright © 2010 John Wiley & Sons, Ltd.     

An enhanced genetic algorithm for structural topology optimization

Genetic algorithms (GAs) have become a popular optimization tool for many areas of research and topology optimization an effective design tool for obtaining efficient and lighter structures. In this paper, a versatile, robust and enhanced GA is proposed for structural topology optimization by using problem‐specific knowledge. The original discrete black‐and‐white (0–1) problem is directly solved by using a bit‐array representation method. To address the related pronounced connectivity issue effectively, the four‐neighbourhood connectivity is used to suppress the occurrence of checkerboard patterns. A simpler version of the perimeter control approach is developed to obtain a well‐posed problem and the total number of hinges of each individual is explicitly penalized to achieve a hinge‐free design. To handle the problem of representation degeneracy effectively, a recessive gene technique is applied to viable topologies while unusable topologies are penalized in a hierarchical manner. An efficient FEM‐based function evaluation method is developed to reduce the computational cost. A dynamic penalty method is presented for the GA to convert the constrained optimization problem into an unconstrained problem without the possible degeneracy. With all these enhancements and appropriate choice of the GA operators, the present GA can achieve significant improvements in evolving into near‐optimum solutions and viable topologies with checkerboard free, mesh independent and hinge‐free characteristics. Numerical results show that the present GA can be more efficient and robust than the conventional GAs in solving the structural topology optimization problems of minimum compliance design, minimum weight design and optimal compliant mechanisms design. It is suggested that the present enhanced GA using problem‐specific knowledge can be a powerful global search tool for structural topology optimization. Copyright © 2005 John Wiley & Sons, Ltd.     

Robust topology optimization of phononic crystals with random field uncertainty

The uncertain spatial variation of material properties can remarkably affect the band gap characteristics of phononic crystals (PnCs). It is necessary to consider this issue when designing and manufacturing PnC materials/structures. This paper investigates a robust topology optimization method for designing the microstructures of PnCs by considering random‐field material properties. Herein, the spatial distribution of the material properties is first represented by a random field and then discretized into uncorrelated stochastic variables with the expansion optimal linear estimation method; stochastic band gap analysis is then conducted with polynomial chaos expansion. Furthermore, a robust topology optimization formulation of PnCs is proposed on the basis of the relative elemental density, where a weighted objective function handles the compromise of the mean value and standard deviation of the PnC band gap. The band gap response is analyzed, employing the finite element method for each sample of polynomial chaos expansion. In this context, the sensitivities of the stochastic band gap behaviors to the design variables are also derived. Numerical examples demonstrate that the proposed method can generate meaningful optimal topologies of PnCs with a relatively large width and less sensitive band gap. Additionally, the effects of the weight factors in the objective function and the variation coefficient of material properties are discussed.     

A first‐order hyperbolic framework for large strain computational solid dynamics: An upwind cell centred Total Lagrangian scheme

This paper builds on recent work developed by the authors for the numerical analysis of large strain solid dynamics, by introducing an upwind cell centred hexahedral finite volume framework implemented within the open source code OpenFOAM [ http://www.openfoam.com/ ]. In Lee, Gil and Bonet (2013), a first‐order hyperbolic system of conservation laws was introduced in terms of the linear momentum and the deformation gradient tensor of the system, leading to excellent behaviour in two‐dimensional bending dominated nearly incompressible scenarios. The main aim of this paper is the extension of this algorithm into three dimensions, its tailor‐made implementation into OpenFOAM and the enhancement of the formulation with three key novelties. First, the introduction of two different strategies in order to ensure the satisfaction of the underlying involutions of the system, that is, that the deformation gradient tensor must be curl‐free throughout the deformation process. Second, the use of a discrete angular momentum projection algorithm and a monolithic Total Variation Diminishing Runge–Kutta time integrator combined in order to guarantee the conservation of angular momentum. Third, and for comparison purposes, an adapted Total Lagrangian version of the hyperelastic‐GLACE nodal scheme of Kluth and Després (2010) is presented. A series of challenging numerical examples are examined in order to assess the robustness and accuracy of the proposed algorithm, benchmarking it against an ample spectrum of alternative numerical strategies developed by the authors in recent publications. Copyright © 2016 John Wiley & Sons, Ltd.