T型微通道内的幂律流体液滴破裂行为的格子Boltzmann方法模拟

采用格子Boltzmann方法(LBM)研究了T型微通道内幂律流体液滴运动行为及其流型相图.主要研究了液滴幂律指数n对液滴破裂时颈部厚度、前端运动距离等形变特性以及流型相图的影响.数值结果表明,幂律流体液滴在T型微通道内存在阻塞破裂、隧道破裂以及不破裂三种流型.在阻塞破裂过程中,液滴颈部厚度随时间逐渐减小,且液滴幂律指数n越大,液滴颈部厚度随时间减小得越慢.同时液滴前端运动距离随时间线性增加,且随着n的增加,液滴破裂时前端运动距离越长.在隧道破裂过程中,液滴颈部厚度也随时间逐渐减小,与阻塞破裂相似,n越大液滴颈部厚度减小得越慢.与阻塞破裂相比,液滴隧道破裂时对应的临界颈部厚度有所增加,且液滴前端运动距离随时间先快速增加,然后再缓慢增加,隧道宽度随时间近似呈对数增长.此外,液滴未破裂时液滴颈部厚度以及液滴前端运动距离出现波动现象.液滴的幂律指数n越大,液滴越容易破裂,但越不容易达到阻塞破裂.根据数值模拟结果得到了各流型相图之间幂函数形式临界分界线的拟合公式,该拟合公式可以预测不同流型.     

T型微通道内液滴在幂律流体中运动机理的格子Boltzmann方法研究北大核心CSCD

采用非Newton不可压两相流格子Boltzmann模型研究了T型微通道内Newton液滴在非Newton幂律流体中的运动过程.研究了非Newton流体幂律指数n、主管道毛细数Ca、两相流量比Q、两相黏度比M以及主管道壁面润湿性θ对液滴在T型微通道内的形成尺寸、形成时间和变形参数(DI)的影响.研究结果表明:首先,主管道流体幂律指数n从0.4增加到1.6时,液滴的形成尺寸近似呈线性减小,而液滴的形成时间和变形参数先快速减小,然后缓慢减小;其次,黏度比对液滴形成尺寸、液滴形成以及变形参数的影响与幂律指数的影响基本一致;再者,随着Ca和主管道壁面润湿性的增加,形成液滴的尺寸近似呈线性减小,形成液滴的时间和变形参数先快速减小然后缓慢减小,且减小趋势随幂律指数的增加而减缓;最后,研究结果还表明主管道和子管道的流量比Q越大,液滴形成时间越长,液滴形成尺寸和变形参数越小.     

任意梯度分布功能梯度圆环的热弹性分析

对沿径向任意变化的材料参数的功能梯度圆环进行了热弹性分析.与以前关于该问题的分析不同,既不需要预先给定具体的梯度变化形式,也不需要对结构进行细分.给出一种新的有效解法将问题转换为求解Fredholm积分方程,从而通过Fredholm积分方程的解给出热应力和位移的分布情况.最后通过算例分析了内外表面受不同温度作用时,材料参数呈现梯度变化对圆环的应力和位移变化的影响,计算结果表明某些特定的材料梯度可有效缓解圆环内的热应力分布.该文得到的结果对功能梯度圆环在结构安全设计方面有重要的理论指导意义.     

含有两个覆盖层的功能梯度磁电弹半空间中表面波的波速方程

研究了含有两个均匀弹性覆盖层的半无限大功能梯度磁电弹材料中的表面波.假设基底为材料性质沿厚度方向指数变化的磁电弹材料,两个覆盖层为不同的均匀弹性材料.通过考察表面及界面的条件,分别给出电磁学开路和短路时,含有两个弹性覆盖层的功能梯度磁电弹半空间中表面波的波速方程.     

功能梯度圆板的轴对称自由振动

基于三维弹性理论,用直接位移法求解了横观各向同性功能梯度圆板的轴对称自由振动,其材料特性沿板厚度方向按指数形式变化,在弹性简支和刚性滑动两种边界条件下得到了三维精确解,即它逐点满足基本方程和边界条件,最后给出了数值算例并与以前的工作进行了对比.该方法也可推广应用于材料特性沿板厚度任意变化情形.     

