一种骨小管中液体流动产生的流量及切应力模型

骨组织受力变形后其内部液体就会流动,同时在其微观结构——骨单元壁中扩散,并进一步产生一系列与骨液流动相关的物理效应,如流体剪切应力、流动电位等,这些物理效应被细胞感知并做出破骨或成骨等反应,来使骨适应外部载荷环境.鉴于骨组织产生的内部液体流动很难实验测定,理论模拟是目前的主要研究手段.基于骨单元的多孔弹性性质建立了骨小管内部液体的流动模型,该模型将骨单元所受的外部载荷与骨小管内部液体的压力、流速、流量和切应力联系起来,并进一步可以研究其力传导与力电传导机制.骨小管模型的建立分别基于中空和考虑哈弗液体的骨单元模型,并考虑了骨单元外壁的弹性约束和刚性位移约束两种边界条件.最终得到骨单元在外部轴向载荷作用下,骨小管内部液体的流量及流体切应力的解析解.结果表明:骨小管中的液体流量与流体切应力都正比于应变载荷幅值和频率,并由载荷的应变率决定.因此应变率可以作为控制流量和流体切应力的一种生理载荷因素.流量随着骨小管半径的增大而非线性增大,而流体切应力则随着骨小管半径的增大而线性增大.此外,在相同的载荷下,含哈弗液体的骨单元的模型中,骨小管中液体的流量和切应力均大于中空骨单元模型.     

粘弹性结合材料界面端奇异场

在电子封装等结构中存在大量的粘弹性界面问题,其破坏一般均始于界面端,但目前尚无关于粘弹性界面端奇异场的解.粘弹性问题在拉普拉斯域内与弹性问题有对应关系,理论上可以利用对应性原理由弹性解经拉氏逆变换得到粘弹性问题的解.但是,对于粘弹性界面端,由于奇异场的奇异指数也是与时间有关的,因此进行严密的拉氏逆变换是非常困难的.本文借鉴弹性界面端奇异场,近似地给出了线性粘弹性体界面端奇异场的具体形式,并通过数值计算验证了近似理论解的有效性.     

J.L.KATZ教授在成都和西安作生物力学学术报告

美国纽约伦塞勒工学院(Rensselaer Polytechnic Institute)生物医学工程中心和物理系主任,生物力学杂志编委,J. Lawrence Katz教授应成都科技大学副校长康振黄教授的邀请,于1982年9月26日至10月14日先后访问了北京、成都和西安。9月28日至10月9日在成都科技大学连续作了八次报告。题目是:1.骨的结构.2.骨的粘弹性理论,3.超声波在骨中的传播(一),4.超声波在骨中的传播(二),5.骨的复合材料模型,弹性和粘弹性,6.骨的声发射,7.与骨直接粘合的生物材料,8.骨的声扫描;回顾与展望。10月12日在西安作了题为“骨的结构与声发射      

关于粘弹性结构和复合结构有限元动力学方程的探讨

对粘弹性结构有限元动力学方程和复合结构有限元动力学方程进行了探讨。对于粘弹性结构有限元动力学方程，研究了表示成MCK形式的可能性，并对几种典型模型进行了分析研究。根据对阻尼的假设，给出了考虑弹性材料阻尼及结构外部粘性阻尼的拉氏域内弹性－－粘弹性复合结构有限元动力学方程。     

关于复合结构振动分析中粘弹性材料本构方程的形式

弹性——粘弹性材料复合结构的振动分析必然要涉及粘弹性材料的本构方程形式。就两种常用的形式:标准形式和分数阶导数形式,本文从粘弹性材料的力学性能曲线拟合和复合结构振动分析两方面比较了这两种本构方程形式的优劣,认为采用标准形式具有更多的优越性。     

加载速率对骨的流动电位的影响

骨受力变形时所引起的压力驱动骨内哈佛氏管或骨小管内的液体流动,是骨内出现流动电位的外部原因.由于人体经常受动态载荷作用,研究骨在动态加载过程中流动电位的变化,更有助于了解骨细胞周围电环境的性质.为此,设计了流动电位测试系统,在骨试样两端施加梯形压力波.实验测量了在不同加载速率下流动电位的变化波形.结果显示,加载速率越高,流动电位有减小的趋势.分析认为这是由于微管中局部产生的湍流引起的.     

