动脉管壁切应力的确定

血管为了适应所处力学环境的变化会产生结构和功能的变化,血液流动作用于血管壁上的切应力在这种血管重建中起着重要的作用.目前直接测量活体血管壁切应力存在着许多技术上的困难.通过分析动脉中血液脉动流的特性,提出一种利用测量管轴上的血流速度,计算血管壁切应力的方法,为定量确定血管壁切应力,进而讨论血管壁切应力对血管重建的影响提供必要的手段.     

弹性动脉血管中血液流动特性的模拟和分析

研究了两个不同的非牛顿血液流动模型:低粘性剪切简单幂律模型和低粘性剪切及粘弹性振荡流的广义Maxwell模型．同时利用这两个非牛顿模型和牛顿模型,研究了磁场中刚性和弹性直血管中血液的正弦型脉动．在生理学条件下,大动脉中血液的弹性对其流动性态似乎并不产生影响,单纯低粘性剪切模型可以逼真地模拟这种血液流动．利用高剪切幂律模型模拟弹性血管中的正弦型脉动流,发现在同一压力梯度下,与牛顿流体相比较,幂律流体的平均流率和流率变化幅度都更小．控制方程用Crank-Niclson方法求解．弹性动脉中血液受磁场作用是产生此结果的直观原因．在主动脉生物流的模拟中,与牛顿流体模型比较,发现在匹配流率曲线上,幂律模型的平均壁面剪切应力增大,峰值壁面剪切应力减小．讨论了弹性血管横切磁场时的血液流动,评估了血管形状和表面不规则等因素的影响．     

主动脉弓及分支血管内非稳态血流分析

运用流体力学中的三维非定常Navier-Stokes方程作为血液流动的控制方程,并采用计算流体力学方法对人体主动脉弓及分支血管内非Newton(牛顿)血液黏度模型下血流进行瞬态数值模拟.分析了一个心动周期内不同时刻血流动力学特征参数的分布对动脉粥样硬化斑块形成的影响,并与Newton血液黏度模型下的血管壁面压力和壁面切应力特征参数进行对比.结果表明:与Newton血液模型相比,非Newton血液模型下血流分布更符合真实血流特性;在心动收缩期,分支血管外侧壁附近存在面积较大的低速涡流区,该区域内血管壁面压力与壁面切应力具有较大的变化量,血液中的血小板、脂质和纤维蛋白等易沉积,血管内壁易疲劳损伤并发生血管重构,促使动脉粥样硬化斑块形成;而在心动舒张期,分支血管内血流速度分布均匀,血管壁面压力与壁面切应力变化量较小,血管壁受到较小的应力作用,对动脉粥样硬化斑块形成的作用较小.     

肾动脉变窄和流-固结构相互作用对非Newton脉动流的影响——一个真实的腹部主动脉和肾动脉模型

研究肾动脉狭窄(RAS)对血液流动和血管壁的影响.根据CT扫描图像,重建腹部主动脉和肾动脉的解剖模型,通过模型的脉动流进行了仿真计算,计算中考虑了流体-固体结构的相互作用(FSI).研究RAS对血管壁剪切应力和位移的影响,RAS使得肾动脉中流量减少,肾素-血管紧缩素系统可能被激活,从而导致严重的高血压.     

锥形血管入口区域内管壁与血液的耦合运动

本文研究了锥形血管入口区域内血管壁与血液间的耦合问题。对具有锥度角的弹性血管入口区域内的管壁运动和血液流动建立的相互耦合作用的数学模型，在满足相应的边界条件下求得了一组血液流动的速度分布公式、压力分布公式以及管壁运动公式，得出了一些重要的结论。     

存在周围粘弹组织约束时动脉管内的脉动流

由于在动脉管外侧牢牢地粘附着结缔组织,而且这些结缔组织是一种粘弹体。为了分析周围粘弹组织对动脉中脉动流的影响,本文将周围结缔组织考虑为Voigt粘弹体,由动脉管内血液流动方程(Navier-Stokes方程)与管壁运动方程(Lamb方程)导出动脉管中脉搏波波速的一般表达式与管内脉动流的速度表达式。Womersley关于不存在周围组织约束的结果[4]与仅考虑周围组织是纯弹性体的结果[6]都做为本文结果的特例而包含在本文的结果之中。     

动脉中脉搏波传播分析

将血管简化为弹性管,并考虑组织对血管壁的约束,利用力学方法建立血液流过血管的力学模型．通过理论分析对脉搏波在血管中的传播规律进行研究,同时分析了血液粘性、血管壁弹性模量、管径对波的传播的影响．通过对考虑血液粘性和不考虑血液粘性的结果比较,发现血液的粘性对脉搏波的传播的影响不能忽略,并且当弹性模量增大时,传播速度增大,血流的压力值增高;血管直径减小时,血流压力也增高,脉搏波速度增大．理论分析得到的结果也有助于利用脉搏波的信息来分析和辅助诊断一些人体疾病的病因．     

