a mixed finite element model for the elastic contact problem the static elastic contact problem is approached using lagrange multipliers leading to a mixed finite element problem a non linear friction law is introduced explicitly and the non local character of the friction phenomena is implicitly assumed in order to avoid stress oscillations near singular points a perturbed lagrangian functional is considered the algorithms herein proposed do not impose nodal dependencies over the contact surfaces allowing for the independent discretization of both bodies the method is able to model simultaneous contact over different regions of any geometrical shape computer code examples and results presented here are restricted to axisymmetrical and bidimensional cases

interface element method iem for a partitioned system with non matching interfaces abstract a novel method for non matching interfaces on the boundaries of finite elements in partitioned domains is presented by introducing interface elements in this paper the interface element method iem satisfies the continuity conditions exactly through interface element domains without recourse to lagrange multiplier technique the interface elements take an interpolation scheme used in meshless methods to join partitioned domains the moving least square mls approximations in the present study are implemented to construct the shape functions of the interface elements the compatibility of displacements on the boundaries of the finite elements and the interface elements is always preserved in this method and the completeness of the shape functions of the interface elements transfers rigid body motions and strain fields correctly moreover alignment of the boundaries of sub domains in the mls approximations and integration domains provides a consistent numerical integration the numerical examples show that the iem is a useful tool for the analysis of a partitioned system and for a global local analysis

an algebraically partitioned feti method for parallel structural analysis algorithm description an algebraically partitioned feti method for the solution of structural engineering problems on parallel computers is presented the present algorithm consists of three attributes an explicit generation of the orthogonal null space matrix associated with the interface nodal forces the floating subdomain rigid body modes computed from the subdomain static equilibrium equation of the classical force method and the identification of redundant interface force constraint operator that emanates when the interface force computations are localized comparisons of the present method with the previously developed differentially partitioned feti method are offered in terms of the saddle point formulations at the end of the paper a companion paper reports implementation details and numerical performance of the proposed algorithm 1997 john wiley sons ltd

a method of finite element tearing and interconnecting and its parallel solution algorithm a novel domain decomposition approach for the parallel finite element solution of equilibrium equations is presented the spatial domain is partitioned into a set of totally disconnected subdomains each assigned to an individual processor lagrange multipliers are introduced to enforce compatibility at the interface nodes in the static case each floating subdomain induces a local singularity that is resolved in two phases first the rigid body modes are eliminated in parallel from each local problem and a direct scheme is applied concurrently to all subdomains in order to recover each partial local solution next the contributions of these modes are related to the lagrange multipliers through an orthogonality condition a parallel conjugate projected gradient algorithm is developed for the solution of the coupled system of local rigid modes components and lagrange multipliers which completes the solution of the problem when implemented on local memory multiprocessors this proposed method of tearing and interconnecting requires less interprocessor communications than the classical method of substructuring it is also suitable for parallel vector computers with shared memory moreover unlike parallel direct solvers it exhibits a degree of parallelism that is not limited by the bandwidth of the finite element system of equations

modelling unilateral frictionless contact using the null space method and cubic b spline interpolation the analysis of unilateral sliding contact in elasticity is equivalent to a minimisation problem subjected to a set of inequality constraints however the presence of boundary discontinuities such as those stemming from the spatial discretisation appears as a major problem to determine the set of active constraints this work introduces a smoothing technique of the master surface resorting to cubic b spline interpolation which is c 1 continuous in contact situations between elastic and rigid bodies and g 1 continuous in elastic elastic contact problems the technique is applied in conjunction with the null space method where the solution is searched in an unconstrained manifold the resulting formulation eases the contact transition along the master surface and recovers the quadratic convergence of the iterative newton raphson process the robustness of the method is demonstrated using 2d and 3d examples

a unified formulation for interface coupling and frictional contact modeling with embedded error estimation summary we present a derivation of a new interface formulation via a merger of continuous galerkin and discontinuous galerkin concepts enhanced by the variational multiscale method developments herein provide treatment for the pure displacement form and mixed form of small deformation elasticity as applied to the solution of two problem classes domain decomposition and contact mechanics with friction the proposed framework seamlessly accommodates merger of different element types within subregions of the computational domain and nonmatching element faces along embedded interfaces these features are retained in the treatment of multibody small deformation contact problems as well where an unbiased treatment of the contact interface stands in contrast to classical master slave constructs numerical results for problems in two and three spatial dimensions illustrate the robustness and versatility of the proposed method copyright 2012 john wiley sons ltd