初应力下含有两个覆盖层的功能梯度磁电弹半空间中表面波研究

假设基底为带有梯度初应力的功能梯度磁电弹材料,材料的性质和初应力都沿厚度方向指数变化.基底上含有两个不同材料的均匀覆盖层,一个是带有常数初应力的磁电弹覆盖层,一个是均匀弹性覆盖层.通过考察不同的电磁学边界条件及界面的连续条件,给出此结构的表面波波速方程.     

辛求解体系下功能梯度材料平面梁的完整解析解

采用辛弹性力学解法,求取弹性模量沿轴向指数变化,而Poisson比保持不变的功能梯度材料平面梁的完整解析解．通过求解被Saint-Venant原理覆盖的一般本征解,建立起完整的解析分析过程,进而给出平面梁位移和应力的精确分布规律．传统的弹性力学分析方法常常忽略被Saint-Venant原理覆盖的解,但这些衰减的本征解对材料的局部效应起着较大的影响作用,可能导致材料或结构的突然失效．采用辛求解方法,充分利用本征向量之间的辛共轭正交关系,得到了功能梯度材料梁的完整解析解．两个数值算例分别将功能梯度材料平面梁的位移和应力分布与相应均匀材料情形的结果进行比较,研究了材料非均匀性对位移和应力解的影响．     

Elastic Analysis of Anisotropic Functionally Graded Rotating Disks With Non-Uniform Thicknesses北大核心CSCD

研究了任意梯度变化的变厚度各向异性转动圆盘的弹性问题.假设圆盘绕刚性轴匀速转动,其材料性能和厚度沿径向任意梯度变化.考虑圆盘在中心转轴处受位移约束,外侧自由,根据各向异性转动圆盘的平衡微分方程,得到关于径向应力的Fredholm积分方程,继而通过对Fredholm积分方程进行数值求解,得到结构的位移场和应力场.对具体梯度变化情况仅需代入相应梯度变化进行求解即可.数值算例部分,通过假设厚度、弹性模量等参数为特殊的幂函数形式,将由Fredholm积分方程求出的数值解与对应的精确解进行对比,以及针对常见的Voigt模型,将由该方法算得的数值解和ANSYS有限元计算结果进行对比,验证了该方法的准确性和精度.其次,针对Voigt模型,重点分析了厚度变化、材料性能梯度参数、各向异性度等对应力场和位移场的影响.提出了针对材料性能和厚度沿径向呈任意梯度变化的圆盘结构弹性分析方法,将为优化功能梯度圆盘的结构和材料参数、有效调整构件应力分布、提高结构安全性,提供强有力的工具;算例分析结果对功能梯度圆盘在复杂条件下的结构安全设计有重要的理论指导意义.     

基于应变梯度理论的粘塑性厚壁圆筒和球壳极限内压分析

基于应变梯度塑性理论,分析了内压作用下厚壁圆筒和球壳的塑性极限荷载．结果表明:圆筒内径在微米量级时,存在尺度效应现象,内径减小,其尺度效应增强;变形越大,影响越大;应变速率敏感指数越大,尺度效应越明显．经典塑性理论结果是当前解的特例．     

弹性边界约束旋转功能梯度圆柱壳结构自由振动行波特性分析

对旋转功能梯度圆柱壳自由振动行波特性及边界约束影响进行了分析研究.将功能梯度材料的物理特性表示成沿壳体厚度方向指数变化的函数,基于Love壳体理论,将圆柱壳3个方向的振动位移场采用改进Fourier(傅立叶)级数方法展开, 进而改善位移函数在边界位置求导连续性,结合旋转圆柱壳结构能量原理描述与Rayleigh Ritz法,推导旋转功能梯度圆柱壳自由振动特征方程.通过将计算结果与现有文献结果对比验证了该文模型的正确性与收敛性.随后,通过算例讨论分析了功能梯度材料特性参数、几何参数、边界条件及约束弹簧刚度对旋转功能梯度圆柱壳自由振动行波振动特性的影响.结果表明:边界条件在环向波数n较小或长径比L/R较小的情况下对行波特性影响较为明显;随着厚径比H/R的增大,边界条件的影响逐渐减小;边界约束弹簧对行波特性影响程度取决于模态阶数情况;功能梯度材料特性参数对前后行波频率的影响随着模态序数的增大而逐渐增大.     

