弯曲动脉壁非线性弹性力学性质的理论分析

血管壁非线性力学特性是探索心血管中的血液流动规律及脉搏波传播现象的一个重要前提．在血管力学研究中,直管动脉壁的本构方程已有相当深人的研究,唯独象主动脉弓那样的弯管动脉壁本构方程;至今还没有建立一个理论模型．文中在已有的直管动脉壁本构方程研究基础上,提出一个理论方法来分析弯曲动脉壁的非线性弹性力学性质,在弯曲动脉壁被模拟为均质、正交各向异性、不可压缩材料的前提下,作者从理论上建立了一个表达弯管动脉壁非线性弹性性质的三维ｅ指数型本构方程文中还探讨了弯曲动脉壁内的残余应变分析．     

旋转,弯曲,扩散流道中流固两相流场的数值模拟及其分析

本文简要介绍了旋转,弯曲,扩散流道中,液－固两相流流数值计算结果,进而从理论上分析了泵叶轮压力面遭受磨损和输送两相流时扬程度降低的原因。     

一种确定均匀动脉壁面切应力的非线性方法

从Ling和Atabek提出的``局部流'理论出发,提出一种利用测量血液黏度、管轴上的血流速度、压力和管径波形计算均匀动脉管壁切应力的非线性方法. 将这种方法与柳兆荣等提出的利用测量血液黏度、管轴上的血流速度和平均管径计算切应力的线性方法比较,结果表明,当管壁脉动幅度较小时,两种方法计算的压力梯度、流速剖面和管壁切应力差别较小;而当管壁脉动幅度增大时,两种方法计算的压力梯度、流速剖面和管壁切应力差别增大.对于小幅脉动均匀动脉,用线性方法计算管壁切应力有较高的精度;而对于大变形均匀动脉,则需要考虑非线性因素对管壁切应力的影响. 由于作为输入量的血液黏度、轴心血流速度、压力波形和管径波形可在活体上通过无损伤或微损伤的检测方法得到,所提出的计算切应力的方法为在体或离体研究切应力与动脉重建的关系提供了方法学基础.     

动脉狭窄对血液流速的影响

为了定量计算动脉局部狭窄对动脉管中血液流动速度的影响，本文分别对狭窄区域内定常流和非定常流动进行了求解，得出了狭窄区域内定常流和脉动流的速度表达式。本文将均匀段的流速形经Ｆｏｕｒｉｅｒ分解成定常和脉动两部分，然后分别计算出狭窄区域内对应的定常和脉动流速，经Ｆｏｕｒｉｅｒ合成还原成流速时域波形，同时针对各种情况将不同狭窄对不同的流速波形的作了分析比较。     

动脉局部狭窄时脉动流的有限元分析

本文利用有限元方法研究动脉局部狭窄下的脉动流流场,重点考查在50％与80％面积狭窄下的速度分布、压力分布、壁面剪应力分布及流动分离情况。几何形状及边界条件均模拟相应的脉动流实验模型。采用测得的随时间变化的速度分布作为入口端条件,并利用罚函数和逆风格式等计算技巧得出了光滑的与实验基本相符的速度、压力波形。本文讨论了不同狭窄下速度、压力、壁面剪应力的分布形态,给出了脉动流中狭窄处局部流动分离的间歇性变化规律,并结合实验与临床应用进行了讨论。     

旋转机械系统多自由度非线性动力学数值分析

首先将转子系统的动力响应问题归结为2n维未知向量v的一阶非线性动力学方程dv／dt=Hv＋f（v,t）,并给出了求解这一方程的一次近似式法和三次多项式迫近法。在非稳态、非线性油膜力等作用下,以刚性Jeffcott转子与112个自由度的汽轮发电机组低压转子系统为例,用上述求解方法分析了它们的动力响应及非线性动力学特性;其问,还将计算结果与Runge—Kutta法、Newmark法的相应结果进行了比较,并深入讨论了数值稳定性问题。汽轮发电机组的算例表明对一些具有较复杂的非线性右端项,、同时规模又较大的问题,如果采用四阶Runge—Kutta法,才算几步就因数值骤然增大而失控;但若用同样步长的一次近似式,由于它是一种显式的无条件稳定算式,则计算过程迅速且结果合理可靠。     