stabilized interface methods for mechanical joints physics based models and variationally consistent embedding abstract this paper presents the application of a new method for interfacial modeling utilizing a merger of continuous galerkin and discontinuous galerkin concepts to simulate the behavior of mechanical joints the interfacial flux terms arising naturally from the discontinuous galerkin treatment provide a mechanism to embed friction models in a variationally consistent fashion due to the unbiased implementation of the interface facilitated by avoiding the master slave concept the deformation of the two interacting surfaces conforms to the local material and geometric attributes of the surfaces this results in a better preservation of physics in interface mechanics additionally the method is incorporated into a variational multiscale framework that comes equipped with a built in error estimation module providing numerical estimation of convergence and distinguishing discretization errors from modeling errors a series of quasi static numerical simulations of a lap joint under fretting conditions are conducted to compare the performance of two friction models i classical coulomb friction model and ii physics based multiscale model hysteresis study of a three dimensional double bolted lap joint for the two friction models is also presented and the computed results are shown to be consistent between conforming and nonconforming meshes

mixed fem of higher order for time dependent contact problems abstract in this paper mixed finite element methods of higher order for time dependent contact problems are discussed the mixed methods are based on resolving the contact conditions by the introduction of lagrange multipliers dynamic signorini problems with and without friction are considered involving thermomechanical and rolling contact rothe s method is used to provide a suitable time and space discretization to discretize in time a stabilized newmark method is applied as an adequate time stepping scheme the space discretization relies on finite elements of higher order in each time step the resulting problems are solved by uzawa s method or alternatively by methods of quadratic programming via a suitable formulation in terms of the lagrange multipliers numerical results are presented towards an application in production engineering the results illustrate the performance of the presented techniques for a variety of problem formulations

mls based variable node elements compatible with quadratic interpolation part i formulation and application for non matching meshes two dimensional variable node elements compatible with quadratic interpolation are developed using the moving least square mls approximation the mapping from the parental domain to the physical element domain is implicitly obtained from mls approximation with the shape functions and their derivatives calculated and saved only at the numerical integration points it is shown that the present mls based variable node elements meet the patch test if a sufficiently large number of integration points are employed for numerical integration the cantilever problem with non matching meshes is chosen to check the feasibility of the present mls based variable node elements and the result is compared with that from the lower order case compatible with linear interpolation copyright 2005 john wiley sons ltd

a mortar segment to segment contact method for large deformation solid mechanics abstract contact modeling is still one of the most difficult aspects of non linear implicit structural analysis most 3d contact algorithms employed today use node on segment approaches for contacting dissimilar meshes two pass node on segment contact approaches have the well known deficiency of locking due to over constraint in this work we develop and demonstrate a segment to segment contact approach based on the mortar method by eliminating the over constraint this method appears to be much more robust than the node on segment approach furthermore node on segment approaches suffer when individual nodes slide out of contact at contact surface boundaries or when contacting nodes slide from facet to facet this causes jumps in the contact forces due to the discrete nature of the constraint enforcement and difficulties in convergence for implicit solution techniques these jumps in forces are avoided in the segment to segment approach since penetration is measured by the occluded weighted volume which varies smoothly with perturbations until now all mortar schemes have been developed for small deformation sliding and have mainly been applied on flat surfaces here we will present the integration scheme and the linearization approach used for handling the contact surface integrals in the event of large sliding the schemes are designed to conserve linear momentum and can be applied to arbitrarily curved 3d surfaces although expensive our numerical examples demonstrate that the approach is far more robust than even smoothed node on segment contact techniques that is problems previously not solvable with node on segment can now be solved with our segment to segment contact

a two step two field hybrid method for the static and dynamic analysis of substructure problems with conforming and non conforming interfaces abstract the need for assembling independent finite element substructure solutions arises in several engineering and scientific problems including the design and analysis of complex structural systems component mode synthesis global local analysis adaptive refinement and parallel processing in this paper we discuss the solution of such problems by a two field hybrid method where the substructures are jointed with low order polynomial or piece wise polynomial lagrange multipliers and present a rayleigh ritz based smoothing procedure for improving the accuracy of the computed coupled solution in the presence of various substructure heterogeneities we consider both conforming and nonconforming substructure meshes and demonstrate the benefits of the proposed two step solution method with several examples from structural mechanics

finite deformation frictional mortar contact using a semi smooth newton method with consistent linearization a two dimensional finite deformation frictional contact formulation with coulomb s law is presented the approach considers multibody contact and is based on a mortar formulation the enforcement of contact constraints is realized with dual lagrange multipliers these alternative multiplier spaces are constructed in a way that the multipliers can easily be eliminated from the global system of equations by static condensation such that the system size does not increase friction kinematic variables are formulated in an objective way and enter non smooth complementarity functions for expressing the contact constraints an active set strategy is derived by applying a semismooth newton method which treats contact nonlinearities material and geometrical nonlinearities in one single iterative scheme by further carrying out a consistent linearization for both normal and frictional contact forces and constraints a robust and highly efficient algorithm for linear and higher order quadratic interpolation is achieved efficiency of the proposed method and quality of results are demonstrated in serveral exemples

