Application of the u-p Finite Element Method to the Study of Articular Cartilage

1991 ◽  
Vol 113 (4) ◽  
pp. 397-403 ◽  
Author(s):  
Jennifer S. Wayne ◽  
Savio L.-Y. Woo ◽  
Michael K. Kwan
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Articular Cartilage ◽  
Numerical Solutions ◽  
Virtual Work ◽  
Fluid Pressure ◽  
Numerical Procedure ◽  
Solid Matrix ◽  
Confined Compression ◽  
Element Method

The finite element method using the principle of virtual work was applied to the biphasic theory to establish a numerical routine for analyses of articular cartilage behavior. The matrix equations that resulted contained displacements of the solid matrix (u) and true fluid pressure (p) as the unknown variables at the element nodes. Both small and large strain conditions were considered. The algorithms and computer code for the analysis of two-dimensional plane strain, plane stress, and axially symmetric cases were developed. The u-p finite element numerical procedure demonstrated excellent agreement with available closed-form and numerical solutions for the configurations of confined compression and unconfined compression under small strains, and for confined compression under large strains. The model was also used to examine the behavior of a repaired articular surface. The differences in material properties between the repair tissue and normal cartilage resulted in significant deformation gradients across the repair interface as well as increased fluid efflux from the tissue.

scholarly journals A Corotational Finite Element Method Combined with Floating Frame Method for Large Steady-State Deformation and Free Vibration Analysis of a Rotating-Inclined Beam

2011 ◽  
Vol 2011 ◽  
pp. 1-29 ◽  
Author(s):  
Ming Hsu Tsai ◽  
Wen Yi Lin ◽  
Yu Chun Zhou ◽  
Kuo Mo Hsiao
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Steady State ◽  
Virtual Work ◽  
Numerical Procedure ◽  
Free Vibration Analysis ◽  
Euler Beam ◽  
Floating Frame ◽  
Frame Method ◽  
Element Method

A corotational finite element method combined with floating frame method and a numerical procedure is proposed to investigate large steady-state deformation and infinitesimal-free vibrationaround the steady-state deformation of a rotating-inclined Euler beam at constant angular velocity. The element nodal forces are derived using the consistent second-order linearization of the nonlinear beam theory, the d'Alembert principle, and the virtual work principle in a current inertia element coordinates, which is coincident with a rotating element coordinate system constructed at the current configuration of the beam element. The governing equations for linear vibration are obtained by the first-order Taylor series expansion of the equation of motion at the position of steady-state deformation. Numerical examples are studied to demonstrate the accuracy and efficiency of the proposed method and to investigate the steady-state deformation and natural frequency of the rotating beam with different inclined angle, angular velocities, radius of the hub, and slenderness ratios.

scholarly journals Influence of Retardation Time on Viscoelastic Behavior of Plane Beams Modeled by the Nonlinear Positional Formulation of the Finite Element Method

2021 ◽  
Vol 2021 ◽  
pp. 1-22
Author(s):  
Juliano dos Santos Becho ◽  
Marcelo Greco
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Rheological Model ◽  
Iterative Process ◽  
Numerical Procedure ◽  
Viscoelastic Behavior ◽  
Retardation Time ◽  
The Finite Element Method ◽  
Divergence Problem ◽  
Element Method

A numerical procedure is presented to avoid the divergence problem during the iterative process in viscoelastic analyses. This problem is observed when the positional formulation of the finite element method is adopted in association with the finite difference method. To do this, the nonlinear positional formulation is presented considering plane frame elements with Bernoulli–Euler kinematics and viscoelastic behavior. The considered geometrical nonlinearity refers to the structural equilibrium analysis in the deformed position using the Newton–Raphson iterative method. However, the considered physical nonlinearity refers to the description of the viscoelastic behavior through the adoption of the stress-strain relation based on the Kelvin–Voigt rheological model. After the presentation of the formulation, a detailed analysis of the divergence problem in the iterative process is performed. Then, an original numerical procedure is presented to avoid the divergence problem based on the retardation time of the adopted rheological model and the penalization of the nodal position correction vector. Based on the developments and the obtained results, it is possible to conclude that the presented formulation is consistent and that the proposed procedure allows for obtaining the equilibrium positions for any time step value adopted without presenting divergence problems during the iterative process and without changing the analysis of the final results.

Hybrid Nodal Integral / Finite Element Method for Time-Dependent Convection Diffusion Equation

2021 ◽  
Author(s):  
Sundar Namala ◽  
Rizwan Uddin
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Diffusion Equation ◽  
Numerical Solutions ◽  
Second Order ◽  
Time Dependent ◽  
Convection Diffusion ◽  
Convection Diffusion Equation ◽  
Integral Methods ◽  
Element Method

Abstract Nodal integral methods (NIM) are a class of efficient coarse mesh methods that use transverse averaging to reduce the governing partial differential equation(s) (PDE) into a set of ordinary differential equations (ODE). The standard application of NIM is restricted to domains that have boundaries parallel to one of the coordinate axes/palnes (in 2D/3D). The hybrid nodal-integral/finite-element method (NI-FEM) reported here has been developed to extend the application of NIM to arbitrary domains. NI-FEM is based on the idea that the interior region and the regions with boundaries parallel to the coordinate axes (2D) or coordinate planes (3D) can be solved using NIM, and the rest of the domain can be discretized and solved using FEM. The crux of the hybrid NI-FEM is in developing interfacial conditions at the common interfaces between the NIM regions and FEM regions. We here report the development of hybrid NI-FEM for the time-dependent convection-diffusion equation (CDE) in arbitrary domains. Resulting hybrid numerical scheme is implemented in a parallel framework in Fortran and solved using PETSc. The preliminary approach to domain decomposition is also discussed. Numerical solutions are compared with exact solutions, and the scheme is shown to be second order accurate in both space and time. The order of approximations used for the development of the scheme are also shown to be second order. The hybrid method is more efficient compared to standalone conventional numerical schemes like FEM.

