PAMM · Proc. Appl. Math. Mech. 7, 4040005–4040006 (2007) / DOI 10.1002/pamm.200700122 Energy-momentum schemes for large deformation contact problems Christian Hesch∗1 und P. Betsch ∗∗1 1 Universit¨at Siegen, Paul-Bonatz-Str. 9-11, 57068 Siegen, Germany. Dynamic contact problems in elasticity are dealt with in the framework of nonlinear finite element methods. A new energy- momentum conserving time-stepping scheme for the mortar contact formulation is presented. The proposed method relies on a reparametrization of the contact constraints in terms of specific invariants. For the time discretisation of the contact forces emanating from the mortar formulation the notion of a discrete gradient is applied. © 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 1 Hamiltonian formulation of discrete elastodynamics We start with the space finite element discretisation of nonlinear elastodynamics, which gives rise to a discrete strain energy function, given by V int(q)= W (C)dV (1) B where q represents a possible discrete configuration and C the discrete version of the deformation tensor (right Cauchy-Green tensor). Furthermore, we assume, that the external forces can be derived from an energy potential ext V (q)=− ρ0b · ϕ d V − t · ϕ d A (2) B ∂Bσ where ρ0 denotes the reference mass density, b the applied body force, ϕ the actual configuration and t the prescribed traction boundary condition. Due to the presence of contact constraints, the equations of motion pertaining to the fully discrete system can be written by using a mid-point-type discretisation of the underlying system of differential algebraic equations: ∆t q − q = (v + v ) n+1 n 2 n n+1 m M(v − v )=−∆t∇ V (q , q ) − ∆t (λ ) ∇ Φ (q , q ) (3) n+1 n q n n+1 l n+1 q l n n+1 l=1 0 = Φ(q ) n+1 Here, the discrete gradient is defined as ∇ V (q , q )=Dπ(q )T ∇ V (π(q ), π(q )) q n n+1 n+q q n n+1 (4) with V (π ) − V (π ) −∇V (π 1 ) · (π − π ) n+1 n π n+ 2 n+1 n ∇ V (π(q ), π(q )) = ∇ V (π 1 )+ (π − π ) (5) q n n+1 π n+ 2 2 n+1 n πn+1 − πn using a reparametrisation with possible invariants π, which have to be members of either the set or ¡ (q)={q · q , 1 ≤ A ≤ B ≤ nnodes} A B (6) ¡ (q)={det([q , q ]), 1 ≤ A ≤ B ≤ n } A B nodes Possible invariants of the strain energy function can be identified as the components of the deformation tensor C.Thesame approach can be applied to the contact constraint functions, as shown next. ∗ e-mail: [email protected], Phone: +49 271 740 2101, Fax: +49 271 740 2436 ∗∗ e-mail: [email protected], Phone: +49 271 740 2224, Fax: +49 271 740 2436 © 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim GAMM Sections 4040006 2 Mortar method The reparametrisation of the mortar constraints can be carried out with at least five invariants, three out of and two out of ¡ π (q )=(x(1) − x(1)) · (x(1) − x(1)) 1 seg 2 1 2 1 π (q )=(x(1) − x(1)) · (x(2) − x(1)) 2 seg 2 1 1 1 π (q )=(x(1) − x(1)) · (x(2) − x(1)) 3 seg 2 1 2 1 (7) π (q )=(x(1) − x(1)) · Λ(−2x(1) + x(2) + x(2)) 4 seg 2 1 1 1 2 π (q )=(x(1) − x(1)) · Λ(x(2) − x(2)) 5 seg 2 1 1 2 x(B) Λ where A denotes the four element nodes, defining a mortar segment (see Betsch & Hesch [1]) and is a constant skew- symmetric matrix with 01 Λ = −10 (8) A straightforward calculation shows, that the contact constraint functions can be recast as 1 1 1 Φseq π(q ) = ξ(1) − ξ(1) π ξ(1) − 1 dη + π ξ(2) − ξ(1)ξ(2) dη (9) 1 seg 16 b a 4 5 −1 −1 and 1 1 1 Φseq π(q ) = ξ(1) − ξ(1) π ξ(2) + ξ(1)ξ(2) dη − π ξ(1) +1 dη (10) 2 seg 16 b a 5 4 −1 −1 again defined for one specific mortar segment, restricted by the local parametrisation ξ. A similar approach can be applied to the node-to-segment contact formulation (see Betsch & Hesch [2]). 3 Numerical example The numerical example deals with the planar model of a bearing depicted in Fig. 1. The bearing consists of two rings (Youngs’s 5 modulus E =10, Poissons’s ratio ν =0.1 and mass density R =0.001), which are discretized by 4-node isoparametric displacement-based plain strain elements. The discretization of the outer ring relies on 10x48 elements, for the inner ring 10x40 have been used. For t ∈ [0, 0.5], a torque acts on the inner ring in form of a hat function over time. Then, for t ∈ (0.5, 2], no external loads are acting on the bearing anymore. Fig. 2 shows that for t ≥ 0.5 the present scheme does indeed conserve the total energy for the frictionless contact problem under consideration. 100 6000 5000 50 4000 0 E 3000 kinetic energy 2000 potential energy −50 total energy 1000 −100 −100 −50 0 50 100 0 0.5 1 1.5 2 t Fig. 1: Discretized bearing Fig. 2: Energy versus time Literatur [1] C. Hesch and P. Betsch, On the energy-momentum conserving integration of large deformation contact problems. In Proceedings of the ECCOMAS Thematic Conference, Mailand, 25–28.06.2007 [2] P. Betsch and C. Hesch, Energy-momentum conserving schemes for frictionless dynamic contact problems. Part I: NTS method IUTAM Bookseries, Springer-Verlag, 3:77–96, 2007. © 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.