consistent treatment of boundaries with mortar contact formulations using dual lagrange multipliers abstract computational modelling of contact problems raises two basic questions which method should be used to enforce the contact conditions and how should this method be discretised the most popular enforcement methods are the lagrange multiplier method the penalty method and combinations of these two a frequently used discretisation method is the so called node to segment approach however this approach might lead to problems like jumps in contact forces loss of convergence or failure to pass the patch test thus in the last few years several segment to segment contact algorithms based on the mortar method were proposed combination of a mortar discretisation with a penalty based enforcement of the contact conditions leads to unphysical penetrations on the other hand a lagrange multiplier mortar method requires additional unknowns hence condensation of the lagrange multipliers is desirable to preserve the initial size of the system of equations this can be achieved by interpolating the lagrange multipliers with so called dual shape functions discretising two contacting bodies leads to opposed contact surface representations of finite element edges called slave and master elements respectively in current versions of dual lagrange multiplier mortar formulations an inconsistency at the boundary appears when only a part of a slave element instead of the entire element belongs to the contact area we present a modified definition of the dual shape functions in such slave elements the basic idea is to construct dual shape functions that fulfill the so called biorthogonality condition within the contact area this leads to consistent mortar matrices also in the boundary region to avoid ill conditioning of the stiffness matrix the modified mortar matrices are weighted with appropriate weighting factors in doing so the corresponding modified lagrange multiplier nodal values are of the same order as the unmodified ones various examples demonstrate the performance of the modified mortar contact algorithm

two dimensional mortar contact methods for large deformation frictional sliding this paper presents a mortar based formulation for the solution of two dimensional frictional contact problems involving finite deformation and large sliding as is widely recognized traditional node to surface contact formulations have several drawbacks in solution of deformable to deformable contact problems including lack of general patch test passage degradation of spatial convergence rates and robustness issues associated with the faceted representation of contacting surfaces the mortar finite element method initially proposed as a technique to join dissimilarly meshed domains has been shown to preserve optimal convergence rates in tied contact problems see discretization methods and iterative solvers based on domain decomposition springer verlag heidelberg 2001 for a recent review and is examined here as an alternative spatial discretization method for large sliding contact in particular a novel description for frictional sliding conditions in large deformation mortar formulations is proposed in this work in recent years the mortar element method has already been successfully implemented to solve frictional contact problems with linearized kinematics see int j numer meth engng 1993 36 3451 however in the presence of large deformations and finite sliding one must face difficulties associated with the definition and linearization of contact virtual work in the case where the mortar projection has a direct dependence on the tangential relative motion along the interface in this paper such a formulation is presented with particular emphasis on key aspects of the linearization procedure and on the robust description of the friction kinematics some novel techniques are proposed to treat the non smoothness in the contact geometry and the searching required to define mortar segments a number of numerical examples illustrate the performance and accuracy of the proposed formulation copyright 2004 john wiley sons ltd

an unbiased nitsche s approximation of the frictional contact between two elastic structures most of the numerical methods dedicated to the contact problem involving two elastic bodies are based on the master slave paradigm it results in important detection difficulties in the case of self contact and multi body contact where it may be impractical if not impossible to a priori nominate a master surface and a slave one in this work we introduce an unbiased finite element method for the finite element approximation of frictional contact between two elastic bodies in the small deformation framework in the proposed method the two bodies expected to come into contact are treated in the same way no master and slave surfaces the key ingredient is a nitsche based formulation of contact conditions as in chouly et al math comput 84 1089 1112 2015 we carry out the numerical analysis of the method and prove its well posedness and optimal convergence in the h 1 norm numerical experiments are performed to illustrate the theoretical results and the performance of the method

a contact smoothing method for arbitrary surface meshes using nagata patches abstract this paper presents a contact surface smoothing method combined with the node to segment discretization technique to solve large deformation frictional contact problems between deformable bodies the nagata patch interpolation is used to smooth the surface mesh providing a master surface with quasi g 1 continuity between patches moreover the local support of the interpolation method allows to deal with surface meshes of arbitrary topology regular and irregular finite element discretizations as well as hybrid meshes the non physical oscillations in the contact force evolution induced by the faceted contact surface representation are reduced using the proposed smoothing method furthermore the smooth representation of the master surface allows a more accurately evaluation of the resulting stresses and forces while providing an important improvement in convergence behaviour four representative numerical examples are used to demonstrate the advantages of the proposed contact smoothing method the results show a significant improvement in the accuracy robustness and performance of the numerical simulations using the smoothing approach when compared with the piecewise faceted contact surface description

