Development of a Finite Element Procedure of Contact Analysis for Articular Cartilage with Large Deformation Based on the Biphasic Theory
|
|
- Branden Cross
- 5 years ago
- Views:
Transcription
1 537 Development of a Finite Element Procedure of Contact Analysis for Articular Cartilage with Large Deformation Based on the Biphasic Theory Xian CHEN,YuanCHEN and Toshiaki HISADA Despite the importance of sliding contact in diarthrodial joints, the contact analysis algorithms presented over the past decade have been limited to cases of infinitesimal deformation and thus cannot reflect the real mechanical behavior of articular cartilage in daily life In this study, a new finite element contact analysis approach allowing a large amount of sliding between articular cartilages is presented based on the biphasic theory, which is an effective model for articular cartilage The geometric constraint condition and the continuity condition of the fluid phase on the contact surfaces are introduced by applying Lagrange multipliers The formulation is carried out by transmitting the contact traction of the tissue and the hydrostatic pressure of the fluid phase equivalently between the contact surfaces by means of integrating virtual work due to contact over the contact area The effectiveness of the proposed algorithm is verified by two numerical examples Key Words: Computational Mechanics, Finite Element Method, Nonlinear Problem, Contact Problem, Biphasic Theory, Articular Cartilage 1 Introduction The extracellular matrix of articular cartilage is composed of proteoglycans and collagen fibers, and protected from excessive deformation and stress by the hydrostatic pressure of the interstitial fluid, which supports about 90% of the external load on the cartilage On the other hand, since mechanical stimulation is regarded a factor that influences the metabolism of chondrocytes, the reduction in deformation and stress of the extracellular matrix may contribute to the maintenance of cartilage function A decrease in the protective function of the interstitial fluid, due to aging or excessive loading, causes destruction of the extracellular matrix and chondrocyte death, thereby leading to joint diseases (1) Therefore, understanding the mechanical behavior of articular cartilage is useful not only for the prevention or treatment of joint diseases but also for the development of artificial tissues Taking the characteristic biphasic structure of cartilage into account, biphasic mixture (2) and poroelastic material (3) mod- Received 23rd May, 2005 (No ) Institute of Environmental Studies, Graduate School of Frontier Sciences, The University of Tokyo, Hongo, Bunkyo-ku, Tokyo , Japan xchen@smlku-tokyoacjp els have been proposed over the past two decades Moreover, since the synovial joint functions as a load transfer mechanism by contact between cartilage surfaces, the finite element method has been applied to contact analyses of articular cartilage based on the biphasic or poroelastic theories (4) (6) However, in these studies, the deformation was assumed to be infinitesimal and the pair of contacting nodes was fixed during the analysis Thus, the sliding that occurs between cartilage surfaces in synovial joints cannot be simulated by these approaches On the other hand, the formulation of a finite element contact analysis for large deformation has been carried out by discretizing the contact surfaces into contact elements and limiting the nodes on the surface of one side (slave surface) to penetrate the surface of the other side (master surface) The so-called one-pass node-to-segment approach (7) has been widely applied due to its ability to reflect a large amount of sliding However, it has been pointed out (8) that this approach may not pass the contact patch test and may lead to errors in the contact traction, since the contact virtual work is evaluated by a concentrated contact force If this approach is applied to a biphasic structure such as cartilage, errors will be produced not only in the contact traction of the tissue but also in the flow and hydrostatic pressure of the interstitial fluid Therefore, JSME International Journal Series C, Vol 48, No 4, 2005
2 538 in order to predict the mechanical behavior of articular cartilage correctly, a new contact analysis algorithm needs to be developed In this study based on the biphasic theory, the impenetrability of the soft tissue and the continuity of the interstitial fluid on the contact surfaces are introduced as constraint conditions by applying a Lagrange multiplier method Furthermore, a new contact algorithm applicable to the contact problems of a biphasic structure with a large amount of sliding is developed, by extending an approach proposed by the authors (9) that can pass a contact patch test developed for single phase structures Furthermore, the validity and effectiveness of the proposed algorithm are verified by two numerical examples 2 Basic Equations of the Biphasic Theory In the biphasic theory, articular cartilage is considered to be a mixture of interstitial fluid as the fluid phase and extracellular matrix as the solid phase Denoting the whole mixture volume as V, the volume fraction of each phase is defined as φ i = dv i /dv i= s, f (1) with superscripts s and f denoting the solid and fluid phases, respectively The saturation condition of cartilage holds as φ f +φ s = 1 (2) For each phase, the apparent density ρ α (mass per unit mixture volume) is related to the true density ρ α T (mass per unit volume of each phase) by ρ α = φ α ρ α T (3) By assuming incompressibility of the solid and fluid phases and combining Eqs (2) and (3) with the mass conservation equation ρ α + (ρ α u α ) = 0 (4) t the following relationship can be obtained (φ f u f +φ s u s ) = (u s +w) = 0 (5) where u s and u f are the velocities of the solid and fluid phases, respectively, w = φ f (u f u s ) gives the relative velocity of the fluid phase with respect to the solid phase Moreover, if the inertial force, volume force and viscosity of the fluid phase are omitted (2), (10), the momentum equation of the mixture and Darcy s law, which governs the motion of the fluid phase in the mixture, are given as follows σ = 0 (6) w = κ p (7) where σ denotes the Cauchy stress tensor of the mixture, the second order tensor κ represents the permeability and p is the hydrostatic pressure of the interstitial fluid By Series C, Vol 48, No 4, 2005 omitting the viscosity of the fluid phase, the Cauchy stress of the mixture can be expressed as σ = pi +σ e (8) where σ e denotes the amount of stress due to elastic deformation of the solid phase and can be determined by applying the constitutive relationship of a hyperelastic material for large deformation problems The kinetic boundary conditions are given by σ T n= ˆt on t (9) p = ˆp on p (10) where ˆt and ˆp are the total external traction of the mixture and the hydrostatic pressure of the external fluid on the portions t and p of the surface, respectively Thus, Eqs (5) (10) and the constitutive relationship of the solid phase compose the boundary value problem of the biphasic theory By multiplying the weight functions δu to Eq (6), δw to Eq (7) and δp to Eq (5), and then integrating them over the volume domain of the solid phaseω, the summation of these integrations leads to δu σdv+ δw ( p+κ 1 w)dv ω ω + δp (u s +w)dv = 0 (11) ω Furthermore, with δu as an admissible displacement field of the solid phase and utilizing Eqs (9) and (10), after integrating Eq (11) by parts followed by use of the divergence theorem and transformation of the volume domain ω and surface domains t, p to the corresponding domains Ω, Γ t and Γ p in the reference configuration, the weak form of the governing equation of the biphasic theory is obtained as (11) δe :(S e pjc 1 )dv + δw κ 1 wjdv Ω Ω ( δw)pjdv + δp (u s +w)jdv Ω Ω = δu tds δw pnds (12) Γ t Γ p where E and C are the Green-Lagrange strain and right Cauchy-Green deformation tensor of the solid phase, respectively, S e is the second Piola-Kirchhoff stress due to deformation of the solid phase, J = det F is the determinant of the solid phase deformation gradient tensor, t and p are the nominal mixture traction and hydrostatic pressure, respectively, referring to the reference configuration, and N is the outward normal unit vector of the surface in the reference configuration 3 Finite Element Formulation for Contact Problems of Articular Cartilage In conventional studies (4) (6), since infinitesimal deformations were assumed or pairs of contacting nodes JSME International Journal
3 539 were fixed in finite element analyses, these approaches cannot be applied to the general case of the large amount of sliding that occurs between cartilage surfaces in synovial joints in daily life On the other hand, although the so-called node-to-segment approach (7), in which sliding is dealt with by considering contact between nodes and contact elements on discretized contact surfaces, is widely adopted due to its simplicity, it has been pointed out that this approach cannot pass the contact patch test and may produce errors of traction in the contact surfaces If this approach is applied to contact analysis of a biphasic structure such as cartilage, it is inevitable that not only errors in the contact traction of the tissue but also errors in the flow and hydrostatic pressure of the interstitial fluid in the contact surfaces will be induced In the present study, in order to overcome this obstacle, a contact analysis algorithm was formulated for articular cartilage with large deformation based on the biphasic theory, by extending an approach proposed by the authors (9) that can pass the contact patch test for a single phase material 3 1 Basic equations for contact analysis Figure 1 shows a model of the contacting surfaces of articular cartilage Considering that the point of cartilage 2 comes into contact with the point of cartilage 1, which is defined as a projection of the point in cartilage 2 to the surface of cartilage 1, the penetration between contact surfaces is defined as follows g = (x 2 x 1 ) n c1 (13) Since the friction coefficient of healthy cartilage surfaces is low, it is assumed that the friction force can be ignored in the formulation By means of the relationship between the outward unit normal vectors in both the contact surfaces n c1 = n c2 (14) the relationship between the contact tractions, t c1 and t c2, of the cartilage surfaces satisfies t c1 = t cn1 n c1, t c2 = t cn2 n c2, t cn t cn1 = t cn2 (15) Therefore, the geometric and kinetic constraint conditions between contacting cartilage surfaces are given as g 0, t cn 0, g t cn = 0 (16) Furthermore, the continuity requirement of an incompressible fluid and the balance condition of hydrostatic pressure between the two sides need (12) Fig 1 Contact surfaces (w c1 w c2 ) n c1 = 0 (17) p c1 = p c2 p c (18) By comparing Eq (17) with Eq (13), it is noted that the penetration defined by Eq (13) specifies the relative location of the contacting cartilage surfaces, while Eq (17) ensures mass balance between the contacting surfaces, which is a characteristic property of a biphasic structure This difference will cause different treatments for these two constraint conditions in the derivation of contact tangent stiffness, as described later By denoting the volume domains of the contacting cartilages in the reference configuration as Ω 1, Ω 2,the surface domains in the reference configuration involving nominal contact tractions t 1, t 2 and nominal hydrostatic pressures p 1, p 2 as Γ t1, Γ t2 and Γ p1, Γ p2, respectively, and the contact boundaries in the current configuration as c1, c2, Eq (12) can be rewritten as δe :(S e pjc 1 )dv + δw κ 1 wjdv Ω 1 Ω 2 Ω 1 Ω 2 ( δw)pjdv + δp (u s +w)jdv Ω 1 Ω 2 Ω 1 Ω 2 = δu tds δw pnds +δw c Γ t1 Γ t2 Γ p1 Γ p2 (19) where δw c = δu t cn1 n c1 ds δw p c1 n c1 ds c1 c1 + δu t cn2 n c2 ds δw p c2 n c2 ds (20) c2 c2 corresponds to the virtual work performed by the contact traction of the mixture and the hydrostatic pressure, and is evaluated in the current configuration for the convenience of introducing continuity condition of interstitial fluid given by Eq (17) 3 2 Discretization of contact surfaces In the context of finite element analysis, the candidate contact articular surfaces are divided into contact elements as shown in Fig 2 In the one-pass approach, the surface of one side (slave surface) is restrained such that it is unable to penetrate the surface of the other side (master surface) Instead of contact between the cartilage surfaces, contact between node x s (slave node) in the slave surface and the element (master element) in the master surface is treated by means of determining the counterpart point x m (master point) of the slave node in the master surface by drawing a line perpendicular to the master element from the slave node Therefore, it becomes possible to treat contact problems with a large amount of sliding between cartilage surfaces by defining the geometric relationship between the node and the element, instead of that between nodes as in the node-to-node approach In Fig 2, ξ 1 and ξ 2 indicate the convected coordinates at the master point x m, while JSME International Journal Series C, Vol 48, No 4, 2005
4 540 Fig 2 Discretization of contact surfaces t m 1, tm 2 and nm are the covariant base vectors and outward unit normal vector, respectively, ie t m 1 = xm ξ 1, tm 2 = xm ξ 2 (21) r m t m 1 tm 2, r m (22) Moreover, the differential area element in the convected coordinate system can be expressed as ds= J m dξ 1 dξ 2 = r m dξ 1 dξ 2 (23) By defining the shape function matrix with shape functions Ni m (i = 1, n) of the contact element as N1 m 0 0 N m n 0 0 [M] = 0 N1 m 0 0 Nn m 0 (24) 0 0 N1 m 0 0 Nn m the weight functions δu, δw and location vector x m at the master point are interpolated as δu m = m] {δu m } (25) δw m = m] {δw m } (26) x m = m] {x m } (27) where {δu m } = { } δu m T 1,1,δum 1,2,δum 1,3,,δum n,1,δum n,2,δum n,3 (28) {δw m } = { } δw m T 1,1,δwm 1,2,δwm 1,3,,δwm n,1,δwm n,2,δwm n,3 (29) {x m } = { x m 1,1, xm 1,2, xm 1,3,, xm n,1, xm n,2, T n,3} xm (30) are vectors consisting of the respective nodal values in the contact element 3 3 Finite element formulation In contrast to evaluating contact virtual work by the concentrated contact force in the conventional one-pass formulation of the finite element contact analysis, the authors have proposed an algorithm (9) that can pass the patch test by performing an integration of the contact virtual work over all the contacting elements The nodes of the master surface are projected onto the slave elements and the contact nodal tractions of the master nodes are transferred from the contact tractions of the projected point, which are obtained by interpolation of the nodal values in the slave elements If this approach is applied to both the contact traction of the tissue and the hydrostatic pressure, nm = rm Series C, Vol 48, No 4, 2005 and Eqs (23) (30) are used, the contact virtual work of Eq (20) in the whole discretized contact surfaces is given by δw c = {δu m } T m m ]T {n m }t m J m dξ 1 dξ 2 +{δu s } T {δw m } T s s ]T {n s }t s J s dξ 1 dξ 2 m m ]T {n m } p m J m dξ 1 dξ 2 {δw s } T s s ]T {n s } p s J s dξ 1 dξ 2 (31) where {δu α }, {δw α } (α = m, s) are defined as the nodal values of the contact element, α similar to Eqs (28) and (29) To perform numerical integration of Eq (31), the contact tractions of the tissue t m, t s and hydrostatic pressures p m, p s at the quadrature point are interpolated from their nodal values and the shape functions of the corresponding contact elements, ie t m = {N m } T {t m }, t s = {N s } T {t s } (32) p m = {N m } T {p m }, p s = {N s } T {p s } (33) {N m } T = {N1 m,,nm n }, {N s } T = {N s 1,,N n} s (34) {t m } = {t1 m,,tm n } T, {t s } = {t s 1,,ts n} T (35) {p m } = {p m 1,, pm n } T, {p s } = {p s 1,, ps n} T (36) On the other hand, by projecting the node of the master element onto the slave element, the nodal contact traction of the tissue ti m and the hydrostatic pressure p m i at node i of the master element are considered to be equal to those at the projected point in the slave element, which can be obtained from interpolation of the nodal values of the slave element, ie ti m = { N s i,1,,n } i,n s p m i = { N s i,1,,n } i,n s t s i,1 t s i,n p s i,1 p s i,n (37) (38) where N s i,l (l = 1, n) is a shape function, and t s i,l and p s i,l are the nodal contact traction of the tissue and the hydrostatic pressure of the projected slave element, respectively With substitution of Eqs (37) and (38) into Eq (31), the contact virtual work for the entire contact surfaces can be expressed as δw c = {δu m } T m m ]T {n m }{N m } T [L]{T s } J m dξ 1 dξ 2 +{δu s } T {δw m } T s s ]T {n s }{N s } T {t s } J s dξ 1 dξ 2 m m ]T {n m }{N m } T [L]{P s } J m dξ 1 dξ 2 JSME International Journal
5 541 {δw s } T s s ]T {n s }{N s } T {p s } J s dξ 1 dξ 2 (39) where N s 1,1,,Ns 1,n [L] = N s n,1,,n n,n s (40) {T s } = {t s 1,1,,ts 1,n ts n,1,,ts n,n} T (41) {P s } = {p s 1,1,, ps 1,n ps n,1,, ps n,n} T (42) denote the shape function matrix, vectors of the nodal contact traction of the tissue and the hydrostatic pressure in all the slave elements related to the nodes in the master surface Furthermore, in this way, only the nodal contact tractions of the tissue and the hydrostatic pressures on the slave surface become the contact nodal variables as Lagrange multipliers In the above formulations, since this research is focused on the patch test problem, the interpolation functions for the contact traction of the tissue and the hydrostatic pressure are chosen to be the same as those of the displacement of the tissue and the relative velocity of fluid phase But it has to be noted that selecting proper interpolations for the contact traction of the tissue and the hydrostatic pressure to satisfy the inf-sup condition is another crucial issue for insuring the stability and optimal convergence of contact analysis of hydrated tissues Since the geometric constraint condition and continuity condition of the fluid phase for the contacting cartilage surfaces need to be satisfied at the contact nodes after discretization, these conditions are discretized by using Eqs (13), (17) and (24) (30) as follows {n m } T ( m ] {x m } [I] { x s }) = 0 (43) {n m } T ( m ] {w m } [I] { w s }) = 0 (44) where { x s } is the location vector, { w s} is the relative velocity vector of the fluid phase of the contact node x s,and[i] is a 3 3 unit matrix To derive the contact tangent stiffness, variations of the unit outward normal vector, the Jacobian transforming the differential area to the convected coordinate system and the covariant base vectors due to deformation are first obtained from Eqs (21) (23) and expressed in matrix form as follows {n α } = 1 J α ( [I] {n α }{n α } T ) {r α } (45) J α = r α = {n α } T {r α } (46) { } t α [ 1 = M α ] {u α } ξ 1 (47) { } t α [ 2 = M α ] {u α ξ 2 } (48) ( [W ] {r α α [ } = 2 M α ] + [ ] W α [ ξ 1 1 M α ] ) {u α } ξ 2 (49) where [M α ]and{u α } are the shape function matrix and the nodal displacement vector of the contact element, respectively, defined similarly to Eqs (24) and (28), and [W1 α] and [W α 2 ] are the matrix expressions of skew-symmetric tensors with the covariant base vectors t α 2 and tα 1 as their axial vectors, ie [ ] 0 t1,3 W α 1 = α t1,2 α t α 1,3 0 t α 1,1 (50) t1,2 α t1,1 α 0 [ ] 0 t2,3 W α 2 = α t2,2 α t2,3 α 0 t2,1 α (51) t α 2,2 t α 2,1 0 In addition, since the nodal contact traction of the tissue and the hydrostatic pressure of the master node are determined from those at the projected point of the slave element, which depends on the relative movement of the contacting surfaces, as long as the variation in the local coordinates of the projected point is computed (13), the contact tangent stiffness can be obtained by differentiating Eq (39) The contributions of the geometric constraint condition and the continuity condition of the fluid phase to the contact tangent stiffness are determined by differentiating Eqs (43) and (44), respectively Since the geometric constraint reflects the relative location of the contacting surfaces, the relationships n m (x m x s ) = 0andn m t m i = 0, (i = 1, 2) considerably simplify the calculation, and differentiation of Eq (43) leads to {n m } T ( m] {u m } [I] { u s } ) = 0 (52) However, for the continuity condition of the fluid phase, similar relationships do not exist, and therefore, considering the variations in the local coordinates and outward normal vector at the master point due to deformation, the differentiation of Eq (44) is derived as {n m } T ( m] {w m } [I] { w s } ) ( +{n m } T m ] ξ 1 + m ] ) ξ 2 {w m ξ 1 ξ 2 } + {n m } T ( m] {w m } [I] { w s } ) = 0 (53) In the above equation, the variations of local coordinates, ξ 1 and ξ 2, have been taken into account, since the master point is defined as the projection of the slave node onto the master element, and is therefore dependent on the relative movement of the contact surfaces 4 Numerical Examples The formulation performed in section 3 is implemented by building the nodal contact traction of the tissue, JSME International Journal Series C, Vol 48, No 4, 2005
6 542 nodal hydrostatic pressure and contact tangent stiffness associated with contacting articular cartilage surfaces into a self-coded biphasic finite element analysis program Next, numerical examples are carried out to assess the validity and effectiveness of the proposed approach using a hexahedral type of element with 8 nodes for displacement and relative fluid velocity and 1 node for hydrostatic pressure In addition, the conventional one-pass approach evaluating the discretized contact virtual work by a concentrated contact force is also performed for comparison 4 1 Contact problem between a cartilage layer and a rigid cylinder To assess the validity of the proposed approach, the contact problem between a cartilage layer and an impermeable rigid cylinder (as shown in Fig 3) is analyzed for comparison with the results (5) of finite element contact analysis of articular cartilage under infinitesimal deformation Since the impermeable rigid cylinder plays the role of imposing a shape on the contact surface, only the part near the surface of the cylinder is modeled and, furthermore, the symmetry is utilized in the modeling All the nodal degrees of freedom of displacement and relative fluid velocity at the bottom of the cartilage layer are constrained by assuming that the cartilage layer is connected to an impermeable rigid bone Furthermore, the nodal degrees of freedom of displacement and relative fluid velocity in the outward direction of the plane are constrained to simulate a 2-D problem In the contact analysis, the impermeable rigid cylinder is set as the master side For comparison with the conventional studies, the geometric nonlinearity is ignored and the solid phase of the cartilage is assumed to be linear elastic with the material constants (5) shown in Fig 3 The load is first applied to the cylinder as a ramp function of time with the maximum value of 1 N per unit length of the cylinder reached within 10e 3 s and then held to induce creep deformation of the cartilage layer The distributions of the contact traction of the tissue and the hydrostatic pressure within the contact area at 100 s after loading obtained by the proposed approach are shown in Fig 4, together with those of the conventional research (5) These results show good agreement, and the small discrepancy near the border of the contact area is considered to be induced by the different finite element meshes of the contact surfaces used in each computation In addition, it is noted that the hydrostatic pressure accounts for about 86% of the contact traction of the tissue, indicating that the hydrostatic pressure supports the majority of the contact load Figure 5 shows the deformation of the cartilage and distribution of the relative velocity of the fluid phase immediately after the load reaches its maximum Since the fluid is prohibited from flowing in and out of the surface of the impermeable rigid cylinder, the continuity condition of the fluid phase on the contact surface results in fluid flow in the tangential direction of the contact surface within the contact area The fluid flow out of the cartilage layer occurs almost at the border of the contact area with a maximum value of 48 µm/s In addition, the contact analysis was also carried out by applying the conventional one-pass approach and, leaving out the details, the results were mostly in agreement with the results of the proposed approach, when the impermeable rigid cylinder was set as the master body in the conventional approach This fact can be explained as follows In the one-pass approach, the nodes in the slave surface, ie the contact surface of the cartilage layer, are constrained such that they are unable to penetrate the master surface Hence, if the master surface is rigid, the nodal displacement in the slave surface is prescribed by the shape Fig 3 Contact model of an impermeable rigid cylinder and a cartilage layer Fig 4 Distributions of the contact and interstitial fluid pressures at 10 s after loading Fig 5 Deformation of cartilage (amplified by 5) and distribution of the relative fluid velocity (arrows) at maximum loading Series C, Vol 48, No 4, 2005 JSME International Journal
7 543 of the master surface In a similar way, the nodal relative fluid velocity in the slave surface is also imposed by means of the continuity condition if the master surface is impermeable Accordingly, the displacement and relative fluid velocity fields, and subsequently the nodal forces in the slave surface, can be evaluated appropriately Consequently, since the conventional one-pass approach involves discretization of the contact virtual work in the slave surface with the concentrated forces at the slave nodes, these appropriately evaluated nodal forces become equivalent to the contact traction of the tissue and hydrostatic pressure in the slave surface in the sense of virtual work However, in real synovial joints, the contact occurs between deformable articular cartilage layers In addition, since it has been pointed out that the conventional one-pass approach for contact problems of a single phase elastic solid cannot pass the patch test due to the nonconforming meshes of the contact surfaces (8), it can be predicted that considerable errors in the contact traction of the tissue and the hydrostatic pressure will be induced by the non-conforming meshes of the contact surfaces From a practical standpoint, the meshes cannot always remain conformable under a large amount of sliding even when they are generated conformably in the initial state Therefore, the effectiveness of the proposed approach for contact problems between deformable bodies with a large amount of sliding will be discussed in the next section 4 2 Contact problem between articular cartilage layers The finite element mesh generated without special attention to the conformability is shown in Fig 6 All the nodal degrees of freedom of displacement and relative fluid velocity at the bottom of the lower cartilage layer are constrained by assuming that the cartilage layer is connected to an impermeable rigid bone The shaded portion of the upper body is also assigned as rigid impermeable bone by specifying a very high stiffness and fixing the degree of freedom of the relative fluid velocity, while the remainder represents another cartilage layer The nodal degrees of freedom of displacement and relative fluid velocity in the outward direction of the plane are constrained to simulate a 2-D problem The solid phase of the two cartilage layers is assumed to be a hyperelastic material of the Saint Venant-Kirchhoff type with the same material constants specified in Fig 6 The analyses are carried out using both the conventional and proposed approaches for comparison by applying a load from 0 1 N within 001 s In the analysis using the conventional approach, due to unevenness of the deformed contact surfaces caused by factors discussed below, certain contact nodes near the edges of the contact elements could not be matched to an unique contact element, ie such contact nodes chattered (14) during the iteration and caused the computation to diverge after the load reached 019 N Therefore, only the results for the two approaches at a load level of 019 N are compared below As shown in Figs 7 and 8, considerable disorder appears in the distributions of the contact traction of the tissue and the hydrostatic pressure using the conventional approach, while smooth distributions of these results are obtained using the proposed approach Consequently, as shown in Fig 9, the conventional approach leads to an unevenly deformed contact surface and a disturbed distribution of the relative fluid velocity, while the proposed approach gives a smooth deformation and regular distribution of the relative fluid velocity From these results, it can be concluded that the proposed approach is applicable to contact problems between deformable cartilages, whereas the conventional approach may cause considerable errors in this case The errors of the conventional approach are consid- Fig 7 Comparison of the total contact pressures obtained by the two approaches at a load level of 019 N Fig 6 FEM model of contacting cartilage layers Fig 8 Comparison of the hydrostatic pressures obtained by the two approaches at a load level of 019 N JSME International Journal Series C, Vol 48, No 4, 2005
8 544 (a) Conventional approach Fig 9 (b) Proposed approach Comparison of the deformations (amplified by 15) and relative fluid velocities (arrows) obtained by the two approaches ered to be contributed by the discretization of the contact virtual work with the concentrated force related to the contact traction of the tissue and the hydrostatic pressure In general, this discretization is not guaranteed to be equivalent to the original virtual work and therefore transfers of the contact traction of the tissue and the hydrostatic pressure between contact surfaces are unable to be carried out appropriately The errors not only influence the deformation of the cartilage but also the elastic stress of the solid phase and the flow of the interstitial fluid, and furthermore, they cannot be reduced by mesh refinement On the other hand, in the proposed approach, the contact virtual work is numerically integrated in both the master and slave elements, and instead of a concentrated force, the nodal contact traction of the tissue and the hydrostatic pressure are transferred between the contact surfaces In this way, the nodal forces with respect to contact can be evaluated equivalently in the sense of virtual work and thus reasonable results are obtained Next, contact analysis of articular cartilage layers with a large amount of sliding is carried out using the proposed approach After the load reaches 1 N at 001 s, a rightward horizontal velocity of 10 mm/s is applied to the upper cartilage to produce relative sliding over the lower cartilage within the zone shown in Fig 6, while the load remains unchanged The distributions of the contact traction of the tissue, hydrostatic pressure and normal stress of the solid phase due to deformation when the upper cartilage moves near the center of the lower cartilage are shown in Fig 10 All the results are obtained smoothly, with the maximum values appearing in the center of the contact area In addition, the compression of the cartilage Series C, Vol 48, No 4, 2005 Fig 10 Distributions of the total contact pressures, hydrostatic pressures and normal elastic stresses on the contact surfaces when the upper cartilage layer slides near the center of the lower cartilage layer due to the contact induces the interstitial fluid to flow toward the outside of the contact area, thereby producing the tensile stress state of the solid phase shown in Fig 10 and the bulge deformation of the tissue shown in Fig 11 near the border of the contact area Moreover, Fig 11 shows the distributions of the vertical displacement of the surface of the lower cartilage when the upper cartilage moves at the left side, the center and the right side of the lower cartilage, respectively Without showing the details of the other locations, smooth deformations are obtained during the whole sliding process In the early stage of sliding (Fig 11 (a)), the leading edge of the contact area is pushed to increase the bulge deformation Since the interstitial fluid flows more easily backward into the area that the upper cartilage has passed, the bulge deformation at the trailing edge of the contact area increases gradually along with the sliding due to the increased amount of fluid flow- JSME International Journal
9 545 Fig 12 Distribution of the relative fluid velocity on the surface of the lower cartilage layer (deformation amplified by 7) (a) Sliding near the left side (time = 01 s) For the prevention and medical treatment of joint diseases, as well as the development of artificial tissues that are close to the living tissues, it is necessary to clarify the mechanical behaviors of articular cartilage under physiological conditions In this study, a finite element contact analysis approach for articular cartilage based on the biphasic theory is proposed by taking into account the large amount of sliding that occurs between cartilage surfaces of synovial joints in daily life The impenetrability of the tissue and the continuity condition of the fluid phase on the contact surfaces are enforced by applying Lagrange multipliers The contact virtual work is integrated in the contact surfaces and the contact forces are transferred between the contact surfaces equivalently in the sense of virtual work The proposed approach is implemented and two numerical examples have been carried out to confirm the validity and effectiveness of the proposed approach References (b) Sliding near the center (time = 05 s) (c) Sliding near the right side (time = 09 s) Fig 11 Vertical nodal displacements on the upper surface of the lower cartilage layer ing backward Eventually, the bulge deformation at the trailing edge of the contact area becomes larger than that at the leading edge (Fig 11 (c)) The distribution of the relative velocity of the fluid phase shown in Fig 12 also demonstrates that a larger relative velocity appears at the trailing edge than at the leading edge Based on the results of the contact analysis for articular cartilage with a large amount of sliding discussed above, the validity and effectiveness of the proposed approach for contact problems of articular cartilage with a large amount of sliding are confirmed, and it is expected that the proposed approach will be a useful tool for actual analyses of joint mechanics 5 Conclusions ( 1 ) Milentijevic, D, Helfet, DL and Torzilli, PA, Influence of Stress Magnitude on Water Loss and Chondrocyte Viability in Impacted Articular Cartilage, Journal of Biomechanical Engineering, Vol125 (2003), pp ( 2 ) Mow, VC, Kuei, SC, Lai, WM and Armstrong, CG, Biphasic Creep and Stress Relaxation of Articular Cartilage in Compression: Theory and Experiments, Journal of Biomechanical Engineering, Vol102 (1980), pp73 84 ( 3 ) Simon, BR, Wu, JSS, Carlton, MW, Evans, JH and Kazarian, LE, Structural Models for Human Spinal Motion Segments Based on a Poroelastic View of the Intervertebral Disk, Journal of Biomechanical Engineering, Vol107 (1985), pp ( 4 ) Macirowski, T, Tepic, S and Mann, RW, Cartilage Stresses in the Human Hip Joint, Journal of Biomechanical Engineering, Vol116 (1994), pp10 18 ( 5 ) Donzelli, PS and Spilker, RL, A Contact Finite Element Formulation for Biological Soft Hydrated Tissues, Computer Methods in Applied Mechanics and Engineering, Vol153 (1998), pp63 79 ( 6 ) Federico, S, Rosa, GL, Herzog, W and Wu, JZ, Effect of Fluid Boundary Conditions on Joint Contact Mechanics and Applications to the Modeling of Osteoarthritic Joints, Journal of Biomechanical Engineering, Vol126 (2004), pp ( 7 ) Wriggers, P, Computational Contact Mechanics, (2002), John Wiley & Sons Ltd ( 8 ) El-Abbasi, N and Bathe, KJ, Stability and Patch Test Performance of Contact Discretizations and a New Solution Algorithm, Computers and Structures, Vol79 (2001), pp ( 9 ) Chen, X and Hisada, T, Development of Finite Element Contact Analysis Algorithm Passing Patch Test, JSME International Journal Series C, Vol 48, No 4, 2005
10 546 Trans Jpn Soc Mech Eng, (in Japanese), (in Press) (10) Ateshian, GA and Wang, H, A Theoretical Solution for the Frictionless Rolling Contact of Cylindrical Biphasic Articular Cartilage Layers, Journal of Biomechanics, Vol28 (1995), pp (11) Chen, X, Chen, Y and Hisada, T, A Study on Mechanical Model of Soft Tissues by Nonlinear Finite Element Analysis Based on Biphasic Theory, Trans Jpn Soc Mech Eng, (in Japanese), Vol70, No697, A (2004), pp (12) Hou, JS, Holmes, MH, Lai, WM and Mow, VC, Boundary Conditions at the Cartilage-Synovial Fluid Interface for Joint Lubrication and Theoretical Verifications, Journal of Biomechanical Engineering, Vol111 (1989), pp78 87 (13) Laursen, TA and Simo, JC, A Continuum-Based Finite Element Formulation for the Implicit Solution of Multibody, Large Deformation Frictional Contact Problems, Int J Numer Methods Eng, Vol36 (1993), pp (14) Puso, MA and Laursen, TA, A 3D Contact Smoothing Method Using Gregory Patches, Int J Numer Eng, Vol54 (2002), pp Series C, Vol 48, No 4, 2005 JSME International Journal
Modelling Anisotropic, Hyperelastic Materials in ABAQUS
Modelling Anisotropic, Hyperelastic Materials in ABAQUS Salvatore Federico and Walter Herzog Human Performance Laboratory, Faculty of Kinesiology, The University of Calgary 2500 University Drive NW, Calgary,
More informationFinite Element Analysis of Permeation Tests on Articular Cartilage under Different Testing Conditions Using COMSOL Multiphysics.
Excerpt from the Proceedings of the COMSOL Conference 2010 Paris Finite Element Analysis of Permeation Tests on Articular Cartilage under Different Testing Conditions Using COMSOL Multiphysics. Grazia
More informationA Comparison Between Mechano-Electrochemical and Biphasic Swelling Theories for Soft Hydrated Tissues
A Comparison Between Mechano-Electrochemical and Biphasic Swelling Theories for Soft Hydrated Tissues W. Wilson C. C. van Donkelaar 1 J. M. Huyghe Department of Biomedical Engineering, Eindhoven University
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationBiomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More informationAn Evaluation of Three Dimensional Diarthrodial Joint Contact Using Penetration Data and the Finite Element Method
An Evaluation of Three Dimensional Diarthrodial Joint Contact Using Penetration Data and the Finite Element Method W. L. Dunbar, Jr.*, K. Ün, P. S. Donzelli and R. L. Spilker Department of Biomedical Engineering
More informationFOR PROOFREADING ONLY
Annals of Biomedical Engineering, Vol. 33, No. 4, May 2005 ( 2005) pp. 492 499 DOI: 10.1007/s10439-005-2506-3 Large Deformation Finite Element Analysis of Micropipette Aspiration to Determine the Mechanical
More informationCOMSOL Used for Simulating Biological Remodelling
COMSOL Used for Simulating Biological Remodelling S. Di Stefano 1*, M. M. Knodel 1, K. Hashlamoun 2, S. Federico, A. Grillo 1 1. Department of Mathematical Sciences G. L. Lagrange, Politecnico di Torino,
More informationTraction on the Retina Induced by Saccadic Eye Movements in the Presence of Posterior Vitreous Detachment
Traction on the Retina Induced by Saccadic Eye Movements in the Presence of Posterior Vitreous Detachment Colangeli E., Repetto R., Tatone A. and Testa A. Grenoble, 24 th October 2007 Table of contents
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationEngineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.
Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More informationMATHEMATICAL MODELLING OF HUMAN LIMB IN MATLAB
MATHEMATICAL MODELLING OF HUMAN LIMB IN MATLAB J. Daněk Department of Mathematics, University of West Bohemia Abstract The contribution deals with mathematical modelling of human limb and numerical simulation
More informationin this web service Cambridge University Press
CONTINUUM MECHANICS This is a modern textbook for courses in continuum mechanics. It provides both the theoretical framework and the numerical methods required to model the behavior of continuous materials.
More information(MPa) compute (a) The traction vector acting on an internal material plane with normal n ( e1 e
EN10: Continuum Mechanics Homework : Kinetics Due 1:00 noon Friday February 4th School of Engineering Brown University 1. For the Cauchy stress tensor with components 100 5 50 0 00 (MPa) compute (a) The
More informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More informationCVEN 7511 Computational Mechanics of Solids and Structures
CVEN 7511 Computational Mechanics of Solids and Structures Instructor: Kaspar J. Willam Original Version of Class Notes Chishen T. Lin Fall 1990 Chapter 1 Fundamentals of Continuum Mechanics Abstract In
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationLectures on. Constitutive Modelling of Arteries. Ray Ogden
Lectures on Constitutive Modelling of Arteries Ray Ogden University of Aberdeen Xi an Jiaotong University April 2011 Overview of the Ingredients of Continuum Mechanics needed in Soft Tissue Biomechanics
More informationEQUIVALENCE BETWEEN INSTANTANEOUS BIPHASIC AND INCOMPRESSIBLE ELASTIC MATERIAL RESPONSE. and Biomedical Engineering. Columbia University.
1 1 1 1 1 1 1 1 0 1 EQUIVALENCE BETWEEN INSTANTANEOUS BIPHASIC AND INCOMPRESSIBLE ELASTIC MATERIAL RESPONSE Gerard A. Ateshian a, Benjamin J. Ellis b, Jeffrey A. Weiss b a Departments of Mechanical Engineering
More informationTransactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN
On the Computation of Elastic Elastic Rolling Contact using Adaptive Finite Element Techniques B. Zastrau^, U. Nackenhorst*,J. Jarewski^ ^Institute of Mechanics and Informatics, Technical University Dresden,
More informationUniversity of Groningen
University of Groningen Nature-inspired microfluidic propulsion using magnetic actuation Khaderi, S. N.; Baltussen, M. G. H. M.; Anderson, P. D.; Ioan, D.; den Toonder, J.M.J.; Onck, Patrick Published
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationPredicting Articular Cartilage Behavior with a Non-Linear Microstructural Model
The Open Mechanics Journal, 27, 1, 11-19 11 Predicting Articular Cartilage Behavior with a Non-Linear Microstructural Model Fulin Lei and Andras Z. Szeri * Center for Biomedical Engineering Research, Department
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More information06 - concept of stress concept of stress concept of stress concept of stress. me338 - syllabus. definition of stress
holzapfel nonlinear solid mechanics [2000], chapter 3, pages 109-129 holzapfel nonlinear solid mechanics [2000], chapter 3, pages 109-129 1 2 me338 - syllabus definition of stress stress [ stres] is a
More informationAnalytical formulation of Modified Upper Bound theorem
CHAPTER 3 Analytical formulation of Modified Upper Bound theorem 3.1 Introduction In the mathematical theory of elasticity, the principles of minimum potential energy and minimum complimentary energy are
More informationNIH Public Access Author Manuscript J Biomech. Author manuscript; available in PMC 2008 January 1.
NIH Public Access Author Manuscript Published in final edited form as: J Biomech. 2007 ; 40(9): 2071 2077. Three-dimensional Inhomogeneous Triphasic Finite Element Analysis of Physical Signals and Solute
More informationARTICLE IN PRESS. Journal of Biomechanics
Journal of Biomechanics 43 (2010) 673 679 Contents lists available at ScienceDirect Journal of Biomechanics journal homepage: www.elsevier.com/locate/jbiomech www.jbiomech.com A linearized formulation
More informationarxiv: v1 [physics.comp-ph] 27 Nov 2018
Monolithic coupling of implicit material point method with finite element method Eugenio Aulisa Giacomo Capodaglio arxiv:1811.11150v1 [physics.comp-ph] 27 Nov 2018 Abstract A monolithic coupling between
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationFinite Deformation Analysis of Dynamic Behavior of Embankment on Liquefiable Sand Deposit Considering Pore Water Flow and Migration
6 th International Conference on Earthquake Geotechnical Engineering 1-4 November 215 Christchurch, New Zealand Finite Deformation Analysis of Dynamic Behavior of Embankment on Liquefiable Sand Deposit
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationUse of Elastic Constitutive Relations in Total Lagrangian Formulation
Topic 15 Use of Elastic Constitutive Relations in Total Lagrangian Formulation Contents: Basic considerations in modeling material response Linear and nonlinear elasticity Isotropic and orthotropic materials
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationFlow and Transport. c(s, t)s ds,
Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationDuality method in limit analysis problem of non-linear elasticity
Duality method in limit analysis problem of non-linear elasticity Igor A. Brigadnov Department of Computer Science North-Western State Technical University Millionnaya 5, St. Petersburg, 191186, Russia
More informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
More informationConstitutive model of brain tissue suitable for finite element analysis of surgical procedures
Journal of Biomechanics 32 (1999 531 537 Technical Note Constitutive model of brain tissue suitable for finite element analysis of surgical procedures Karol Miller* Department of Mechanical and Materials
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More information1 Nonlinear deformation
NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected
More informationPrediction of dynamic behavior of workpieces in ultrasonic plastic welding
Prediction of dynamic behavior of workpieces in ultrasonic plastic welding Takao HIRAI** Fumiyasu KURATANI** Tatsuya YOSHIDA** and Saiji WASHIO*** **Department of Mechanical Engineering, Graduate School
More informationA multiscale framework for lubrication analysis of bearings with textured surface
A multiscale framework for lubrication analysis of bearings with textured surface *Leiming Gao 1), Gregory de Boer 2) and Rob Hewson 3) 1), 3) Aeronautics Department, Imperial College London, London, SW7
More informationFormulation and numerical implementation of micro-scale boundary conditions for particle aggregates
Granular Matter (2017) 19:72 DOI 10.1007/s10035-017-0757-y ORIGINAL PAPER Formulation and numerical implementation of micro-scale boundary conditions for particle aggregates J. Liu 1 E. Bosco 1 A. S. J.
More informationA FINITE ELEMENT FORMULATION FOR NONLINEAR 3D CONTACT PROBLEMS
A FINITE ELEMENT FORMULATION FOR NONLINEAR 3D CONTACT PROBLEMS Federico J. Cavalieri, Alberto Cardona, Víctor D. Fachinotti and José Risso Centro Internacional de Métodos Computacionales en Ingeniería
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationUnconfined Compression of Articular Cartilage: Nonlinear Behavior and Comparison With a Fibril-Reinforced Biphasic Model
M. Fortin Institute of Biomedical Engineering, Ecole Polytechnique, Montreal, Quebec, Canada J. Soulhat A. Shirazi-Adl Department of Mechanical Engineering, Ecole Polytechnique, Montreal, Quebec, Canada
More informationUNIVERSITY OF CALGARY. Nonlinear Elasticity, Fluid Flow and Remodelling in Biological Tissues. Aleksandar Tomic A THESIS
UNIVERSITY OF CALGARY Nonlinear Elasticity, Fluid Flow and Remodelling in Biological Tissues by Aleksandar Tomic A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT OF THE REQUIREMENTS
More informationPerformance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis
Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan
More informationGeometrically exact beam dynamics, with and without rotational degree of freedom
ICCM2014 28-30 th July, Cambridge, England Geometrically exact beam dynamics, with and without rotational degree of freedom *Tien Long Nguyen¹, Carlo Sansour 2, and Mohammed Hjiaj 1 1 Department of Civil
More informationNonlinear analysis in ADINA Structures
Nonlinear analysis in ADINA Structures Theodore Sussman, Ph.D. ADINA R&D, Inc, 2016 1 Topics presented Types of nonlinearities Materially nonlinear only Geometrically nonlinear analysis Deformation-dependent
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationA Phenomenological Theory and Numerical Procedure for Chemo-Mechanical Coupling Behavior of Hydrogel
Copyright 2009 Tech Science Press CMC, vol.12, no.1, pp.39-55, 2009 A Phenomenological Theory and Numerical Procedure for Chemo-Mechanical Coupling Behavior of Hydrogel Q. S. Yang 1, B. S. Liu and L. T.
More informationMulti-Point Constraints
Multi-Point Constraints Multi-Point Constraints Multi-Point Constraints Single point constraint examples Multi-Point constraint examples linear, homogeneous linear, non-homogeneous linear, homogeneous
More informationArbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact
Arbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact Philip P. Garland 1 and Robert J. Rogers 2 1 School of Biomedical Engineering, Dalhousie University, Canada 2 Department
More informationMathematical Background
CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous
More informationMechanics of Materials and Structures
Journal of Mechanics of Materials and Structures A THREE DIMENSIONAL CONTACT MODEL FOR SOIL-PIPE INTERACTION Nelly Piedad Rubio, Deane Roehl and Celso Romanel Volume 2, Nº 8 October 2007 mathematical sciences
More information4 NON-LINEAR ANALYSIS
4 NON-INEAR ANAYSIS arge displacement elasticity theory, principle of virtual work arge displacement FEA with solid, thin slab, and bar models Virtual work density of internal forces revisited 4-1 SOURCES
More informationUsing Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses
Using Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses Brian Baillargeon, Ramesh Marrey, Randy Grishaber 1, and David B. Woyak 2 1 Cordis Corporation,
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationarxiv:q-bio/ v1 [q-bio.qm] 30 Nov 2003
Material forces in the context of biotissue remodelling arxiv:q-bio/0312002v1 [q-bio.qm] 30 Nov 2003 K. Garikipati, H. Narayanan, E. M. Arruda, K. Grosh, S. Calve University of Michigan, Ann Arbor, Michigan
More informationPrediction of the bilinear stress-strain curve of engineering material by nanoindentation test
Prediction of the bilinear stress-strain curve of engineering material by nanoindentation test T.S. Yang, T.H. Fang, C.T. Kawn, G.L. Ke, S.Y. Chang Institute of Mechanical & Electro-Mechanical Engineering,
More informationFinite deformation effects in soft elastohydrodynamic lubrication problems
Finite deformation effects in soft elastohydrodynamic lubrication problems Stanis law Stupkiewicz a,, Jakub Lengiewicz a, Przemys law Sadowski a, Stanis law Kucharski a a Institute of Fundamental Technological
More informationDownloaded from ascelibrary.org by Colorado University at Boulder on 06/14/18. Copyright ASCE. For personal use only; all rights reserved.
onlinear, Large Deformation Finite-Element Beam/Column Formulation for the Study of the Human Spine: Investigation of the Role of Muscle on Spine Stability Franck J. Vernerey 1 and Brian Moran 2 Downloaded
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationSolution of contact problems in linear elasticity using a feasible interior point algorithm for nonlinear complementarity problems
7 th World Congress on Structural and Multidisciplinary Optimization COEX Seoul, May - 5 May 007, Korea Solution of contact problems in linear elasticity using a feasible interior point algorithm for nonlinear
More informationRigid bodies - general theory
Rigid bodies - general theory Kinetic Energy: based on FW-26 Consider a system on N particles with all their relative separations fixed: it has 3 translational and 3 rotational degrees of freedom. Motion
More informationPorous Fibre-Reinforced Materials Under Large Deformations
Porous Fibre-Reinforced Materials Under Large Deformations Salvatore Federico Department of Mechanical and Manufacturing Engineering, The University of Calgary 2500 University Drive NW, Calgary, Alberta,
More information3-dimensional joint torque calculation of compression sportswear using 3D-CG human model
3-dimensional joint torque calculation of compression sportswear using 3D-CG human model Akihiro Matsuda, University of Tsukuba Hirokazu Tanaka, University of Tsukuba Hitoshi Aoki, University of Tsukuba
More informationDeformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera
Karel October Tůma 24, Simulation 2018 of a bovine eye 1/19 Deformation of bovine eye fluid structure interaction between viscoelastic vitreous, non-linear elastic lens and sclera Karel Tůma 1 joint work
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015 Institute of Structural Engineering Method of Finite Elements II 1 Constitutive
More informationLecture 15: Revisiting bars and beams.
3.10 Potential Energy Approach to Derive Bar Element Equations. We place assumptions on the stresses developed inside the bar. The spatial stress and strain matrices become very sparse. We add (ad-hoc)
More informationFinal Project: Indentation Simulation Mohak Patel ENGN-2340 Fall 13
Final Project: Indentation Simulation Mohak Patel ENGN-2340 Fall 13 Aim The project requires a simulation of rigid spherical indenter indenting into a flat block of viscoelastic material. The results from
More informationICES REPORT March R. Kruse, N. Nguyen-Thanh, L. De Lorenzis, and T.J.R. Hughes
ICES REPORT 15-08 March 2015 Isogeometric collocation for large deformation elasticity and frictional contact problems by R. Kruse,. guyen-thanh, L. De Lorenzis, and T.J.R. Hughes The Institute for Computational
More informationReconsideration on the Use of Elastic Models to Predict the. Instantaneous Load Response of the Knee Joint
Prepared for Journal of Engineering in Medicine Submitted on February 3, 011 Revision submitted on April 18, 011 Reconsideration on the Use of Elastic Models to Predict the Instantaneous Load Response
More informationPoroelastic Finite Element Analysis of a Heterogeneous Articular Cartilage Explant Under Dynamic Compression in ABAQUS
Poroelastic Finite Element Analysis of a Heterogeneous Articular Cartilage Explant Under Dynamic Compression in ABAQUS A Thesis Presented to the Faculty of the College of Engineering California Polytechnic
More informationUNIVERSITY OF CALGARY. Viscoelasticity of Articular Cartilage and Ligament: Constitutive Modeling and Experiments. Sahand Ahsanizadeh A THESIS
UNIVERSITY OF CALGARY Viscoelasticity of Articular Cartilage and Ligament: Constitutive Modeling and Experiments by Sahand Ahsanizadeh A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT
More informationLinear Constitutive Relations in Isotropic Finite Viscoelasticity
Journal of Elasticity 55: 73 77, 1999. 1999 Kluwer Academic Publishers. Printed in the Netherlands. 73 Linear Constitutive Relations in Isotropic Finite Viscoelasticity R.C. BATRA and JANG-HORNG YU Department
More informationDebonding process in composites using BEM
Boundary Elements XXVII 331 Debonding process in composites using BEM P. Prochazka & M. Valek Czech Technical University, Prague, Czech Republic Abstract The paper deals with the debonding fiber-matrix
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
More informationA Penetration-Based Finite Element Method for Hyperelastic Three-Dimensional Biphasic Tissues in Contact. Kerem Ün
A Penetration-Based Finite Element Method for Hyperelastic Three-Dimensional Biphasic Tissues in Contact by Kerem Ün A Thesis Submitted to the Graduate Faculty of Rensselaer Polytechnic Institute in Partial
More informationMultiscale modeling and simulation of shock wave propagation
University of Iowa Iowa Research Online Theses and Dissertations Spring 2013 Multiscale modeling and simulation of shock wave propagation Cory Nelsen University of Iowa Copyright 2013 Cory Nelsen This
More informationFETI domain decomposition method to solution of contact problems with large displacements
FETI domain decomposition method to solution of contact problems with large displacements Vít Vondrák 1, Zdeněk Dostál 1, Jiří Dobiáš 2, and Svatopluk Pták 2 1 Dept. of Appl. Math., Technical University
More informationComputational non-linear structural dynamics and energy-momentum integration schemes
icccbe 2010 Nottingham University Press Proceedings of the International Conference on Computing in Civil and Building Engineering W Tizani (Editor) Computational non-linear structural dynamics and energy-momentum
More informationNonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess
Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable
More informationInternational Journal of Pure and Applied Mathematics Volume 58 No ,
International Journal of Pure and Applied Mathematics Volume 58 No. 2 2010, 195-208 A NOTE ON THE LINEARIZED FINITE THEORY OF ELASTICITY Maria Luisa Tonon Department of Mathematics University of Turin
More informationA Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements
A Numerical Study of Finite Element Calculations for Incompressible Materials under Applied Boundary Displacements A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment
More informationContinuum Mechanics Fundamentals
Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are
More informationGlossary. Glossary of Symbols. Glossary of Roman Symbols Glossary of Greek Symbols. Contents:
Glossary Glossary of Symbols Contents: Glossary of Roman Symbols Glossary of Greek Symbols Glossary G-l Glossary of Roman Symbols The Euclidean norm or "two-norm." For a vector a The Mooney-Rivlin material
More informationA surface Cauchy Born model for nanoscale materials
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 26; 68:172 195 Published online 9 May 26 in Wiley InterScience (www.interscience.wiley.com). DOI: 1.12/nme.1754 A surface
More informationStress intensity factor analysis for an interface crack between dissimilar isotropic materials
Stress intensity factor analysis for an interface crack between dissimilar isotropic materials under thermal stress T. Ikeda* & C. T. Sun* I Chemical Engineering Group, Department of Materials Process
More information