The coupled thermo-mechanical strain gradient plasticity theory that accounts for microstructure based size effects is outlined within this work. This incorporates spatial gradients of selected micro-structural fields...
详细信息
The coupled thermo-mechanical strain gradient plasticity theory that accounts for microstructure based size effects is outlined within this work. This incorporates spatial gradients of selected micro-structural fields based on length-scales that describe the evolving dissipative mechanisms. In the mechanical part, the model problem of von Mises plasticity with gradient-extended hardening/softening response is considered as discussed in Miehe et al. (2013, 2014a). In the thermal part, we follow the investigations of Simo and Miehe (1992) that demonstrate the effect of temperature on the mechanical fields resulting in a thermal expansion. To this end, two classes of solution schemes for the coupled problem are considered: (i) Global product formula algorithm arising from operator split which leads to a two step solution procedure, and (ii) an implicit coupled algorithm which employs simultaneous solution of the coupled system of equations. In the product formula algorithm, the mechanical and thermal problems are solved separately, resulting in a symmetric problem. However, in the implicit coupled algorithm, a simultaneous solution of the coupled system of equations for gradient thermo-plasticity is employed. A noteworthy drawback of this solution scheme arises from the high computational efforts in comparison with the product formula algorithm. From the computational viewpoint, the standard Galerkin finite element method fails in the context of isochoric plastic flow due to the over constrained pressure field. To circumvent these difficulties, we extend the well-known Q1P0-type and MINI-type mixed finite elements design of gradient plasticity to account for thermal effects. The performance of the formulation is demonstrated by means of some representative examples. (C) 2017 Elsevier Ltd. All rights reserved.
We propose numerically efficient solutions for the open- and closed-loop time-optimal soft-constrained control of a linear system representing a large flexible space structure. The open-loop solution is expressed in t...
详细信息
We propose numerically efficient solutions for the open- and closed-loop time-optimal soft-constrained control of a linear system representing a large flexible space structure. The open-loop solution is expressed in terms of the controllability Grammian matrix, for which we have obtained a closed-form expression for the undamped system. The qualitative dependence of the control on the initial state and the existence of many solutions satisfying the necessary conditions are shown. A nominal closed-loop control policy is subsequently formulated, but it is shown to be numerically expensive due to the nonuniqueness of extremal solutions. A continuation-based algorithm is proposed to alleviate the computational problem. Finally, the open- and closed-loop controls are shown to exhibit a saturation property reminiscent of the hard-constrained problem.
The implementation of the projected algorithm and of the consistent tangent tensor for general isotropic three-invariant elastoplastic models under plane stress conditions discussed in Part I of this paper [Valoroso, ...
详细信息
The implementation of the projected algorithm and of the consistent tangent tensor for general isotropic three-invariant elastoplastic models under plane stress conditions discussed in Part I of this paper [Valoroso, N., Rosati, L., 2008. Consistent derivation of the constitutive algorithm for plane stress isotropic plasticity. Part 1: Theoretical formulation. International journal of Solids and Structures, doi: 10.1016/***.2008.08.012.] is addressed. The connections between the general three-dimensional case and the plane stress problem are analyzed in detail and an algorithmic treatment taking full advantage of the isotropic properties of the model is presented. In particular, intrinsic (matrix-free) expressions are provided for all steps of the stress computation scheme that allow one to carry out the numerical implementation in a way that is completely independent from the matrix representations. The numerical performances of the present solution scheme are evaluated through representative numerical examples. (c) 2008 Elsevier Ltd. All rights reserved.
Compaction of crystallographic texture data is highly desirable in crystal plasticity simulations because the computational time involved in such calculations scales linearly with the number of crystal orientations. I...
详细信息
Compaction of crystallographic texture data is highly desirable in crystal plasticity simulations because the computational time involved in such calculations scales linearly with the number of crystal orientations. In a recent publication, we have reported a rigorous procedure for reducing large datasets of crystal orientations for cubic-orthotropic and hexagonal-orthotropic polycrystalline metals using symmetrized generalized spherical harmonics (GSH) functions. The procedure relies on a quantitative description of crystallographic texture using an orientation distribution function (ODF) and its series representation using GSH. The core procedure consists of matching the spectral representation of a fullsize ODF containing any number of crystal orientations with that of an ODF containing a compact set of orientations. In this paper, we generalize the procedure to any crystal structure with no restrictions to sample symmetry. These major extensions are accompanied by dealing with significantly more dimensions as well as imaginary terms. Two approaches for generating an initial set of orientations in the compact ODF are explored, one based on binning of a given fundamental zone in the Bunge-Euler orientation space and another that takes advantage of MTEX to maximize the compaction. The overall procedure has been successfully applied to compaction of large ODFs for cubic, hexagonal, and orthorhombic polycrystalline metals with orthotropic and no sample symmetry. It is quantitatively demonstrated that texture evolution, twin volume fraction evolution, stress-strain response, and geometrical changes of samples can be accurately simulated to large plastic strains with compact ODFs using crystal plasticity finite element models. (C) 2018 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.
We examine the effectiveness of absorbing layers as nonreflecting computational boundaries for the Euler equations. The absorbing-layer equations are simply obtained by splitting the governing equations in the coordin...
详细信息
We examine the effectiveness of absorbing layers as nonreflecting computational boundaries for the Euler equations. The absorbing-layer equations are simply obtained by splitting the governing equations in the coordinate directions and introducing absorption coefficients in each split equation. This methodology is similar to that used by Berenger for the numerical solutions of Maxwell's equations. Specifically, we apply this methodology to three physical problems-shock-vortex interactions, a plane free shear how, and an axisymmetric jet-with emphasis on acoustic wave propagation. Our numerical results indicate that the use of absorbing layers effectively minimizes numerical reflection in all three problems considered.
A simple dynamic analysis algorithm is presented in this paper for both forward and backward calculations of a pavement system consisting of an asphalt concrete layer, underlain by a uniform subgrade to a depth H wher...
详细信息
A simple dynamic analysis algorithm is presented in this paper for both forward and backward calculations of a pavement system consisting of an asphalt concrete layer, underlain by a uniform subgrade to a depth H wherein the bedrock is located. The subgrade soil is represented by a higher-order continuum model-the modified Vlasov model. The asphalt concrete layer is represented by a three-parameter complex compliance function in a frequency domain. The governing equations of the dynamic pavement system, along with the solution algorithms for both forward and backward computations, are presented in detail. A numerical example is provided to illustrate the importance of considering dynamic effect in predicting pavement response under dynamic load. In addition, numerical examples are given to demonstrate the use of nondestructive testing data to back calculate the material properties, such as the modulus, damping, creep compliance and fatigue cracking speed for an asphalt concrete layer and the modulus damping for the subgrade layer. (C) 1998 Elsevier Science Ltd. All rights reserved.
The impact of computational fluid dynamics (CFD) on the traditional role of aeronautical ground test facilities over the next 15 years is assessed via a National Research Council study. More powerful scientific comput...
详细信息
This is the story of exactly how the Pentium floating-point division problem was discovered, and what you need to know about the mathematics and computer engineering involved before deciding whether to replace the chi...
详细信息
This is the story of exactly how the Pentium floating-point division problem was discovered, and what you need to know about the mathematics and computer engineering involved before deciding whether to replace the chip, install the workaround provided here, or do nothing.
We construct filter bank transforms that are adapted to bounded domains. The transforms are constructed with the aim of solving PDEs together with finite difference methods, and the properties important for this kind ...
详细信息
We construct filter bank transforms that are adapted to bounded domains. The transforms are constructed with the aim of solving PDEs together with finite difference methods, and the properties important for this kind of application are analyzed. As the filter banks do not need to correspond to wavelets, short filters can be used. The price paid is extra growth factors in the number of coefficients and in the error, compared to the multiresolution analysis approach. However, these factors are small. We give upper bounds of the growth factors, and show numerical examples for an interval, a three-dimensional box, and a triangle. We also show examples of fast differentiation and other operations on thresholded expansions. (C) 2000 IMACS. Published by Elsevier Science B.V. All rights reserved.
A flow stress model describing precipitate hardening in the nickel based alloy Inconel (R) 718 following thermal treatment is presented. The interactions between precipitates and dislocations are included in a disloca...
详细信息
A flow stress model describing precipitate hardening in the nickel based alloy Inconel (R) 718 following thermal treatment is presented. The interactions between precipitates and dislocations are included in a dislocation density based material model. Compression tests have been performed using solution annealed, fully-aged and half-aged material. Models were calibrated using data for solution annealed and fully-aged material, and validated using data from half-aged material. Agreement between experimental data and model predictions is good. (C) 2013 Elsevier B.V. All rights reserved.
暂无评论