a mortar approach for fluid structure interaction problems immersed strategies for deformable and rigid bodies abstract this paper proposes a new fluid structure interaction immersed computational framework the coupling between an underlying incompressible fluid and an embedded solid is formulated by means of the overlapping domain decomposition method in conjunction with a mortar approach leading to a variationally consistent scheme which is capable of unifying a range of methodologies currently available in the literature this novel framework provides great flexibility and enables the modelling of immersed deformable solids compressible and incompressible as well as rigid bodies through the use of a weak director based formulation a novel null space reduction scheme is employed in order to enhance the conditioning of the resulting system of equations and reduce the computational cost an implicit structure preserving time integration algorithm is used to yield extra stability and robustness and the use of a segmentation technique near the boundary between fluid and solid also leads to enhanced accuracy the methodology is benchmarked against results obtained by using alternative boundary fitted methodologies

quotients exactness and nuclearity in the operator system category abstract we continue our study of tensor products in the operator system category we define operator system quotients and exactness in this setting and refine the notion of nuclearity by studying operator systems that preserve various pairs of tensor products one of our main goals is to relate these refinements of nuclearity to the kirchberg conjecture in particular we prove that the kirchberg conjecture is equivalent to the statement that every operator system that is min er nuclear is also el c nuclear we show that operator system quotients are not always equal to the corresponding operator space quotients and then study exactness of various operator system tensor products for the operator system quotient we prove that an operator system is exact for the min tensor product if and only if it is min el nuclear we give many characterizations of operator systems that are min er nuclear el c nuclear min el nuclear and el max nuclear these characterizations involve operator system analogues of various properties from the theory of c algebras and operator spaces including the wep and llp

a formulation based on localized lagrange multipliers for bem fem coupling in contact problems this paper presents a unified formulation for the combination of the finite element method fem and the boundary element method bem in 3d frictional contact problems that is based on the use of localized lagrange multipliers llms resolution methods for the contact problem between non matching meshes have traditionally been based on a direct coupling of the contacting solids using classical lagrange multipliers these methods tend to generate strongly coupled systems that require a deep knowledge of the discretization characteristics on each side of the contact zone complicating the process of mixing different numerical techniques in this work a displacement contact frame is inserted between the fe and be interface meshes discretized and finally connected to the contacting substructures using llms collocated at the mesh interface nodes this methodology will provide a partitioned formulation which preserves software modularity and facilitates the connection of non matching fe and be meshes

algorithms for the analysis of 3d finite strain contact problems from the constraint imposition aspects in 3d to friction regularization various ideas are exposed in this paper a variation of the rockafellar lagrangian is proposed which results in continuous second order derivatives if lagrange multiplier estimates are greater or equal than one this fact allows the adoption of a full second order i e lagrange newton method avoiding sequential unconstrained minimization techniques an algorithm for global and local contact detection is presented which is developed for dealing with large step sizes typical of implicit methods a modified constraint definition to deal with non smooth situations is presented aspects of friction implementation including a regularization scheme which ensures stepwise objectivity are detailed finally several illustrative examples are carried out with success

a point to segment contact formulation for isogeometric nurbs based finite elements formulation of isogeometric finite elements has received a great deal of attention in the recent past the present study deals with the treatment of problems in structural mechanics including large elastic deformations and contact one decisive difference between isogeometric finite elements based on nurbs functions and standard finite elements fe based on lagrange polynomials is higher inter element continuity this is a promising characteristics to model contact problems as all issues associated with kinks between elements are naturally avoided particularly in cases with large sliding this has the potential to be an attractive feature we present a bilateral isogeometric collocation contact formulation for geometrically non linear two dimensional problems using greville and botella points to collocate the contact integrals different methods to obtain accurate and physically meaningful stress distributions are investigated and compared the results show that the higher inter element continuity in context of non smooth problems may imply undesired effects in the numerical solution such as unphysical stress oscillations similar effects have been observed in standard p fem numerical experiments indicate that these oscillations may be avoided if the basis functions of contact and non contact zones are separated by knot relocation and knot repetition

a monolithic smoothing gap algorithm for contact impact based on the signed distance function a new algorithm has been developed for smoothing the surfaces in finite element formulations of contact impact a key feature of this method is that the smoothing is done implicitly by constructing smooth signed distance functions for the bodies these functions are then employed for the computation of the gap and other variables needed for implementation of contact impact the smoothed signed distance functions are constructed by a moving least squares approximation with a polynomial basis results show that when nodes are placed on a surface the surface can be reproduced with an error of about one per cent or less with either a quadratic or a linear basis with a quadratic basis the method exactly reproduces a circle or a sphere even for coarse meshes results are presented for contact problems involving the contact of circular bodies copyright c 2002 john wiley sons ltd

