含液体松质骨的压缩实验分析

在本文中,假设不含液体松质骨为红弹性体,以及松质骨固液两相结构特性,建立了含液体松质骨的单向压缩本构方程,其形式与三参量线粘弹性体的本构方程相同。并通过拟合实验结果,得到本构方程中有关参数的值。     

Finite Element Simulation of the Fatigue Behaviour of Cancellous Bone*

Meccanica - Fatigue failure of bone has been implicated in a number of clinical failure scenarios. At the material level, the fatigue behaviour of cancellous bone is poorly understood. At the...     

Fluid Flow Analysis of Integrated Porous Bone Scaffold and Cancellous Bone at Different Skeletal Sites: In Silico Study

The dynamic characteristic of bone is its ability to remodel itself through mechanobiological responses. Bone regeneration is triggered by mechanical cues from physiological activities that generate structural strain and cause bone marrow movement. This phenomenon is crucial for bone scaffold when implanted in the cancellous bone as host tissue. Often, the fluid movement of bone scaffold and cancellous bone is studied separately, which does not represent the actual environment once implanted. In the present study, the fluid flow analysis properties of bone scaffold integrated into the cancellous bone at different skeletal sites are investigated. Three types of porous bone scaffolds categorized based on pore size configurations: 1 mm, 0.8 mm and hybrid (0.8 mm interlaced with 0.5 mm) were used. Three different skeletal sites of femoral bone were selected: neck, lateral condyle and medial condyle. Computational fluid dynamics was utilized to analyze the fluid flow properties of bone scaffold integrated cancellous bone. The results of this study reveal that the localization and maximum value of shear stress in an independent bone scaffold are significantly different compared to the bone scaffold integrated with cancellous bone by about 160% to 448% percentage difference. Low shear stress and high permeability were found across models that have higher Tb.Sp (trabecular separation). Specimen C and femoral lateral condyle showed the highest permeability in their respective category.

     

The evolution of Hooke's law due to texture development in FCC polycrystals

Initially isotropic aggregates of crystalline grains show a texture-induced anisotropy of both their inelastic and elastic behavior when submitted to large inelastic deformations. The latter, however, is normally neglected, although experiments as well as numerical simulations clearly show a strong alteration of the elastic properties for certain materials. The main purpose of the work is to formulate a phenomenological model for the evolution of the elastic properties of cubic crystal aggregates. The effective elastic properties are determined by orientation averages of the local elasticity tensors. Arithmetic, geometric, and harmonic averages are compared. It can be shown that for cubic crystal aggregates all of these averages depend on the same irreducible fourth-order tensor, which represents the purely anisotropic portion of the effective elasticity tensor. Coupled equations for the flow rule and the evolution of the anisotropic part of the elasticity tensor are formulated. The flow rule is based on an anisotropic norm of the stress deviator defined by means of the elastic anisotropy. In the evolution equation for the anisotropic part of the elasticity tensor the direction of the rate of change depends only on the inelastic rate of deformation. The evolution equation is derived according to the theory of isotropic tensor functions. The transition from an elastically isotropic initial state to a (path-dependent) final anisotropic state is discussed for polycrystalline copper. The predictions of the model are compared with micro–macro simulations based on the Taylor–Lin model and experimental data.     

Generalized Darcy's Law and the Structure of the Phase and Relative Phase Permeabilities for Two-Phase Flows through Anisotropic Porous Media

For two-phase immiscible fluid flows a generalized Darcy's law is written in invariant tensor form for crystallographic point symmetry groups and anisotropic textures. The representation of the phase permeability coefficient tensors and the structure of the expressions for the relative phase permeabilities are analyzed for all symmetry groups. The relation between the phase and absolute permeability coefficient tensors is specified by a fourth-rank tensor with the external symmetry coinciding with external symmetry of the phase permeability tensors. It is shown that the external symmetry of the phase permeability coefficient tensors can differ from the external symmetry of the absolute permeability tensor. For triclinic and monoclinic symmetry groups it is shown that the phase permeability coefficient tensors may not be coaxial with each other and with the absolute permeability tensor; moreover, the directions of the principal axes of the phase permeability coefficient tensors can depend on the saturation.     

