Computational Wear and Contact Modeling for Fretting Analysis with Isogeometric Dual Mortar Methods

2016 ◽  
Vol 681 ◽  
pp. 1-18 ◽  
Author(s):  
Philipp Farah ◽  
Markus Gitterle ◽  
Wolfgang A. Wall ◽  
Alexander Popp
Keyword(s):  
Finite Element ◽  
Three Dimensions ◽  
Fretting Wear ◽  
Additional Contribution ◽  
Shape Functions ◽  
Element Discretization ◽  
B Splines ◽  
Mortar Methods ◽  
Deformation Regime ◽  
Incremental Scheme

A finite element framework based on dual mortar methods is presented for simulating fretting wear effects in the finite deformation regime. The mortar finite element discretization is realized with Lagrangean shape functions as well as isogeometric elements based on non-uniform rational B-splines (NURBS) in two and three dimensions. Fretting wear effects are modeled in an incremental scheme with the help of Archard’s law and the worn material is considered as additional contribution to the gap function. Numerical examples demonstrate the robustness and accuracy of the presented algorithm.

Development of a Nonlinear, C1-Beam Finite-Element Code for Actuator Design of Slender Flexible Robots

2017 ◽  
Author(s):  
Hidenori Murakami ◽  
Oscar Rios ◽  
Takeyuki Ono
Keyword(s):  
Finite Element ◽  
Large Deformation ◽  
Cross Section ◽  
Beam Model ◽  
Finite Element Code ◽  
Shape Functions ◽  
Lagrange Equations ◽  
Element Discretization ◽  
Flexible Robots ◽  
Actuator Design

For actuator design and motion simulations of slender flexible robots, planar C1-beam elements are developed for Reissner’s large deformation, shear-deformable, curbed-beam model. Internal actuation is mechanically modeled by a rate-form of beam constitutive relation, where actuation curvature is prescribed at each time. Geometrically, a curbed beam is modeled as a frame bundle, whereby at each point on beam’s curve of centroids a moving orthonormal frame is attached to a cross section. After a finite element discretization, a curve of centroids is modeled as a C1-curve, employing cubic shape functions for both planar coordinates with an arc-parameter. The cubic shape functions have already been utilized in linear Euler-Bernoulli beams for the interpolation of transverse displacement. To define the rotation angle of each cross section or the attitude of the moving frame, quadratic shape functions are used introducing a middle node, resulting in three angular nodal displacements. As a result, each beam element has total eleven nodal coordinates. The implementation of a nonlinear finite element code is facilitated by the principle of virtual work, which yields Reissner’s large deformation curbed beam model as the Euler-Lagrange equations. For time integration, the Newmark method is utilized. Finally, as applications of the code, a few inchworm motions induced by different actuation curvature fields are presented.

Hybrid Isogeometric-Finite Element Discretization Applied to Stress Concentration Problems

2018 ◽  
Vol 10 (08) ◽  
pp. 1850081 ◽  
Author(s):  
Saeed Maleki-Jebeli ◽  
Mahmoud Mosavi-Mashhadi ◽  
Mostafa Baghani
Keyword(s):  
Finite Element ◽  
Stress Concentration ◽  
Finite Elements ◽  
Field Variable ◽  
Higher Order ◽  
Stable Convergence ◽  
Order Continuity ◽  
Element Discretization ◽  
B Splines ◽  
Computer Methods

Isogeometric analysis (IGA) employs non-uniform rational B-splines (NURBS) or other B-spline-based variants to represent both the geometry and the field variable. Exact geometry representation and higher order global continuity (at least [Formula: see text] even on elements’ boundaries) are two favorable properties that would make IGA an appropriate discretization technique in problems with responses associated with the derivatives of the primary field variable. As a category of these problems, in this paper, 2D elastostatic problems involving stress concentration sites are analyzed with a hybrid isogeometric-finite element (IG-FE) discretization. To exploit higher order continuity of NURBS basis functions, IGA discretization is applied selectively at pre-identified locations of high displacement gradients where the stress concentration occurs. In addition, considering computational efficiency, the rest of problem domain is discretized by means of linear Lagrangian finite elements. The connection of NURBS and Lagrangian domain is carried out through employing specially devised elements [Corbett, C. J. and Sauer, R. A. [2014] “NURBS-enriched contact finite elements”, Computer Methods in Applied Mechanics and Engineering 275, 55–75]. The methodology is applied in some 2D elastostatic examples. Increasing the number of DOFs and comparing convergence of the concentrated stress value using different discretizations, it is shown that the hybrid IG-FE discretizations generally have faster and more stable convergence response compared with pure FE discretizations especially at lower DOFs.