unilateral non linear dynamic contact of thin walled structures using a primal dual active set strategy the efficient modelling of three dimensional contact problems is still a challenge in non linear implicit structural analysis we use a primal dual active set strategy siam j optim 2003 13 865 888 based on dual lagrange multipliers siam j numer anal 2000 38 989 1012 to handle the non linearity of the contact conditions this allows us to enforce the contact constraints in a weak integral sense without any additional parameter due to the biorthogonality condition of the basis functions the lagrange multipliers can be locally eliminated we perform a static condensation to achieve a reduced system for the displacements the lagrange multipliers representing the contact pressure can be easily recovered from the displacements in a variationally consistent way for the application to thin walled structures we adapt a three dimensional non linear shell formulation including the thickness stretch of the shell to contact problems a reparametrization of the geometric description of the shell body gives us a surface oriented shell element which allows the application of contact conditions directly to nodes lying on the contact surface shell typical locking phenomena are treated with the enhanced assumed strain method and the assumed natural strain method the discretization in time is done with the implicit generalized method j appl mech 1993 60 371 375 and the generalized energy momentum method comp methods appl mech eng 1999 178 343 366 to compare the development of energies within a frictionless contact description in order to conserve the total energy within the discretized frictionless contact framework we follow an approach from laursen and love int j numer methods eng 2002 53 245 274 who introduced a discrete contact velocity to update the velocity field in a post processing step various examples show the good performance of the primal dual active set strategy applied to the implicit dynamic analysis of thin walled structures

an interior point method for isogeometric contact abstract the interior point method is applied to frictionless contact mechanics problems and is shown to be a viable alternative to the augmented lagrangian approach the method is derived from a mixed formulation which induces a contact discretization scheme in the spirit of the mortar method and naturally delivers slack variables that help constrain the solution to the feasible region the derivation of the algorithm as well as its robustness benefits from the contact interface description that is induced by nurbs based isogeometric volume discretizations various interior point algorithms are discussed including a primal dual approach that satisfies the unilateral contact constraints exactly in addition to two primal approaches that retain an arbitrary barrier parameter the developed algorithms can easily be pursued starting from an augmented lagrangian implementation numerical investigations on benchmark problems demonstrate the efficiency and the robustness of the framework but also highlight current limitations that suggest paths for future research overall the results indicate that the interior point method can challenge the augmented lagrangian method in contact mechanics even displaying potential for higher efficiency and robustness

resolution of fourth order problems by the mortar element method we present a resolution of fourth order problems by the mortar element method in the case of spectral discretizations we prove that this approach leads to a well posed discrete problem and the numerical analysis gives an optimal error estimate

a new computational approach to contact mechanics using variable node finite elements in this paper a new computational strategy for two dimensional contact problems is developed with the aid of variable node finite elements within the range of infinitesimal deformations the variable node elements which are among mls moving least square based finite elements enable us to transform node to surface contact problems into node to node contact problems this contact formulation with variable node elements leads to an accurate and effective solution procedure needless to mention that the contact patch test is passed without any additional treatment through several numerical examples we demonstrate its simplicity and the effectiveness of the proposed scheme

completely bounded module maps and the haagerup tensor product abstract we obtain a representation theorem for module maps defined on the algebra k h of compact operators this is used to show that the haagerup tensor product of two operator algebras enjoys strong versions of tomita s commutant theorem and the slice map property we also give a general theorem concerning the automatic complete boundedness of module maps

a formulation of conserving impact system based on localized lagrange multipliers this paper proposes a new finite element method of frictional impact of elastic bodies the formulation introduces a contact frame that is placed in between contacting bodies and represents the contact surface the nonpenetration condition and the slip stick condition are defined between the contacting body and the contact frame with the aid of the independent localized lagrange multipliers representing the contact forces the position of the contact frame and the local coordinate of the contacting node along the contact frame are also treated as the independent variable which enables the exact satisfaction of the constraint conditions without the deficiency or redundant constraint the energy and momentum conservation algorithm is applied to the proposed impact system for the case of frictional impact the linear momentum is exactly conserved and the angular momentum is approximately conserved with negligible error copyright 2006 john wiley sons ltd

a large deformation mortar formulation of self contact with finite sliding this paper presents a new numerical method in which self contact phenomena associated with a body undergoing large deformations and sliding can be described in particular the approach relies on a particular extension of the mortar approach appropriate for this class of problems a bounding volume hierarchy organized as a binary tree is built for the self contact surface based on the geometry and the mesh connectivity of the surface a curvature criterion using a new algorithm to detect subsurface adjacency is used to accelerate the self contact searching procedure to ensure that the mortar traction fields are properly defined on contiguous surface patches a novel facet sorting algorithm is also proposed based on the mesh connectivity of the contact element pairs found by the self contact searching algorithm several two and three dimensional numerical examples show the new self contact mortar formulation to be very efficient and also demonstrate that it can be combined with multi body contact algorithms to simulate a very general class of contact problems

an improved interface element with variable nodes for non matching finite element meshes extensive improvements of the interface element method iem are proposed for the efficient treatment of non matching finite element meshes our approach enables us to establish the master element via the moving least square mls approximation and so to remove the cumbersome process of constructing interface elements the values of shape functions and their derivatives are therefore mapped from the master element as in the conventional finite element method for the assurance of convergence and compatibility condition a patch test is demonstrated through several examples of 2d linear elasticity the convergence rate is compared between the present interface element and the previous version