The Anisotropic Lagrangian Averaged Euler and Navier-Stokes Equations

 The purpose of this paper is twofold. First, we give a derivation of the Lagrangian averaged Euler (LAE-α) and Navier-Stokes (LANS-α) equations. This theory involves a spatial scale α and the equations are designed to accurately capture the dynamics of the Euler and Navier-Stokes equations at length scales larger than α, while averaging the motion at scales smaller than α. The derivation involves an averaging procedure that combines ideas from both the material (Lagrangian) and spatial (Eulerian) viewpoints. This framework allows the use of a variant of G. I. Taylor's ``frozen turbulence' hypothesis as the foundation for the model equations; more precisely, the derivation is based on the strong physical assumption that fluctutations are frozen into the mean flow. In this article, we use this hypothesis to derive the averaged Lagrangian for the theory, and all the terms up to and including order α² are accounted for. The equations come in both an isotropic and anisotropic version. The anisotropic equations are a coupled system of PDEs (partial differential equations) for the mean velocity field and the Lagrangian covariance tensor. In earlier works by Foias, Holm & Titi [10], and ourselves [16], an analysis of the isotropic equations has been given. In the second part of this paper, we establish local in time well-posedness of the anisotropic LANS-α equations using quasilinear PDE type methods. (Accepted September 2, 2002) Published online November 26, 2002 Dedicated to Stuart Antman on the occasion of his 60th birthday Communicated by S. MüLLER     

Anisotropic damage mechanics for viscoelastic ice

We present a formulation of continuum damage in glacier ice that incorporates the induced anisotropy of the damage effects but restricts these formally to orthotropy. Damage is modeled by a symmetric second rank tensor that structurally plays the role of an internal variable. It may be interpreted as a texture measure that quantifies the effective specific areas over which internal stresses can be transmitted. The evolution equation for the damage tensor is motivated in the reference configuration and pushed forward to the present configuration. A spatially objective constitutive form of the evolution equation for the damage tensor is obtained. The rheology of the damaged ice presumes no volume conservation. Its constitutive relations are derived from the free enthalpy and a dissipation potential, and extends the classical isotropic power law by elastic and damage tensor dependent terms. All constitutive relations are in conformity with the second law of thermodynamics.PACS 83.60.Df, 62.20.Mk     

圆柱型各向异性弹性力学平面问题

本文对圆柱型各向异性弹性力学平面问题的基本方程进行了改写。在此基础上，导出了应力函数Ｇ和位移函数φ，它们满足相同的控制方程，比文〔１〕的应力函数Ｆ的控制方程要简单，便于求得特解，并有Ｆ＝ｒＧ的关系。还对若干经典问题进行了求解。     

复变形式的各向异性板弯曲问题的基本解

提出了求解各向异性板弯曲问题基本解的新方法。得到的基本解简捷明了，相应的法向弯矩和相当剪力的表达式易求，故便于应用在一般边界条件的各向异性板弯曲问题的边界积分方程。     

Non-linear Model for Wood Saturation

The paper presents a non-linear model of saturation with fluid of anisotropic capillary porous bodies and the results of experimental investigations of wood saturation with methacrylate. The obtained experimental curves illustrating the distribution of methacrylate in wood samples allow the estimation of material coefficients and verification of the theoretical model. The theoretical model is developed based on the balance equations of mass, momentum, and energy, and the thermodynamics of irreversible processes. The non-linear differential equation, describing the distribution of methacrylate content in wood and its evolution during the saturation is solved numerically. The theoretical curves obtained on the basis of both linear and non-linear models are compared with the experimental data and better agreement between them for the non-linear model is stated.     

On Backus Average for Generally Anisotropic Layers

In this paper, following the Backus (in J. Geophys. Res. 67(11):4427–4440, 1962) approach, we examine expressions for elasticity parameters of a homogeneous generally anisotropic medium that is long-wave-equivalent to a stack of thin generally anisotropic layers. These expressions reduce to the results of Backus (1962) for the case of isotropic and transversely isotropic layers.In the over half-a-century since the publications of Backus (1962) there have been numerous publications applying and extending that formulation. However, neither George Backus nor the authors of the present paper are aware of further examinations of the mathematical underpinnings of the original formulation; hence this paper.We prove that—within the long-wave approximation—if the thin layers obey stability conditions, then so does the equivalent medium. We examine—within the Backus-average context—the approximation of the average of a product as the product of averages, which underlies the averaging process.In the presented examination we use the expression of Hooke’s law as a tensor equation; in other words, we use Kelvin’s—as opposed to Voigt’s—notation. In general, the tensorial notation allows us to conveniently examine effects due to rotations of coordinate systems.     

Sharp Interface Limits for Non-Local Anisotropic Interactions

We prove a convexity property of the surface tension corresponding to a non-local, anisotropic free-energy functional of van der Waals type which implies that the Wulff shape is strictly convex and smooth. We also prove that the transport coefficients of the limiting anisotropic motion by mean curvature obtained in [33] are strictly positive and equal to the stiffness parameters determined by the surface tension.     

关于泡沫塑料各向异性模型的修正

本文研究了开孔泡沫塑料的力学各向异性问题。在已有模型的基础上，结合扫描电镜的分析结果提出了对Ｋａｎａｋｋａｎａｔ模型和Ｈｕｂｅｒ－Ｇｉｂｓｏｎ模型的修正。这种修正反映了开孔泡沫塑料支柱的横向尺寸比对泡沫塑料各向异性性质的影响     