Calculating Energy Release Rate as a Function of Crack Length Using a Multiple-Step Crack Closure Technique in Tire Finite Element Models

2018 ◽  
Vol 46 (3) ◽  
pp. 130-152
Author(s):  
Dennis S. Kelliher
Keyword(s):  
Finite Element ◽  
Energy Release Rate ◽  
Crack Length ◽  
Energy Release ◽  
Crack Closure ◽  
Release Rate ◽  
Fatigue Behavior ◽  
Three Dimensions ◽  
Quantitative Prediction ◽  
Closure Technique

ABSTRACT When performing predictive durability analyses on tires using finite element methods, it is generally recognized that energy release rate (ERR) is the best measure by which to characterize the fatigue behavior of rubber. By addressing actual cracks in a simulation geometry, ERR provides a more appropriate durability criterion than the strain energy density (SED) of geometries without cracks. If determined as a function of crack length and loading history, and augmented with material crack growth properties, ERR allows for a quantitative prediction of fatigue life. Complications arise, however, from extra steps required to implement the calculation of ERR within the analysis process. This article presents an overview and some details of a method to perform such analyses. The method involves a preprocessing step that automates the creation of a ribbon crack within an axisymmetric-geometry finite element model at a predetermined location. After inflating and expanding to three dimensions to fully load the tire against a surface, full ribbon sections of the crack are then incrementally closed through multiple solution steps, finally achieving complete closure. A postprocessing step is developed to determine ERR as a function of crack length from this enforced crack closure technique. This includes an innovative approach to calculating ERR as the crack length approaches zero.

scholarly journals Higher Order Multiscale Finite Element Method for Heat Transfer Modeling

Materials ◽  
2021 ◽  
Vol 14 (14) ◽  
pp. 3827
Author(s):  
Marek Klimczak ◽  
Witold Cecot
Keyword(s):  
Heat Transfer ◽  
Finite Element Method ◽  
Finite Element ◽  
Degrees Of Freedom ◽  
Higher Order ◽  
Shape Functions ◽  
Multiscale Finite Element Method ◽  
Multiscale Finite Element ◽  
Steady State Heat ◽  
Steady State Heat Transfer

In this paper, we present a new approach to model the steady-state heat transfer in heterogeneous materials. The multiscale finite element method (MsFEM) is improved and used to solve this problem. MsFEM is a fast and flexible method for upscaling. Its numerical efficiency is based on the natural parallelization of the main computations and their further simplifications due to the numerical nature of the problem. The approach does not require the distinct separation of scales, which makes its applicability to the numerical modeling of the composites very broad. Our novelty relies on modifications to the standard higher-order shape functions, which are then applied to the steady-state heat transfer problem. To the best of our knowledge, MsFEM (based on the special shape function assessment) has not been previously used for an approximation order higher than p = 2, with the hierarchical shape functions applied and non-periodic domains, in this problem. Some numerical results are presented and compared with the standard direct finite-element solutions. The first test shows the performance of higher-order MsFEM for the asphalt concrete sample which is subject to heating. The second test is the challenging problem of metal foam analysis. The thermal conductivity of air and aluminum differ by several orders of magnitude, which is typically very difficult for the upscaling methods. A very good agreement between our upscaled and reference results was observed, together with a significant reduction in the number of degrees of freedom. The error analysis and the p-convergence of the method are also presented. The latter is studied in terms of both the number of degrees of freedom and the computational time.

A Full-System Approach of the Elastohydrodynamic Line/Point Contact Problem