三点弯曲型试件Ⅰ/Ⅱ复合型断裂的数值与实验研究

三点弯曲试件常被用来研究材料的纯Ⅰ型及Ⅰ/Ⅱ复合型断裂性能。使用Abaqus软件通过有限元法获得三点弯曲试件不同跨长比2S/L、相对裂纹长度a和裂纹倾角β下的Ⅰ型与Ⅱ型无量纲应力强度因子。由有限元结果可知,在相对裂纹长度一定时,跨长比越小,裂纹倾角越小,Ⅰ型无量纲应力强度因子可取得0,即越容易得到纯Ⅱ型断裂。当跨长比相同时,相对裂纹长度越大越容易得到纯Ⅱ型断裂。使用6组共24个三点弯曲试件,研究砂岩的纯Ⅰ型、纯Ⅱ型以及复合型加载时的断裂韧度,得到该种砂岩的纯Ⅰ型断裂韧度为0.786MPa·m0.5,纯Ⅱ型断裂韧度为0.344MPa·m0.5。实验结果与最大周向应力(MTS)准则和广义最大周向应力(GMTS)准则的预测结果进行对比,比较结果显示GMTS准则更能准确地预测实验结果。     

柱面扁壳的非线性数值分析

本文研究了柱面扁壳在均布荷载作用下的大变形弯曲问题。首先导出了位移型非线性控制微分方程,然后结合两种边界情况(简支与固定)给出了用正交配点法进行解算的详细公式。     

血流动力学数值模拟与动脉粥样硬化研究进展

血流动力学因素被认为与动脉粥样硬化等病理改变密切相关。目前血流动力学数值模拟的对象,主要集中于分支动脉、弯曲动脉以及因血管内膜增生而导致的局部狭窄动脉,这些都是动脉粥样硬化多发的病灶部位。精确的血流动力学数值模拟,必须依赖于解剖精确的血管几何模型和生理真实的血流与管壁有限变形的非线性瞬态流－固耦合。只有在“虚拟血液流动”的基础上,综合考虑血管内的壁面剪应力、粒子滞留时间和氧气的跨血管壁传输等多种因素,血流动力学的数值模拟才能真正有助于人们理解动脉粥样硬化的血流动力学机理,才有可能应用于有关动脉疾病的外科手术规划中。      

钨合金三点弯曲破坏的细观数值模拟

运用有限元方法模拟了钨合金材料的三点弯曲实验;采用单胞单元计算模型,分析了断口形貌与材料力学性能之间的关系;建立了宏观断裂形式与微观结构参数之间的关联.结果表明:裂纹主要产生在 W-W 界面、W-M 界面,且在基体中向前发展,钨合金材料的破坏可归结为延性材料的破坏;裂纹的尖端始终处于应力集中的状态,随着裂纹的扩展,拉应力集中不断释放,并且转移到新的裂纹尖端;裂纹尖端的单元为拉伸破坏,拉应力是导致裂纹扩展的主要因素.     

Numerical modelling of shear-dependent mass transfer in large arteries

A numerical scheme for the simulation of blood flow and transport processes in large arteries is presented. Blood flow is described by the unsteady 3D incompressible Navier–Stokes equations for Newtonian fluids; solute transport is modelled by the advection–diffusion equation. The resistance of the arterial wall to transmural transport is described by a shear-dependent wall permeability model. The finite element formulation of the Navier–Stokes equations is based on an operator-splitting method and implicit time discretization. The streamline upwind/Petrov–Galerkin (SUPG) method is applied for stabilization of the advective terms in the transport equation and in the flow equations. A numerical simulation is carried out for pulsatile mass transport in a 3D arterial bend to demonstrate the influence of arterial flow patterns on wall permeability characteristics and transmural mass transfer. The main result is a substantial wall flux reduction at the inner side of the curved region. © 1997 John Wiley & Sons, Ltd.     