structurally balanced controller order reduction with guaranteed closed loop performance abstract this paper considers controller order reduction with the consideration of closed loop stability and performance a closed loop structurally balancing method is proposed to reduce the controller orders for a class of systems and controllers the resulting reduced order controllers if they exist are guaranteed to stabilize the corresponding systems error bounds are also derived for the closed loop performance degradation

a spectral element methodology tuned to parallel implementations abstract we present in this paper a modification of the conforming spectral element method for the approximation of the poisson equation this new procedure allows the implementation on a parallel machine in an easier and faster way the resulting method is nonconforming but the numerical analysis shows that the error between the exact solution and the discrete solution is of the same asymptotic order as the error of the original conforming method

on some aspects for contact with rigid surfaces surface to rigid surface and curves to rigid surface algorithms abstract special algorithms allowing a simplified description of contact between deformable body and rigid surfaces are developed based on the geometrically exact covariant description of contact a special attention is given to various geometric combinations where the contact can be represented as a contact between surfaces and b contact between a curve and a surface for contact between surfaces leading to the segment to analytical surface stas approach two algorithms can be distinguished based on the selection of a coordinate system for the closest point projection cpp procedure a rigid surface is a slave surface and b rigid surface is a master surface a special combination of both contact kinematics for the surface to surface and for the curve to curve approaches is employed for the contact between a curve and a surface leading to the curve to rigid surface ctrs approach the last algorithm is verified with the well known euler formula for the rope cylinder interaction as well as with a new derived generalization into a 3d spiral rope on a cylinder the developed algorithms can be straightforwardly implemented within an iso geometric approach as well as within the conventional finite elements where rigid surfaces are given by cad patches any type of elements can be employed for the contacting deformable surface curve because the algorithms are formulated in a covariant form

a contact formulation based on localized lagrange multipliers formulation and application to two dimensional problems the non penetration condition in contact problems is traditionally based on the classical lagrange multiplier method this method makes extensive use of modelling details of the contacting bodies for contact enforcement as the contact surface meshes are in general non matching to deal with this problem we introduce a novel element in the lagrange multiplier approach of contact modelling namely a contact frame placed in between contacting bodies it acts as a medium through which contact forces are transferred without violating equilibrium in the contact domain for discrete contact models only nodal information of the contacting bodies is required which makes the proposed contact enforcement generic the contact frame has its own independent freedoms which allows the formulation to pass contact patch tests by design copyright 2002 john wiley sons ltd

a segment to segment mortar contact method for quadratic elements and large deformations abstract this paper describes a new extension of mortar based segment to segment algorithms for large deformation contact analysis permitting their use when quadratic interpolations are used in the underlying finite element method the work builds upon results previously described by puso and laursen m a puso t a laursen a mortar segment to segment contact method for large deformation solid mechanics comput methods appl mech engrg 193 2004 601 629 m a puso t a laursen a mortar segment to segment frictional contact method for large deformations comput methods appl mech engrg 193 2004 601 629 where methods were developed for quasistatic and implicit dynamic analysis in the context of linear element interpolations here we see that many algorithmic ideas from the linear case can be adopted for use with quadratic elements along with different alternative interpolations of the mortar contact multipliers the result is a class of methods appropriate for use for either mesh tying true large sliding simulation frictionless or frictional and readily admits mixing of element types several numerical examples are given to demonstrate the effectiveness of the approach

a mortar method for energy momentum conserving schemes in frictionless dynamic contact problems summary in the present work the mortar method is applied to planar large deformation contact problems without friction in particular the proposed form of the mortar contact constraints is invariant under translations and rotations these invariance properties lay the foundation for the design of energy momentum timestepping schemes for contact impact problems the iterative solution procedure is embedded into an active set algorithm lagrange multipliers are used to enforce the mortar contact constraints the solution of generalized saddle point systems is circumvented by applying the discrete null space method numerical examples demonstrate the robustness and enhanced numerical stability of the newly developed energymomentum scheme copyright q 2008 john wiley sons ltd

a finite deformation mortar contact formulation using a primal dual active set strategy in recent years nonconforming domain decomposition techniques and in particular the mortar method have become popular in developing new contact algorithms here we present an approach for two dimensional frictionless multibody contact based on a mortar formulation and using a primal dual active set strategy for contact constraint enforcement we consider linear and higher order quadratic interpolations throughout this work so called dual lagrange multipliers are introduced for the contact pressure but can be eliminated from the global system of equations by static condensation thus avoiding an increase in system size for this type of contact formulation we provide a full linearization of both contact forces and normal non penetration and tangential frictionless sliding contact constraints in the finite deformation frame the necessity of such a linearization in order to obtain a consistent newton scheme is demonstrated by further interpreting the active set search as a semi smooth newton method contact nonlinearity and geometrical and material nonlinearity can be resolved within one single iterative scheme this yields a robust and highly efficient algorithm for frictionless finite deformation contact problems numerical examples illustrate the efficiency of our method and the high quality of results

