Extension of Generalized Plasticity Model for Thermocyclic Loading

Author(s):  
J. C. Sobotka ◽  
R. H. Dodds

This work extends the generalized plasticity model for structural metals under cyclic loading proposed by Lubliner et al. (1993, “A New Model of Generalized Plasticity and its Numerical Implementation,” Int. J. Solids Struct., 22, pp. 3171–3184) to incorporate temperature-dependence into the elastic-plastic response. Proposed flow equations satisfy the Clausius–Duhem inequality through a thermodynamically consistent energy functional and retain key aspects of conventional plasticity models: Mises yield surface, normal plastic flow, and additive decomposition of strain. Uniaxial specialization of the 3D rate equations leads to a simple graphical method to estimate model properties. The 3D integration scheme based on backward Euler discretization leads to a scalar quadratic expression to determine the plastic strain rate multiplier and has a symmetric algorithmic tangent matrix. Both properties of the integration lead to a computationally efficient implementation especially suited to large-scale, finite element analyses. In comparison studies using experimental data from a Cottrell–Stokes test, the modified rate equations for the generalized plasticity model capture a thermally activated increase in the flow stress.

2020 ◽  
Vol 48 (4) ◽  
pp. 987-1003
Author(s):  
Hans Georg Bock ◽  
Jürgen Gutekunst ◽  
Andreas Potschka ◽  
María Elena Suaréz Garcés

AbstractJust as the damped Newton method for the numerical solution of nonlinear algebraic problems can be interpreted as a forward Euler timestepping on the Newton flow equations, the damped Gauß–Newton method for nonlinear least squares problems is equivalent to forward Euler timestepping on the corresponding Gauß–Newton flow equations. We highlight the advantages of the Gauß–Newton flow and the Gauß–Newton method from a statistical and a numerical perspective in comparison with the Newton method, steepest descent, and the Levenberg–Marquardt method, which are respectively equivalent to Newton flow forward Euler, gradient flow forward Euler, and gradient flow backward Euler. We finally show an unconditional descent property for a generalized Gauß–Newton flow, which is linked to Krylov–Gauß–Newton methods for large-scale nonlinear least squares problems. We provide numerical results for large-scale problems: An academic generalized Rosenbrock function and a real-world bundle adjustment problem from 3D reconstruction based on 2D images.


2006 ◽  
Vol 129 (2) ◽  
pp. 255-264 ◽  
Author(s):  
Xin Lei ◽  
Cliff J. Lissenden

Discontinuously reinforced aluminum (DRA) is currently used where design considerations include specific stiffness, tailorable coefficient of thermal expansion, or wear resistance. Plastic deformation plays a role in failures due to low cycle fatigue or simple ductile overload. DRA is known to exhibit pressure dependent yielding. Plastic deformation in metals is widely regarded to be incompressible, or very nearly so. A continuum plasticity model is developed that includes a Drucker–Prager pressure dependent yield function, plastic incompressibility via a nonassociative Prandtl–Reuss flow rule, and a generalized Armstrong–Frederick kinematic hardening law. The model is implemented using a return mapping algorithm with backward Euler integration for stability and the Newton method to determine the plastic multiplier. Material parameters are characterized from uniaxial tension and uniaxial compression experimental results. Model predictions are compared to experimental results for a nonproportional compression–shear load path. The tangent stiffness tensor is nonsymmetric because the flow rule is not associated with the yield function, which means that the commonly used algorithms that require symmetric matrices cannot be used with this material model. Model correlations with tension and compression loadings are excellent. Model predictions of shear and nonproportional compression–shear loadings are reasonably good. The nonassociative flow rule could not be validated by comparison of the plastic strain rate direction with the yield function and the flow potential due to scatter in the experimental results. The model is capable of predicting the material response obtained in the experiments, but additional validation is necessary for the condition of high hydrostatic pressure.


Author(s):  
B. Aparna ◽  
S. Madhavi ◽  
G. Mounika ◽  
P. Avinash ◽  
S. Chakravarthi