树脂基复合材料板的粘弹性损伤本构关系

1 引言一般树脂基复合材料板具有相当强烈的各向异性和非均匀性,受载后很容易发生基体裂纹群等损伤.同时,即使在常温下,这类材料也显示粘弹性.而且,损伤与粘弹性都是各向异性的.因此,复合材料力学响应的分析比均匀无损的粘弹性材料要复杂得多,困难得多.实际上,在复合材料的结构强度和尺寸稳定性设计中,它的时间相关性和存在损伤是两个不能回避的重要问题.为此,特别需要复合材料粘弹性损伤本构关系的知识.最近一个时期,复合材料的粘弹性本构关系已得到一定的研究.作者曾提出适用于复合材料分析的弹脆性损伤模型以及考虑损伤的粘弹性本构关系.在此基础     

G/E材料动态力学性能的应力波实时测定

线性粘弹性复合材料的动态力学性能测定不同于弹性材料静态力学性能的测定,而应力波方法又是不同于其它测定粘弹性本构的方法.本文用纵波沿 G/E长杆中的传播,测出了它的拉压复模量,同时也引建了一套完善的应力波实时测定系统.     

切应力对牛骨蠕变变形的影响研究

骨尤其是湿骨,在恒定载荷作用下会发生蠕变变形。为了确定切应力是否影响骨的蠕变变形,采用对骨薄板试样分别施加集中载荷和均布载荷的方式,测量试样挠度实时的变化曲线。结果显示,在载荷恒定时,骨试样的挠度随时间不断增加,体现了典型的蠕变特性。集中荷载下骨的蠕变变形远大于均布荷载下骨的蠕变变形,湿骨的蠕变位移比干骨高近7倍。分析认为,对试样粘弹性性质的影响不仅有正应力的作用,也有切应力的作用;切应力产生的蠕变变形约为正应力所产生蠕变变形的0.85倍。     

非线性粘弹性平面问题的数学模型及其失记效应

本文以位移为基本未知量,利用非线性粘弹性力学中的Leaderman本构关系和线性几何假设,建立了非线性粘弹性平面问题的数学模型;在粘弹性泊松比为常数的情况下,利用Titchmarsh定理和Laplace变换法证明了非线性粘弹性平面问题与非线性弹性平面问题之间存在着某些对应关系,对应关系为粘弹性问题的求解提供了一种新的思路,利用这些关系可直接从相应弹性问题获得粘弹性平面问题的部分响应,与传统的时域有限差分法相比,计算时间明显缩短,另外,对应关系也揭示了粘弹性结构的失记效应,即结构的部分响应仅与外部输入的现时值有关,而与其历史无关。     

Poroelastic behaviors of the osteon: A comparison of two theoretical osteon models

In the paper, two theoretical poroelastic osteon models are presented to compare their poroelastic behaviors, one is the hollow osteon model (Haversian fluid is neglected) and the other is the osteon model with Haversian fluid considered. They both have the same two types of impermeable exterior boundary conditions, one is elastic restraint and the other is displacement constrained, which can be used for analyzing other experiments performed on similarly shaped poroelastic specimens. The obtained analytical pressure and velocity solutions demonstrate the effects of the loading factors and the material parameters, which may have a significant stimulus to the mechanotransduction of bone remodeling signals. Model comparisons indicate: (1) The Haversian fluid can enhance the whole osteonal fluid pressure and velocity fields. (2) In the hollow model, the key loading factor governing the poroelastic behavior of the osteon is strain rate, while in the model with Haversian fluid considered, the strain rate governs only the velocity. (3) The pressure amplitude is proportional to the loading frequency in the hollow model, while in the model with Haversian fluid considered, the loading frequency has little effect on the pressure amplitude.     

MECHANLCAL FAILURE IN HAVERSIAN SYSTEMS ON THE BASIS OF MICROSTRUCTURE

