Microplane Constitutive Model and Computational Framework for Blood Vessel Tissue

2005 ◽  
Vol 128 (3) ◽  
pp. 419-427 ◽  
Author(s):  
Ferhun C. Caner ◽  
Ignacio Carol
Keyword(s):  
Finite Element ◽  
Soft Tissue ◽  
Constitutive Modeling ◽  
Arterial Wall ◽  
Convergence Rates ◽  
Equations Of Motion ◽  
Small Time ◽  
Dynamic Relaxation ◽  
Explicit Integration ◽  
Biological Soft Tissue

This paper presents a nonlinearly elastic anisotropic microplane formulation in 3D for computational constitutive modeling of arterial soft tissue in the passive regime. The constitutive modeling of arterial (and other biological) soft tissue is crucial for accurate finite element calculations, which in turn are essential for design of implants, surgical procedures, bioartificial tissue, as well as determination of effect of progressive diseases on tissues and implants. The model presented is defined at a lower scale (mesoscale) than the conventional macroscale and it incorporates the effect of all the (collagen) fibers which are anisotropic structural components distributed in all directions within the tissue material in addition to that of isotropic bulk tissue. It is shown that the proposed model not only reproduces Holzapfel’s recent model but also improves on it by accounting for the actual three-dimensional distribution of fiber orientation in the arterial wall, which endows the model with advanced capabilities in simulation of remodeling of soft tissue. The formulation is flexible so that its parameters could be adjusted to represent the arterial wall either as a single material or a material composed of several layers in finite element analyses of arteries. Explicit algorithms for both the material subroutine and the explicit integration with dynamic relaxation of equations of motion using finite element method are given. To circumvent the slow convergence of the standard dynamic relaxation and small time steps dictated by the stability of the explicit integrator, an adaptive dynamic relaxation technique that ensures stability and fastest possible convergence rates is developed. Incompressibility is enforced using penalty method with an updated penalty parameter. The model is used to simulate experimental data from the literature demonstrating that the model response is in excellent agreement with the data. An experimental procedure to determine the distribution of fiber directions in 3D for biological soft tissue is suggested in accordance with the microplane concept. It is also argued that this microplane formulation could be modified or extended to model many other phenomena of interest in biomechanics.

scholarly journals Inverse Nonlinear Finite Element Methods for Surgery Simulation and Image Guidance

2008 ◽  
Author(s):  
Philip Pratt
Keyword(s):  
Finite Element ◽  
Finite Element Methods ◽  
Equations Of Motion ◽  
Nonlinear Finite Element ◽  
Weighted Sums ◽  
Global Coordinate System ◽  
Surgery Simulation ◽  
Explicit Integration ◽  
Time Step ◽  
Nonlinear Finite Element Methods

Nonlinear finite element methods are described in which cyclic organ motion is implied from 4D scan data. The equations of motion corresponding to an explicit integration of the total Lagrangian formulation are reversed, such that the sequence of node forces which produces known changes in displacement is recovered. The forces are resolved from the global coordinate system into systems local to each element, and at every simulation time step are expressed as weighted sums of edge vectors. In the presence of large deformations and rotations, this facilitates the combination of external forces, such as tool-tissue interactions, and also positional constraints. Applications in the areas of surgery simulation and minimally invasive robotic interventions are discussed, and the methods are illustrated using CT images of a pneumatically-operated beating heart phantom.

scholarly journals Triply Coupled Bending-Bending-Torsion Vibrational Analysis of Rotating Beams Using Fem and Dynamic Finite Element Formulation

2021 ◽  
Author(s):  
Mohammad Shavezipur
Keyword(s):  
Finite Element ◽  
Natural Frequencies ◽  
Convergence Rates ◽  
Equations Of Motion ◽  
Vibrational Analysis ◽  
Element Formulation ◽  
Shape Functions ◽  
Weighted Residual Method ◽  
Torsion Beam ◽  
Dynamic Finite Element