血液动力学中血管流激波与驻波的相互作用

血液动力学问题是生物力学心血管系统中的重要研究课题.血管内斑块处,血管截面和血管壁的材质发生变化,对血液流动产生重要影响.血液流动中基本波及其相互作用对探究血液流动的规律、生理学意义及与疾病的关系有着重要的意义.本文研究血液动力学血液流动简化数学模型的基本波的相互作用.血管流模型是3×3非严格双曲型方程组.构造性地得到了初值为三段常状态时,血管流问题的解,即解决了激波与驻波的相互作用问题.特别地,给出四种后前激波与驻波的相互作用的结果.     

低频交变磁场对于刚性直圆管中脉动流的影响*

本文通过求解灾变磁场作用下刚性直圆管脉动流的运动方程,得到了它的分析解.计算了流速分布及阻抗.计算结果对于深入了解低频磁场对于血液动力学的影响以及它的临床应用具有一定参考价值.     

计及血浆层的局部狭窄血管内血液脉动流分析

本文研究了血浆层的存在对局部狭窄血管内脉动流动特性的影响。分析表明,在狭窄区域内,血浆层的存在明显地改变了血液流动的速度分布、纵向阻抗和压力梯度等特性,但几乎不改变纵向阻抗和压力梯度的幅角;而且这种影响与Womersley数、狭窄程度等因素有密切关系。分析血浆层影响下动脉狭窄流的特性,对于认识动脉粥样硬化的病理特性有重要意义。     

Dynamic response of pulsatile flow of blood in a stenosed tapered artery

The aim of this paper is to throw some light on the rheological study of pulsatile blood flow in a stenosed tapered arterial segment. Arterial wall is considered to be rigid and flexible separately for improving the similarity to the in vivo situation. The streaming blood is considered to be Newtonian. The governing nonlinear equations of motion are sought using the well‐known stream function‐vorticity method and are solved numerically by finite difference technique. Important rheological parameters, such as axial velocity component, wall shear stress, and flow separation region are estimated in the neighborhood of the stenosis. Effects of stenosis height, vessel tapering, and wall flexibility on the blood flow are investigated properly and are explained in detail through their graphical representations.     

精确计算Bingham流体圆管层流压降及速度分布规律的新方法

Bingham(宾汉)模型情况下,多采用通用公式进行圆管层流压降的解析计算,即将Bingham模型本构方程代入粘性流体圆管层流流动通用公式进行计算,仅能得到压降的解析解.新方法结合Bingham流体本构方程与运动方程,建立有关力学平衡方程,并运用代数方程的根式解理论对圆管层流流动时的非线性方程进行求解,可直接求得Bingham流体圆管层流压降及速度流核区半径的解析解,进一步可求得圆管层流速度解析解;Bingham流体圆管层流速度的直接影响因素为流量、塑性粘度和屈服值,研究发现速度流核宽度与屈服值成正比,与流量及塑性粘度成反比,且流核的宽度越大,流核区的速度越小.     

Initial-boundary-value problem of hyperbolic equations for blood flow in a vessel

In this paper, we study an initial-boundary-value problem of a system of hyperbolic, partial-differential equations that models blood flow in a vessel. The one-spatial-dimensional model assumes that blood flow in the vessel is an incompressible, homogeneous, Newtonian fluid which has a small Womersley number. Boundary conditions with either the pressure or the flow rate at each end of the vessel are considered, and the existence of the global solution is obtained using a form of Glimm's finite-difference scheme.     

血液流动与血管壁运动

本文讨论了哺乳动物循环系统的血液流动与血管壁运动之间的相互作用问题.在假定流动处于稳定的振荡流动情况下,导得了一组血液流动速度分布公式,压力分布公式以及约束应力公式,管壁位移公式.把Kuchar的公式从定常流动情况推广到非定常的振荡流动情况.文中还讨论了动脉血管壁的弹性效应问题.     

Compliance of a Biocomposite Vascular Tissue in Longitudinal and Circumferential Directions as a Basis for Creating Artificial Substitutes

Experimental studies on the mechanical behavior of 11 human common carotid arteries at different values of internal pressure and axial force were performed on a device allowing us to measure the internal pressure, axial force, and circumferential and longitudinal deformations of the vessel. The persons age ranged from 20 to 28 years. Two types of experiments were carried out. In the first series, cylindrical samples were gradually loaded by an internal pressure from 0 to 200 mmHg at different longitudinal stretch ratios. The second series included axial extension of the same samples at different circumferential stretch ratios. The undulation level of elastin membranes in longitudinal and circumferential directions for samples fixated at different values of internal pressure and longitudinal stretch ratio were determined from histological data.     

