Partitioned formulation of frictional contact problems using localized Lagrange multipliers
|
|
- Rosaline Wells
- 6 years ago
- Views:
Transcription
1 COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2000; 00:1 6 Preared using cnmauth.cls [Version: 2002/09/18 v1.02] Partitioned formulation of frictional contact roblems using localized Lagrange multiliers José A. González 1, K. C. Park 2,, and Carlos A. Felia 2 1 Escuela Suerior de Ingenieros de Sevilla, Avda. de los Descubrimientos s/n, Sevilla, Sain 2 Deartment of Aerosace Engineering Sciences, Center for Aerosace Structures, University of Colorado, Boulder, CO , USA SUMMARY Interface treatment methods for the contact roblem between non-matching meshes have traditionally been based on a direct couling of the contacting solids emloying a master-slave strategy or classical Lagrange multiliers. These methods tend to generate strongly couled systems that is deendent on the discretization characteristics on each side of the contact zone. In this work a dislacement frame is intercalated between the interface meshes. The frame is then discretized so that the discrete frame nodes are connected to the contacting substructures using the localized Lagrange multiliers collocated at the interface nodes. The resulting methodology alleviates the need for master-slave book-keeing and rovides a artitioned formulation which reserves software modularity, facilitates non-matching mesh treatment and asses the contact atch test. Frictional contact roblems are used to demonstrate the salient features of the roosed method. Coyright c 2000 John Wiley & Sons, Ltd. key words: contact; friction; finite element method; localized Lagrange multiliers 1. INTRODUCTION In comutational structural mechanics, roblems resenting frictional contact surfaces combined with non-matching meshes remain as a difficult task. The source of this difficulty comes not only from the strong non-linearity of the frictional contact law, involving multivalued relationshis between kinematic and static variables, but also from the severe discontinuity forced by the different meshes used to model the contact interface. Most existing interface treatment methods for the contact roblem between non-matching meshes have traditionally been based on a direct couling of the contacting solids emloying a masterslave strategy, classical Lagrange multiliers or mortar-like methods [1, 2, 3, 4]. Corresondence to: K. C. Park, Center for Aerosace Structures and Deartment of Aerosace Engineering Sciences, University of Colorado at Boulder, Camus Box 429, Boulder, CO , U.S.A. kcark@colorado.edu Contract/grant sonsor: Consejería de Educación y Ciencia, Junta de Andalucía, Sain; contract/grant number: Convocatoria 2004 de concesión de ayudas ara el erfeccionamiento de investigadores en centros de investigación de fuera de Andalucía Received 17 January 2000 Coyright c 2000 John Wiley & Sons, Ltd. Revised 18 Setember 2002
2 2 J. A. GONZÁLEZ, K. C. PARK AND C. A. FELIPPA The resent work is an extension of the formulation roosed by Rebel, Park and Felia [5] to solve the contact roblem with non-matching meshes introducing an intermediate contact surface, or contact frame, endowed with indeendent degrees of freedom and treated with a FEM discretization. This frame is connected to the contacting substructures using localized Lagrange multiliers [6] and can be constructed in order to reserve the constant-stress interface atch test [7]. In the work of Rebel, Park and Felia a general formulation of the roblem is derived using the variational framework roosed by Park and Felia [8, 9] and then is articularized to the two dimensional case. However, they considered the contact zone to be a riori known after alying a two stages redictor-corrector algorithm and decided the contact oint states using a trial and error based algorithm. In the resent work, we extend [5] so that the contact frame may searate from the substructures in all directions, that is, the frame now is a comletely free system that moves while maintaining the contact zone as an unknown. Another new feature is in the way of finding contact states; in the resent work contact conditions are imosed mathematically using the augmented Lagrangian formulation and rojection functions, making the contact search to be art of the overall contact algorithm. 2. THE CONTACT FRAME Let us suose two solids in contact and denote their domains Ω and Ω. To formulate the contact roblem, instead of considering the direct interaction between the two bodies during the contact rocess, we will insert a deforming non-hysical surface Λ between them and reformulate the contact roblem in terms of interaction of the two solids with this new element, called the contact frame, using localized Lagrange multiliers on each side of the frame. The contact tractions acting on the frame are reresented in the exloded view of Figure 1 where the localized Lagrange multiliers connecting solid Ω with the frame are named using the vector quantity λ = (λ n, λ t1, λ t2 ) t and the multiliers connecting solid Ω are named λ = ( λ n, λ t1, λ t2 ) t. These forces are exressed using two locally orthonormal base systems connected to the frame; B = [n, a 1, a 2 ] used to describe λ and B = [ n, ā 1, ā 2 ] used for λ. These local systems for describing the frame kinematics are defined in the following way; a 1 and a 2 are the orthogonal vectors contained in the frame tangent lane at the considered oint and vector n oints towards solid Ω. The barred base system B at the same osition, will be defined in oosite direction to B. In order to describe the motion of the two solids we use the dislacement fields u and ū defined on Ω and Ω resectively, which, when added to the reference configurations X and X, rovide us with the current ositions x and x: x = X + u x = X + ū The motion of the contact frame is described using its dislacements v from its initial configuration Y, roviding a current osition (1) y = Y + v (2) however, this motion will be restricted to ermanently maintain the frame just in the middle between the two contact interfaces. To do that, the relative frame-contact interface distance Preared using cnmauth.cls
3 PARTITIONED FORMULATION OF FRICTIONAL CONTACT PROBLEMS 3 Ω Ω Contact Zone Γ C λ λ n Λ a 2 Ω Contact Frame Λ a 1 a 2 λ n a 1 Ω λ Figure 1. Left: Contact interfaces with intercalated frame. Rigth: Exloded view with localized Lagrange multiliers. vectors k = (k n, k t1, k t2 ) t and k = ( k n, k t1, k t2 ) t are exressed in a local coordinate system of the frame by k = B t (x y) k = B t ( x y) k k (3) = 0 Previous definitions[5] make the contact frame to coincide with the contact zone where k n is equal to zero. Also k t will give us the direction of the relative motion of the contact interfaces, a variable needed to formulate the frictional behavior. 3. ENFORCEMENT OF THE CONTACT CONDITIONS The behavior of the contact interface is governed by the non-enetration condition and the Coulomb friction law. To formulate the contact conditions, let s introduce the augmented Lagrange multilier variable λ(r) = λ + rk with a enalty arameter r > 0, and the Coulomb disk of radius g rojection oerator P Cg (.) : R 2 R 2 { [x, y] t if x 2 + y 2 g 2 P Cg (x, y) = g [x, y]t otherwise (4) 2 x2 +y then, the comlete fulfilment of contact conditions can be assured using the cone rojection oerator P (.) : R 3 R 3 alied to the augmented Lagrangian multiliers in the following Preared using cnmauth.cls
4 4 J. A. GONZÁLEZ, K. C. PARK AND C. A. FELIPPA way [ P (λ(r)) = max(0, λ n (r)) P Cµ max(0,λn) (λ t (r)) This definition allows us to exress the contact conditions with the final exression ] (5) P (λ(r)) = λ (6) 4. VARIATIONAL FORMULATION To derive the equilibrium equations of the constrained system we use the variational formulation roosed by Park and Felia [8, 9] where the roblem is treated as if all bodies were entirely free, formulating the virtual work by summing u the contributions of each body. The localized constraint conditions are constructed and then multilied by unknown coefficients, denoted herein as the localized Lagrange multiliers. The resulting constraint functional is added to the virtual work of the free substructures to yield the total work of the system. The variational functional that reresents the total energy of the system δπ, is then comosed by the energy of the two comletely free substructures lus the interface constraint functional associated with the contact henomena where the contact interface otential δπ total i δπ = δπ free + δ π free + δπ total i (7) consists of contributions from both sides δπ total i = δπ i + δ π i (8) which will be derived in this section. To do that, let us decomose each one of the two interface functionals into two terms δπ i = δπ k + δπ c δ π i = δ π k + δ π c (9) the first one is related with the kinematic ositioning of the frame, equation (3), that is enforced in a weak sense using the variational form δπ k = (δ{λ [B t (x y) k]})d (10) and the second one reresents the virtual work of the contact forces, contribution to the weak form that can be exressed in the following way δπ c = (λ δk)d (11) where the contact forces λ have to satisfy the unilateral contact law and the frictional law. These restrictions can be automatically satisfied relacing the Lagrange multiliers by the rojection oerator (5) obtaining δπ c = (P (λ (r)) δk)d (12) Preared using cnmauth.cls
5 PARTITIONED FORMULATION OF FRICTIONAL CONTACT PROBLEMS 5 Substituting (12) and (10) into the first of (9) rovides the exression for the total variation of the interface otential at the non barred side δπ i = (λ δ{b t (x y)} + δλ {B t (x y) k} + δk { λ + P (λ (r))})d (13) Similarly, one can obtain δ π i as δ π i = ( λ δ{ B t ( x y)} + δ λ { B t ) ( x y) k} + δk { λ + P ( λ (r) })dγc (14) for the barred side. To deal with the first terms of equations (13) and (14) we decomose them in the following way λ δ{b t (x y) = (δu δv) {Bλ} + λ n δn (x y) + λ tα δa α (x y) (15) λ δ{ B t ( x y) = (δū δv) { B λ} + λ n δ n ( x y) + λ tα δā α ( x y) (16) where the variation of the normal and tangential frame unitary vectors can be written, see [5], as δa α = Q α δv,α, δā α = Q α δv,α (17) δn = e βα (a β I)Q α δv,α, δ n = e βα (ā β I) Q α δv,α (18) Substituting into (13) and (14) leads to the final exression for the interface otentials: δπ i = (δλ {B t (x y) k} + (δu δv) {Bλ} + δv,α {Q α Φ α (x y)} + δk { λ + P (λ (r))})d (19) with Φ α = λ tα I + λ n e αβ (a β I), and for the barred: side δ π i = (δ λ { B t ( x y) k} + (δū δv) { B λ} + δv,α { Q α Φ α ( x y)} + δk { λ + P ( λ (r) ) })d (20) with Φ α = λ tα I λ n e αβ (a β I). It is imortant to mention that receding interface constraint functional will not lead to a standard minimization roblem, for the Coulomb disk C g inside the cone rojection oerator (5) is a function of the normal contact through the friction limit g which deends on the solution u. To overcome this difficulty and following to Alart and Curnier [10] we have obtained a articular form of quasi-variational functional by substituting g by the convex set µ max(0, λ n ); for this reason, the minimization roblem associated with (7) is considered as a quasi-variational roblem. 5. DISCRETE EQUATIONS The discrete contact roblem will be defined in terms of contact airs, i.e. coules formed by a set of interface nodes and its associated frame element. Those contact airs are established Preared using cnmauth.cls
6 6 J. A. GONZÁLEZ, K. C. PARK AND C. A. FELIPPA Solid Ω ( v, k ) Λ ( v, k ) ξ a 2 1 ξ 1 ξ n ( u, λ ) 2 a 2 3 ( v, k ) 2 2 ( v, k ) 3 3 Figure 2. Contact airs between the barred contact interface node and the frame. Contact airs between the non-barred contact interface node and the frame element is indeendently defined. before starting each time ste, calculating for every otentially contacting boundary node its nearest frame element and rojecting geometrically on it, see Figure 2. The frame dislacements v(ξ 1, ξ 2 ) and its average distance to the solids k(ξ 1, ξ 2 ) are interolated using isoarametric finite elements in the following form v(ξ 1, ξ 2 ) = N(ξ 1, ξ 2 ) v 1 :, k(ξ 1, ξ 2 ) = N(ξ 1, ξ 2 ) k 1 : (21) where n f is the number of nodes in the frame element, and N N nf 0 0 N(ξ 1, ξ 2 ) = 0 N N nf 0 v nf k nf (22) 0 0 N N nf is the shae functions aroximation matrix. The localized Lagrange multiliers are collocated with the contacting interface nodes so that they can be exressed in terms of Dirac s delta functions λ(ξ 1, ξ 2 ) = λ δ(ξ ξ ) (23) with ξ = (ξ 1, ξ 2 ) and ξ the frame coordinates of the node rojection, Figure 2. This definition of contact forces will reduce integrations over the contact zone to summations over the contact airs, i.e. n λ(ξ 1, ξ 2 ) f(ξ 1, ξ 2 )d = λ f(ξ ) (24) Preared using cnmauth.cls =1
7 PARTITIONED FORMULATION OF FRICTIONAL CONTACT PROBLEMS 7 with n the total number of contact airs. Exression (24) is useful in order to maintain the contact interface generic, leading to modular couling software. To manage the interface discrete variables we introduce the substructural interface nodal indicator L, the well known Boolean finite element assembling oerator defined in the following way u = L u u, λ = L λ λ (25) v = L v v, k = L k k where L is used to extract the variable associated with a contact interface node from the global unknowns vector. Using the receding definitions, the discrete variational form of the total energy functional can be finally written δπ = δu f + δū f + δλ g + δ λ ḡ δv (q + q) + δk ( + ) (26) where f = n =1 n =1 L t uf, g = n =1 n L t λ g q = L t vq, = L t k =1 with the corresonding exressions for the barred quantities. For each contact airs, the following equations hold (27) f = {f int f ext + Bλ}, f = { f int f ext + B λ} (28) g = {B t (x y) Nk}, ḡ = { B t ( x y) Nk} (29) q = { N t Bλ + N ț αq α Φ α (x y)}, q = { N t B λ + N ț α Q α Φ α ( x y)} (30) = {N t ( λ + P (λ(r)))}, = { N t ( λ + P ( λ(r) ) )} (31) where equation (28) reresents the equilibrium equation of each substructure (term Bλ are the contact forces exressed in the global system), equation (29) governs the relative motion between the frame and the substructures, equation (30) evaluates the forces acting on the frame and equation (31) imoses the fulfilment of the contact conditions. The stationary of the above variational equation, viz., δπ = 0, is obtained by solving for the unknown z = ( u, ū, λ, λ, v, k ) the following B-differentiable system F(z) = f f g ḡ q q + = 0 (32) where non F-differentiability occurs because equation + can resent non-linear directional derivatives. The artitioned equations of motion for the frame-based contact roblem are obtained from the Jacobian of system (32), that can be exressed as Preared using cnmauth.cls
8 8 J. A. GONZÁLEZ, K. C. PARK AND C. A. FELIPPA K 0 B 0 L v 0 u f 0 K 0 B L v 0 ū f B t L b N k λ 0 B t 0 0 L b N k λ = g ḡ L t v L t v L t b L t b D v D v 0 v q + q 0 0 P λ P λ 0 P k + P k k (33) where only non-differentiable terms P λ, P λ, P k and P k need a secial treatment. For examle, comonents of the receding equation for the non barred quantities are obtained by assembling contributions of each contact airs : B = n =1 n =1 L t ub L λ, L v = n =1 n L t ul v L v L b = L t λ L bl v, N k = L t λ N kl k P λ = n =1 D v = n =1 =1 L t vd v L v L t k P λl λ, P k = n =1 L t k P kl k (34) with sums extended to the n contact airs from that side and where each term comes from differentiation of equations (28), (29), (30) and (31) L v = { Φ t αq α } N,α (35) (x y) t (a 2 I) (x y) t (a 1 I) L b = Bt N (x y) t Q 1 N,1 0 Q 2 N,2 0 (x y) t D v = +N ț α (36) N k = {N} (37) N t [Φ t αq α ]N,α + N ț α[q α Φ α ]N 1 y,α {Qα Φ α (x y) a α + a α Q α Φ α (x y) +[Φ α (x y)] a α Q α }N,α } Exressions for the barred side can be similarly obtained. On the other hand, terms P λ and P λ are decomosed into their normal and tangential arts: P λ = P λn + P λt and P k = P kn + P λt (39) with the following definition for the normal art: (38) 1. Case λ n (r) < 0 P λn = Nt Preared using cnmauth.cls
9 PARTITIONED FORMULATION OF FRICTIONAL CONTACT PROBLEMS 9 2. Case λ n (r) 0 P kn = Nt r N In the receding equations one can see that instead of calculating a comlicated non-linear directional derivative for the non-differentiable case λ n (r) = 0, it has been relaced by a more simle linear derivative coming from the right side. A similar technique is used for nondifferentiable oints of the tangential contribution: 1. Case λ n 0 2. Case λ n > 0 (a) When λ t (r) µλ n P λt = Nt P kt = Nt r 0 N 0 0 r (b) otherwise, if λ t (r) > µλ n with α = P λt = P kt = µλ n, β = λ 2 t1 (r)+λ 2 t (r) 2 and R = λ t (r) λ t (r). Nt Nt α λt 1 (r) λ n Ψ 11 Ψ 12 α λt 2 (r) λ n Ψ 21 Ψ rθ 11 rθ 12 N 0 rθ 21 rθ 22 µλ n (λ 2 t 1 (r)+λ 2 t 2 (r)) 3 2, Ψ = (α 1)I βr, Θ = αi βr In the above exressions we are using a linearized substitute for the non-linear directional derivative aearing when any of these equalities hold λ n = 0, λ n (r) = 0 or λ t (r) = µλ n, although the function can be exected to be normally differentiable in the large majority of ractical cases. 6. SOLVING THE NON-LINEAR SYSTEM To solve system (32) the Generalized Newton s Method with Line Search (GNMLS) has been used. GNMLS is an effective extension of the Newton s Method for B-differentiable functions roosed by Pang [11] in a general context and articularized by Alart [12] and Christensen [13] for the contact case. This method is based on the comutation of the non-linear directional Preared using cnmauth.cls
10 10 J. A. GONZÁLEZ, K. C. PARK AND C. A. FELIPPA derivative of the objective function; however, it is well known that in contact roblems this non-linear directional derivative rarely needs to be comuted and can be substituted by a linearized version without affecting the algorithm convergence; this is the aroach adoted herein in obtaining (39). If we define the scalars β (0, 1), σ (0, 1/2), and ε > 0 but small, the alication of GNMLS algorithm to solve the non-linear equations F(z) = 0, can be summarized in the following stes: 1. Time integration t: Solve for z t+ t with the known state variable vector z t. (a) Inner GNMLS iterations, loo k. (b) Use GMRES to solve for z t k in the system F(zt k ; zt k ) = F(zt k ). (c) Obtain first integer m = 1, 2,... that fulfills H(z t k + βm z t k ) (1 2σβm )H(z t k ) with H(z) = 1 2 Ft (z)f(z). (d) Udate the solution z t k+1 = zt k + τ k z t k with τ k = β m. (e) If H(z t k+1 ) ɛ continue, else comute new GNMLS iteration k k Make z t+ t = z t k+1 and solve for next time ste t t + t. The algorithm uses two nested loos, the external one t cares for time marching and the internal one k does the Newton subiterations. On each one of this subiterations a linear system has to be solved to comute the search direction, where we use a sarse matrix storage scheme combined with the GMRES solver. When a direction z t k is obtained, it is scaled by a factor of τ k obtained from a decreasing error condition. 7. APPLICATIONS The roblem considered is the Hertzian contact of two geometrically identical elastic sheres of radius R 1 = R 2 = 8 mm indented alying an external load P = 1 N normal to the contact zone. The material roerties for the first shere are E 1 = 100 GPa, ν 1 = 0.2 and E 2 = 100 GPa, ν 2 = 0.4 for the second. This difference in the Poisson s ratio, together with a friction coefficient µ = 0.25 will roduce a distribution of tangential stresses in the contact zone that will not affect considerably the normal ressure, which will be very close to the Hertz solution that redicts a contact zone radius of mm. The meshes for both sheres are comletely identical, generated using 1280 hexahedral elements and 1569 nodes for each one of them, and refining the otentially contact zone (a square of 1mm 2 ) with 100 quadrilateral elements and 121 contact nodes, see Figure 3. The contact frame is a lane of 1mm 2 that uses the same discretization as the otential contact zone and the normal load is alied using ten time stes.it should be noted that because Poisson s ratios for the two contacting sheres are different, the initially matching interface meshes will lead to non-matching meshes. The convergence rate for one ste solution is resented in Figure 4 in terms of relative error H(z k )/H(z 1 ) for each GNMLS subiteration k together with the scaling arameter τ k. On the right side of the same Figure is shown the normal distance from the frame to the sheres in the deformed state where the contact zone is defined by k n 0 and with the mesh resolution used its radius is estimated to be around 0.35mm. Preared using cnmauth.cls
11 PARTITIONED FORMULATION OF FRICTIONAL CONTACT PROBLEMS 11 Figure 3. Mesh and frame used for the Hert roblem. -0.5mm Relative error H(z k )/H(z 1 ) E τ 1E-4 k E-5 H(z k )/H(z 1 ) 0.3 1E E E mm Iteration number (k) -0.5mm Solution advance a=0.35mm 0.5mm Disacement (mm) k n Figure 4. Convergence rate of the GNMLS algorithm and normal distance k n between the deformed sheres and the frame. Preared using cnmauth.cls
12 12 J. A. GONZÁLEZ, K. C. PARK AND C. A. FELIPPA FEM Hertz λ n λ n y=0.0 y= y= x(mm) y(mm) Contact zone X (mm) Figure 5. Distribution of normal multiliers in the contact zone and comarison with Hertz solution. λt λt y=0.0 y=0.3 y=0.2 FEM y(mm) x(mm) Contact zone X (mm) Figure 6. Distribution of tangential multiliers with µ = 0.25 and the normal load N alied with ten time stes. The normal Lagrange multiliers are resented in Figure 5 comared with the Hertz solution (right) for three cross sections along the x direction. The maximum normal Lagrange multilier λ n max =0.035 N is obtained at the center of the contact zone and the differences between the numerical results and the Hertz solution can be attributed to a coarse discretization combined with the use of linear elements. It is well known that the couling between the normal and tangential stresses for this roblem is very low, meaning that normal tractions are not considerably affected by the frictional henomena. It was observed also that the number of stes used to aly the normal load didn t modify the normal solution but affected considerably the tangential multiliers inside the central stick zone. Finally, the tangential Lagrange multiliers in x direction can be seen in Figure 6 obtaining the same solution rotated 90 o for the orthogonal direction y. It reresents a comlex shae that have to satisfy the sli condition at the external annular region making the tangential forces equal to the friction limit. Preared using cnmauth.cls
13 PARTITIONED FORMULATION OF FRICTIONAL CONTACT PROBLEMS CONCLUSION A formulation of the frictional contact roblem using localized Lagrange multiliers to connect the contacting substructures to an adatative contact frame is resented. The contact conditions are mathematically formulated using rojection functions and the contact frame is allowed to move freely between the substructures maintaining the contact zone as an unknown. The GNMLS is alied to solve the non-smooth system of equations reresenting the equations of motion and different regularization techniques to imrove convergence are roosed and tested. The algorithm and introduced formulations rove to be very robust and efficient when solving 3D non-matching contact roblems. Suggested methodology simlifies the rocedure to solve this kind of roblems, minimizing geometrical knowledge of the contacting substructures needed to formulate the contact interface behavior. This fact can facilitate further extensions for contact roblems, like iterative arallel comutations or connection between different numerical techniques like the FEM and the Boundary Element Method. REFERENCES 1. Sim, J. C., Wriggers, P. and Taylor, R. L. A erturbed lagrangian formulation for the finite element solution of contact roblems. Comuter Methods in Alied Mechanics and Engineering (50), (1985). 2. Paadooulos, P. and Taylor, R. L. A mixed formulation for the finite element solution of contact roblems. Comuter Methods in Alied Mechanics and Engineering (94), (1992). 3. Sim, J. C. and Laursen, T.A. An augmented lagragian treatment of contact roblems involving friction. Comuters & Structures 42, (1992). 4. Puso, M. A. and Laursen, T. A. A mortar segment-to-segment contact method for large deformation solid mechanics. Comuter Methods in Alied Mechanics and Engineering (193), (2004). 5. Rebel G, Park KC, Felia CA. A contact formulation based on localized Lagrange multiliers: formulation and alication to two-dimensional roblems. International Journal for Numerical Methods in Engineering 2002; 54(2): Park KC, Felia CA, Gumaste UA. A localized version of the metod of Lagrange multiliers and its alications. Comutational Mechanics 2000; 24: Park KC, Felia CA, Rebel G. A simle algorithm for localized construction of non-matching structural interfaces. International Journal for Numerical Methods in Engineering 2002; 53: Park KC, Felia CA. A variational framework for solution method develoments in structural mechanics. Journal of Alied Mechanics 1998; 65: Park KC, Felia CA. A localized version of the metod of Lagrange multiliers and its alications. Comutational Mechanics 2000; 24: Alart P, Curnier A. A mixed formulation for frictional contact roblems rone to Newton like solution methods. Com. Al. Mech. Engng. 1991; 32: Pang JS. Newton s method for B-differentiable equations. Math. Oer. Research 1990; 15(2): Alart P. Méthode de Newton généralisée en mécanique du contact. J. Math. Pures Al. 1997; 76: Christensen PW, Klarbring A, Pang JS, Strömberg N. Algorithms for frictional contact roblems based on mathematical rogramming. Linköing Univ., Sweden 1998; Preared using cnmauth.cls
Partitioned Formulation for Solving 3D Frictional Contact Problems with BEM using Localized Lagrange Multipliers
Copyright c 2007 ICCES ICCES, vol.2, no.1, pp.21-27, 2007 Partitioned Formulation for Solving 3D Frictional Contact Problems with BEM using Localized Lagrange Multipliers L. Rodríguez-Tembleque 1, J.A.
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 informationFE 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 informationCode_Aster. Connection Harlequin 3D Beam
Titre : Raccord Arlequin 3D Poutre Date : 24/07/2014 Page : 1/9 Connection Harlequin 3D Beam Summary: This document exlains the method Harlequin develoed in Code_Aster to connect a modeling continuous
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 information5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell:
5. PRESSURE AND VELOCITY SPRING 2019 5.1 The momentum equation 5.2 Pressure-velocity couling 5.3 Pressure-correction methods Summary References Examles 5.1 The Momentum Equation Each comonent of momentum
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 informationarxiv: v1 [quant-ph] 22 Apr 2017
Quaternionic Quantum Particles SERGIO GIARDINO Institute of Science and Technology, Federal University of São Paulo (Unifes) Avenida Cesare G. M. Lattes 101, 147-014 São José dos Camos, SP, Brazil arxiv:1704.06848v1
More informationCalculation of eigenvalue and eigenvector derivatives with the improved Kron s substructuring method
Structural Engineering and Mechanics, Vol. 36, No. 1 (21) 37-55 37 Calculation of eigenvalue and eigenvector derivatives with the imroved Kron s substructuring method *Yong Xia 1, Shun Weng 1a, You-Lin
More informationRecursive Estimation of the Preisach Density function for a Smart Actuator
Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationSYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY
SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in
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 informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationPreconditioning techniques for Newton s method for the incompressible Navier Stokes equations
Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College
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 informationTopology Optimization of Three Dimensional Structures under Self-weight and Inertial Forces
6 th World Congresses of Structural and Multidiscilinary Otimization Rio de Janeiro, 30 May - 03 June 2005, Brazil Toology Otimization of Three Dimensional Structures under Self-weight and Inertial Forces
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 informationIsogeometric analysis based on scaled boundary finite element method
IOP Conference Series: Materials Science and Engineering Isogeometric analysis based on scaled boundary finite element method To cite this article: Y Zhang et al IOP Conf. Ser.: Mater. Sci. Eng. 37 View
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 informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationContact Dynamics with Lagrange multipliers
Contact Dynamics with Lagrange multiliers Stehan Brunßen 1, Stefan Hüeber 1, and Barbara Wohlmuth 1 IANS, University of Stuttgart, Pfaffenwaldring 57, 70569 Stuttgart, Germany {brunssen,hueeber,wohlmuth}@ians.uni-stuttgart.de
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
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 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 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 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 informationOn Two-Finger Grasping of Deformable Planar Objects
2011 IEEE International Conference on Robotics and Automation Shanghai International Conference Center May 9-13, 2011, Shanghai, China On Two-Finger Grasing of Deformable Planar Objects Yan-Bin Jia Feng
More informationOptimization of Gear Design and Manufacture. Vilmos SIMON *
7 International Conference on Mechanical and Mechatronics Engineering (ICMME 7) ISBN: 978--6595-44- timization of Gear Design and Manufacture Vilmos SIMN * Budaest Universit of Technolog and Economics,
More information22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 13: PF Design II The Coil Solver
.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 13: PF Design II The Coil Solver Introduction 1. Let us assume that we have successfully solved the Grad Shafranov equation for a fixed boundary
More informationMODULAR LINEAR TRANSVERSE FLUX RELUCTANCE MOTORS
MODULAR LINEAR TRANSVERSE FLUX RELUCTANCE MOTORS Dan-Cristian POPA, Vasile IANCU, Loránd SZABÓ, Deartment of Electrical Machines, Technical University of Cluj-Naoca RO-400020 Cluj-Naoca, Romania; e-mail:
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 informationNONLINEAR OPTIMIZATION WITH CONVEX CONSTRAINTS. The Goldstein-Levitin-Polyak algorithm
- (23) NLP - NONLINEAR OPTIMIZATION WITH CONVEX CONSTRAINTS The Goldstein-Levitin-Polya algorithm We consider an algorithm for solving the otimization roblem under convex constraints. Although the convexity
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 informationOn the relationship between sound intensity and wave impedance
Buenos Aires 5 to 9 Setember, 16 Acoustics for the 1 st Century PROCEEDINGS of the nd International Congress on Acoustics Sound Intensity and Inverse Methods in Acoustics: Paer ICA16-198 On the relationshi
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
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 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 informationSTRONG TYPE INEQUALITIES AND AN ALMOST-ORTHOGONALITY PRINCIPLE FOR FAMILIES OF MAXIMAL OPERATORS ALONG DIRECTIONS IN R 2
STRONG TYPE INEQUALITIES AND AN ALMOST-ORTHOGONALITY PRINCIPLE FOR FAMILIES OF MAXIMAL OPERATORS ALONG DIRECTIONS IN R 2 ANGELES ALFONSECA Abstract In this aer we rove an almost-orthogonality rincile for
More informationNotes on duality in second order and -order cone otimization E. D. Andersen Λ, C. Roos y, and T. Terlaky z Aril 6, 000 Abstract Recently, the so-calle
McMaster University Advanced Otimization Laboratory Title: Notes on duality in second order and -order cone otimization Author: Erling D. Andersen, Cornelis Roos and Tamás Terlaky AdvOl-Reort No. 000/8
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 informationICES REPORT Contact Treatment in Isogeometric Analysis with NURBS
ICES REPORT - July Contact Treatment in Isogeometric Analysis with NURBS by I. Temizer, P. Wriggers, T.J.R. Hughes The Institute for Comutational Engineering and Sciences The University of Texas at Austin
More informationOn a class of Rellich inequalities
On a class of Rellich inequalities G. Barbatis A. Tertikas Dedicated to Professor E.B. Davies on the occasion of his 60th birthday Abstract We rove Rellich and imroved Rellich inequalities that involve
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 informationThe Fekete Szegő theorem with splitting conditions: Part I
ACTA ARITHMETICA XCIII.2 (2000) The Fekete Szegő theorem with slitting conditions: Part I by Robert Rumely (Athens, GA) A classical theorem of Fekete and Szegő [4] says that if E is a comact set in the
More informationA Model for Randomly Correlated Deposition
A Model for Randomly Correlated Deosition B. Karadjov and A. Proykova Faculty of Physics, University of Sofia, 5 J. Bourchier Blvd. Sofia-116, Bulgaria ana@hys.uni-sofia.bg Abstract: A simle, discrete,
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 informationMulti-Operation Multi-Machine Scheduling
Multi-Oeration Multi-Machine Scheduling Weizhen Mao he College of William and Mary, Williamsburg VA 3185, USA Abstract. In the multi-oeration scheduling that arises in industrial engineering, each job
More informationAsymptotically Optimal Simulation Allocation under Dependent Sampling
Asymtotically Otimal Simulation Allocation under Deendent Samling Xiaoing Xiong The Robert H. Smith School of Business, University of Maryland, College Park, MD 20742-1815, USA, xiaoingx@yahoo.com Sandee
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 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 informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationDeformation Effect Simulation and Optimization for Double Front Axle Steering Mechanism
0 4th International Conference on Comuter Modeling and Simulation (ICCMS 0) IPCSIT vol. (0) (0) IACSIT Press, Singaore Deformation Effect Simulation and Otimization for Double Front Axle Steering Mechanism
More informationFor q 0; 1; : : : ; `? 1, we have m 0; 1; : : : ; q? 1. The set fh j(x) : j 0; 1; ; : : : ; `? 1g forms a basis for the tness functions dened on the i
Comuting with Haar Functions Sami Khuri Deartment of Mathematics and Comuter Science San Jose State University One Washington Square San Jose, CA 9519-0103, USA khuri@juiter.sjsu.edu Fax: (40)94-500 Keywords:
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 informationAn Algorithm for Total Variation Minimization and Applications
Journal of Mathematical Imaging and Vision 0: 89 97, 004 c 004 Kluwer Academic Publishers. Manufactured in The Netherlands. An Algorithm for Total Variation Minimization and Alications ANTONIN CHAMBOLLE
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 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 informationA filter-based computational homogenization method for handling non-separated scales problems
A filter-based comutational homogenization method for handling non-searated scales roblems Julien Yvonnet, Amen Tognevi, Guy Bonnet, Mohamed Guerich To cite this version: Julien Yvonnet, Amen Tognevi,
More informationSTABILITY ANALYSIS AND CONTROL OF STOCHASTIC DYNAMIC SYSTEMS USING POLYNOMIAL CHAOS. A Dissertation JAMES ROBERT FISHER
STABILITY ANALYSIS AND CONTROL OF STOCHASTIC DYNAMIC SYSTEMS USING POLYNOMIAL CHAOS A Dissertation by JAMES ROBERT FISHER Submitted to the Office of Graduate Studies of Texas A&M University in artial fulfillment
More informationAn effective blob approach to unsteady thermal buoyant flows
Comutational Methods and Exerimental Measurements XII 685 An effective blob aroach to unsteady thermal buoyant flows C. Golia & B. Buonomo Deartment of Aerosace and Mechanical Engineering, Second University
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationNode-voltage method using virtual current sources technique for special cases
Node-oltage method using irtual current sources technique for secial cases George E. Chatzarakis and Marina D. Tortoreli Electrical and Electronics Engineering Deartments, School of Pedagogical and Technological
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 informationA Closed-Form Solution to the Minimum V 2
Celestial Mechanics and Dynamical Astronomy manuscrit No. (will be inserted by the editor) Martín Avendaño Daniele Mortari A Closed-Form Solution to the Minimum V tot Lambert s Problem Received: Month
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationComputer arithmetic. Intensive Computation. Annalisa Massini 2017/2018
Comuter arithmetic Intensive Comutation Annalisa Massini 7/8 Intensive Comutation - 7/8 References Comuter Architecture - A Quantitative Aroach Hennessy Patterson Aendix J Intensive Comutation - 7/8 3
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 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 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 informationSolution sheet ξi ξ < ξ i+1 0 otherwise ξ ξ i N i,p 1 (ξ) + where 0 0
Advanced Finite Elements MA5337 - WS7/8 Solution sheet This exercise sheets deals with B-slines and NURBS, which are the basis of isogeometric analysis as they will later relace the olynomial ansatz-functions
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 informationDevelopment of self-adaptively loading for planetary roller traction-drive transmission
Available online www.jocr.com Journal of Chemical and Pharmaceutical Research, 013, 5(9):498-506 Research Article ISSN : 0975-7384 CODEN(USA) : JCPRC5 Develoment of self-adatively loading for lanetary
More informationStrong Matching of Points with Geometric Shapes
Strong Matching of Points with Geometric Shaes Ahmad Biniaz Anil Maheshwari Michiel Smid School of Comuter Science, Carleton University, Ottawa, Canada December 9, 05 In memory of Ferran Hurtado. Abstract
More informationCorrespondence Between Fractal-Wavelet. Transforms and Iterated Function Systems. With Grey Level Maps. F. Mendivil and E.R.
1 Corresondence Between Fractal-Wavelet Transforms and Iterated Function Systems With Grey Level Mas F. Mendivil and E.R. Vrscay Deartment of Alied Mathematics Faculty of Mathematics University of Waterloo
More informationOn Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm
On Line Parameter Estimation of Electric Systems using the Bacterial Foraging Algorithm Gabriel Noriega, José Restreo, Víctor Guzmán, Maribel Giménez and José Aller Universidad Simón Bolívar Valle de Sartenejas,
More informationMODELING AND SIMULATION OF A SATELLITE PROPULSION SUBSYSTEM BY PHYSICAL AND SIGNAL FLOWS. Leonardo Leite Oliva. Marcelo Lopes de Oliveira e Souza
Satellite Proulsion Subsystem MODELING AND SIMULATION OF A SATELLITE PROPULSION SUBSYSTEM BY PHYSICAL AND SIGNAL FLOWS Leonardo Leite Oliva National Institute for Sace Research, INPE Av. dos Astronautas,
More informationSolution of time domain electric field integral equation for arbitrarily shaped dielectric bodies using an unconditionally stable methodology
RADIO SCIENCE, VOL. 38, NO. 3, 148, doi:1.129/22rs2759, 23 Solution of time domain electric field integral equation for arbitrarily shaed dielectric bodies using an unconditionally stable methodology Young-seek
More informationMultiplicative group law on the folium of Descartes
Multilicative grou law on the folium of Descartes Steluţa Pricoie and Constantin Udrişte Abstract. The folium of Descartes is still studied and understood today. Not only did it rovide for the roof of
More informationAn Inverse Problem for Two Spectra of Complex Finite Jacobi Matrices
Coyright 202 Tech Science Press CMES, vol.86, no.4,.30-39, 202 An Inverse Problem for Two Sectra of Comlex Finite Jacobi Matrices Gusein Sh. Guseinov Abstract: This aer deals with the inverse sectral roblem
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 informationc Copyright by Helen J. Elwood December, 2011
c Coyright by Helen J. Elwood December, 2011 CONSTRUCTING COMPLEX EQUIANGULAR PARSEVAL FRAMES A Dissertation Presented to the Faculty of the Deartment of Mathematics University of Houston In Partial Fulfillment
More informationA Numerical Method for Critical Buckling Load for a Beam Supported on Elastic Foundation
A Numerical Method for Critical Buckling Load for a Beam Suorted on Elastic Foundation Guo-ing Xia Institute of Bridge Engineering, Dalian University of Technology, Dalian, Liaoning Province, P. R. China
More informationEvaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models
Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Ketan N. Patel, Igor L. Markov and John P. Hayes University of Michigan, Ann Arbor 48109-2122 {knatel,imarkov,jhayes}@eecs.umich.edu
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 informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationLINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL
LINEAR SYSTEMS WITH POLYNOMIAL UNCERTAINTY STRUCTURE: STABILITY MARGINS AND CONTROL Mohammad Bozorg Deatment of Mechanical Engineering University of Yazd P. O. Box 89195-741 Yazd Iran Fax: +98-351-750110
More informationMath 751 Lecture Notes Week 3
Math 751 Lecture Notes Week 3 Setember 25, 2014 1 Fundamental grou of a circle Theorem 1. Let φ : Z π 1 (S 1 ) be given by n [ω n ], where ω n : I S 1 R 2 is the loo ω n (s) = (cos(2πns), sin(2πns)). Then
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 informationA PRIMAL-DUAL ACTIVE SET ALGORITHM FOR THREE-DIMENSIONAL CONTACT PROBLEMS WITH COULOMB FRICTION
A PRIMAL-DUAL ACTIVE SET ALGORITHM FOR THREE-DIMENSIONAL CONTACT PROBLEMS WITH COULOMB FRICTION S. HÜEBER, G. STADLER, AND B.I. WOHLMUTH Abstract. In tis aer, efficient algoritms for contact roblems wit
More informationKinetics of Protein Adsorption and Desorption on Surfaces with Grafted Polymers
1516 Biohysical Journal Volume 89 Setember 2005 1516 1533 Kinetics of Protein Adsortion and Desortion on Surfaces with Grafted Polymers Fang Fang,* Javier Satulovsky, y and Igal Szleifer* *Deartment of
More informationChapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method
CIVL 7/87 Chater - The Stiffness Method / Chater Introdction to the Stiffness (Dislacement) Method Learning Objectives To define the stiffness matrix To derive the stiffness matrix for a sring element
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 informationThickness and refractive index measurements using multiple beam interference fringes (FECO)
Journal of Colloid and Interface Science 264 2003 548 553 Note www.elsevier.com/locate/jcis Thickness and refractive index measurements using multile beam interference fringes FECO Rafael Tadmor, 1 Nianhuan
More informationBEST CONSTANT IN POINCARÉ INEQUALITIES WITH TRACES: A FREE DISCONTINUITY APPROACH
BEST CONSTANT IN POINCARÉ INEQUALITIES WITH TRACES: A FREE DISCONTINUITY APPROACH DORIN BUCUR, ALESSANDRO GIACOMINI, AND PAOLA TREBESCHI Abstract For Ω R N oen bounded and with a Lischitz boundary, and
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 informationRe-entry Protocols for Seismically Active Mines Using Statistical Analysis of Aftershock Sequences
Re-entry Protocols for Seismically Active Mines Using Statistical Analysis of Aftershock Sequences J.A. Vallejos & S.M. McKinnon Queen s University, Kingston, ON, Canada ABSTRACT: Re-entry rotocols are
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 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 informationDucted Wind/Water Turbines and Propellers Revisited By Michael, J. Werle, PhD 1 and Walter M. Presz, Jr., PhD 2 FLODESIGN, INC. WILBRAHAM, MA.
Introduction Ducted Wind/Water Turbines and roellers Revisited By Michael, J. Werle, hd and Walter M. resz, Jr., hd FLODEIGN, IN. WILBRAHAM, MA. 0095 There has been considerable effort and discussion in
More information