The aim of this paper is to observe variations in the mechanical properties of Haversian systems including Haversian lamellae and canals,by means of the compressive test of five bone specimens subjected to different loads corresponding to 2KN,3KN,dKN,5KN and 6KN respectively. From a series of related microscopic pictures it is clearly seen that:(1) in the lower level of loading (2KN) ,the Haversian canals much the same as the one unsubjected to loads,and conditions of Harversian lamellae are also the same,but,the individual Harversian canal was slightly bent; (2) in the case of the higher levels of loading,it is first deformed in the weaked portion of the Haversian lamellae; (3) with increasing loads,the fracture shapes of the Haversian system are very complex,but the fractures always take place in the cement line between osteons.At the same time,deformed dimensions of the Harversian lamellae and canals were measured. The observed phenomena were qualitatively interpreted by the theory of the linear viscoelas-ticity.     

Transverse isotropic poroelastic osteon model under cyclic loading

The poroelastic problem associated with a hollow cylinder under cyclic loading is solved. This cylinder models an osteon, basic unit of cortical bone. Both fluid and solid phases are supposed compressible. Solid matrix is modeled as an elastic transverse isotropic material. An explicit close-form solution for the steady state is obtained. Fluid flow distribution as a function of poroelastic properties and cyclic loading is discussed as it could influence bone remodeling. Strain rate of loading is shown to play a significant role in mass flux in the porous material.     

微振激励下黏弹性阻尼器微观链结构力学模型

减小微振动对高精密仪器至关重要,利用黏弹性阻尼器进行微振动抑制是一个新兴而又具有挑战性的课题.本文采用分子链网络模型方法分析了黏弹性材料的微观分子链结构,综合考虑材料分子链结构中的网络链和自由链对黏弹性材料力学性能的影响,提出一种基于材料微观分子链结构的微振激励下黏弹性阻尼器力学模型.模型分别采用标准线性固体模型和Maxwell模型来描述网络链和自由链中单个链的力学性能,并分别采用8链网络模型和3链网络模型考虑两种类型分子链的综合效应,引入温频等效原理描述温度对微振激励下黏弹性阻尼器力学性能的影响.该模型能够描述温度和频率对黏弹性阻尼器动态力学性能的影响,并能够反映黏弹性材料的微观结构与材料力学性能的关系.为验证所提模型的有效性及考察黏弹性阻尼器在微振激励下的耗能能力和动态力学性能,在微振条件下对黏弹性阻尼器进行了动态力学性能试验.研究结果表明黏弹性阻尼器具有较好的微振耗能能力,其动态力学性能受温度和频率影响较大,所提的力学模型能够精确地描述微振激励下黏弹性阻尼器动态力学性能随温度和频率的变化关系.      

骨陷窝-骨细胞形状和方向对骨单元内液体流动行为的影响

骨组织内的流体流动不仅为骨细胞的生存提供了充足营养供应及代谢物排放途径,也在骨重建过程中起到关键作用. 为了更精确地阐明骨内液体流动的具体形式,这项研究利用骨陷窝-骨细胞的密度,形态和方向等参数来计算骨单元内液体的流动行为. 首先,计算出不同形状和方向的骨陷窝周围骨小管的数量及分布情况,其次利用算出的参数以及骨组织其他微结构数据来估计骨组织的渗透率和孔隙率等参数,最后根据计算所得的参数建立骨单元的多孔弹性力学有限元模型,并分析了在轴向位移载荷作用下骨陷窝形状和方向对骨单元内液体渗流行为的影响. 结果表明,在所研究的参数范围内不同骨单元模型的相同区域上,骨陷窝形状影响下的骨单元最大压力和流速比最小的分别增加了86%和18%;骨陷窝方向影响下的最大压力和流速比最小的分别增加了125%和56%. 伸长形骨陷窝对单个骨单元局部压力的影响远大于扁平形和圆形骨陷窝. 骨陷窝从0°绕$x$轴旋转到90°过程中压力是逐渐降低的,且30°,45°和60°的模型对骨单元内局部流速有显著影响. 该模型表示骨陷窝的形状和方向以及骨小管的三维分布对骨单元内液体压力和流速幅值及沿不同方向的流动差异有显著的影响. 这项研究将有助于精确量化描述骨内液体的流体行为.      

Drained and Undrained Poroelastic Properties of Healthy and Pathological Bone: A Poro-Micromechanical Investigation

