H. T. Xiao¹, Z. Q. Yue²

CMES-Computer Modeling in Engineering & Sciences, Vol.52, No.2, pp. 159-180, 2009, DOI:10.3970/cmes.2009.052.159

Abstract This paper presents a three-dimensional boundary element method for contact problems of an elastic indenter on the surface of functionally graded materials (FGMs). The FGM elastic properties can have any irregular variations with depth. The indenter is subjected to the loading normal to the flat contact surface. The classical Kelvin solution is used for the mathematical formulation of the homogeneous elastic indenter. The generalized Kelvin solution is used for the mathematical formulation of the FGM base. The contact variables are defined with respect to each of the surfaces using local coordinate systems. The corresponding contact More >

Peter J. Ireman, Anders Klarbring¹, Niclas Strömberg

CMES-Computer Modeling in Engineering & Sciences, Vol.52, No.2, pp. 125-158, 2009, DOI:10.3970/cmes.2009.052.125

Abstract In this paper finite element approaches for fretting fatigue are proposed on the basis of a non-local model of continuum damage coupled to friction and wear. The model is formulated in the frame-work of a standard material. In a previous paper this was done in the spirit of Maugin, where an extra entropy flux is introduced in the second law in order to include the gradient of the internal variable in a proper manner. In this paper we follow instead the ideas of Frémond and others, where this extra entropy flux is no longer needed,… More >

J. Ma², H. Lu², B. Wang², R. Hornung³, A. Wissink³, R. Komanduri^2,*

CMES-Computer Modeling in Engineering & Sciences, Vol.14, No.2, pp. 101-118, 2006, DOI:10.3970/cmes.2006.014.101

Abstract A new method for multiscale simulation bridging two scales, namely, the continuum scale using the generalized interpolation material point (GIMP) method and the atomistic scale using the molecular dynamics (MD), is presented and verified in 2D. The atomistic strain from the molecular dynamics simulation is determined through interpolation of the displacement field into an Eulerian background grid using the same generalized interpolation functions as that in the GIMP method. The atomistic strain is consistent with that determined from the virial theorem for interior points but provides more accurate values at the boundary of the MD… More >

S. Hartmann¹, R. Weyler², J. Oliver¹, J.C. Cante², J.A. Hernández¹

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.3, pp. 211-270, 2010, DOI:10.3970/cmes.2010.055.211

Abstract This work describes a three-dimensional contact domain method for large deformation frictionless contact problems. Theoretical basis and numerical aspects of this specific contact method are given in [Oliver, Hartmann, Cante, Weyler and Hernández (2009)] and [Hartmann, Oliver, Weyler, Cante and Hernández (2009)] for two-dimensional, large deformation frictional contact problems. In this method, in contrast to many other contact formulations, the necessary contact constraints are formulated on a so-called contact domain, which can be interpreted as a fictive intermediate region connecting the potential contact surfaces of the deformable bodies. This contact domain has the same dimension More >

A. Colombo¹, U. Galvanetto²

CMES-Computer Modeling in Engineering & Sciences, Vol.53, No.3, pp. 235-254, 2009, DOI:10.3970/cmes.2009.053.235

Abstract The paper addresses the problem of computing the stable manifolds of equilibria and limit cycles of saddle type in piecewise smooth dynamical systems. All singular points that are generically present along one-dimensional or two-dimensional manifolds are classified and such a classification is then used to define a method for the numerical computation of the stable manifolds. Finally the proposed method is applied to the case of a stick-slip oscillator. More >

S. Tangaramvong, F. Tin-Loi

CMES-Computer Modeling in Engineering & Sciences, Vol.49, No.3, pp. 269-296, 2009, DOI:10.3970/cmes.2009.049.269

Abstract We consider the optimal synthesis, namely minimum weight design, of rigid perfectly-plastic structures for which some supports involve unilateral frictional contact. This problem is of interest as it is not only encountered in practice but it also involves, in the general friction case, a nonassociative complementarity condition that makes it theoretically and numerically challenging. For simplicity of exposition, we focus on the class of bar structures for which yielding is governed by either pure bending or by combined axial and flexural forces. In view of possible multiplicity of solutions due to nonassociativity, a direct optimization More >

K.Willner ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.48, No.3, pp. 303-336, 2009, DOI:10.3970/cmes.2009.048.303

Abstract The dynamic behavior of structures with joints is strongly influenced by the constitutive behavior within the contact areas. In this paper the influence of an elaborate constitutive contact model based on a rough surface model is investigated. The contact model is able to describe several effects like pressure dependent contact stiffness in normal and tangential direction as well as microslip effects. The corresponding constitutive contact laws are implemented in a finite element code. Numerical simulations are compared to experimental results of a clamped double-beam experiment. More >

G. Narbona-Reina¹, J.D.D. Zabsonré², E.D. Fernández-Nieto¹, D. Bresch³

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.1, pp. 27-72, 2009, DOI:10.3970/cmes.2009.043.027

Abstract In this work we present a new two-dimensional bilayer Shallow-Water model including viscosity and friction effects on the bottom and interface level. It is obtained following [Gerbeau and Perthame (2001)] from an asymptotic analysis of non-dimensional and incompressible Navier-Stokes equations with hydrostatic approximation. In order to obtain the viscosity effects into the model we must have into account a second order approximation. To evaluate this model we perform two numerical tests consisting of an internal dam-break problem for both, one and two dimensional cases. In the first one we make a comparison between the model More >

V.V. Zozulya¹

CMES-Computer Modeling in Engineering & Sciences, Vol.42, No.3, pp. 187-216, 2009, DOI:10.3970/cmes.2009.042.187

Abstract The mathematical statement for contact problem with unilateral restrictions and friction is done in classical and weak forms. Different variational formulation of unilateral contact problems with friction based on principles of virtual displacements and virtual stresses are considered. Especially boundary variational functionals that are used with boundary integral equations have been established. Nonsmooth optimization algorithms of Udzawa type for solution of unilateral contact problem with friction have been developed. Some theoretical results of existence and uniqueness in elastostatic unilateral contact problem with friction are outlined. More >

A.N. Guz¹, V.V. Zozulya²

CMES-Computer Modeling in Engineering & Sciences, Vol.22, No.2, pp. 119-128, 2007, DOI:10.3970/cmes.2007.022.119

Abstract The frictional contact interaction of the edges of a finite plane crack is studied for the case of normal incidence of a harmonic SH-shear wave which produces antiplane deformation. The forces of contact interaction and displacement discontinuity are analyzed. Influence of the wave frequency on the stress intensity factor for different coefficients of friction is studied here. More >