Mathematical analysis of the quasilinear effects in a hyperbolic model blood flow through compliant axi‐symmetric vessels

In this paper, we present a mathematical analysis of the quasilinear effects arising in a hyperbolic system of partial differential equations modelling blood flow through large compliant vessels. The equations are derived using asymptotic reduction of the incompressible Navier–Stokes equations in narrow, long channels. To guarantee strict hyperbolicity we first derive the estimates on the initial and boundary data which imply strict hyperbolicity in the region of smooth flow. We then prove a general theorem which provides conditions under which an initial–boundary value problem for a quasilinear hyperbolic system admits a smooth solution. Using this result we show that pulsatile flow boundary data always give rise to shock formation (high gradients in the velocity and inner vessel radius). We estimate the time and the location of the first shock formation and show that in a healthy individual, shocks form well outside the physiologically interesting region (2.8m downstream from the inlet boundary). In the end we present a study of the influence of vessel tapering on shock formation. We obtain a surprising result: vessel tapering postpones shock formation. We provide an explanation for why this is the case. Copyright © 2003 John Wiley & Sons, Ltd.     

Mathematical modelling of blood flow through an overlapping arterial stenosis

Of concern in the paper is a theoretical study of blood flow in an arterial segment in the presence of a time-dependent overlapping stenosis using an appropriate mathematical model. A remarkably new shape of the stenosis in the realm of the formation of the arterial narrowing caused by atheroma is constructed mathematically. The artery is simulated as an elastic (moving wall) cylindrical tube containing a viscoelastic fluid representing blood. The unsteady flow mechanism of the present investigation is subjected to a pulsatile pressure gradient arising from the normal functioning of the heart. The equations governing the motion of the system are sought in the Laplace transform space and their relevant solutions supplemented by the suitable boundary conditions are obtained numerically in the transformed domain through the use of an appropriate finite difference technique. Laplace inversion is also carried out by employing numerical techniques. A thorough quantitative analysis is performed at the end of the paper for the flow velocity, the flux, the resistive impedances, and the wall shear stresses together with their variations with the time, the pressure gradient, and the severity of the stenosis in order to illustrate the applicability of the present mathematical model under consideration.     

Axisymmetric deformation of cylinders made of homogeneous and fiber-reinforced elastic materials in butt-end torsion

The equations of axisymmetric deformation of a circular cylinder made of homogeneous and weakly fiber-reinforced materials in a macroscopically two-dimensional problem statement are presented. The fibers are placed in the axial, circumferential, and radial directions and along spirals in the cylinder wall. They also can lie along opposite spirals in adjacent cylindrical layers, forming an angle-ply reinforcement scheme. At a low content of fibers, it is assumed that they are in the uniaxial stress state, determined by the axial tension or constraint compression in the deformed cylinder. Based on this mathematical model, the deformation of homogeneous and fiber-reinforced elastic cylinders of various length, with free outer surfaces, in butt-end torsion is investigated. The effect of length of the cylinders on their deformation character in torsion is evaluated. The torsion of the cylinders under the conditions of sliding fit along the inner surface and under the action of pressure on the surface, when the initial and deformed configurations are practically congruent, is considered.     

The limit equilibrium of an anisotropic medium under the general plasticity condition

A theory of the limit equilibrium of an anisotropic medium under the general plasticity condition in the plane strain state is developed. The proposed yield criterion (the limit equilibrium condition) is obtained by combining the von Mises–Hill yield criterion of an ideally plastic anisotropic material and Prandtl's limit equilibrium condition for a medium under the general plasticity law. It is shown that the problem is statically determinate, i.e., if the boundary conditions are specified in stresses, the stress state in plastic region can only be obtained using equilibrium equations. It is established that the equations describing the stress state are hyperbolic and have two families of characteristic curves that intersect at variable angles. In deriving the equations describing the velocity field, the material is assumed to be rigid plastic, and the associated law of flow is applied. It is shown that the equations for the velocities are also hyperbolic, and their characteristic curves are identical with those of the equations for stresses. However, the directions of the principal values of the stress and strain rate tensors are different due to the anisotropy of the material. The characteristic directions differ from the isotropic case in that the normal and tangential components of the stress tensor do not satisfy the limit conditions. It is established that the equations obtained allow of partial solutions, and in this case, at least one family of characteristic curves consists of straight lines. The conditions along the lines of discontinuity of the velocity are investigated, and it is shown that, as in the isotropic case, these are characteristic curves of the system of governing equations. In the anisotropic formulation, the well-known Rankine problem of the limit state of a ponderable layer is solved. From an analysis of the velocity field it is shown that plastic flow of the entire layer is possible only for a slope angle equal to the angle of internal friction. For slope angles less than the angle of internal friction, the solutions obtained are solutions of problems of the pressure of the medium on the retaining walls. The change in this pressure as a function of the parameters of anisotropy is investigated, and turns out to be significant.