An analysis model of pulsatile blood flow in arteries

IntroductionTheperiodicallypulsatilebloodflowinthearterycausesthecircumferentialandaxialmotionoftheelasticbloodvesselandinturntheoscillationofthevesselaffectsthatofthebloodflow .Womersley[1]resolvedsuccessfullythisfluid_solidcouplingproblembysolvingbothlinearNavier_Stokesequationsandthemotionequationsofthethin_walledelastictubeandgainedtheexpressionsofthebloodflowvelocitiesandthevasculardisplacements.Histheoryhasbeenthebasisforthequantitativeanalysisoftherelationshipofthearterialstructureandi…     

Non-isothermal flow of polymer melts in a curved tube

The results of a numerical study (using finite differences) of heat transfer in polymer melt flow is presented. The rheological behaviour of the melt is described by a temperature-dependent power-law model. The curved tube wall is assumed to be at constant temperature. Convective and viscous dissipation terms are included in the energy equation. Velocity, temperature and viscosity profiles, Nusselt numbers, bulk temperatures, etc. are presented for a variety of flow conditions. Br — Brinkman number - c specific heat, J/kg K - De — Dean number - E dimensionless apparent viscosity, eq. (14d) - G dimensionless shear rate, eq. (19) - k parameter of the power-law model, °C^–1, eq. (7) - mass flow rate, kg/s - m ₀ parameter of the power-law model, Pa · s ⁿ , eq. (7) - n parameter of the power-law model, eq. (7) - Nu 2r _p/ — Nusselt number, eqs. (28,31) - p pressure, Pa - Pe — Péclet number - P —(p/)/r _c — pressure gradient, Pa/m - dissipated energy, W, eq. (29) - total energy, W, eq. (30) - r radial coordinate, m - r _c radius of tube-curvature, m, fig. 1 - r _p radius of tube, m, fig. 1 - r _t variable, m, eq. (6) - R dimensionless radial coordinate, eq. (14a) - R _c dimensionlessr _c, eq. (14a) - R _t dimensionlessr _t, eq. (14a) - t temperature, °C - bulk temperature, °C, eq. (27) - t ₀ inlet temperature of the melt, °C - t _w tube wall temperature, °C - T dimensionless temperature, eq. (14c) - T _w dimensionless tube wall temperature - T dimensionless bulk temperature - u ₁ variable, s^–1, eq. (4) - u ₂ variable, s^–1, eq. (5) - U ₁ dimensionlessu ₁, eq. (18) - U ₂ dimensionlessu ₂, eq. (18) - v velocity in-direction, m/s - average velocity of the melt, m/s - V dimensionlessv, eq. (14b) - dimensionless , eq. (15) - z r _c — centre length of the tube, m - Z dimensionlessz, eq. (14e) - heat transfer coefficient, W/m² K - shear rate, s^–1, eq. (8) - — shear rate, s^–1 - apparent viscosity, Pa · s, eq. (7) - ₀ — apparent viscosity, Pa · s - angular coordinate, rad, fig. 1 - thermal conductivity, W/m K - melt density, kg/m³ - axial coordinate, rad, fig. 1 - rate of strain tensor, s^–1, eq. (8) - (—p) pressure drop, Pa     

棱柱绕流流动的数值模拟

利用数值方法对长宽比为1／3, 1和3的棱柱绕流在雷诺数为100的非稳态流动特性进行了分析和研究。采用有限体积法对棱柱绕流的二维流动N-S方程进行离散求解,分析和研究了非稳态的棱柱绕流流场,升力系数,阻力系数和涡动特性,数值模拟的结果与相关文献的数据比较吻合。通过上述研究能够为了解棱柱绕流的非稳态流动特性提供有力的帮助。而对棱柱三维流动的模拟分析和对雷诺数的变化对棱柱流动特性的影响进行研究,将为掌握棱柱绕流的工程特性打下基础。     

Bifurcation of flow and mass transport in a curved blood vessel

