On the numerical algorithm for isotropic kinematic hardening with the Armstrong Frederick evolution of the back stress
|
|
- Ruth Henderson
- 5 years ago
- Views:
Transcription
1 Comut. Methods Al. Mech. Engrg. 191 (2002) On the numerical algorithm for isotroic kinematic hardening with the Armstrong Frederick evolution of the back stress Vlado A. Lubarda *, David J. Benson Deartment of Mechanical and Aerosace Engineering, University of California, 9500 Gilman Drive, San Diego, La Jolla, CA , USA Received 20 Setember 2001 Abstract The algorithmic consistency conditions are derived which yield the lastic loading index for the combined isotroic kinematic hardening lasticity with the Armstrong Frederick evolution equation for the back stress. Rate-indeendent elastic lastic, and rate-deendent viscoelastic lastic and elastic viscolastic material models are all encomassed by the analysis. The directions of the lastic and viscolastic increments of strain are determined from a develoed numerical algorithm, in conjunction with the radial return method. Numerical evaluations are resented to comare the redictions of various models under cyclic loading in simle shear. The effects of different material arameters on stress resonse are discussed. Ó 2002 Elsevier Science B.V. All rights reserved. Keywords: Algorithmic consistency; Radial return; Viscolasticity; Anisotroic hardening 1. Introduction Numerical algorithms for treating the governing equations of lastic resonse have been well develoed for various rate-indeendent and rate-deendent constitutive models. Particular attention was given to both continuum henomenological models and hysical models based on the single crystal lasticity. The recent books by Simo and Hughes [1], and Belytschko et al. [2] can be consulted for the review and assessment of the roosed methods and their origin. The urose of this aer is to elaborate on the algorithmic consistency conditions and numerical algorithms for the combined isotroic kinematic hardening lasticity with the Armstrong Frederick evolution equation for the back stress. The rate-indeendent elastic lastic, and the rate-deendent viscoelastic lastic, and elastic viscolastic constitutive models are encomassed by the analysis. The directions of the lastic and viscolastic increments of strain are determined from a develoed algorithm using the method of radial return and a numerical recie based * Corresonding author. Tel.: ; fax: address: vlubarda@ucsd.edu (V.A. Lubarda) /02/$ - see front matter Ó 2002 Elsevier Science B.V. All rights reserved. PII: S (02)
2 3584 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) on the generalized midoint rule. The obtained results extend the analysis of Simo and Hughes [1], which is restricted to combined isotroic kinematic hardening with the Prager s tye evolution equation for the back stress, and of Zienkiewicz and Cormeau [3] and Hornberger and Stamm [4]. Numerical evaluations are resented to comare the redictions of various models under cyclic forward and backward loading in simle shear. The effects of hardening and viscous arameters on stress resonse are then discussed. 2. Analytical background In the model of combined isotroic kinematic hardening J 2 lasticity the yield condition at an arbitrary state of deformation is secified by 1 ffiffiffi ks ak ¼Kð#Þ: ð1þ 2 The norm kk acts on a second-order tensor A by the square root trace oeration kak ¼ðA : AÞ 1=2 : ð2þ The deviatoric art of the Cauchy stress r in Eq. (1) is S ¼ r ðtrrþi=3. The rectangular comonents of the second-order unit tensor I are the Kronecker deltas d ij. The center of the yield circle in the deviatoric lane trr ¼ 0 is the back stress a. This is a traceless tensor, whose counterart in the lane trr ¼ constant is b ¼ a þ 1trr 3 I: ð3þ The current radius of the yield circle is K, which is assumed to deend on the equivalent lastic strain # ¼ ffiffiffi Z t 2 kde k: ð4þ 0 The integral is taken along the deformation ath from the initial to current instant of time on some quasistatic scale t. The lastic art of the strain increment de is denoted by de. The rate-indeendent materials are considered which comly with the normality rule, so that de is codirectional with the outward normal to the current yield surface in stress sace, i.e., de ¼ dc S a ; dc > 0: ð5þ ks ak The loading index dc is related to the increment of equivalent lastic strain by d# ¼ ffiffi 2 dc; ð6þ since kde k¼dc. According to the Armstrong Frederick [5] model, the evolution equation for the back stress a is da ¼ c 1 de c 2 akde k: ð7þ The arameters c 1 and c 2 can be either constants or functions of #. Ifc 2 ¼ 0, the Prager [6,7] hardening model is obtained in which the instantaneous translation of the yield surface is in the direction of the lastic increment of strain de, thus in the direction of the stress difference S a. The addition of the second term (the so-called recall term) in Eq. (7), roortional to a, modifies the direction of this translation and gives rise to hardening moduli for the reversed lastic loading that are in better agreement with exerimental data. Indeed, from the consistency condition the loading index can be exressed as dc ¼ 1 ðs aþ : ds ; ð8þ H ks ak
3 where H is the hardening modulus given by H ¼ 2 dk d# þ c ðs aþ : a 1 c 2 ffiffi : ð9þ 2 K This clearly deends not only on # and Kð#Þ, but also exlicitly on the current state of stress S and a. For examle, in uniaxial tension r along the direction 1 we have r ij ¼ rd i1 d j1 ; b ij ¼ b 1 d i1 d j1 þ b 2 ðd i2 d j2 þ d i3 d j3 Þ; ð10þ and S ij ¼ rd i1 d j1 1d 3 ij ; ð11þ a ij ¼ ad i1 d j1 1aðd 2 i2d j2 þ d i3 d j3 Þ¼ 3ad 2 i1d j1 1d 3 ij ; ð12þ with the connections b 1 ¼ a þ 1r; 3 b 2 ¼ 1a þ 1 r: 2 3 ð13þ The stress deendent term in Eq. (9) therefore becomes ðs aþ : a ¼ r 3a a ¼ ffiffiffi 3 2 Ka: ð14þ The minus sign corresonds to lastic loading in the reversed direction (at the stress level r ¼ 3a=2 ffiffi 3 K). In the numerical evaluations it will be assumed that the radius of the yield surface, initially equal to K 0, saturates to the value K 1 at large # (saturation kinematic hardening), such that Kð#Þ ¼K 0 þðk 1 K 0 Þ 1 ex b 0# : ð15þ K 1 K 0 The arameter b 0 is an aroriate constant whose hysical interretation is available from an easily verified roerty dk b 0 ¼ lim #!0 d# : ð16þ The total increment of strain is the sum of elastic and lastic arts. For infinitesimal elastic comonent of deformation, and for an isotroic material, we write de ¼ 1 2l ds þ de : ð17þ The deviatoric art of the strain increment de is denoted by de, and l is the elastic shear modulus. Thus, ds ¼ 2lde 2lde : ð18þ Here, we identify ds e ¼ 2lde V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) as the elastic art of the deviatoric stress increment ds, while ð19þ ds ¼ 2l de ¼ 2ldc S a ks ak ð20þ reresents its lastic art. The lastic loading condition is ds : de < 0; i:e:; ðs aþ : de > 0: ð21þ
4 3586 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) Viscoelastic lastic resonse A simle viscoelastic lastic model of material resonse is obtained by adding a viscous art of the strain increment to its elastic and lastic contributions [8]. If viscous resonse is governed by the Newton s viscosity law, we can write de ¼ 1 Sdt þ 1 S a ds þ dc g e 2l ks ak ; ð22þ where dt now reresents the actual time increment during which ds and de take lace. The viscosity coefficient is denoted by g e. The consistency condition for the continuing lastic deformation gives the lastic loading index dc according to Eq. (8), and the hardening modulus H according to Eq. (9). It can be readily shown that the stress increment ds is exressed in terms of the increments de and Sdt by " 1 ds ¼ 2l J 1 þ H=2l ðs aþðs aþ ks ak 2 # : de 1 Sdt : ð23þ g e The outer tensor roduct is denoted by, and J stands for the symmetric fourth-order unit tensor with the rectangular comonents ðd ik d jl þ d il d jk Þ=2. The lastic loading condition is ðs aþ : de 1 Sdt > 0: ð24þ g e The material time arameter s e will be conveniently introduced as s e ¼ g e 2l ; such that dt ¼ 1 dt : g e 2l s e ð25þ ð26þ 2.2. Elastic viscolastic resonse In an elastic viscolastic material model the strain increment is the sum of elastic and viscolastic arts, such that de ¼ 1 ds þ dev ð27þ 2l Using a viscolasticity model based on the concet of overstress [9 12], with anisotroic hardening accounted for via the back stress a, we may write de v ¼ dt D ks ak ffiffiffi E S a 2 Kð#Þ g ks ak : ð28þ The viscosity arameter is g, and the Macauley brackets hi are used, such that hwi ¼ðwþjwjÞ=2, for any scalar w. Consequently, the stress increment is ds ¼ 2l de dt s D ks ak ffiffiffi E S a 2 Kð#Þ ks ak : ð29þ
5 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) The material time arameter s is introduced as s ¼ g 2l : ð30þ More involved constitutive equations of viscolastic resonse could also be used, such as those considered by Chaboche [13,14], Bammann [15], McDowell [16], M uller-hoee and Wriggers [17], Freed and Walker [18], and Lubarda and Benson [19]. 3. Algorithmic consistency condition for rate-indeendent resonse Fig. 1 shows a trace of the yield surface in the deviatoric lane at the two nearby instants t n and t nþ1. The center of the first yield surface is a n, and its radius is Kð# n Þ. The index n designates the number of deformation increments through which the current state has been reached. The center of the subsequent yield surface, after an increment of strain de nþ1 has been alied along the lastic loading branch, is denoted by a nþ1, and its radius by Kð# nþ1 Þ. We have # nþ1 ¼ # n þ d# nþ1 ; ð31þ and a nþ1 ¼ a n þ da nþ1 : ð32þ The radial return method with the backward Eulerian scheme will be emloyed. The udated stress state on the new yield surface is defined by S nþ1 ¼ S nþ1 2lde nþ1 : ð33þ The elastic redictor state of stress is S nþ1 ¼ S n þ ds e ¼ S n þ 2lde; while the remaining term on the right-hand side of Eq. (33) is the lastic corrector term [1,20], ð34þ ds nþ1 ¼ 2lde nþ1 ¼ 2ldc nþ1 k k : ð35þ Fig. 1. The yield surface with the center at the back stress a n, and the subsequent yield surface with the center at a nþ1 ¼ a n þ da nþ1. The elastic redictor state of stress S nþ1 is brought to the stress state S nþ1 on the yield surface by the lastic corrector 2lde nþ1.
6 3588 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) Consequently, and S nþ1 ¼ S n þ 2lde 2l dc nþ1 k k ð36þ ¼ S n a n þ 2lde 2l dc nþ1 k k da nþ1: ð37þ The increment of the back stress will be comuted from the Armstrong Frederick evolution equation, and the algorithmic recie da nþ1 ¼ c 1 de nþ1 c 2kde nþ1 k½ha n þð1 hþa nþ1 Š; 0 6 h 6 1: ð38þ The generalized midoint rule is used for the recall term (the value h ¼ 0 and 1 corresonding to the forward and backward Euler rules, resectively, and h ¼ 1=2 to the midoint rule). Substitution of Eq. (5) into Eq. (38) gives da nþ1 ¼ c 1 dc nþ1 k k c 2 dc nþ1 ½ha n þð1 hþa nþ1 Š; ð39þ having in mind that kde nþ1 k¼dc nþ1. Eq. (39) can be rewritten as da nþ1 ¼ a nþ1 dc nþ1 c 2 a n ; ð40þ k k c 1 where c 1 a nþ1 ¼ : ð41þ 1 þ c 2 ð1 hþdc nþ1 When Eq. (40) is inserted into Eq. (37), we obtain þ ð2l þa nþ1 Þdc nþ1 k k ¼ B n; ð42þ with B n ¼ S n a n þ 2l de þ b nþ1 dc nþ1 a n ð43þ and b nþ1 ¼ c 2 c 2 a nþ1 ¼ : ð44þ c 1 1 þ c 2 ð1 hþdc nþ1 To derive an equation for the loading arameter dc nþ1 we take a trace roduct of Eq. (42) with itself, which gives ks nþ1 a nþ1 kþð2l þ a nþ1 Þdc nþ1 2 ¼kS n a n k 2 þk2lde þ b nþ1 dc nþ1 a n k 2 þ 2ðS n a n Þ : ð2lde þ b nþ1 dc nþ1 a n Þ: ð45þ In view of Eq. (1), alied at the instants t n and t nþ1, Eq. (45) is rewritten as ffiffiffi 2 Kð#n þ ffiffi h 2 dcnþ1 Þþð2lþa nþ1 Þdc nþ1 ¼ 2K 2 ð# n Þ þk2ldeþb nþ1 dc nþ1 a n k 2 þ 2ðS n a n Þ : ð2l de þ b nþ1 dc nþ1 a n Þi 1=2: ð46þ
7 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) This is the key equation of the numerical method. It reresents an algorithmic consistency condition for the considered hardening model with the Armstrong Frederick evolution of the back stress. Being in the form of a nonlinear equation for the loading arameter dc nþ1, its solution is sought by numerical means. If the Newton s iterative method is emloyed, the function F ¼ F ðdc nþ1 Þ is defined by F ¼ ffiffi 2 Kð#n þ ffiffiffi 2 dcnþ1 Þþð2lþa nþ1 Þdc nþ1 h 1=2: 2K 2 ð# n Þ þk2ldeþb nþ1 dc nþ1 a n k 2 þ 2ðS n a n Þ : ð2lde þ b nþ1 dc nþ1 a n Þi ð47þ whose gradient with resect to dc nþ1 is F 0 ¼ 2 dk d# ð# n þ ffiffiffi 2 dcnþ1 Þþa 0 nþ1 dc nþ1 þ 2l þ a nþ1 h 2K 2 ð# n Þ þk2ldeþb nþ1 dc nþ1 a n k 2 þ 2ðS n a n Þ : ð2lde þ b nþ1 dc nþ1 a n Þ ðb 0 nþ1 dc nþ1 þ b nþ1 ÞðS n a n þ 2lde þ b nþ1 dc nþ1 a n Þ : a n : i 1=2 ð48þ The iterative rocedure to find the zero value of F is then based on the recie dc iþ1 nþ1 ¼ dci nþ1 F ðdci nþ1 Þ F 0 ðdc i nþ1 Þ : ð49þ The iterations are ended when a desired accuracy has been reached, i.e., when jf ðdc iþ1 nþ1þj falls to within a rescribed error tolerance. The convergence is guaranteed because F is a convex function of dc nþ1. In the iterative rocess we use # i nþ1 ¼ # n þ ffiffi 2 dc i ð50þ nþ1 : It is assumed in Eq. (38) that c 1 and c 2 are either constants or deendent on # n only. If the generalization is made so that c i ¼ w i c i ð# n Þþð1 w i Þc i ð# nþ1 Þði ¼ 1; 2Þ, only the rates a 0 nþ1 and b0 nþ1 are affected. If c 2 ¼ 0, the numerical algorithm corresonding to Prager s evolution of the back stress is obtained. In this case, Eq. (46) simlifies to ffiffi 2 Kð#n þ ffiffi 2 dcnþ1 Þþð2lþc 1 Þdc nþ1 ¼kS nþ1 a nk; in accord with the result from Simo and Hughes [1]. It remains to secify the direction of de nþ1, i.e., to determine the normalized tensor ð51þ k k : To that goal, we rewrite Eq. (42) as ½k kþð2lþa nþ1 Þdc nþ1 Š k k ¼ B n: By taking a norm there follows: k kþð2lþa nþ1 Þdc nþ1 ¼kB n k: Consequently, Eq. (54) gives k k ¼ B n kb n k : ð52þ ð53þ ð54þ ð55þ
8 3590 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) This is a desired result, because, with the reviously determined loading arameter dc nþ1, the tensor B n is comletely secified. The back stress can now be udated by using Eq. (40). 4. Algorithmic consistency conditions for rate-deendent resonse 4.1. Viscoelastic lastic resonse The udated stress state on the new yield surface is given according to the backward Eulerian scheme by S nþ1 ¼ S nþ1 2l de nþ1 : ð56þ The viscoelastic redictor state of stress is from Eq. (22) S nþ1 ¼ S n þ 2lde ½uS n þð1 uþs nþ1 Š dt ; ð57þ s e where the generalized midoint rule with 0 6 u 6 1 is used for the viscous contribution, and s e ¼ g e =2l. The remaining term on the right-hand side of Eq. (56) is the lastic corrector term ds nþ1 ¼ 2lde nþ1 ¼ 2ldc nþ1 k k : ð58þ Introducing the arameter q ¼ 1 þð1 uþ dt 1 ; ð59þ s e it is readily found that dt S nþ1 ¼ qs n þ q 2lde us n 2lqdc s nþ1 ð60þ e k k and dt ¼ qðs n a n Þþq 2lde us n 2lqdc s nþ1 e k k þðq 1Þa n da nþ1 : ð61þ If the increment of back stress is comuted from the Armstrong Frederick evolution equation, and the algorithmic recie of Eqs. (38) and (40), there follows: þð2lq þ a nþ1 Þdc nþ1 k k ¼ B n; ð62þ where dt B n ¼ qðs n a n Þþq 2lde us n þ A n ð63þ s e and A n ¼ðb nþ1 dc nþ1 þ q 1Þa n : ð64þ The equation for the loading arameter dc nþ1 can be deduced by taking a trace roduct of Eq. (62) with itself. This gives
9 ks nþ1 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) a nþ1 kþð2lq þ a nþ1 Þdc nþ1 ¼ q 2 ks n a n k 2 þ dt 2 q 2l de us n þ A n s e dt þ 2qðS n a n Þ : q 2l de us n s e þ A n : ð65þ When k k¼ ffiffiffi 2 Kð#n þ ffiffiffi 2 dcnþ1 Þ; ks n a n k¼ ffiffi 2 Kð#n Þ; ð66þ are substituted into Eq. (65) we obtain the algorithmic consistency condition in the form of a nonlinear equation for the lastic loading index dc nþ1. The equation can be solved by the Newton s iterative method along the same lines as described in the revious section. If g e!1, i.e., dt=s e! 0, then q ¼ 1, and Eq. (65) reduces to Eq. (45). It can be readily verified that k k ¼ B n kb n k ; which secifies the direction of the lastic increment of strain de nþ1, i.e., B n de nþ1 ¼ dc nþ1 kb n k : ð67þ ð68þ The udated stress state follows from Eq. (60) Elastic viscolastic resonse The radial return method with the backward Eulerian scheme will again be emloyed. The udated stress state is defined by S nþ1 ¼ S nþ1 2ldev nþ1 : ð69þ The elastic redictor state of stress is (see Fig. 2) S nþ1 ¼ S n þ 2lde; ð70þ while the remaining term on the right-hand side of Eq. (69) is the viscolastic corrector term. The viscolastic art of the strain increment is defined by Eq. (28) as de v nþ1 ¼ dc nþ1 k k ; ð71þ where 2ldc nþ1 ¼ dt h k k ffiffi 2 Kð#nþ1 Þi ; ð72þ s rovided that the term within the square brackets is ositive. In the case of the Armstrong Frederick evolution of the back stress, we can write from Eq. (45) k k¼ ð2lþa nþ1 Þdc nþ1 h 1=2: þ ks n a n k 2 þk2lde þ b nþ1 dc nþ1 a n k 2 þ 2ðS n a n Þ : ð2lde þ b nþ1 dc nþ1 a n Þi ð73þ
10 3592 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) Fig. 2. The static yield surface with the center at the back stress a n. The static state of stress S 0 n is on the yield surface, while the corresonding dynamic state of stress S n is outside the yield surface. The subsequent static yield surface has the center at a nþ1 ¼ a n þ da nþ1, with the static and dynamic states of stress S 0 nþ1 and S nþ1. The elastic redictor state of stress S nþ1 is brought to the stress state S nþ1 by the viscolastic corrector 2lde v nþ1. As in Section 3, the arameters are introduced c 1 a nþ1 ¼ ; ð74þ 1 þ c 2 ð1 hþdc nþ1 b nþ1 ¼ c 2 c 2 a nþ1 ¼ : ð75þ c 1 1 þ c 2 ð1 hþdc nþ1 The substitution of Eq. (73) into Eq. (72) yields a nonlinear equation for the viscolastic loading index dc nþ1. This is 2l 1 þ a nþ1 2l þ s dc dt nþ1 þ ffiffi 2 Kð#n þ ffiffi 2 dcnþ1 Þ h 1=2: ¼ ks n a n k 2 þk2lde þ b nþ1 dc nþ1 a n k 2 þ 2ðS n a n Þ : ð2lde þ b nþ1 dc nþ1 a n Þi ð76þ With a determined dc nþ1, the direction of the viscolastic increment of strain de v nþ1 k k ¼ B n kb n k ; is secified by where B n ¼ S n a n þ 2l de þ b nþ1 dc nþ1 a n : ð78þ The viscolastic model corresonding to the Prager s evolution of the back stress is obtained by setting c 2 ¼ 0 and a nþ1 ¼ c 1. The resulting equation for dc nþ1 is 2l 1 þ c 1 2l þ s dc dt nþ1 þ ffiffiffi 2 Kð#n þ ffiffi 2 dcnþ1 Þ¼kS nþ1 a nk: ð79þ In the limit s =dt! 0, we recover the equations for the loading index of the rate-indeendent lasticity from Section 3. In articular, Eq. (76) reduces to Eq. (45), while Eq. (79) reduces to Eq. (51). ð77þ
11 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) Numerical results and conclusions To illustrate the alication of the derived formulas and comare the redictions of various constitutive models, we consider a strain history shown in Fig. 3. The shear strain of amount is first alied at constant strain-rate of 100 s 1, followed by a segment of constant strain during s, and the reversal of strain to zero value at the strain-rate of negative 100 s 1. The material is subsequently subjected to an oosite ulse of shear strain. The calculated time history of the corresonding shear stress in the case of viscoelastic lastic model is shown in Fig. 4. The solid curve corresonds to Prager hardening, and the dashed curve to Armstrong Frederick hardening. The utilized material arameters were: l ¼ 80 GPa, K 0 ¼ 100 MPa, K 1 ¼ 200 MPa, b 0 ¼ 2 GPa, and g e ¼ 4 MPa s. In the case of Prager hardening the constant c 1 ¼ 10 GPa, while for the Armstrong Frederick hardening we took c 1 ¼ 10 GPa and c 2 ¼ 100 MPa. The stress relaxations in Fig. 4 occur during the time intervals of constant strain. Fig. 5 shows the stress strain curves obtained in the two cases. The difference in the hardening rates between the two models is clear. The hardening rates can be further adjusted by varying the material arameters, as required by a articular exerimental data. For examle, the effect of the arameter c 2 is illustrated in Fig. 6, which shows that the increase of this arameter decreases the hardening rate and makes the stress increase more gradual. The stress time variation corresonding to shear strain inut from Fig. 3 in the case of elastic viscolastic model is shown in Fig. 7. The two curves again corresond to two different hardening models. Fig. 3. The time history of alied forward and backward cycles of shear strain. Fig. 4. The time history of the shear stress corresonding to shear strain cycles from Fig. 3, in the case of a viscoelastic lastic material model. The darker curve is for the Armstrong Frederick hardening, and the lighter curve for the Prager kinematic hardening.
12 3594 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) Fig. 5. The shear stress shear strain curves for a viscoelastic lastic material model corresonding to time variations of strain and stress deicted in Figs. 3 and 4. Fig. 6. The effect of the Armstrong Frederick arameter c 2 on the stress resonse of a viscoelastic lastic material. The darker curve is for c 2 ¼ 50 MPa, and the lighter curve for c 2 ¼ 100 MPa. Fig. 7. The time history of the shear stress corresonding to shear strain cycles from Fig. 3, in the case of an elastic viscolastic material model. The darker curve is for the Perzyna Armstrong Frederick model, and the lighter curve for the Perzyna Prager model. The utilized material arameters were as in the revious case, excet that the viscosity arameter was g ¼ 2 MPa s. Fig. 8 shows the corresonding stress strain curve and the hysteretic loos obtained in the two cases. The effect of the Armstrong Frederick arameter c 2 on the stress resonse is illustrated in Fig. 9.
13 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) Fig. 8. The shear stress shear strain curves for an elastic viscolastic material model corresonding to time variations of strain and stress deicted in Figs. 3 and 7. Fig. 9. The effect of the Armstrong Frederick arameter c 2 on the stress resonse of an elastic viscolastic material. The darker curve is for c 2 ¼ 50 MPa, and the lighter curve for c 2 ¼ 100 MPa. Fig. 10. The effect of the viscosity arameter g on the stress resonse of an elastic viscolastic material. The darker curve is for g ¼ 4 MPa s, and the lighter curve for g ¼ 2 MPa s. More gradual change of stress is achieved by higher values of this arameter. Fig. 10 shows the effect of the viscosity arameter g. The resonse in the straining under constant rate of strain is stiffer, and the stress relaxation under constant strain less raid for higher values of g.
14 3596 V.A. Lubarda, D.J. Benson / Comut. Methods Al. Mech. Engrg. 191 (2002) Acknowledgements Research funding rovided by the Los Alamos National Laboratories is kindly acknowledged. References [1] J.C. Simo, T.J.R. Hughes, Comutational Inelasticity, Sringer, New York, [2] T. Belytschko, W.K. Liu, B. Moran, Nonlinear Finite Elements for Continua and Structures, Wiley, New York, [3] O.C. Zienkiewicz, I.C. Cormeau, Visco-lasticity lasticity and cree in elastic solids a unified numerical solution aroach, Int. J. Numer. Meth. Engrg. 8 (1974) [4] K. Hornberger, H. Stamm, An imlicit integration algorithm with a rojection method for viscolastic constitutive equations, Int. J. Numer. Meth. Engrg. 28 (1989) [5] P.J. Armstrong, C.O. Frederick, A mathematical reresentation of the multiaxial Bauschinger effect, GEGB Re. RD/B/N 731, [6] W. Prager, The theory of lasticity: a survey of recent achievements (James Clayton Lecture), Proc. Inst. Mech. Engrg. 169 (1955) [7] W. Prager, A new method of analyzing stresses and strains in work-hardening lastic solids, J. Al. Mech. 23 (1956) [8] F.J. Zerilli, R.W. Armstrong, Thermal activation based constitutive equations for olymers, J. Phys. IV 10 (2000) [9] V.V. Sokolovskii, Proagation of elastic viscolastic waves in bars, Prikl. Mat. Mekh. 12 (1948) , in Russian. [10] L.E. Malvern, The roagation of longitudinal waves of lastic deformation in a bar of material exhibiting a strain-rate effect, J. Al. Mech. 18 (1951) [11] P. Perzyna, The constitutive equations for rate sensitive lastic materials, Q. Al. Math. 20 (1963) [12] P. Perzyna, Fundamental roblems in viscolasticity, Adv. Al. Mech. 9 (1966) [13] J.-L. Chaboche, Constitutive equations for cyclic lasticity and cyclic viscolasticity, Int. J. Plasticity 5 (1989) [14] J.-L. Chaboche, Unified cyclic viscolastic constitutive equations develoment, caabilities, and thermodynamic framework, in: A.S. Krausz, K. Krausz (Eds.), Unified Constitutive Laws for Plastic Deformation, Academic Press, San Diego, 1996, [15] D.J. Bammann, Modeling temerature and strain rate deendent large deformations of metals, Al. Mech. Rev. 43 (5/2) (1990) S312 S319. [16] D.L. McDowell, A nonlinear kinematic hardening theory for cyclic thermolasticity and thermoviscolasticity, Int. J. Plasticity 8 (1992) [17] N. M uller-hoee, P. Wriggers, Numerical integration of viscolastic material laws at large strains with alication to imact roblems, Comut. Struct. 45 (1992) [18] A.D. Freed, K.P. Walker, Viscolasticity with cree and lasticity bounds, Int. J. Plasticity 9 (1993) [19] V.A. Lubarda, D.J. Benson, On the evolution equation for the rest stress in rate-deendent lasticity, Int. J. Plasticity 18 (2002) [20] R.D. Krieg, D.B. Krieg, Accuracies of numerical solution methods for the elastic-erfectly lastic model, J. Pressure Vessel Techn. 99 (1977)
FE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationThe effect of dynamic bending moments on the ratchetting behavior of stainless steel pressurized piping elbows
International Journal of echanical Engineering and Alications 2014; 2(2): 31-37 Published online ay 30, 2014 (htt://www.scienceublishinggrou.com/j/ijmea) doi: 10.11648/j.ijmea.20140202.12 The effect of
More informationCharacterization of Material Parameters
Proceedings of the World Congress on Engineering 29 Vol II WCE 29, July 1-3, 29, London, U.K. Characterization of Material Parameters S. M. Humayun Kabir, Tae-In Yeo, Sang-Ho Kim Abstract The resent work
More informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More information8.7 Associated and Non-associated Flow Rules
8.7 Associated and Non-associated Flow Rules Recall the Levy-Mises flow rule, Eqn. 8.4., d ds (8.7.) The lastic multilier can be determined from the hardening rule. Given the hardening rule one can more
More informationLower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
More informationFinite Element Analysis of V-Bending of Polypropylene Using Hydrostatic-Pressure-Dependent Plastic Constitutive Equation*
Materials Transactions, Vol. 48, No. 1 (7). 6 to 664 #7 The Jaan Society for Technology of Plasticity Finite Element Analysis of V-Bending of Polyroylene Using Hydrostatic-Pressure-Deendent Plastic onstitutive
More informationANALYSIS OF ULTRA LOW CYCLE FATIGUE PROBLEMS WITH THE BARCELONA PLASTIC DAMAGE MODEL
XII International Conerence on Comutational Plasticity. Fundamentals and Alications COMPLAS XII E. Oñate, D.R.J. Owen, D. Peric and B. Suárez (Eds) ANALYSIS OF ULTRA LOW CYCLE FATIGUE PROBLEMS WITH THE
More informationExact Solutions in Finite Compressible Elasticity via the Complementary Energy Function
Exact Solutions in Finite Comressible Elasticity via the Comlementary Energy Function Francis Rooney Deartment of Mathematics University of Wisconsin Madison, USA Sean Eberhard Mathematical Institute,
More informationPressure-sensitivity Effects on Toughness Measurements of Compact Tension Specimens for Strain-hardening Solids
American Journal of Alied Sciences (9): 19-195, 5 ISSN 1546-939 5 Science Publications Pressure-sensitivity Effects on Toughness Measurements of Comact Tension Secimens for Strain-hardening Solids Abdulhamid
More informationComputational Solid Mechanics Computational Plasticity
Comutational Solid Mechanics Comutational Plasticity Chater 3. J Plasticity Models C. Agelet de Saracibar ETS Ingenieros de Caminos, Canales y Puertos, Universidad Politécnica de Cataluña (UPC), Barcelona,
More informationA compression line for soils with evolving particle and pore size distributions due to particle crushing
Russell, A. R. (2011) Géotechnique Letters 1, 5 9, htt://dx.doi.org/10.1680/geolett.10.00003 A comression line for soils with evolving article and ore size distributions due to article crushing A. R. RUSSELL*
More informationChurilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics
Churilova Maria Saint-Petersburg State Polytechnical University Deartment of Alied Mathematics Technology of EHIS (staming) alied to roduction of automotive arts The roblem described in this reort originated
More informationSession 5: Review of Classical Astrodynamics
Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve
More informationAdiabatic Shear Bands in Simple and Dipolar Plastic Materials
Adiabatic Shear Bands in Simle and Diolar Plastic Materials T W \-1RIGHT us Army Ballistic Research Laboratory Aberdeen Proving Ground, MD 215 R C BATRA University of Missouri-Rolla Rolla, Missouri 6541
More informationImplementing the New First and Second Differentiation of a General Yield Surface in Explicit and Implicit Rate- Independent Plasticity
Journal of Solid Mechanics Vol. 6, No. 3 (2014). 310-321 Imlementing the New First and Second Differentiation of a General Yield Surface in Exlicit and Imlicit Rate- Indeendent Plasticity F. Moayyedian,
More informationA Simple And Efficient FEM-Implementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars
Aalborg Universitet A Simle And Efficient FEM-Imlementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars Published in: Proceedings of the 9th Nordic Seminar on Comutational
More informationSiping Road 1239, , Shanghai, P.R. China
COMPARISON BETWEEN LINEAR AND NON-LINEAR KINEMATIC HARDENING MODELS TO PREDICT THE MULTIAXIAL BAUSCHINGER EFFECT M.A. Meggiolaro 1), J.T.P. Castro 1), H. Wu 2) 1) Department of Mechanical Engineering,
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationTHE INFLUENCE OF DISLOCATION DENSITY ON THE BEHAVIOUR OF CRYSTALLINE MATERIALS
THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume 7, Number 4/6, 358 365 THE INFLUENCE OF DISLOCATION DENSITY ON THE BEHAVIOUR OF CRYSTALLINE MATERIALS
More informationMHA042 - Material mechanics: Duggafrågor
MHA042 - Material mechanics: Duggafrågor 1) For a static uniaxial bar problem at isothermal (Θ const.) conditions, state principle of energy conservation (first law of thermodynamics). On the basis of
More informationFeedback-error control
Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller
More informationAdam Paweł Zaborski. 8 Plasticity. reloading. 1. Bauschinger s effect. 2. unchanged yielding limit. 3. isotropic hardening
8 lasticity Introduction Definitions loading/unloading/reloading words commonly used in lasticity lastic strain a ermanent strain that doesn t vanish after unloading, lastically assive rocess a rocess
More informationKeywords: pile, liquefaction, lateral spreading, analysis ABSTRACT
Key arameters in seudo-static analysis of iles in liquefying sand Misko Cubrinovski Deartment of Civil Engineering, University of Canterbury, Christchurch 814, New Zealand Keywords: ile, liquefaction,
More informationChapter 1 Fundamentals
Chater Fundamentals. Overview of Thermodynamics Industrial Revolution brought in large scale automation of many tedious tasks which were earlier being erformed through manual or animal labour. Inventors
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationA General Damage Initiation and Evolution Model (DIEM) in LS-DYNA
9th Euroean LS-YNA Conference 23 A General amage Initiation and Evolution Model (IEM) in LS-YNA Thomas Borrvall, Thomas Johansson and Mikael Schill, YNAmore Nordic AB Johan Jergéus, Volvo Car Cororation
More informationQuasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight Channel with a Square Cross Section
Article Nihon Reoroji Gakkaishi Vol.34, No.2, 105~113 (Journal of the Society of Rheology, Jaan) 2006 The Society of Rheology, Jaan Quasi-Three-Dimensional Simulation of Viscoelastic Flow through a Straight
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationFlow liquefaction instability prediction using finite elements
Acta Geotechnica (215) 1:83 1 DOI 1.17/s1144-14-342-z RESEARCH PAPER Flow instability rediction using finite elements Toktam Mohammadnejad José E. Andrade Received: 14 November 213 / Acceted: 6 June 214
More informationFinite Element Analysis of Concrete Structures Using PlasticDamage Model in 3-D Implementation
Finite Element Analysis of Concrete Structures Using PlasticDamage Model in 3-D Imlementation O. Omidi1 and V. Lotfi,* Received: October 009 Acceted: August 010 Abstract: Neither damage mechanics model
More informationSAMARIS. Draft report
SAMARIS Work Package 5 - Performance based secifications Selection and evaluation of models for rediction of ermanent deformations of unbound granular materials in road avements Draft reort Pierre Hornych
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationOn the elasticity of transverse isotropic soft tissues (L)
J_ID: JAS DOI: 10.1121/1.3559681 Date: 17-March-11 Stage: Page: 1 Total Pages: 5 ID: 3b2server Time: 12:24 I Path: //xinchnasjn/aip/3b2/jas#/vol00000/110099/appfile/ai-jas#110099 1 2 3 4 5 6 7 AQ18 9 10
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationEngineering Fracture Mechanics
Engineering Fracture Mechanics 90 (2012) 172 179 Contents lists available at SciVerse ScienceDirect Engineering Fracture Mechanics journal homeage: www.elsevier.com/locate/engfracmech Discussion of: Mode
More informationREFINED STRAIN ENERGY OF THE SHELL
REFINED STRAIN ENERGY OF THE SHELL Ryszard A. Walentyński Deartment of Building Structures Theory, Silesian University of Technology, Gliwice, PL44-11, Poland ABSTRACT The aer rovides information on evaluation
More informationGeo-E2010 Advanced Soil Mechanics L Wojciech Sołowski. 07 March 2017
Geo-E2010 Advanced Soil Mechanics L Wojciech Sołowski 07 March 2017 Soil modeling: critical state soil mechanics and Modified Cam Clay model Outline 1. Refresh of the theory of lasticity 2. Critical state
More informationMECHANICAL BEHAVIOR AND CONSTITUTIVE MODELING OF METAL CORES
JOURNL OF MECHNICS OF MTERILS ND STRUCTURES Vol., No., MECHNICL EHVIOR ND CONSTITUTIVE MODELING OF METL CORES SHKN VZIRI ND ZHENYU XUE Studying the mechanical behavior of metal cores rovides insight into
More informationImplementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis
CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty
More informationNumerical and experimental investigation on shot-peening induced deformation. Application to sheet metal forming.
Coyright JCPDS-International Centre for Diffraction Data 29 ISSN 197-2 511 Numerical and exerimental investigation on shot-eening induced deformation. Alication to sheet metal forming. Florent Cochennec
More informationA note on variational representation for singular values of matrix
Alied Mathematics and Comutation 43 (2003) 559 563 www.elsevier.com/locate/amc A note on variational reresentation for singular values of matrix Zhi-Hao Cao *, Li-Hong Feng Deartment of Mathematics and
More informationConference Paper A Finite Element Approach for the Elastic-Plastic Behavior of a Steel Pipe Used to Transport Natural Gas
Conference Paers in Energy, Article ID 267095, 10 ages htt://dx.doi.org/10.1155/2013/267095 Conference Paer A Finite Element Aroach for the Elastic-Plastic Behavior of a Steel Pie Used to Transort Natural
More information97.398*, Physical Electronics, Lecture 8. Diode Operation
97.398*, Physical Electronics, Lecture 8 Diode Oeration Lecture Outline Have looked at basic diode rocessing and structures Goal is now to understand and model the behavior of the device under bias First
More informationFDTD solutions for the distribution of radiation from dipoles embedded in dielectric particles
Journal of Quantitative Sectroscoy & Radiative Transfer 106 (007) 57 61 www.elsevier.com/locate/jqsrt FDTD solutions for the distribution of radiation from dioles embedded in dielectric articles Changhui
More informationANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE
ANALYTICAL MODEL FOR THE BYPASS ALE IN A LOOP HEAT PIPE Michel Seetjens & Camilo Rindt Laboratory for Energy Technology Mechanical Engineering Deartment Eindhoven University of Technology The Netherlands
More informationNonlinear Static Analysis of Cable Net Structures by Using Newton-Raphson Method
Nonlinear Static Analysis of Cable Net Structures by Using Newton-Rahson Method Sayed Mahdi Hazheer Deartment of Civil Engineering University Selangor (UNISEL) Selangor, Malaysia hazheer.ma@gmail.com Abstract
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More informationVarious Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various Families of R n Norms and Some Open Problems
Int. J. Oen Problems Comt. Math., Vol. 3, No. 2, June 2010 ISSN 1998-6262; Coyright c ICSRS Publication, 2010 www.i-csrs.org Various Proofs for the Decrease Monotonicity of the Schatten s Power Norm, Various
More informationA NUMERICAL MODEL OF FINITE DIFFERENCE (F.D) FOR DYNAMIC PILE DRIVING
A NUMERICAL MODEL OF FINITE DIFFERENCE (F.D) FOR DYNAMIC PILE DRIVING *Amir Taban 1, Masoud Mirmohammad Sadeghi 2 and Mohammad Ali Roshan Zamir 3 1 Deartment of Civil Engineering, Isfahan (Khorasgan) Branch,
More informationInternational Journalof Fatigue
International Journal of Fatigue 29 (2007) 1634 1646 International Journalof Fatigue www.elsevier.com/locate/ijfatigue Variable amlitude fatigue crack growth, exerimental results and modeling R. Hamam
More informationFlying characteristics in the free molecular region (influence of accommodation coefficients)
Microsyst Technol (25) 11: 85 811 DOI 1.17/s542-5-538- TECHNICAL PAPER Shigehisa Fukui Æ Hidekazu Shimada Æ Kiyomi Yamane Hiroshige Matsuoka Flying characteristics in the free molecular region (influence
More informationA PIEZOELECTRIC BERNOULLI-EULER BEAM THEORY CONSIDERING MODERATELY CONDUCTIVE AND INDUCTIVE ELECTRODES
Proceedings of the 6th International Conference on Mechanics and Materials in Design, Editors: J.F. Silva Gomes & S.A. Meguid, P.Delgada/Aores, 26-30 July 2015 PAPER REF: 5513 A PIEZOELECTRIC BERNOULLI-EULER
More informationIdentification of the source of the thermoelastic response from orthotropic laminated composites
Identification of the source of the thermoelastic resonse from orthotroic laminated comosites S. Sambasivam, S. Quinn and J.M. Dulieu-Barton School of Engineering Sciences, University of Southamton, Highfield,
More informationThe. Consortium. Continuum Mechanics. Original notes by Professor Mike Gunn, South Bank University, London, UK Produced by the CRISP Consortium Ltd
The C R I S P Consortium Continuum Mechanics Original notes b Professor Mike Gunn, South Bank Universit, London, UK Produced b the CRISP Consortium Ltd THOR OF STRSSS In a three dimensional loaded bod,
More information1 University of Edinburgh, 2 British Geological Survey, 3 China University of Petroleum
Estimation of fluid mobility from frequency deendent azimuthal AVO a synthetic model study Yingrui Ren 1*, Xiaoyang Wu 2, Mark Chaman 1 and Xiangyang Li 2,3 1 University of Edinburgh, 2 British Geological
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationSpatiotemporal Gradient Analysis of Differential Microphone Arrays*
Satiotemoral Gradient Analysis of Differential Microhone Arrays* MIHAILO KOLUNDŽIJA, CHRISTOF FALLER, AND MARTIN VETTERLI** (mihailo.kolundzija@efl.ch) (christof.faller@illusonic.com) (martin.vetterli@efl.ch)
More informationNumerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model
Numerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model Q. Sun, Y.J. Liu, R.H, Jin School of Aeronautics, Northwestern Polytechnical University, Xi an 710072,
More informationSHAPE OPTOMIZATION OF H-BEAM FLANGE FOR MAXIMUM PLASTIC ENERGY DISSIPATION
The Fourth China-Jaan-Korea Joint Symosium on Otimization of Structural and Mechanical Systems Kunming, Nov. 6-9, 2006, China SHAPE OPTOMIZATION OF H-BEAM FLANGE FOR MAXIMUM PLASTIC ENERGY DISSIPATION
More informationAn elasto-plastic model to describe the undrained cyclic behavior of saturated sand with initial static shear
University of Wollongong Research Online Faculty of Engineering - Paers (Archive) Faculty of Engineering and Information Sciences 211 An elasto-lastic model to describe the undrained cyclic behavior of
More informationDomain Dynamics in a Ferroelastic Epilayer on a Paraelastic Substrate
Y. F. Gao Z. Suo Mechanical and Aerosace Engineering Deartment and Princeton Materials Institute, Princeton University, Princeton, NJ 08544 Domain Dynamics in a Ferroelastic Eilayer on a Paraelastic Substrate
More informationThe Properties of Pure Diagonal Bilinear Models
American Journal of Mathematics and Statistics 016, 6(4): 139-144 DOI: 10.593/j.ajms.0160604.01 The roerties of ure Diagonal Bilinear Models C. O. Omekara Deartment of Statistics, Michael Okara University
More informationDesign of Isolated Bridges from the Viewpoint of Collapse under Extreme Earthquakes
Design of Isolated Bridges from the Viewoint of Collase under Extreme Earthquakes D.W. Chang, Y.T. Lin, C.H. Peng, C.Y. Liou CECI Engineering Consultants, Inc., Taiwan T.Y. Lee National Central University,
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Jul 1998
arxiv:cond-mat/98773v1 [cond-mat.stat-mech] 5 Jul 1998 Floy modes and the free energy: Rigidity and connectivity ercolation on Bethe Lattices P.M. Duxbury, D.J. Jacobs, M.F. Thore Deartment of Physics
More informationINTRODUCING THE SHEAR-CAP MATERIAL CRITERION TO AN ICE RUBBLE LOAD MODEL
Symosium on Ice (26) INTRODUCING THE SHEAR-CAP MATERIAL CRITERION TO AN ICE RUBBLE LOAD MODEL Mohamed O. ElSeify and Thomas G. Brown University of Calgary, Calgary, Canada ABSTRACT Current ice rubble load
More informationDesign of NARMA L-2 Control of Nonlinear Inverted Pendulum
International Research Journal of Alied and Basic Sciences 016 Available online at www.irjabs.com ISSN 51-838X / Vol, 10 (6): 679-684 Science Exlorer Publications Design of NARMA L- Control of Nonlinear
More informationSELF-SIMILAR FLOW UNDER THE ACTION OF MONOCHROMATIC RADIATION BEHIND A STRONG CYLINDRICAL SHOCK WAVE IN A NON-IDEAL GAS
SELF-SIMILAR FLOW UNDER THE ACTION OF MONOCHROMATIC RADIATION BEHIND A STRONG CYLINDRICAL SHOCK WAVE IN A NON-IDEAL GAS *J. P. Vishwakarma and Vijay Kumar Pandey Deartment of Mathematics & Statistics,
More informationA Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition
A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino
More informationTemperature, current and doping dependence of non-ideality factor for pnp and npn punch-through structures
Indian Journal of Pure & Alied Physics Vol. 44, December 2006,. 953-958 Temerature, current and doing deendence of non-ideality factor for n and nn unch-through structures Khurshed Ahmad Shah & S S Islam
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationCREEP AND SHRINKAGE EFFECT ON THE DYNAMIC ANALYSIS OF REINFORCED CONCRETE SLAB-AND-BEAM STRUCTURES
ECCM 99 Euroean Conference on Comutational Mechanics August 31 Setember 3 München, Germany CREEP AND SHRINKAGE EFFECT ON THE DYNAMIC ANALYSIS OF REINFORCED CONCRETE SLABANDBEAM STRUCTURES Evangelos J.
More informationModeling Volume Changes in Porous Electrodes
Journal of The Electrochemical Society, 53 A79-A86 2006 003-465/2005/53/A79/8/$20.00 The Electrochemical Society, Inc. Modeling olume Changes in Porous Electrodes Parthasarathy M. Gomadam*,a,z John W.
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationTowards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK
Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)
More informationA constitutive model for spring-back prediction in which the change of Young s modulus with plastic deformation is considered
International Journal of Machine Tools & Manufacture 47 (27) 1791 1797 Short communication A constitutive model for spring-back prediction in which the change of Young s modulus with plastic deformation
More informationApplication of Automated Ball Indentation for Property Measurement of Degraded Zr2.5Nb
Journal of Minerals & Materials Characterization & Engineering, Vol. 10, No.7,.661-669, 011 jmmce.org Printed in the USA. All rights reserved Alication of Automated Ball Indentation for Proerty Measurement
More informationINFLUENCE OF POROSITY AND MOISTURE ON MECHANICAL PROPERTIES OF CERAMIC BRICK ANALYTICAL HOMOGENISATION APPROACH
53 Proc. Int. Sym. Brittle Matrix Comosites 11 A.M.Brandt, J.Olek, M.A.Glinicki, C.K.Y.Leung, J.Lis, eds. Warsaw, Setember 28-30, 2015 Institute of Fundamental Technological Research INFLUNC OF POROSITY
More informationModelling a Partly Filled Road Tanker during an Emergency Braking
Proceedings of the World Congress on Engineering and Comuter Science 217 Vol II, October 25-27, 217, San Francisco, USA Modelling a Partly Filled Road Tanker during an Emergency Braking Frank Otremba,
More informationResponse of a sandwich panel subject to fire or elevated temperature on one of the surfaces
Comosites: Part A 37 (26) 981 988 www.elsevier.com/locate/comositesa Resonse of a sandwich anel subject to fire or elevated temerature on one of the surfaces V. Birman a, *, G.A. Kardomateas b, G.J. Simitses
More informationChapter 5. Transient Conduction. Islamic Azad University
Chater 5 Transient Conduction Islamic Azad University Karaj Branch 1 Transient Conduction Many heat transfer roblems are time deendent Changes in oerating conditions in a system cause temerature variation
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
More informationMeshless Methods for Scientific Computing Final Project
Meshless Methods for Scientific Comuting Final Project D0051008 洪啟耀 Introduction Floating island becomes an imortant study in recent years, because the lands we can use are limit, so eole start thinking
More informationMODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
More informationCOMPARISON OF VARIOUS OPTIMIZATION TECHNIQUES FOR DESIGN FIR DIGITAL FILTERS
NCCI 1 -National Conference on Comutational Instrumentation CSIO Chandigarh, INDIA, 19- March 1 COMPARISON OF VARIOUS OPIMIZAION ECHNIQUES FOR DESIGN FIR DIGIAL FILERS Amanjeet Panghal 1, Nitin Mittal,Devender
More informationAnalysis of cold rolling a more accurate method
Analysis of cold rolling a more accurate method 1.1 Rolling of stri more accurate slab analysis The revious lecture considered an aroximate analysis of the stri rolling. However, the deformation zone in
More informationNumerical Methods for Particle Tracing in Vector Fields
On-Line Visualization Notes Numerical Methods for Particle Tracing in Vector Fields Kenneth I. Joy Visualization and Grahics Research Laboratory Deartment of Comuter Science University of California, Davis
More informationTheory of magnetic oscillations in Josephson-junction stacks
Physica C 437 438 (2006) 157 161 www.elsevier.com/locate/hysc Theory of magnetic oscillations in Josehson-junction stacks A.E. Koshelev Materials Science Division, Argonne National Laboratory, Argonne,
More informationClosed-form solution for the size of plastic zone in an edge-cracked strip
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Publications from the Deartment of Engineering Mechanics Mechanical & Materials Engineering, Deartment of May 2002
More informationWeek 8 lectures. ρ t +u ρ+ρ u = 0. where µ and λ are viscosity and second viscosity coefficients, respectively and S is the strain tensor:
Week 8 lectures. Equations for motion of fluid without incomressible assumtions Recall from week notes, the equations for conservation of mass and momentum, derived generally without any incomressibility
More informationGuidelines for Constructing Strain Gradient Plasticity Theories
. A. Fleck Deartment of Engineering, Cambridge University, Trumington Street, Cambridge CB2 1Z, UK J. W. Hutchinson School of Engineering and Alied Sciences, Harvard University, Cambridge, MA 2138 J. R.
More informationUB Electronic Packaging Laboratory, Department of Civil Engineering, 102 Ketter Hall, University at Buffalo, Buffalo, NY 14260, USA
Mechanics of Materials 38 (2006) 585 598 www.elsevier.com/locate/mechmat Damage mechanics constitutive model for Pb/Sn solder joints incorporating nonlinear kinematic hardening and rate dependent effects
More informationarxiv: v1 [cond-mat.stat-mech] 1 Dec 2017
Suscetibility Proagation by Using Diagonal Consistency Muneki Yasuda Graduate School of Science Engineering, Yamagata University Kazuyuki Tanaka Graduate School of Information Sciences, Tohoku University
More informationThree-dimensional modeling of strain localization in Cosserat continuum theory
hree-dimensional modeling of strain localization in Cosserat continuum theory A.R. Khoei, S. Yadegari and M. Anahid Downloaded from ijce.iust.ac.ir at 1:42 IRS on Saturday October 13th 2018 Center of Excellence
More informationSimple extension of Lemaitre's elastoplastic damage model to account for hydrolytic degradation
884 Simle extension of Lemaie's elastolastic damage model to account for hydrolytic degradation Absact Being motivated by the technological alications of bioabsorbable olymeric materials in the fields
More informationVIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES
Journal of Sound and Vibration (998) 22(5), 78 85 VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES Acoustics and Dynamics Laboratory, Deartment of Mechanical Engineering, The
More informationFEM simulation of a crack propagation in a round bar under combined tension and torsion fatigue loading
FEM simulation of a crack roagation in a round bar under combined tension and torsion fatigue loading R.Citarella, M.Leore Det. of Industrial Engineering University of Salerno - Fisciano (SA), Italy. rcitarella@unisa.it
More informationFlexible Pipes in Trenches with Stiff Clay Walls
Flexible Pies in Trenches with Stiff Clay Walls D. A. Cameron University of South Australia, South Australia, Australia J. P. Carter University of Sydney, New South Wales, Australia Keywords: flexible
More informationKEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS
4 th International Conference on Earthquake Geotechnical Engineering June 2-28, 27 KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS Misko CUBRINOVSKI 1, Hayden BOWEN 1 ABSTRACT Two methods for analysis
More informationCONCRETE MATERIAL MODELING IN EXPLICIT COMPUTATIONS ABSTRACT
Worksho on Recent Advances in Comutational Structural Dynamics and High Performance Comuting USAE Waterways Exeriment Station Aril 24-26 1996 CONCRETE MATERIAL MODELING IN EXPLICIT COMPUTATIONS L. Javier
More information