a high order generalized fem for through the thickness branched cracks this paper presents high order implementations of a generalized finite element method for through the thickness three dimensional branched cracks this approach can accurately represent discontinuities such as triple joints in polycrystalline materials and branched cracks independently of the background finite element mesh representative problems are investigated to illustrate the accuracy of the method in combination with various discretizations and refinement strategies the combination of local refinement at crack fronts and high order continuous and discontinuous enrichments proves to be an excellent combination which can deliver convergence rates close to that of problems with smooth solutions copyright 2007 john wiley sons ltd

a non mortar mixed finite element method for elliptic problems on non matching multiblock grids abstract we consider the approximation of second order elliptic equations on domains that can be described as a union of sub domains or blocks we assume that a grid is defined on each block independently so that the resulting grid over the entire domain need not be conforming i e match across the block boundaries several techniques have been developed to approximate elliptic equations on multiblock grids that utilize a mortar finite element space defined on the block boundary interface itself we define a mixed finite element method that does not use such a mortar space the method has an advantage in the case where adaptive local refinement techniques will be used in that there is no mortar grid to refine as is typical of mixed methods our method is locally conservative element by element it is also globally conservative across the block boundaries theoretical results show that the approximate solution converges at the optimal rate to the true solution we present computational results to illustrate and confirm the theory

injectivity of operator spaces we study the structure of injective operator spaces and the existence and uniqueness of the injective envelopes of operator spaces we give an easy example of an injective operator space which is not completely isometric to any c algebra this answers a question of wittstock 23 furthermore we show that an operator space e is injective if and only if there exists an injective c algebra a and two projections p and q in a such that e is completely isometric to paq

a special focus on 2d formulations for contact problems using a covariant description a fully covariant description based on the consideration of contact conditions especially for the 2d case is proposed the description is based on a reconsideration of contact kinematics and all necessary operations such as derivatives in a specially chosen curvilinear coordinate system based on a curved geometry in plane in addition details of the finite element implementation are presented for the simple linear contact element special cases requiring the update of history variables as well as their careful transfer over the element boundaries are illustrated by numerical examples with these procedures artificial jumps in the contact forces can be avoided copyright 2005 john wiley sons ltd

nonconforming mortart element methods for the spectral discretization of two dimensional fourth order problems this paper deals with two versions of the mortar element method applied to the spectral discretization of two dimensional fourth order problems in each case the analysis of the discrete problem is carried out the key results are the uniform ellipticity of the discrete bilinear form leading to error estimates which are optimal with respect to the discretization parameter and the regularity of the exact solution some numerical results are given to illustrate various aspects of the method

a perturbed lagrangian formulation for the finite element solution of contact problems making use of a perturbed lagrangian formulation a finite element procedure for contact problems is developed for the general case in which node to node contact no longer holds the proposed procedure leads naturally to a discretization of the contact interface into contact segments within the context of a bilinear interpolation for the displacement field a mixed finite element approximation is introduced by assuming discont naous contact pressure constant on the contact segment because of this piece wise constant approximation the gap function enters into the formulation in an average sense instead of through a point wise definition numerical examples are presented that illustrate the performance of the proposed procedure current finite element formulations for contact problems based on either the classical lagrange parameter procedure l 3 12 203 or the penalty function method 4 6 11 are characterized by a point wise enforcement of the contact constraint condition in the sense that penetration of the bodies is established on a nodal basis moreover in this methods the recovery of the contact pressure over the element from the contact nodal forces generally requires an additional procedure within the framework of classical lagrange multiplier methods the contact condition is exactly satisfied by transforming the constrained problem into an unconstrained one with the introduction of additional variables lagrange multipliers these extra variables add computational effort to the solution process which often requires special procedures to handle the presence of zero diagonal terms penalty methods on the other hand enable one to transform the constrained problem into an unconstrained one without introducing additional variables the constraint condition is now satisfied only approximately for finite values of the penalty parameter the main difhculty associated with these methods however lies in the poor conditioning of the problem as the penalty is increased to more accurately enforce the constraint condition this is a well understood phenomenon particularly in the context of the incompressible and nearly incompressible problem in solid and fluid mechanics e g see l 22 25 for a review recently augmented lagrangian procedures have been proposed as a promising way to partially overcome these difficulties and regularize the penalty formulation e g see the survey in 7 8 within the framework of zinearited e u c it is possible to restrict the finite element formulation of contact problems by assuming that node to node contact occurs this is in fact the case often considered in the literature l 2 4 6 12 20 211 in the general context of fully