Anisotropic slope limiting for discontinuous Galerkin methods

In this paper, we present an anisotropic version of a vertex‐based slope limiter for discontinuous Galerkin methods. The limiting procedure is carried out locally on each mesh element utilizing the bounds defined at each vertex by the largest and smallest mean value from all elements containing the vertex. The application of this slope limiter guarantees the preservation of monotonicity. Unnecessary limiting of smooth directional derivatives is prevented by constraining the x and y components of the gradient separately. As an inexpensive alternative to optimization‐based methods based on solving small linear programming problems, we propose a simple operator splitting technique for calculating the correction factors for the x and y derivatives. We also provide the necessary generalizations for using the anisotropic limiting strategy in an arbitrary rotated frame of reference and in the vicinity of exterior boundaries with no Dirichlet information. The limiting procedure can be extended to elements of arbitrary polygonal shape and three dimensions in a straightforward fashion. The performance of the new anisotropic slope limiter is illustrated by two‐dimensional numerical examples that employ piecewise linear discontinuous Galerkin approximations. Copyright © 2017 John Wiley & Sons, Ltd.     

CAVITATION BIFURCATION FOR COMPRESSIBLE ANISOTROPIC HYPERELASTIC MATERIALS

I. INTRODUCTION The cavitation bifurcation problem, sudden formation and growth of voids in solid materials, haslong attracted much attention because of the fundamental role it plays in the local failure and fractureof materials. For hyperelastic materi…     

Anisotropic model for granulated sea ice dynamics

A continuum model describing sea ice as a layer of granulated thick ice, consisting of many rigid, brittle floes, intersected by long and narrow regions of thinner ice, known as leads, is developed. We consider the evolution of mesoscale leads, formed under extension, whose lengths span many floes, so that the surrounding ice is treated as a granular plastic. The leads are sufficiently small with respect to basin scales of sea ice deformation that they may be modelled using a continuum approach. The model includes evolution equations for the orientational distribution of leads, their thickness and width expressed through second-rank tensors and terms requiring closures. The closing assumptions are constructed for the case of negligibly small lead ice thickness and the canonical deformation types of pure and simple shear, pure divergence and pure convergence. We present a new continuum-scale sea ice rheology that depends upon the isotropic, material rheology of sea ice, the orientational distribution of lead properties and the thick ice thickness. A new model of lead and thick ice interaction is presented that successfully describes a number of effects: (i) because of its brittle nature, thick ice does not thin under extension and (ii) the consideration of the thick sea ice as a granular material determines finite lead opening under pure shear, when granular dilation is unimportant.     

Modelling Slight Compressibility for Hyperelastic Anisotropic Materials

In order to avoid the numerical difficulties in locally enforcing the incompressibility constraint using the displacement formulation of the Finite Element Method, slight compressibility is typically assumed when simulating the mechanical response of nonlinearly elastic materials. The current standard method of accounting for slight compressibility of hyperelastic materials assumes an additive decomposition of the strain-energy function into a volumetric and an isochoric part. A new proof is given to show that this is equivalent to assuming that the hydrostatic stress is a function of the volume change only and that uniform dilatation is a possible solution to the hydrostatic stress boundary value problem, with therefore no anisotropic contribution to the mechanical response. An alternative formulation of slight compressibility is proposed, one that does not suffer from this defect. This new model generalises the standard model by including a mixed term in the volume change and isochoric response. Specific models of slight compressibility are given for isotropic, transversely isotropic and orthotropic materials.     

Sharp Sobolev Asymptotics for Critical Anisotropic Equations

We investigate blow-up theory and prove sharp Sobolev asymptotics for a general class of anisotropic critical equations in bounded domains of .     

Anisotropic unstructured mesh adaption for flow simulations

Three new ideas for anisotropic adaption of unstructured triangular grids are presented, with particular emphasis on fluid flow computations. © 1997 John Wiley & Sons, Ltd.     

Anisotropic Elasticity and Multi-Material Singularities

Multi-material wedges associated with convergence of geometrical and material discontinuity lines generally show singular stress fields around the vertex of the wedge. In this paper, the eigenvalue problem for a multi-material wedge composed of several anisotropic elastic sectors is formulated in a completely generally manner, including the cases of degenerate and extra-degenerate material sectors, and various types of edge conditions for both open and closed wedges. General representation of the elasticity solution in a degenerate or extra-degenerate anisotropic sector requires higher-order eigenmodes (generalized eigenfunctions) in addition to zeroth-order eigenmodes. Such higher-order eigenmodes are obtained from appropriate analytical expressions of the zeroth-order eigenmode by using the derivative rule. The analysis is applied to one bisector wedge and one trisector wedge in a three-layer cracked composite model to obtain accurate elasticity solutions of the singular stress fields. These solutions were determined using the traction data generated on a circular collocation path by a conventional finite element analysis.