This research presents the numerical analysis of the triply coupled flap-wise, cord-wise and torsional vibrations of flexible rotating blades. Euler-Bernoulli bending and St. Venant torsion beam theories are considered to derive the governing differential equations of motion. Based on Finite Element Methodology (FEM), the cubic "Hermite" shape functions are implemented where the solution of the equations results in a linear engine problem. Then, the Dynamic (frequency dependent) Trigonometric Shape Functions (DISF's) for beam's uncoupled displacements are derived. The application of the Dynamic Finite Element (DFE) approach to the solution of the governing equations is then presented. The DFE formulation, based on the weighted residual method and the DTSF's results in a nonlinear engine problem representing eigenvalues and engine modes of the system. The applicability of the DFE method is then demonstrated by illustrative examples, where a Wittrick-Williams root counting technique is used to find the system's natural frequencies. The DFE approach, an intermediate method between FEM and "Exact" formulation, is characterized by higher convergence rates, and can be advantageously used when multiple natural frequencies and/or higher modes of beam-like structures are to be evaluated.

scholarly journals Triply Coupled Bending-Bending-Torsion Vibrational Analysis of Rotating Beams Using Fem and Dynamic Finite Element Formulation

2021 ◽  
Author(s):  
Mohammad Shavezipur
Keyword(s):  
Finite Element ◽  
Natural Frequencies ◽  
Convergence Rates ◽  
Equations Of Motion ◽  
Vibrational Analysis ◽  
Element Formulation ◽  
Shape Functions ◽  
Weighted Residual Method ◽  
Torsion Beam ◽  
Dynamic Finite Element

This research presents the numerical analysis of the triply coupled flap-wise, cord-wise and torsional vibrations of flexible rotating blades. Euler-Bernoulli bending and St. Venant torsion beam theories are considered to derive the governing differential equations of motion. Based on Finite Element Methodology (FEM), the cubic "Hermite" shape functions are implemented where the solution of the equations results in a linear engine problem. Then, the Dynamic (frequency dependent) Trigonometric Shape Functions (DISF's) for beam's uncoupled displacements are derived. The application of the Dynamic Finite Element (DFE) approach to the solution of the governing equations is then presented. The DFE formulation, based on the weighted residual method and the DTSF's results in a nonlinear engine problem representing eigenvalues and engine modes of the system. The applicability of the DFE method is then demonstrated by illustrative examples, where a Wittrick-Williams root counting technique is used to find the system's natural frequencies. The DFE approach, an intermediate method between FEM and "Exact" formulation, is characterized by higher convergence rates, and can be advantageously used when multiple natural frequencies and/or higher modes of beam-like structures are to be evaluated.

The Influence of Loading Position in A Priori High Stress Detection using Mode Superposition

2020 ◽  
Vol 1 (1) ◽  
pp. 93-102
Author(s):  
Carsten Strzalka ◽  
◽  
Manfred Zehn ◽  
Keyword(s):  
Finite Element ◽  
Degrees Of Freedom ◽  
Equations Of Motion ◽  
A Priori ◽  
High Stress ◽  
Von Mises Stress ◽  
Detailed Knowledge ◽  
Variable Loading ◽  
Numerical Effort ◽  
Von Mises

For the analysis of structural components, the finite element method (FEM) has become the most widely applied tool for numerical stress- and subsequent durability analyses. In industrial application advanced FE-models result in high numbers of degrees of freedom, making dynamic analyses time-consuming and expensive. As detailed finite element models are necessary for accurate stress results, the resulting data and connected numerical effort from dynamic stress analysis can be high. For the reduction of that effort, sophisticated methods have been developed to limit numerical calculations and processing of data to only small fractions of the global model. Therefore, detailed knowledge of the position of a component’s highly stressed areas is of great advantage for any present or subsequent analysis steps. In this paper an efficient method for the a priori detection of highly stressed areas of force-excited components is presented, based on modal stress superposition. As the component’s dynamic response and corresponding stress is always a function of its excitation, special attention is paid to the influence of the loading position. Based on the frequency domain solution of the modally decoupled equations of motion, a coefficient for a priori weighted superposition of modal von Mises stress fields is developed and validated on a simply supported cantilever beam structure with variable loading positions. The proposed approach is then applied to a simplified industrial model of a twist beam rear axle.

A New Approach to the Rigid Finite Element Method in Modeling Spatial Slender Systems

2018 ◽  
Vol 18 (02) ◽  
pp. 1850017 ◽  
Author(s):  
Iwona Adamiec-Wójcik ◽  
Łukasz Drąg ◽  
Stanisław Wojciech
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Dynamic Analysis ◽  
Equations Of Motion ◽  
Nonlinear Dynamic Analysis ◽  
New Approach ◽  
Static And Dynamic Analysis ◽  
Rigid Finite Element Method ◽  
Longitudinal Flexibility ◽  
Element Method

The static and dynamic analysis of slender systems, which in this paper comprise lines and flexible links of manipulators, requires large deformations to be taken into consideration. This paper presents a modification of the rigid finite element method which enables modeling of such systems to include bending, torsional and longitudinal flexibility. In the formulation used, the elements into which the link is divided have seven DOFs. These describe the position of a chosen point, the extension of the element, and its orientation by means of the Euler angles Z[Formula: see text]Y[Formula: see text]X[Formula: see text]. Elements are connected by means of geometrical constraint equations. A compact algorithm for formulating and integrating the equations of motion is given. Models and programs are verified by comparing the results to those obtained by analytical solution and those from the finite element method. Finally, they are used to solve a benchmark problem encountered in nonlinear dynamic analysis of multibody systems.

Sign in / Sign up

Export Citation Format

Share Document