consistent mesh tying methods for topologically distinct discretized surfaces in non linear solid mechanics this work presents a new technique for joining dissimilar meshes in non linear finite element analysis with the novel result that the adjoining surfaces do not need to be spatially coincident for patch tests to be passed the work builds in some sense on results previously presented by dohrmann et al international journal for numerical methods in engineering 1991 40 1205 but considerably generalizes upon the results presented there most importantly it is shown that consistency of the tied finite element formulation may be achieved such that implementation of element formulations is not affected through the use of a three field variational principle applied over the contacting region a key conclusion of the work is that although the approach advocated bears some resemblance to mortar element schemes for contact analysis see for example mcdevitt and laursen international journal for numerical methods in engineering 2000 48 1525 a straightforward application of a mortaring scheme is not sufficient to guarantee patch test passage when dissimilar meshing of curved surfaces leads to local open volumes or overlaps between tied surfaces in their reference configurations the proposed algorithm is demonstrated through a number of numerical examples performed within a matrix free non linear equation solving framework published in 2003 by john wiley sons ltd

a consistent tangent stiffness matrix for three dimensional non linear contact analysis a consistent tangent stiffness matrix for the analysis of non linear contact problems is presented the associated element has three or four nodes and establishes contact between three dimensional structures like solids and shells it accounts for the non linear kinematics of large deformation analysis and guarantees a quadratic convergence rate two formulations the penalty method and the lagrange multiplier method are investigated

an unbiased nitsche s formulation of large deformation frictional contact and self contact in this paper we propose an extension of nitsche s method for frictional contact in large elastic deformations in fact we develop an unbiased strategy in which no asymmetry is considered between contact surfaces conversely to the master slave strategy this enables to take into account within the same formalism contact of two elastic bodies as well as self contact we provide a numerical study of the performance and robustness of the method

node to node scheme for three dimensional contact mechanics using polyhedral type variable node elements abstract a node to node contact scheme which is applicable for three dimensional contact analyses involving large deformations is developed with the aid of polyhedral elements the key issue is to transform nonmatching meshes into matching meshes in a seamless manner because polyhedral elements are allowed to have arbitrary numbers of polygonal faces and nodes they can be used as transition elements for coupling nonmatching meshes in this paper the polyhedral elements make it possible to always maintain node to node contact during the contact deformation the present approach guarantees that the patch test is passed and the nonpenetration condition is satisfied and hence it yields smoother contact pressure with faster convergence than the conventional node to surface or surface to surface contact scheme

a stable and optimally convergent latin cutfem algorithm for multiple unilateral contact problems in this paper we propose a novel unfitted finite element method for the simulation of multiple body contact the computational mesh is generated independently of the geometry of the interacting solids which can be arbitrarily complex the key novelty of the approach is the combination of elements of the cutfem technology namely the enrichment of the solution field via the definition of overlapping fictitious domains with a dedicated penalty type regularisation of discrete operators and the latin hybrid mixed formulation of complex interface conditions furthermore the novel p1 p1 discretisation scheme that we propose for the unfitted latin solver is shown to be stable robust and optimally convergent with mesh refinement finally the paper introduces a high performance 3d level set cutfem framework for the versatile and robust solution of contact problems involving multiple bodies of complex geometries with more than two bodies interacting at a single point

a mixed formulation of mortar based frictionless contact abstract a class of mortar based frictionless contact formulations is derived based on a classical three field mixed variational framework within a penalty regularization complemented by uzawa augmentations discrete mortar constraints are naturally induced by the variational setting major aspects of earlier mortar approaches are obtained through constrained lumped or unconstrained recovery procedures for the mixed kinematic and kinetic mortar quantities from their projected counterparts two and three dimensional examples at the infinitesimal and finite deformation regimes highlight the local and global quality of the contact interactions

finite elements in space and time for dynamic contact problems dynamic obstacle and signorini problems are discretized by continuous and discontinuous finite elements in space and time the resulting discrete problems are attributed to a standard implicit time stepping scheme through relaxation of impact phenomena and suitable numerical integration the method can cope with dynamic contact problems which is shown by an analysis of a model problem moreover numerical examples demonstrate that it is actually able to approximate the solution of dynamic contact problems which are not fully covered by the theory copyright 2008 john wiley sons ltd

a constrained optimization methodology for the detection phase in contact mechanics simulations the detection phase in computational contact mechanics can be subdivided into a global search and a local detection when potential contact is detected by the former a rigorous local detection determines which surface elements come or may come in contact in the current increment we first introduce a rigorous definition of the closest point for non differentiable lower dimensional manifolds we then simplify the detection by formulating an optimization problem subject to inequality constraints the formulation is then solved using different techniques from the field of mathematical optimization for both linear and quadratic finite element meshes the resulting general and robust detection scheme is tested on a set of problems and compared with other techniques commonly used in computational geometry