A numerical analysis of flow and concentration fields of macromolecules in a, slightly curved blood vessel was carried out. Based on these results, the effect of the bifurcation of a flow on the mass transport in a curved blood vessel was discussed. The macromolecules turned out to be easier to deposit in the inner part of the curved blood vessel near the critical Dean number. Once the Dean number is higher than the critical number, the bifurcation of the flow appears. This bifurcation can prevent macromolecules from concentrating in the inner part of the curved blood vessel. This result is helpful for understanding the possible correlations between the blood dynamics and atherosclerosis. The project supported by National Natural Science Foundation of China (10002003), JSPS Postdoctoral Fellowship for Foreign Researcher and Foundation for University Teachers, the Ministry of Education     

曲线桥分析的精细传递矩阵法

将精细积分与传递矩阵法相结合,提出一种新的精细传递矩阵格式,应用于曲线桥的分析中。与传统的传递矩阵法相比,无需对微分方程组进行求解,只需迭代即可得到所需要的传递矩阵。根据边界条件,得到结构的内力及变形。算例表明,该方法正确有效。     

A numerical study of an unsteady laminar flow in a doubly constricted 3D vessel

Unsteady flow dynamics in doubly constricted 3D vessels have been investigated under pulsatile flow conditions for a full cycle of period T. The coupled non‐linear partial differential equations governing the mass and momentum of a viscous incompressible fluid has been numerically analyzed by a time accurate Finite Volume Scheme in an implicit Euler time marching setting. Roe's flux difference splitting of non‐linear terms and the pseudo‐compressibility technique employed in the current numerical scheme makes it robust both in space and time. Computational experiments are carried out to assess the influence of Reynolds' number and the spacing between two mild constrictions on the pressure drop across the constrictions. The study reveals that the pressure drop across a series of mild constrictions can get physiologically critical and is also found to be sensitive both to the spacing between the constrictions and the oscillatory nature of the inflow profile. The flow separation zone on the downstream constriction is seen to detach from the diverging wall of the constriction leading to vortex shedding with 3D features earlier than that on the wall in the spacing between the two constrictions. Copyright © 2002 John Wiley & Sons, Ltd.     

Geometric nonlinear formulation for curved beams with varying curvature

Based on exact Green strain of spatial curved beam, the nonlinear strain-displacement relation for plane curved beam with varying curvature is derived. Instead of using the previous straight beam elements, curved beam elements are used to approximate the curved beam with varying curvature. Based on virtual work principle, rigid-flexible coupling dynamic equations are obtained. Physical experiments were carried out to capture the large overall motion and the strain of curved beam to verify the present rigid-flexible coupling formulation for curved beam based on curved beam element. Numerical results obtained from simulations were compared with those results from the physical experiments. In order to illustrate the effectiveness of the curved beam element methodology, the simulation results of present curved beam elements are compared with those obtained by previous straight beam elements. The dynamic behavior of a slider-crank mechanism with an initially curved elastic connecting rod is investigated. The advantage of employing generalized-α method is pointed out and the special nonlinear dynamic characteristics of the curved beam are concluded.     

BLOOD FLOW AND MACROMOLECULAR TRANSPORT IN CURVED BLOOD VESSELS

A numerical analysis of the steady/pulsatile flow and macromolecular (such as LDL and Albumin) transport in curved blood vessels was carried out. The computational results predict that the vortex of the secondary flow is time-dependent in the aortic arch. The concentration of macromolecule concentrates at the region of sharp curve, and the wall concentration at the outer part is higher than that at the inner part. Atherosclerosis and thrombosis are prone to develop in such regions with sharp flow.     

Numerical analysis of critical jet penetration

In this paper, a computational program for two dimensional unsteady compressible, elastic-plastic flow is used to study low speed jet penetration into a steel target. Based on the computational results, the dynamic process of penetration can be distinguished into two stages, cratering and steady penetrating. This is in agreement with the jet penetration mechanism developed by Cheng Che Min in 1972^[1,2]. The consumption of the jet energy is discussed and the value of the critical jet velocity, which represents the strength resistance of the target, is estimated.