2008 ◽  
Vol 130 (2) ◽  
Author(s):  
W. Habchi ◽  
D. Eyheramendy ◽  
P. Vergne ◽  
G. Morales-Espejel
Keyword(s):  
Finite Element ◽  
Nonlinear System ◽  
System Approach ◽  
Jacobian Matrix ◽  
Point Contact ◽  
System Of Equations ◽  
Nonlinear System Of Equations ◽  
Full System ◽  
Element Discretization ◽  
Fully Coupled

The solution of the elastohydrodynamic lubrication (EHL) problem involves the simultaneous resolution of the hydrodynamic (Reynolds equation) and elastic problems (elastic deformation of the contacting surfaces). Up to now, most of the numerical works dealing with the modeling of the isothermal EHL problem were based on a weak coupling resolution of the Reynolds and elasticity equations (semi-system approach). The latter were solved separately using iterative schemes and a finite difference discretization. Very few authors attempted to solve the problem in a fully coupled way, thus solving both equations simultaneously (full-system approach). These attempts suffered from a major drawback which is the almost full Jacobian matrix of the nonlinear system of equations. This work presents a new approach for solving the fully coupled isothermal elastohydrodynamic problem using a finite element discretization of the corresponding equations. The use of the finite element method allows the use of variable unstructured meshing and different types of elements within the same model which leads to a reduced size of the problem. The nonlinear system of equations is solved using a Newton procedure which provides faster convergence rates. Suitable stabilization techniques are used to extend the solution to the case of highly loaded contacts. The complexity is the same as for classical algorithms, but an improved convergence rate, a reduced size of the problem and a sparse Jacobian matrix are obtained. Thus, the computational effort, time and memory usage are considerably reduced.

Rayleigh-Ritz Analysis of Vibrating Plates Based on a Class of Eigenfunctions

2009 ◽  
Author(s):  
Giuseppe Catania ◽  
Silvio Sorrentino
Keyword(s):  
Finite Element Analysis ◽  
Finite Element ◽  
Boundary Conditions ◽  
Linear Combination ◽  
Equation Of Motion ◽  
Extensive Literature ◽  
Shape Functions ◽  
Element Analysis ◽  
Vibrating Plates ◽  
General Basis

In the Rayleigh-Ritz condensation method the solution of the equation of motion is approximated by a linear combination of shape-functions selected among appropriate sets. Extensive literature dealing with the choice of appropriate basis of shape functions exists, the selection depending on the particular boundary conditions of the structure considered. This paper is aimed at investigating the possibility of adopting a set of eigenfunctions evaluated from a simple stucture as a general basis for the analysis of arbitrary-shaped plates. The results are compared to those available in the literature and using standard finite element analysis.

An Angular Finite Element Method for the Solution of the Even-Parity Equation

2009 ◽  
Author(s):  
R. Becker ◽  
R. Koch ◽  
M. F. Modest ◽  
H.-J. Bauer
Keyword(s):  
Finite Element ◽  
Finite Element Approximation ◽  
Basis Functions ◽  
New Method ◽  
Discrete Ordinates ◽  
Test Case ◽  
Element Approximation ◽  
Promising Alternative ◽  
Element Discretization ◽  
Spatial Coordinates

The present article introduces a new method to solve the radiative transfer equation (RTE). First, a finite element discretization of the solid angle dependence is derived, wherein the coefficients of the finite element approximation are functions of the spatial coordinates. The angular basis functions are defined according to finite element principles on subdivisions of the octahedron. In a second step, these spatially dependent coefficients are discretized by spatial finite elements. This approach is very attractive, since it provides a concise derivation for approximations of the angular dependence with an arbitrary number of angular nodes. In addition, the usage of high-order angular basis functions is straightforward. In the current paper the governing equations are first derived independently of the actual angular approximation. Then, the design principles for the angular mesh are discussed and the parameterization of the piecewise angular basis functions is derived. In the following, the method is applied to two-dimensional test cases which are commonly used for the validation of approximation methods of the RTE. The results reveal that the proposed method is a promising alternative to the well-established practices like the Discrete Ordinates Method (DOM) and provides highly accurate approximations. A test case known to exhibit the ray effect in the DOM verifies the ability of the new method to avoid ray effects.

Sign in / Sign up

Export Citation Format

Share Document