Computational Methods and Applications in Parachute Maneuvers
|
|
- Marjory Mosley
- 5 years ago
- Views:
Transcription
1 COMPUTATIONAL MECHANICS WCCM VI in conjunction with APCOM 04, Sept. 5-10, 2004, Beijing, China c 2004 Tsinghua University Press & Springer-Verlag Computational Methods and Applications in Parachute Maneuvers Keith R. Stein 1, Tayfun E. Tezduyar 2, Sunil V. Sathe 2, Masayoshi Senga 2, Richard J. Benney 3 and Richard D. Charles 3 1 Department of Physics, Bethel University, St. Paul, Minnesota, U.S.A. 2 Mechanical Engineering, Rice University - MS 321, Houston, Texas, U.S.A. 3 U.S. Army Soldier Systems Center, Natick, Massachusetts, U.S.A. k-stein@bethel.edu, tezduyar@rice.edu, sathe@rice.edu, msenga@rice.edu, richard.benney@natick.army.mil, richard.charles@natick.army.mil Abstract: Fluid structure interactions (FSI) are experienced during all stages of parachute operations. During special parachute maneuvers, and even during terminal descent, parachute canopies experience significant dynamics and undergo large shape changes. The FSI involved in the parachute operations introduce a variety of challenges that must be properly addressed in computational modeling efforts [1]. These challenges provide motivation for using methods that are based on interface-tracking techniques. Accordingly, we use the Deforming-Spatial- Domain/Stabilized Space Time (DSD/SST) [2] formulation for our flow computations and a Lagrangian finite element formulation for the parachute deformation computations [3], with advanced mesh update methods. Here, we present our methods for simulation of parachute FSI behavior and address some of the modeling challenges involved. Results are presented for round parachutes, with focus on maneuvers of the steerable kind. Key words: airdrop systems, fluid structure interaction INTRODUCTION New parachute systems are being developed for aerial delivery of supplies with greater precision than has been achieved in the past. One type of parachute system that offers increased precision is the steerable round parachute, which allows for limited maneuverability and glide to be achieved using a traditional round parachute system. When equipped with a guidance system, this type of parachute can offer automated maneuvering, can increase drop zone accuracy, and can allow for in-flight corrections during parachute descent to overcome wind-induced error. Non-steerable round parachutes are very much under the control of the winds after leaving the aircraft and, for this reason, are often referred to as dumb parachutes. Additionally, non-steerable round parachutes have no capability for limiting drop zone dispersion between the first and last parachutes exiting the aircraft. For this reason, there is interest in developing a class of precision parachute systems, with steerable round parachutes being an economical solution for many airdrop situations. In this paper we focus on the behavior of a steerable round parachute during maneuvers. Specifically, we study the limited glide that can be achieved with a round parachute by deflecting a portion of the canopy through riser control. Each riser connects the payload to multiple suspension lines, which are attached to a section of the canopy skirt. Thus, each riser
2 controls some section of the parachute canopy. Here, initial gliding performance of a G 12 cargo parachute is simulated through the deflection of one half of the canopy by retraction of two risers. Simulation of parachute maneuvers involves a variety of challenges that must be addressed in accurate modeling. Parachute canopies have a strong response to the surrounding flow fields, experiencing significant dynamics and undergoing large shape changes, even during uncontrolled descent. During maneuver (i.e., riser controls), canopy dynamics in fluid structure interactions (FSI) is increased due to the forced deflections in the canopy and the subsequent response of the parachute. Following the initiation of the riser control, the parachute pitches towards (or away from) the retracted (or extended) section of the canopy and achieves a limited amount of glide. The pitching of the parachute leads into gliding descent, which continues as long as the riser control is held. For the simulations described in this paper, the parachute operates at sufficiently low speed and, hence, the aerodynamics is governed by the Navier Stokes equations of incompressible flows. The parachute is treated as a cable membrane tension structure, governed by the linear momentum balance equation. Parachute FSI during control line maneuvers involve canopies that undergo shape changes, with the spatial domain occupied by the fluid varying (i.e. deforming) with respect to time. To accomodate canopy shape changes, we use methods which are based on interface-tracking techniques. We use the Deforming-Spatial- Domain/Stabilized Space Time (DSD/SST) [2] formulation for the flow field computations and a Lagrangian finite element formulation [3] for the parachute deformation computations. In addition to the fluid and structural computations, we use advanced mesh update methods that we have developed to handle parachute deformations within the fluid mesh. FSI behavior is modeled with the quasi-direct coupling method introduced in [4]. We use a six degree-of-freedom (6-dof) modeling approach to simulate the glide behavior of the parachute after the parachute has had some time to respond to the riser control. Fluid object coupling for the 6-dof approach is achieved in a block-iterative fashion, similar to the approach we have widely used in earlier FSI computations [5, 6]. The 6-dof approach is more economical than a coupled FSI computation and allows for faster computation of the glide performance of the parachute with deflected canopy. We assume here that glide performance is driven primarily by the large canopy deflections imposed by the riser controls, and less by the smaller canopy motions due to the FSI that continue throughout the descent. In the following sections we describe our computational methods and demonstrate these methods for a steerable round cargo parachute. We begin in the next section by describing the governing equations and finite element methods, including the DSD/SST finite element formulation. A brief description of our mesh update methods is included. Following that, we describe the quasi-direct coupling used for FSI computations and the block-iterative approach used for 6-dof computations. After this, we present results from the FSI and 6-dof computation involving a steerable round cargo parachute system. Finally, we end with some concluding remarks. FINITE ELEMENT MODEL Fluid Dynamics Let Ω t IR n sd be the spatial fluid mechanics domain with boundary Γt at time t (0, T), where the subscript t indicates the time-dependence of the spatial domain and its boundary.
3 The Navier Stokes equations of incompressible flows can be written on Ω t and t (0, T) as ( ) u ρ t + u u f σ = 0, (1) u = 0, (2) where ρ, u and f are the density, velocity and external force, respectively. The stress tensor σ is defined as σ(p,u) = pi + 2µε(u). Here p, I and µ are the pressure, identity tensor and viscosity, respectively. The strain rate tensor is defined as ε(u) = 1 ( ) ( u) + ( u) T. (4) 2 Both Dirichlet- and Neumann-type boundary conditions are accounted for: u = g on (Γ t ) g, n σ = h on (Γ t ) h. (5) Here (Γ t ) g and (Γ t ) h are complementary subsets of the boundary Γ t, n is the unit normal vector at the boundary, and g and h are given functions. A divergence-free velocity field is specified as the initial condition. The DSD/SST finite element method is suitable for handling spatial domains for the fluid that change in time. In this method, the finite element formulation of the governing equations is written over a sequence of N space time slabs Q n, where Q n is the slice of the space time domain between the time levels t n and t n+1. At each time step, the integrations involved in the finite element formulation are performed over Q n. The space time finite element interpolation functions are continuous within a space time slab, but discontinuous from one space time slab to another. The notation ( ) n and ( ) + n denotes the function values at t n as approached from below and above. Each Q n is decomposed into space time elements Q e n, where e = 1, 2,..., (n el ) n. The subscript n used with n el is to account for the general case in which the number of space time elements may change from one space time slab to another. The Dirichlet- and Neumann-type boundary conditions are enforced over (P n ) g and (P n ) h, the complementary subsets of the lateral boundary of the space time slab. The finite element trial function spaces (Su) h n for velocity and (Sp h ) n for pressure, and the test function spaces (Vu h) n and (Vp h) n = (Sp h) n are defined by using, over Q n, first-order polynomials in both space and time. The DSD/SST formulation is written as follows: given (u h ) n, find u h (Su) h n and p h (Sp h) n such that w h (Vu h) n and q h (Vp h) n: ( ) u w h h ρ Q n t + uh u h f h dq + ε(w h ) : σ(p h,u h )dq Q n w h h h dp + q h u h dq + (w h ) + n ρ ( ) (u h ) + n (u h ) n dω (P n) h Q n Ω n + + [ ( ) ] 1 w h τ SUPG ρ + u h w h + τ PSPG q h Q e n ρ t ( ) ] u h ρ t + uh u h σ(p h,u h ) ρf h dq ν LSIC w h ρ u h dq = 0, (6) (n el ) n e=1 [ (n el ) n e=1 Q e n (3)
4 where τ SUPG (streamline-upwind/petrov-galerkin), τ PSPG (pressure-stabilizing/petrov-galerkin), and ν LSIC (least-squares on incompressibility constraint) are the stabilization parameters (see [7]). This formulation is applied to all space time slabs Q 0, Q 1, Q 2,...,Q N 1, starting with (u h ) 0 = u 0. For an earlier, detailed reference on this stabilized formulation see [2]. For the FSI computations in this paper, w h / t has been dropped from the SUPG stabilization term in Eqn. (6). Structural Dynamics Let Ω s t IRn xd be the spatial domain bounded by Γ s t, where n xd = 2 for membranes and n xd = 1 for cables. The boundary Γ s t is composed of (Γ s t) g and (Γ s t) h. Here, the superscript s corresponds to the structure. The equations of motion for the structural system are: ( d ρ s 2 ) y dt + ηdy 2 dt fs σ s = 0, (7) where, y is the displacement, ρ s is the material density, f s is the external body force, σ s is the Cauchy stress tensor, and η is the mass-proportional damping coefficient. The damping provides additional stability and is used in inflating (i.e., prestressing) the parachute, but not during the FSI computations presented in this paper. Our numerical method uses total Lagrangian formulation of the problem. Thus, stresses are expressed in terms of the 2nd Piola Kirchoff stress tensor S, which is related to the Cauchy stress tensor through a kinematic transformation. Under the assumption of large displacements and rotations, small strains, and no material damping, the membranes and cables are treated as Hookean materials with linear elastic properties. For membranes, under the assumption of plane stress, S becomes S ij = ( λm G ij G kl + µ m [ G il G jk + G ik G jl]) E kl, (8) where λ m = 2λ mµ m (λ m + 2µ m ). (9) Here, E kl are the components of the Cauchy-Green strain tensor, G ij are the components of the contravariant metric tensor in the original configuration, and λ m and µ m are Lamé constants. For cables, under the assumption of uniaxial tension, S becomes S 11 = E c G 11 G 11 E 11, (10) where E c is the cable Young s modulus. The semi-discrete finite element formulation for the structural dynamics is based on the principle of virtual work: Ω s 0 w h ρ sd2 y h dt 2 dωs + w h ηρ sdyh Ω s dt dωs + δe : S h dω s = w h (t + ρ s f s )dω s. (11) 0 Ω s 0 Ω s t Here the weighting function w h is also the virtual displacement. The air pressure force on the canopy surface is represented by vector t, a geometrically nonlinear term that increases the overall nonlinearity of the formulation. The left-hand-side terms of Eq. (11) are referred to in the original configuration and the right-hand-side terms for the deformed configuration at time t.
5 Mesh Update Appropriate mesh update methods are needed to handle spatial domains for the fluid that change in time. Selection of these methods depends on several factors, such as the complexity of the moving boundary or interface, how unsteady the moving boundary or interface is, and how the initial mesh is generated. In general, the mesh update consists of moving the mesh in a way that limits element distortion, with full or partial remeshing when the element distortion becomes excessive. The only condition that mesh motion must satisfy is that at the moving boundary or interface the normal velocity of the mesh has to match the normal velocity of the fluid. Beyond that, we have flexibility in how to move the mesh, with the main objective being to reduce the frequency of remeshing. In computation of parachute FSI reported here we use an automatic mesh moving scheme [8] based on the equations of linear elasticity, with selective treatment of mesh elements and also the deformation modes in terms of shape and volume changes so that smaller elements enjoy more protection from mesh deformation. The motion of the internal nodes is determined by solving these additional equations, with the boundary conditions for these mesh motion equations specified in such a way that they match the normal velocity of the fluid at the interface. Several additional features have been developed and added to our central mesh moving strategy to reduce the frequency of remeshing [9]. These features have proven effective for a variety of test cases [10, 11]. COUPLING Fluid Structure Interactions Parachute FSI involve situations where the structure is very light, with structural response very sensitive to small changes in the fluid dynamics. In these types of problems, achieving acceptable levels of convergence in the nonlinear equation systems for the FSI can be difficult. We use the quasi-direct coupling method, which has been shown to be effective for FSI involving light structures [4]. Spatial and temporal discretizations of the finite element formulations for the fluid and structure (see Eqs. (6) and (11)) lead to a fully-coupled, nonlinear system of equations that can be symbolically written as N 1 (d 1,d 2 ) = F 1, N 2 (d 1,d 2 ) = F 2, (12) (13) where d 1 and d 2 are the vectors of nodal unknowns associated with the fluid unknown function (u 1 ) and structure unknown function (u 2 ), respectively. Solution of these equations with Newton Raphson method would require the solution in each iteration of a linearized system of the form A 11 x 1 + A 12 x 2 = b 1, A 21 x 1 + A 22 x 2 = b 2, (14) where b 1 and b 2 are the residuals of the nonlinear system, x 1 and x 2 are incremental updates for d 1 and d 2, and A βγ = N β d γ. (15) In our situation, the off-diagonal blocks A 12 and A 21 represent the fluid structure coupling matrices. Retaining these coupling matrices in the current form would require accounting
6 for the dependence of A 12 on mesh motion. In this paper, we use the quasi-direct coupling approach, which partially includes the coupling blocks, without accounting for the dependence of A 12 on mesh motion. This approach provides a much tighter fluid structure coupling than the block-iterative coupling and, therefore, delivers significant improvement in convergence. Fluid Object Interactions In certain situations it may be desirable to model a parachute as a rigid body with the dynamics governed by the force and moment balance equations. This 6-dof approach is more ecomonical than a coupled FSI computation and can be used to simulate parachute glide after the canopy is deflected with a riser control. For the 6-dof simulations presented in the following section we use a block-iterative method for the fluid-object coupling. Here, the parachute is treated as a rigid object with the deformed parachute geometry resulting from a preceding FSI computation. The block-iterative coupling technique can be viewed as an approximate Newton Raphson method, where the fluid object coupling matrices are ignored in the solution of the linear equation system each Newton Raphson iteration. This is equivalent to dropping A 12 and A 21 from Eq. (14), where b 2, x 2, and d 2 now come from the system of equations for the 6-dof motion. MANEUVER OF A G 12 PARACHUTE WITH RISER CONTROLS Computations are carried out to study the response of a G 12 cargo parachute to a riser control involving the retraction of two risers, deflecting one half of the canopy. The G 12 is a 64 ft diameter cargo parachute designed to deliver payloads of 2200 lb at descent speeds of 28 ft/s. The G 12 is constructed with 64 suspension lines each 51.2 ft that extend from the canopy to four risers each ft. The risers are each connected to 16 suspension lines and merge to a single confluence point, which in turn is connected to four cables each ft that hold the payload. A schematic with dimensions for the G 12 model is shown in Figure 1 (left). The structural model is composed of membranes, cables, and concentrated masses. The canopy is modeled with triangular membrane elements. Linear cable elements are used to model the suspension lines, radial reinforcements along the canopy, risers, and payload support cables. The payload is modeled with eight concentrated masses that are interconnected with a set of truss elements. Material properties are selected to be representative of the G 12. Our model for the G 12 in its constructed configuration is shown in Figure 1 (middle). Several measures are taken in preparations for the G 12 riser control simulations. First, a stand alone structural deformation computation is carried out to determine the inflated (i.e., prestressed) shape of the G 12 (See Figure 1, right). Secondly, a stand alone flow field computation is carried out to obtain a developed flow about the inflated G 12. Finally, a FSI computation is carried out prior to the riser control to remove any mismatch between the fluid forces and the parachute shape. The simulations presented in the following two subsections are carried out in two stages. In the first stage, canopy shape changes and the onset of glide that results from riser retraction is simulated with a FSI computation. In the second stage, glide behavior is simulated with a 6-dof computation. Here, we assume that glide performance is driven primarily by the significant canopy shape changes that result from the retracted risers, and less by smaller canopy motions from continued FSI behavior.
7 Figure 1: G 12 Parachute: Dimensions (left); Constructed geometry (middle); Prestressed geometry (right). Fluid Structure Interaction of G 12 Parachute with Retraction of Two Risers The response of the G 12 to a two-riser retraction is simulated using the quasi-direct method described. The FSI for the riser control is carried out in two computations. The first computation is carried out to model the dynamic retraction of two risers by 3.84 f t in 0.29 s. Riser retraction is modeled by prescribing that the cables defining the retractable risers shrink in time. This is accomplished by decreasing the natural length of the cables over the duration of the retraction. These reductions of cable length are accounted for in the calculation of cable stress terms, whereas inertia terms are always based on the initial cable lengths. In the second computation the retracted risers are held, allowing the G 12 to respond to the retraction for an additional 0.22 s and to begin to enter into a state of glide. Results for the FSI of the two-riser retraction are shown in Figures 2-3. The G 12 geometry is shown in Figure 2 prior to the retraction (left) and at the end of the FSI computations (right). The shortening of the two retracted risers and the deflection of the right half of the canopy are evident. Although the canopy has undergone a significant change in shape, the resulting pitch and glide of the parachute are just being initiated at this point in time. The pitching and gliding become more apparent during the following 6-dof computation. The flow field surrounding the G 12 at four stages of the FSI computation are shown in Figure 3, with colors representing the pressure field surrounding the canopy. The pressure field from the FSI results in a horizontal component of aerodynamic force in the direction of anticipated glide. As expected, retraction of the two risers results in a horizontal force which will tend to drive the parachute towards the retracted risers (i.e., the +x-direction). The transition to glide is simulated in the 6-dof computations, as described in the following subsection. Fluid Object Interaction for 6-Dof Behavior of G 12 Parachute with Two Retracted Risers The solution at the end of the FSI computation is used as an initial condition for an ensuing 6-dof simulation to study the gliding performance of the G 12 with the deflected canopy due to the retracted risers. Here, the G 12 parachute is treated as a rigid body and the dynamics is determined based on the force and moment balance equations, with forces
8 Figure 2: G 12 geometry prior to (left) and after (right) the two-riser retraction. and moments arising from the unsteady fluid pressure distribution on the canopy surface. The simulation is carried out for a physical time of approximately 5.4 s, during which the canopy pitches towards the deflected section of the canopy and transitions to glide. During transition, the canopy pitches from a nearly vertical orientation due to a moment induced by the deflected canopy geometry. Eventually, the canopy pitches beyond equilibrium and the moment becomes negative. In time, the G 12 will approach an equilibrium glide orientation, although some oscillations will continue due to the unsteady nature of the flow field. To handle the gliding parachute, we specify boundary conditions for the mesh motion equations such that the outer boundaries of the mesh move with the center of mass of the parachute canopy. Implementation of this mesh boundary condition assures that the parachute canopy interface remains centered in the fluid mesh. Results from the 6-dof computations for the gliding performance of the G 12 with deflected risers are shown in Figures 4-6. The trajectory of the payload and the velocity for the center of mass are shown in Figure 4. For our simulations, the two retracted risers are on the +x side of the parachute and will induce a preferred glide in the x z plane. Hence, the trajectory plot is displayed in the x z plane, with the understanding that induced motion in the y- direction is minimal. The trajectory indicates a transition from vertical descent to glide in the +x-direction and is shown along with the anticipated trajectory without the maneuver. The components of velocity in the G 12 glide plane are shown for the center of mass (Figure 4, right). As expected, induced glide is primarily in the +x-direction, with the glide speed building up as the simulation progresses. It is apparent from the figure that the glide speed and descent speed are of the same order of magnitude by the end of the computation. The time history for the horizontal force in the +x-direction is shown in Figure 5. This is the glide-inducing component the total aerodynamic force, and is shown along with the weight of the G 12 parachute payload system. The build up of the x-component of force, which initiates the pitching motion of the G 12 about the y-axis, is evident. The time history
9 Figure 3: Sequence of G-12 pictures (left to right and top to bottom) showing the velocity vectors colored by pressure, during and immediately after two-riser retraction. for the glide-inducing component of the moment (i.e., about the y-axis) acting on the G 12 canopy is shown in Figure 6. This component of the total moment is dominant, and remains positive initially and becomes negative at about 3.25 s. The moment causes the G 12 to begin pitching forward around the y-axis with increasing angular speed for the first 3.25 s and then to continue pitching forward with decreasing angular speed. Eventually, the G 12 will pitch back and settle into an equilibrium glide orientation. CONCLUDING REMARKS We described our computational methods for the simulation of parachute fluid structure interaction that are applicable to situations involving maneuvers for round parachutes. These problems involve parachute canopies that interact strongly with the surrounding air flow and undergo large shape changes in time. Hence, our methods are based on interfacetracking techniques. We have used the Deforming-Spatial-Domain/Stabilized Space Time (DSD/SST) formulation, with advanced mesh update methods, which is effective in solving the fluid dynamics for this class of problems. The DSD/SST formulation and a Lagrangian finite element formulation for parachute deformations were coupled in a quasi-direct fashion for the FSI studies presented in this paper. For the 6-dof studies, we coupled the DSD/SST method to the parachute force and moment balance equations in a block-iterative fashion. Simulations were carried out to study maneuvers for a round cargo parachute. Specifically, we focused on the glide performance of a G 12 parachute due to the retraction of risers controlling one half of the canopy. Results were presented for the dynamic response resulting from the retraction of two risers, the transition to glide, and the gliding performance. FSI computations based on the quasi-direct coupling approach were presented to study the initial behaviors, for which FSI play a most dominant role. These computations were carried out to study the dynamic response of a G 12 parachute to the surrounding flow field during the
10 Trajectory of payload Trajectory of payload without maneuver Glide speed Descent speed Vertical Position (ft) Speed (ft/s) Horizontal Position (ft) Time (s) Figure 4: Trajectory of G 12 payload (left) and velocity of center of mass (right) during 6-dof motion Force in glide direction Weight of parachute system Force (lb) Time (s) Figure 5: Horizontal force in the glide direction and G 12 parachute system weight. retraction of two risers. Follow on 6-dof computations were carried out for the G 12 with the deflected canopy to study the transition to glide for the parachute. The simulations presented here have demonstrated how the quasi-direct coupling method for FSI can be used along with a more efficient 6-dof fluid-solid interaction model to study relevant problems involving parachute systems. Acknowledgments: This work was supported by the Army Natick Soldier Center.
11 Moment (lb-ft) Time (s) Figure 6: Moment acting on a G-12 about an anticipated pitching axis during 6-dof motion. REFERENCES [1] K. Stein, T. Tezduyar, and R. Benney, Applications in airdrop systems: Fluid-structure interaction modeling, in Proceedings of the Fifth World Congress on Computational Mechanics (Web Site), On-line publication: Paper-ID: 81545, Vienna, Austria, (2002). [2] T.E. Tezduyar, Stabilized finite element formulations for incompressible flow computations, Advances in Applied Mechanics, 28 (1992) [3] M.L. Accorsi, J.W. Leonard, R. Benney, and K. Stein, Structural modeling of parachute dynamics, AIAA Journal, 38 (2000) [4] T.E. Tezduyar, S. Sathe, R. Keedy, and K. Stein, Space time techniques for finite element computation of flows with moving boundaries and interfaces, in S. Gallegos, I. Herrera, S. Botello, F. Zarate, and G. Ayala, editors, Proceedings of the III International Congress on Numerical Methods in Engineering and Applied Science, CD-ROM, [5] K. Stein, R. Benney, V. Kalro, T.E. Tezduyar, J. Leonard, and M. Accorsi, Parachute fluid structure interactions: 3-D Computation, Computer Methods in Applied Mechanics and Engineering, 190 (2000) [6] K.R. Stein, R.J. Benney, T.E. Tezduyar, J.W. Leonard, and M.L. Accorsi, Fluid structure interactions of a round parachute: Modeling and simulation techniques, Journal of Aircraft, 38 (2001) [7] T. Tezduyar and Y. Osawa, Methods for parallel computation of complex flow problems, Parallel Computing, 25 (1999) [8] T.E. Tezduyar, M. Behr, S. Mittal, and A.A. Johnson, Computation of unsteady incompressible flows with the finite element methods space time formulations, iterative strategies and massively parallel implementations, in New Methods in Transient Analysis, PVP-Vol.246/AMD-Vol.143, ASME, New York, (1992) 7 24.
12 [9] T. Tezduyar, Finite element interface-tracking and interface-capturing techniques for flows with moving boundaries and interfaces, in Proceedings of the ASME Symposium on Fluid-Physics and Heat Transfer for Macro- and Micro-Scale Gas-Liquid and Phase- Change Flows (CD-ROM), ASME Paper IMECE2001/HTD-24206, ASME, New York, New York, (2001). [10] K. Stein, T. Tezduyar, and R. Benney, Mesh moving techniques for fluid structure interactions with large displacements, Journal of Applied Mechanics, 70 (2003) [11] K. Stein and T. Tezduyar, Advanced mesh update techniques for problems involving large displacements, in Proceedings of the Fifth World Congress on Computational Mechanics, On-line publication: Paper-ID: 81489, Vienna, Austria, (2002).
T.E. Tezduyar, S. Sathe, R. Keedy. K. Stein
III Congreso Internacional sobre Metodos Numericos en Ingenieria y Ciencias Aplicadas S. Gallegos, I. Herrera, S. Botello, F. Zarate, y G. Ayala Editores) c ITESM, Monterrey 004 CIMNE, Barcelona 004 SPACE
More informationAerodynamic Interactions Between Parachute Canopies
K. Stein Department of Physics, Bethel College, St. Paul, MN 55112 T. Tezduyar V. Kumar S. Sathe Mechanical Engineering, Rice University, MS 321, Houston, TX 77005 R. Benney Natick Soldier Center, Natick,
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationSUPG STABILIZATION PARAMETERS CALCULATED FROM THE QUADRATURE-POINT COMPONENTS OF THE ELEMENT-LEVEL MATRICES
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 24 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) W. Rodi and P. Le Quéré (assoc. eds.) Jyväskylä,
More informationADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS
ECCOMAS Computational Fluid Dynamics Conference 2001 Swansea, Wales, UK, 4-7 September 2001 c ECCOMAS ADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS Tayfun E. Tezduyar Team for Advanced
More informationParachute uid±structure interactions: 3-D computation
Comput. Methods Appl. Mech. Engrg. 190 (2000) 373±386 www.elsevier.com/locate/cma Parachute uid±structure interactions: 3-D computation Keith Stein a, *, Richard Benney a, Vinay Kalro b, Tayfun E. Tezduyar
More informationSolution techniques for the fully discretized equations in computation of fluid structure interactions with the space time formulations
Comput. Methods Appl. Mech. Engrg. 95 (2006) 5743 5753 www.elsevier.com/locate/cma Solution techniques for the fully discretized equations in computation of fluid structure interactions with the space
More informationStabilized Formulations and Smagorinsky Turbulence Model for Incompressible Flows
WCCM V Fifth World Congress on Computational Mechanics July 7 2, 22, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Stabilized Formulations and Smagorinsky Turbulence Model for
More informationFINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS
FINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS Lucia Catabriga Department of Computer Science, Federal University of Espírito Santo (UFES) Av. Fernando Ferrari,
More informationStabilization and shock-capturing parameters in SUPG formulation of compressible flows
Comput. Methods Appl. Mech. Engrg. 195 (006) 161 163 www.elsevier.com/locate/cma Stabilization and shock-capturing parameters in SUPG formulation of compressible flows Tayfun E. Tezduyar *, Masayoshi Senga
More informationComputation of Cardiovascular Fluid-Structure Interactions with the DSD/SST Method
COMPUTATIONAL MECHANICS WCCM VI in conjunction with APCOM 04, Sept. 5-10, 2004, Beijing, China c 2004 Tsinghua University Press & Springer-Verlag Computation of Cardiovascular Fluid-Structure Interactions
More informationBiomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationNonlinear analysis in ADINA Structures
Nonlinear analysis in ADINA Structures Theodore Sussman, Ph.D. ADINA R&D, Inc, 2016 1 Topics presented Types of nonlinearities Materially nonlinear only Geometrically nonlinear analysis Deformation-dependent
More informationCoupling techniques for partitioned fluid structure interaction simulations with black box solvers
Coupling techniques for partitioned fluid structure interaction simulations with black box solvers Joris Degroote 1, Robby Haelterman 2, Sebastiaan Annerel 1, Jan Vierendeels 1 1 Ghent University, Department
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationA NOTE ON RELATIONSHIP BETWEEN FIXED-POLE AND MOVING-POLE APPROACHES IN STATIC AND DYNAMIC ANALYSIS OF NON-LINEAR SPATIAL BEAM STRUCTURES
European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 212) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 1-14, 212 A NOTE ON RELATIONSHIP BETWEEN FIXED-POLE
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationStabilized and Coupled FEM/EFG Approximations for Fluid Problems
COMPUTATIONAL MECHANICS WCCM VI in conjunction with APCOM4, Sept. 5-, 24, Beijing, China c 24 Tsinghua University Press & Springer-Verlag Stabilized and Coupled FEM/EFG Approximations for Fluid Problems
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationA Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics
A Harmonic Balance Approach for Large-Scale Problems in Nonlinear Structural Dynamics Allen R, PhD Candidate Peter J Attar, Assistant Professor University of Oklahoma Aerospace and Mechanical Engineering
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationAn arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem
An arbitrary lagrangian eulerian discontinuous galerkin approach to fluid-structure interaction and its application to cardiovascular problem Yifan Wang University of Houston, Department of Mathematics
More informationImplementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible Flapping Wings within Overture
10 th Symposimum on Overset Composite Grids and Solution Technology, NASA Ames Research Center Moffett Field, California, USA 1 Implementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationENGN 2290: Plasticity Computational plasticity in Abaqus
ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationDiese Arbeit wurde vorgelegt am Lehrstuhl für computergestützte Analyse technischer Systeme
Diese Arbeit wurde vorgelegt am Lehrstuhl für computergestützte Analyse technischer Systeme Dreidimensionale Isogeometrische Elastische Festkörper für Fluid-Struktur Interaktionen Three-Dimensional Isogeometric
More informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationOn Verification and Validation of Spring Fabric Model
On Verification and Validation of Spring Fabric Model Zheng Gao, Qiangqiang Shi, Yiyang Yang, Bernard Moore, and Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University Stony
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationUnified Lagrangian formulation for solid and fluid mechanics and FSI problems
Unified Lagrangian formulation for solid and fluid mechanics and FSI problems Alessandro Franci 1,2, Eugenio Oñate 1,2, Josep Maria Carbonell 1,2 1 International Center for Numerical Methods in Engineering
More informationThe Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011
The Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011 Derivation of Navier-Stokes Equation 1 The total stress tensor
More informationSimulation of Aeroelastic System with Aerodynamic Nonlinearity
Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationLarge Scale Fluid-Structure Interaction by coupling OpenFOAM with external codes
Large Scale Fluid-Structure Interaction by coupling OpenFOAM with external codes Thomas Gallinger * Alexander Kupzok Roland Wüchner Kai-Uwe Bletzinger Lehrstuhl für Statik Technische Universität München
More informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationComputer modeling of cardiovascular fluid structure interactions with the deforming-spatial-domain/stabilized space time formulation
Comput. Methods Appl. Mech. Engrg. 195 (2006) 1885 1895 www.elsevier.com/locate/cma Computer modeling of cardiovascular fluid structure interactions with the deforming-spatial-domain/stabilized space time
More informationTable of Contents. II. PDE classification II.1. Motivation and Examples. II.2. Classification. II.3. Well-posedness according to Hadamard
Table of Contents II. PDE classification II.. Motivation and Examples II.2. Classification II.3. Well-posedness according to Hadamard Chapter II (ContentChapterII) Crashtest: Reality Simulation http:www.ara.comprojectssvocrownvic.htm
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationNew Representation of Bearings in LS-DYNA
13 th International LS-DYNA Users Conference Session: Aerospace New Representation of Bearings in LS-DYNA Kelly S. Carney Samuel A. Howard NASA Glenn Research Center, Cleveland, OH 44135 Brad A. Miller
More informationMODELLING MIXED-MODE RATE-DEPENDENT DELAMINATION IN LAYERED STRUCTURES USING GEOMETRICALLY NONLINEAR BEAM FINITE ELEMENTS
PROCEEDINGS Proceedings of the 25 th UKACM Conference on Computational Mechanics 12-13 April 217, University of Birmingham Birmingham, United Kingdom MODELLING MIXED-MODE RATE-DEPENDENT DELAMINATION IN
More informationApplication of Chimera Grids in Rotational Flow
CES Seminar Write-up Application of Chimera Grids in Rotational Flow Marc Schwalbach 292414 marc.schwalbach@rwth-aachen.de Supervisors: Dr. Anil Nemili & Emre Özkaya, M.Sc. MATHCCES RWTH Aachen University
More informationSingle Curved Fiber Sedimentation Under Gravity. Xiaoying Rong, Dewei Qi Western Michigan University
Single Curved Fiber Sedimentation Under Gravity Xiaoying Rong, Dewei Qi Western Michigan University JunYong Zhu, Tim Scott USDA Forest Products Laboratory ABSTRACT Dynamics of single curved fiber sedimentation
More informationDEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD R. A. Amaro
More informationEngineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.
Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems
More informationElmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties. SIF file : phasechange solid-solid
Elmer :Heat transfert with phase change solid-solid in transient problem Application to silicon properties 3 6 1. Tb=1750 [K] 2 & 5. q=-10000 [W/m²] 0,1 1 Ω1 4 Ω2 7 3 & 6. α=15 [W/(m²K)] Text=300 [K] 4.
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationNONLINEAR STRUCTURAL DYNAMICS USING FE METHODS
NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.
More informationFLEXIBILITY METHOD FOR INDETERMINATE FRAMES
UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These
More informationNonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess
Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationA multi-solver quasi-newton method for the partitioned simulation of fluid-structure interaction
A multi-solver quasi-newton method for the partitioned simulation of fluid-structure interaction J Degroote, S Annerel and J Vierendeels Department of Flow, Heat and Combustion Mechanics, Ghent University,
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationStructural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).
Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free
More informationSlow Velocity Flow Fields in Composite Materials
Slow Velocity Flow Fields in Composite Materials A Coupled Problem by the Homogenization Method Noboru Kikuchi and His Associates The University of Michigan Ann Arbor, MI 48109, USA Major Contributors
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationFor an imposed stress history consisting of a rapidly applied step-function jump in
Problem 2 (20 points) MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0239 2.002 MECHANICS AND MATERIALS II SOLUTION for QUIZ NO. October 5, 2003 For
More informationNonlinear bending analysis of laminated composite stiffened plates
Nonlinear bending analysis of laminated composite stiffened plates * S.N.Patel 1) 1) Dept. of Civi Engineering, BITS Pilani, Pilani Campus, Pilani-333031, (Raj), India 1) shuvendu@pilani.bits-pilani.ac.in
More information1 Nonlinear deformation
NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected
More informationAdvantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations
The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University
More informationFinite element analysis of membrane wrinkling
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2001; 50:1017 1038 Finite element analysis of membrane wrinkling K. Lu, M. Accorsi ; and J. Leonard Department of Civil
More informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
More informationAEROELASTIC ANALYSIS OF SPHERICAL SHELLS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationChapter 3. Forces, Momentum & Stress. 3.1 Newtonian mechanics: a very brief résumé
Chapter 3 Forces, Momentum & Stress 3.1 Newtonian mechanics: a very brief résumé In classical Newtonian particle mechanics, particles (lumps of matter) only experience acceleration when acted on by external
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationAnalytical formulation of Modified Upper Bound theorem
CHAPTER 3 Analytical formulation of Modified Upper Bound theorem 3.1 Introduction In the mathematical theory of elasticity, the principles of minimum potential energy and minimum complimentary energy are
More informationBENCHMARK LINEAR FINITE ELEMENT ANALYSIS OF LATERALLY LOADED SINGLE PILE USING OPENSEES & COMPARISON WITH ANALYTICAL SOLUTION
BENCHMARK LINEAR FINITE ELEMENT ANALYSIS OF LATERALLY LOADED SINGLE PILE USING OPENSEES & COMPARISON WITH ANALYTICAL SOLUTION Ahmed Elgamal and Jinchi Lu October 07 Introduction In this study: I) The response
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 informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationA consistent dynamic finite element formulation for a pipe using Euler parameters
111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,
More informationComputational Stiffness Method
Computational Stiffness Method Hand calculations are central in the classical stiffness method. In that approach, the stiffness matrix is established column-by-column by setting the degrees of freedom
More informationOPTIMIZATION OF TAPERED WING STRUCTURE WITH AEROELASTIC CONSTRAINT
ICAS 000 CONGRESS OPTIMIZATION OF TAPERED WING STRUCTURE WITH AEROELASTIC CONSTRAINT Harijono Djojodihardjo, I Wayan Tjatra, Ismojo Harjanto 3 Computational Aeroelastic Group, Aerospace Engineering Department
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationNumerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions
TECHNISCHE MECHANIK, 33, 2, (2013), 97 103 submitted: December 11, 2012 Numerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions R. Glüge, M. Weber
More informationFLOW-INDUCED VIBRATIONS OF A LIGHT CIRCULAR CYLINDER AT REYNOLDS NUMBERS 10 3 TO 10 4
Journal of Sound and
More informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
More informationAN ALGORITHM FOR TOPOLOGY OPTIMIZATION
AN ALGORITHM FOR TOPOLOGY OPTIMIZATION OF MULTIBODY SYSTEMS TOHEED GHANDRIZ Master s thesis 2014:E52 Faculty of Science Centre for Mathematical Sciences Numerical Analysis CENTRUM SCIENTIARUM MATHEMATICARUM
More informationNewton Method: General Control and Variants
23 Newton Method: General Control and Variants 23 1 Chapter 23: NEWTON METHOD: GENERAL CONTROL AND VARIANTS TABLE OF CONTENTS Page 23.1 Introduction..................... 23 3 23.2 Newton Iteration as Dynamical
More informationVirtual Work and Variational Principles
Virtual Work and Principles Mathematically, the structural analysis problem is a boundary value problem (BVP). Forces, displacements, stresses, and strains are connected and computed within the framework
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationarxiv: v1 [physics.comp-ph] 27 Nov 2018
Monolithic coupling of implicit material point method with finite element method Eugenio Aulisa Giacomo Capodaglio arxiv:1811.11150v1 [physics.comp-ph] 27 Nov 2018 Abstract A monolithic coupling between
More informationarxiv: v1 [physics.flu-dyn] 6 Nov 2017
A Variational Projection Scheme for Nonmatching Surface-to-Line Coupling arxiv:1711.01699v1 [physics.flu-dyn] 6 Nov 2017 between 3D Flexible Multibody System and Incompressible Turbulent Flow P. S. Gurugubelli
More informationThe Application of a New Material Porosity Algorithm for Parachute Analysis
9 th International LS-DYNA Users Conference Fluid/Structure The Application of a New Material Porosity Algorithm for Parachute Analysis Benjamin Tutt Irvin Aerospace Inc Santa Ana CA 92704 Abstract This
More informationParachute Dynamic Stability and the Effects of Apparent Inertia
Parachute Dynamic Stability and the Effects of Apparent Inertia Jason M. Ginn Georgia Institute of Technology, Atlanta, GA, 333 Ian G. Clark Jet Propulsion Laboratory, California Institute of Technology,
More informationConservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body
Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference
More informationCHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS
61 CHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS 4.1 INTRODUCTION The analysis of cantilever beams of small dimensions taking into the effect of fringing fields is studied and
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationDUCTILITY BEHAVIOR OF A STEEL PLATE SHEAR WALL BY EXPLICIT DYNAMIC ANALYZING
The 4 th World Conference on arthquake ngineering October -7, 008, Beijing, China ABSTRACT : DCTILITY BHAVIOR OF A STL PLAT SHAR WALL BY XPLICIT DYNAMIC ANALYZING P. Memarzadeh Faculty of Civil ngineering,
More information