We propose a new design for large-scale multimedia content protection systems. Our design leverages cloud infrastructures to provide cost efficiency, rapid deployment, scalability, and elasticity to accommodate varying workloads. The proposed system can be used to protect different multimedia content types, including videos, images, audio clips, songs, and music clips. The system can be deployed on private and/or public clouds. Our system has two novel components: (i) method to create signatures of videos, and (ii) distributed matching engine for multimedia objects. The signature method creates robust and representative signatures of videos that capture the depth signals in these videos and it is computationally efficient to compute and compare as well as it requires small storage. The distributed matching engine achieves high scalability and it is designed to support different multimedia objects. We implemented the proposed system and deployed it on two clouds: Amazon cloud and our private cloud. Our experiments with more than 11,000 videos and 1 million images show the high accuracy and scalability of the proposed system. In addition, we compared our system to the protection system used by YouTube and our results show that the YouTube protection system fails to detect most copies of videos, while our system detects more than 98% of them.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Hongyi Zhang ◽  
Xiaowei Zhan ◽  
Bo Li

AbstractSimilarity in T-cell receptor (TCR) sequences implies shared antigen specificity between receptors, and could be used to discover novel therapeutic targets. However, existing methods that cluster T-cell receptor sequences by similarity are computationally inefficient, making them impractical to use on the ever-expanding datasets of the immune repertoire. Here, we developed GIANA (Geometric Isometry-based TCR AligNment Algorithm) a computationally efficient tool for this task that provides the same level of clustering specificity as TCRdist at 600 times its speed, and without sacrificing accuracy. GIANA also allows the rapid query of large reference cohorts within minutes. Using GIANA to cluster large-scale TCR datasets provides candidate disease-specific receptors, and provides a new solution to repertoire classification. Querying unseen TCR-seq samples against an existing reference differentiates samples from patients across various cohorts associated with cancer, infectious and autoimmune disease. Our results demonstrate how GIANA could be used as the basis for a TCR-based non-invasive multi-disease diagnostic platform.


1992 ◽  
Vol 114 (4) ◽  
pp. 847-857 ◽  
Author(s):  
J. H. Wagner ◽  
B. V. Johnson ◽  
R. A. Graziani ◽  
F. C. Yeh

Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large-scale, multipass, heat transfer model with both radially inward and outward flow. Trip strips on the leading and trailing surfaces of the radial coolant passages were used to produce the rough walls. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages: coolant-to-wall temperature ratio, Rossby number, Reynolds number, and radius-to-passage hydraulic diameter ratio. The first three of these four parameters were varied over ranges that are typical of advanced gas turbine engine operating conditions. Results were correlated and compared to previous results from stationary and rotating similar models with trip strips. The heat transfer coefficients on surfaces, where the heat transfer increased with rotation and buoyancy, varied by as much as a factor of four. Maximum values of the heat transfer coefficients with high rotation were only slightly above the highest levels obtained with the smooth wall model. The heat transfer coefficients on surfaces where the heat transfer decreased with rotation, varied by as much as a factor of three due to rotation and buoyancy. It was concluded that both Coriolis and buoyancy effects must be considered in turbine blade cooling designs with trip strips and that the effects of rotation were markedly different depending upon the flow direction.


Author(s):  
Mahdi Esmaily Moghadam ◽  
Yuri Bazilevs ◽  
Tain-Yen Hsia ◽  
Alison Marsden

A closed-loop lumped parameter network (LPN) coupled to a 3D domain is a powerful tool that can be used to model the global dynamics of the circulatory system. Coupling a 0D LPN to a 3D CFD domain is a numerically challenging problem, often associated with instabilities, extra computational cost, and loss of modularity. A computationally efficient finite element framework has been recently proposed that achieves numerical stability without sacrificing modularity [1]. This type of coupling introduces new challenges in the linear algebraic equation solver (LS), producing an strong coupling between flow and pressure that leads to an ill-conditioned tangent matrix. In this paper we exploit this strong coupling to obtain a novel and efficient algorithm for the linear solver (LS). We illustrate the efficiency of this method on several large-scale cardiovascular blood flow simulation problems.


Author(s):  
S. V. Subramanian ◽  
R. Bozzola ◽  
Louis A. Povinelli