Linearization of the finite element method for gradient flows by Newton’s method

2020 ◽  
Author(s):  
Georgios Akrivis ◽  
Buyang Li
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Newton’S Method ◽  
Newton's Method ◽  
Numerical Solutions ◽  
Resolvent Estimate ◽  
Optimal Order ◽  
Gradient Flows ◽  
The Finite Element Method ◽  
Element Method

Abstract The implicit Euler scheme for nonlinear partial differential equations of gradient flows is linearized by Newton’s method, discretized in space by the finite element method. With two Newton iterations at each time level, almost optimal order convergence of the numerical solutions is established in both the $L^q(\varOmega )$ and $W^{1,q}(\varOmega )$ norms. The proof is based on techniques utilizing the resolvent estimate of elliptic operators on $L^q(\varOmega )$ and the maximal $L^p$-regularity of fully discrete finite element solutions on $W^{-1,q}(\varOmega )$.

A Size-Dependent Coupled Symplectic and Finite Element Method for Steady-State Forced Vibration of Built-Up Nanobeam Systems

2019 ◽  
Vol 19 (07) ◽  
pp. 1950081 ◽  
Author(s):  
Zhenhuan Zhou ◽  
Junhai Fan ◽  
C. W. Lim ◽  
Dalun Rong ◽  
Xinsheng Xu
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Steady State ◽  
Forced Vibration ◽  
Numerical Solutions ◽  
New Approach ◽  
Size Dependent ◽  
Numerical Result ◽  
Proper Comparison ◽  
Element Method

A novel size-dependent coupled symplectic and finite element method (FEM) is proposed to study the steady-state forced vibration of built-up nanobeam system resting on elastic foundations. The overall system is modeled as a combination of nonlocal Timoshenko beams. A new analytical subsystem modeling with formulation and another numerical subsystem modeling are developed and discussed. In the analytical subsystem model, the uniform nanobeams are modeled and solved by a new approach based on a series of analytical symplectic eigensolutions. The numerical subsystem model applies a nonlocal FEM to solve nonuniform nanobeams. Analytical and numerical solutions are presented, and a proper comparison between the two approaches is established. Comprehensive and accurate numerical result is subsequently presented to illustrate the accuracy and reliability of the coupled method. The new results established are expected to have reference values for future studies.

scholarly journals NURSE-2 DoF Device for Arm Motion Guidance: Kinematic, Dynamic, and FEM Analysis

2020 ◽  
Vol 10 (6) ◽  
pp. 2139
Author(s):  
Betsy D. M. Chaparro-Rico ◽  
Daniele Cafolla ◽  
Marco Ceccarelli ◽  
Eduardo Castillo-Castaneda
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Upper Limb ◽  
Numerical Procedure ◽  
Fem Analysis ◽  
Arm Movement ◽  
Assistive Device ◽  
Arm Motion ◽  
Element Method

Patients with neurological or orthopedic lesions require assistance during therapies with repetitive movements. NURSE (cassiNo-qUeretaro uppeR-limb aSsistive dEvice) is an arm movement aid device for both right- and left-upper limb. The device has a big workspace to conduct physical therapy or training on individuals including kids and elderly individuals, of any age and size. This paper describes the mechanism design of NURSE and presents a numerical procedure for testing the mechanism feasibility that includes a kinematic, dynamic, and FEM (Finite Element Method) analysis. The kinematic demonstrated that a big workspace is available in the device to reproduce therapeutic movements. The dynamic analysis shows that commercial motors for low power consumption can achieve the needed displacement, acceleration, speed, and torque. Finite Element Method showed that the mechanism can afford the upper limb weight with light-bars for a tiny design. This work has led to the construction of a NURSE prototype with a light structure of 2.6 kg fitting into a box of 35 × 45 × 30 cm. The latter facilitates portability as well as rehabilitation at home with a proper follow-up. The prototype presented a repeatability of ±1.3 cm that has been considered satisfactory for a device having components manufactured with 3D rapid prototyping technology.

scholarly journals A GALERKIN FINITE ELEMENT METHOD TO SOLVE FRACTIONAL DIFFUSION AND FRACTIONAL DIFFUSION-WAVE EQUATIONS

2013 ◽  
Vol 18 (2) ◽  
pp. 260-273 ◽  
Author(s):  
Alaattin Esen ◽  
Yusuf Ucar ◽  
Nuri Yagmurlu ◽  
Orkun Tasbozan
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Wave Equations ◽  
Numerical Solutions ◽  
Fractional Diffusion ◽  
Galerkin Finite Element Method ◽  
Diffusion Wave ◽  
Galerkin Finite Element ◽  
Base Functions ◽  
Element Method

In the present study, numerical solutions of the fractional diffusion and fractional diffusion-wave equations where fractional derivatives are considered in the Caputo sense have been obtained by a Galerkin finite element method using quadratic B-spline base functions. For the fractional diffusion equation, the L1 discretizaton formula is applied, whereas the L2 discretizaton formula is applied for the fractional diffusion-wave equation. The error norms L 2 and L ∞ are computed to test the accuracy of the proposed method. It is shown that the present scheme is unconditionally stable by applying a stability analysis to the approximation obtained by the proposed scheme.

Sign in / Sign up

Export Citation Format

Share Document