Lattice Boltzmann Study on the Motion of Dual Droplets in Microchannels With Contact Angle Hysteresis北大核心CSCD

接触角滞后表现为流体在非理想固体表面上运动时前进接触角和后退接触角不同,是两相流体在润湿表面上流动的重要现象.该文采用改进的伪势格子Boltzmann(LB)多组分模型,并与几何润湿边界条件相结合,研究了两个液滴在具有接触角滞后性微通道表面上的运动行为,主要研究了通道内特征数、通道表面性质以及液滴初始参数的影响.研究结果表明：毛细数的增大有助于液滴的移动,然而并不利于液滴的排出,且毛细数的增加对上游液滴的影响大于其对下游液滴的影响;另一方面,接触角滞后性窗口越大,液滴运动和形变更迟缓,但形变程度更明显,两液滴更早地发生合并,但更晚地排出管道;液滴间距的增加使液滴的运动行为在不同阶段表现为不同的模式,但都导致通道中残留小液滴,使得液滴排出通道的时间增加.研究结果还表明：上游液滴和下游液滴的相对尺寸差距越大,越不利于液滴排出管道.     

液气界面张力垂直分量引起的基底弹性变形

Young方程是毛细理论和润湿的重要方程之一．但是,该方程只描述了3个界面张力的水平分量之间的平衡与接触角的关系,而对液气界面张力垂直分量未作任何描述．现在,随着软材料的广泛应用,该垂直分量将引起基底的表面变形,并在微流体系统的制造过程中起到重要作用,这已是该研究领域的共识．综述了关于表面变形这一问题在理论分析,实验研究和数值模拟等方面取得的进展．而且,还讨论了由垂直分量引起的表面变形对液滴润湿和铺展行为、微悬臂梁的弯曲、弹性毛细现象、电弹性毛细现象等的影响．不仅对该问题的历史发展和目前的研究进展进行了简单的综述,并且也针对后续的研究提出了几点建议．     

3D elasticity analytical solution for bending of FG micro/nanoplates resting on elastic foundation using modified couple stress theory

This paper addresses a 3D elasticity analytical solution for static deformation of a simply-supported rectangular micro/nanoplate made of both homogeneous and functionally graded (FG) material within the framework of modified couple stress theory. The plate is assumed to be resting on a Winkler–Pasternak elastic foundation, and its modulus of elasticity is assumed to vary exponentially along thickness. By expanding displacement components in double Fourier series along in-plane coordinates and imposing relevant boundary conditions, the boundary value problem (BVP) of plate system, including its governing partial differential equations (PDEs) of equilibrium are reduced to BVP consisting only ordinary ones (ODEs). Parametric studies are conducted among displacement and stress components developed in the plate and FG material gradient index, length scale parameter, and foundation stiffnesses. From the numerical results, it is concluded that the out-of-plane shear stresses are not necessarily zero at the top and bottom surfaces of plate. The results of this investigation may serve as a benchmark to verify further bending analyses of either homogeneous or FG micro/nanoplates on elastic foundation.     

Patient-specific model of lung deformation using spatially dependent constitutive parameters

Breathing-induced spatially dependent lung deformation is predicted using patient-specific elastic properties with the contact–impact analysis model. The lung geometry is derived from 4D CT scan data of real patients. The spatially varying Young’s modulus for the patient is obtained from a previous study that used inverse deformation of the lung. The compact–impact analysis is implemented using the finite element method. The predicted lung deformation is compared with the results based on linear elasticity. The results are consistent with physiology, indicating large deformations near the diaphragm and smaller values at remote locations on the lobe. The effect of non-linearity of elastic property is most significant at the remote locations where the diaphragm-induced deformation is significantly attenuated.     