The performance of a three dimensional computer code developed for predicting the flowfield in stationary and rotating turbomachinery blade rows is described in this study. The four stage Runge-Kutta numerical integration scheme is used for solving the governing flow equations and yields solution to the full, three dimensional, unsteady Euler equations in cylindrical coordinates. This method is fully explicit and uses the finite volume, time marching procedure. In order to demonstrate the accuracy and efficiency of the code, steady solutions were obtained for several cascade geometries under widely varying flow conditions. Computed flowfield results are presented for a fully subsonic turbine stator and a low aspect ratio, transonic compressor rotor blade under maximum flow and peak efficiency design conditions. Comparisons with Laser Anemometer measurements and other numerical predictions are also provided to illustrate that the present method predicts important flow features with good accuracy and can be used for cost effective aerodynamic design studies.


1998 ◽  
Vol 120 (1) ◽  
pp. 24-34 ◽  
Author(s):  
C. Fu ◽  
D. L. McDowell ◽  
I. C. Ume

A finite element procedure using a semi-implicit time-integration scheme has been developed for a cyclic thermoviscoplastic constitutive model for Pb-Sn solder and OFHC copper, two common metallic constituents in electronic packaging applications. The scheme has been implemented in the commercial finite element (FE) code ABAQUS (1995) via the user-defined material subroutine, UMAT. Several single-element simulations are conducted to compare with previous test results, which include monotonic tensile tests, creep tests, and a two-step ratchetting test for 62Sn36Pb2Ag solder; a nonproportional axial-torsional test and a thermomechanical fatigue (TMF) test for OFHC copper. At the constitutive level, we also provide an adaptive time stepping algorithm, which can be used to improve the overall computation efficiency and accuracy especially in large-scale FE analyses. We also compare the computational efforts of fully backward Euler and the proposed methods. The implementation of the FE procedure provides a guideline to apply user-defined material constitutive relations in FE analyses and to perform more sophisticated thermomechanical simulations. Such work can facilitate enhanced understanding thermomechanical reliability issue of solder and copper interconnects in electronic packaging applications.


Water ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 785
Author(s):  
Arman Rokhzadi ◽  
Musandji Fuamba

This paper studies the air pressurization problem caused by a partially pressurized transient flow in a reservoir-pipe system. The purpose of this study is to analyze the performance of the rigid column model in predicting the attenuation of the air pressure distribution. In this regard, an analytic formula for the amplitude and frequency will be derived, in which the influential parameters, particularly, the driving pressure and the air and water lengths, on the damping can be seen. The direct effect of the driving pressure and inverse effect of the product of the air and water lengths on the damping will be numerically examined. In addition, these numerical observations will be examined by solving different test cases and by comparing to available experimental data to show that the rigid column model is able to predict the damping. However, due to simplified assumptions associated with the rigid column model, the energy dissipation, as well as the damping, is underestimated. In this regard, using the backward Euler implicit time integration scheme, instead of the classical fourth order explicit Runge–Kutta scheme, will be proposed so that the numerical dissipation of the backward Euler implicit scheme represents the physical dissipation. In addition, a formula will be derived to calculate the appropriate time step size, by which the dissipation of the heat transfer can be compensated.


Fluids ◽  
2021 ◽  
Vol 6 (4) ◽  
pp. 148
Author(s):  
Seyed Amin Nabavizadeh ◽  
Himel Barua ◽  
Mohsen Eshraghi ◽  
Sergio D. Felicelli

A multi-distribution lattice Boltzmann Bhatnagar–Gross–Krook (BGK) model with a multiple-grid lattice Boltzmann (MGLB) model is proposed to efficiently simulate natural convection over a wide range of Prandtl numbers. In this method, different grid sizes and time steps for heat transfer and fluid flow equations are chosen. The model is validated against natural convection in a square cavity, since extensive benchmark solutions are available for that problem. The proposed method can resolve the computational difficulty in simulating problems with very different time scales, in particular, when using extremely low or high Prandtl numbers. The technique can also enhance computational speed and stability while keeping the simplicity of the BGK method. Compared with the conventional lattice Boltzmann method, the simulation time can be reduced up to one-tenth of the time while maintaining the accuracy in an acceptable range. The proposed model can be extended to other lattice Boltzmann collision models and three-dimensional cases, making it a great candidate for large-scale simulations.


Sign in / Sign up

Export Citation Format

Share Document