Poro-micromechanics allows for the quantification of poroelastic properties such as the Biot and Skempton coefficients, once a continuum micromechanics model for the material under consideration has been developed and validated. Employing such a model for the transversely isotropic elasticity of cortical and trabecular bone, we determine the tensors of Biot and Skempton coefficients as functions of the volume fractions of mineral, collagen, and the micropore space (Haversian and Volkmann canals, and the inter-trabecular space). Increase of microporosity, as experienced in osteoporosis, as well as decrease of mineral content, as experienced in osteomalacia, lead to an increase of Biot and Skempton coefficients, i. e. to magnification of the mechanical role of the marrow filling the micropore space. For quantification of the marrow pressure rise upon downfall, undrained conditions are appropriate, as can be shown by model predictions of non-destructive impact experiments.     

Cylindre transverse isotrope poroélastique chargé cycliquement : application à un ostéonCyclic loading of a transverse isotropic poroelastic cylinder: a model for the osteon

The poroelastic problem associated with a hollow cylinder under cyclic loading is solved. Both fluid and solid phases are supposed compressible. Solid matrix is modeled as an elastic transverse isotropic material. An explicit close-form solution for the steady state is obtained. This cylinder is considered as a model for an osteon, the basic unit of cortical bone. The fluid flow distribution as a function of poroelastic properties and cyclic loading is discussed, as this could influence bone remodeling. To cite this article: A. Rémond, S. Naili, C. R. Mecanique 332 (2004).     

Modeling and measuring visco-elastic properties: From collagen molecules to collagen fibrils

Collagen is the main structural protein in vertebrate biology, determining the mechanical behavior of connective tissues such as tendon, bone and skin. Although extensive efforts in the study of the origin of collagen exceptional mechanical properties, a deep knowledge of the relationship between molecular structure and mechanical properties remains elusive, hindered by the complex hierarchical structure of collagen-based tissues. Understanding the viscoelastic behavior of collagenous tissues requires knowledge of the properties at each structural level. Whole tissues have been studied extensively, but less is known about the mechanical behavior at the submicron, fibrillar and molecular level. Hence, we investigate the viscoelastic properties at the molecular level by using an atomistic modeling approach, performing in silico creep tests of a collagen-like peptide. The results are compared with creep and relaxation tests at the level of isolated collagen fibrils performed previously using a micro-electro-mechanical systems platform. Individual collagen molecules present a non-linear viscoelastic behavior, with a Young's modulus increasing from 6 to 16 GPa (for strains up to 20%), a viscosity of 3.84±0.38 Pa s, and a relaxation time in the range of 0.24–0.64 ns. At the fibrils level, stress–strain–time data indicate that isolated fibrils exhibit viscoelastic behavior that could be fitted using the Maxwell–Weichert model. The fibrils showed an elastic modulus of 123±46 MPa. The time-dependent behavior was well fit using the two-time-constant Maxwell–Weichert model with a fast time response of 7±2 s and a slow time response of 102±5 s.     

黏弹性材料等效分数阶微观结构标准线性固体模型

从黏弹性材料微观链结构出发,以橡胶基黏弹性材料超弹性理论分子网链高斯(Gauss)统计模型和黏滞流动理论为基础,研究黏弹性材料的微观结构、填料等对黏弹性性能的影响.用温频等效原理描述温度对黏弹性材料力学性能的影响,建立了可以有效描述黏弹性材料耗能特性的等效分数阶微观结构标准线性固体模型.采用动态热机械分析仪(DMA)对高聚物黏弹性材料力学性能、耗能能力进行测试.试验表明:在低温区域,储能模量较大,随着温度的升高,储能模量下降显著;能量损耗因子在高温和低温区域数值较小,在玻璃化转变温度附近数值较高.根据测试数据对所提等效分数阶微观结构标准线性固体模型进行验证,该力学模型能够较好地描述黏弹性材料储能模量和能量损耗因子随温度的变化趋势.用9050A和ZN22黏弹性材料对模型的有效性进一步验证,结果表明:9050A和ZN22黏弹性材料具有较好的耗能能力,所提出的等效分数阶微观结构标准线性固体模型能够准确地描述微观结构和填料对黏弹性材料宏观性能的影响,能够准确地描述黏弹性材料在不同温度和频率下的动态力学性能.