Effect of gravity on subject-specific human lung deformation

A biomechanical model of human lung is developed and used to investigate the effect of gravity on lung deformation. The lung is assumed to behave as a poro-elastic medium with spatially dependent elastic property. Finite element analysis is performed on a three-dimensional (3D) lung geometry reconstructed from a four-dimensional Computed Tomography (4DCT) scan dataset of human patient. The spatially dependent Young’s modulus (YM) values are estimated using inverse analysis from a linear elastic deformation model. The predicted deformation of selected landmarks is monitored with and without gravity, and compared with data obtained from 4DCT registration. The results show that gravity indeed significantly affects the magnitude and distribution of lung deformation with the maximum displacement enhanced by 54% in the direction of gravity, for the conditions investigated. In summary, the accuracy of predicted deformation is improved through incorporation of gravity in the biomechanical model of lung.     

Identification of the elastic modulus of polymeric materials by using thin-walled cylindrical specimens

A method for determining the elastic modulus of polymeric materials from deformation diagrams of thin-walled circular cylindrical shells in compression and tension in the region of geometrical nonlinearity has been elaborated. A numerical solution is found by the finite-element method (ANSYS.) The existence of a unified deformation diagram in generalized coordinates is established, from which the elastic modulus is determined. To validate the method, the eigenfrequencies of cylindrical specimens were found experimentally. The results obtained are compared with FEM calculations.     

Stability of strain gradient elastic bars in tension

Equilibrium of a bar under uniaxial tension is considered as optimization problem of the total potential energy. Uniaxial deformations are considered for a material with linear constitutive law of strain second gradient elasticity. Applying tension on an elastic bar, necking is shown up in high strains. That means the axial strain forms two homogeneously deformed sections in the ends of the bars and a section in the middle with high variable strain. The interactions of the intrinsic (material) lengths with the non linear strain displacement relations develop critical states of bifurcation with continuous Fourier’s spectrum. Critical conditions and post-critical deformations are defined with the help of multiple scales perturbation method. An erratum to this article can be found at     

Buckling and postcritical behaviour of the elastic infinite plate strip resting on linear elastic foundation

In this paper the von Kármán model for thin, elastic, infinite plate strip resting on a linear elastic foundation of Winkler type is studied. The infinite plate strip is simply-supported and subjected to evenly distributed compressive loads. The critical values of bifurcation parameters and buckling modes for given frequency of longitudinal waves are found on the basis of investigation of linearized problem. The mathematical nonlinear model is reduced to operator equation with Fredholm type operator of index 0 depending on parameters defined in corresponding Hölder spaces. The Lyapunov-Schmidt reduction and the Crandall-Rabinowitz bifurcation theorem (gradient case) are used to examine the postcritical behaviour of the plate. It is proved that there exists maximal frequency of longitudinal waves depending on the compressive load and the stiffness modulus of foundation.     

Contact angle effects on microdroplet deformation using CFD

In this paper we use computational fluid dynamics (CFD) to study the effect of contact angle on droplet shape as it moves through a contraction. A new non-dimensional number is proposed in order to predict situations where the deformed droplet will form a slug in the contraction and thus have the opportunity to interact with the channel wall. It is proposed that droplet flow into a contraction is a useful method to ensure that a droplet will wet a channel surface without a trapped lubrication film, and thus help ensure that a slug will remain attached to the wall downstream of the contraction. We demonstrate that when a droplet is larger than a contraction, capillary and Reynolds numbers, and fluid properties may not be sufficient to fully describe the droplet dynamics through a contraction. We show that, with everything else constant, droplet shape and breakup can be controlled simply by changing the wetting properties of the channel wall. CFD simulations with contact angles ranging from 30° to 150° show that lower contact angles can induce droplet breakup while higher contact angles can form slugs with contact angle dependent shape.