细胞趋硬性迁移的研究进展

细胞力学微环境可以调控许多细胞生理功能.特别地,在细胞力学微环境各种信号梯度的作用下,细胞可以定向地迁移.这些定向迁移可以显著影响伤口愈合、癌细胞转移和组织形貌发育等生理过程.目前为止,细胞的定向迁移主要包括:在化学药物梯度作用下的趋药性迁移,在黏附分子梯度作用下的趋触性迁移,以及在细胞外基质硬度梯度作用下的趋硬性迁移等.虽然细胞趋药性和趋触性迁移的力学机理得到了很好的研究,但是关于细胞趋硬性迁移的机制和作用还不清楚.该文重点介绍了细胞趋硬性的相关实验和理论研究进展,分析了不同研究间的联系与区别,讨论了细胞趋硬性迁移的潜在力学机制,提出尚存在的问题和未来可能的研究方向.     

细胞骨架与细胞外基质的力学建模与分析

细胞和生物组织需要适应人体内复杂的力学生物学环境,一方面要被动地承受外部环境中的机械力,另一方面在组织生长、修复等生理过程中要积极主动地产生机械力调整自身的结构和形态.细胞和生物组织的力学性质主要由细胞骨架和细胞外基质决定,它们在微观上都是生物聚合物交联形成的复杂的、各向异性的三维网络结构.这方面早期的力学研究主要集中在通过各种网络模型,理解其普遍存在的非线性响应和硬化行为.近年来随着实验方法、理论建模和计算机模拟技术的大幅进步,这些生命介质的力学性质及其潜在的力学机理得到了更深入的理解.该文回顾了近些年细胞骨架和细胞外基质研究方面取得的部分进展,主要侧重动态交联属性、生物聚合物力学化学耦合赋予的主动材料属性、交联网络塑性和断裂,以及力学训练引发的自适应网络重构.发展细胞骨架与细胞外基质的力学模型与计算方法,分析该类生命介质的复杂力学行为,理解这些力学行为的潜在机制,可以加深我们对细胞和组织的力学生物学认识,并为人造生物材料和细胞组织工程提供基础和参考.     

分子刚度对膜锚定受体-配体结合动力学影响的理论与模拟研究

人体内大部分生物学过程都离不开细胞黏附.细胞黏附行为主要由锚定于细胞膜上的特异性分子（又称受体和配体）的结合动力学关系来决定.已有研究表明，特异性分子的结合关系受外力及细胞膜波动等多种因素影响.然而，特异性分子刚度对细胞膜锚定受体 配体结合关系的影响机制仍不清楚.近期关于新冠病毒强传染力的研究表明，特异性黏附分子刚度对病毒与细胞结合具有重要影响.该文通过建立生物膜黏附的粗粒度模型，借助分子模拟和理论分析来研究分子刚度在黏附中的作用.结果表明，始终存在一个最佳膜间距及最佳分子刚度值，使得黏附分子亲和力和结合动力学参数达到最大值.这项研究不仅能加深人们对细胞黏附的认知，还有助于指导药物设计、疫苗研发等.     

圆形均布荷载作用下分数导数型黏弹性土基变形分析

土体的蠕变特性是影响工后沉降和工程安全的重要因素.基于半空间弹性土基受圆形均布荷载作用弹性理论解,根据弹性与黏弹性理论的对应原理,建立了分数导数型黏弹性土基在竖向圆形均布荷载作用下的地表位移与分数阶导数等参数的关系,并分析了不同分数阶下地表变形的时效特性.结果表明,与经典黏弹性本构模型相比,分数导数黏弹性模型能够在较宽的范围内描述黏弹性土基变形的特性,采用分数导数Kelvin黏弹性本构模型计算的地表沉降较经典的Kelvin黏弹性模型小,土基的蠕变特性与分数导数的阶数有关,具有更为广泛的适用性和应用前景.     

计及时变演化特征的硅泡沫垫层非线性黏弹性模型研究

基于黏弹性基本理论,引入材料非线性特征,考虑了材料加载、保载应力松弛历史、老化效应以及黏弹性模型各运动单元退化的差异性,并从两种老化机制出发,获得了老化硅泡沫垫层力学模型以及长时应力松弛硅泡沫垫层接续加载力学模型.模型机理清晰,能够反映材料服役历史信息及其对力学效应的影响.     

基于弹性壳的三维群体细胞动力学模型

群体细胞迁移常见于胚胎发育、伤口愈合和肿瘤侵袭等各种生理和病理过程中,关于其动力学的研究对于揭示群体细胞迁移机理、深刻理解有关生物过程十分重要.该文构建了群体细胞的三维可变形壳状模型,提出了一种考虑细胞弹性形变和细胞间接触与黏附相互作用的群体细胞动力学理论,并发展了相应的数值算法.基于所发展的动力学模型与算法,对多细胞在嚢腔里的受限旋转运动进行了模拟,复现了相关实验现象,分析了细胞极性、细胞形变、胞间相互作用等因素对多细胞三维动力学的影响规律.     

黏弹性圆柱壳计及切变形和转动惯量时的动力学稳定性

根据修正的Timoshenko理论,在几何非线性中考虑了剪切变形和转动惯量,对黏弹性圆柱壳的动力稳定性进行了研究.根据Bubnov-Galerkin法,结合基于求积公式的数值方法,将问题简化为求解具有松弛奇异核的非线性积分-微分方程的问题.针对物理-力学和几何参数在大范围内的变化,研究壳体的动力特性,显示了材料的黏弹性对圆柱壳动力稳定性的影响.最后,比较了通过不同的理论得到的结果.     

受约束高分子链的拉伸

高分子生物材料微观力学性质的定量刻画,以及先进生物微纳米技术与器件的发展均需要定量分析生物大分子等高分子链在复杂微环境中的统计热力学性质与行为.在实现这一目标的过程中,连续介质力学与统计热力学的交叉研究扮演着很重要的角色.针对这一领域的力学问题,该综述先从DNA分子的受力拉伸出发,通过引入描述高分子链统计热力学性质的几类理论模型,指出了蠕虫链相较其他理想随机链模型在描述半柔性高分子链力与构型变化关系时具有较为显著的优势,从而使得人们对高分子在复杂微环境下,统计热力学性质与行为的定性与定量认识在很大程度上取决于基于蠕虫链模型的相关研究进展.根据这一事实,通过回顾与梳理空间几何约束对高分子链随机构象分布影响、高分子链在拉力与约束同时作用时的统计热力学建模、以及基于高性能计算机的高分子链统计物理性质仿真等各方面研究的现状,系统总结了蠕虫链在不同约束与受力微环境下,其统计热力学性质与行为研究的最新进展和依旧存在的挑战性难题.最后,通过总结分析,指出了蠕虫链在复杂微环境下的统计热力学研究是从分子与细胞尺度理解生命现象、发展先进微纳米技术以及构建软物质本构关系的重要基础,目前已成为极富挑战性的力学交叉科学前沿课题.     

抗血管生成药物Endostantin作用下实体肿瘤血管生成的数值模拟:考虑基质力学环境及血管生成抑素的影响

为研究抗血管生成药物Endostatin作用下,肿瘤血管生成过程中基质力学环境及血管生成抑素的影响,考虑内皮细胞(EC)和细胞外基质(ECM)两相,耦合抗血管生成药物Endostatin和血管生成抑素Angiostatin的抑制效应,建立肿瘤内外血管生成的二维数值模型.抗血管生成因子Angiostatin和药物Endostatin耦合作用时,可明显降低肿瘤组织内的微血管密度,对肿瘤快速生长起到一定的抑制作用.所给出的模型,可以较好模拟基质力学环境影响下,肿瘤抗血管生成因子对内皮细胞迁移和增殖的抑制作用.     

二维黏弹性力学问题的无网格自然单元法

基于无网格自然单元法,建立了求解二维黏弹性力学问题的一条新途径.基于弹性 黏弹性对应原理和Laplace(拉普拉斯)变换技术,首先将黏弹性问题转换成Laplace域内与弹性力学问题相同的形式,然后推导出基于自然单元法分析黏弹性问题的基本公式.作为一种新兴的无网格数值计算方法,自然单元法的实质是一种基于自然邻近插值的Galerkin(伽辽金)法.相对于其他无网格法,自然单元法的形函数具有插值性和支持域各向异性等特点.算例结果证明了所提分析方法的有效性.     

A general stability result for the semilinear viscoelastic wave model under localized effects

Our main goal in the present work is to address an integro-differential model under localized viscoelastic and frictional effects arising in the Boltzmann theory of viscoelasticity. More precisely, we consider a general version in the history context of the pioneer localized viscoelastic problem approached by Cavalcanti and Oquendo (2003) in the null history scenario, and more recently by Cavalcanti et al. (2018) in the history framework. By means of a new observability inequality, we prove a general stability result to the model under a weaker assumption on the localized frictional damping and a slower condition on the decreasing memory kernel (of polynomial type) than the previously mentioned works. To achieve such stability results, we still work in a general setting by removing the assumption on complementary damping mechanisms and show, in some reasonable situations concerning the density coefficient, that the localized viscoelastic effect is enough to ensure the general stability (of polynomial type) to the problem.     

Linear and nonlinear vibrations of fractional viscoelastic Timoshenko nanobeams considering surface energy effects

In this paper, the linear and nonlinear vibrations of fractional viscoelastic Timoshenko nanobeams are studied based on the Gurtin–Murdoch surface stress theory. Firstly, the constitutive equations of fractional viscoelasticity theory are considered, and based on the Gurtin–Murdoch model, stress components on the surface of the nanobeam are incorporated into the axial stress tensor. Afterward, using Hamilton's principle, equations governing the two-dimensional vibrations of fractional viscoelastic nanobeams are derived. Finally, two solution procedures are utilized to describe the time responses of nanobeams. In the first method, which is fully numerical, the generalized differential quadrature and finite difference methods are used to discretize the linear part of the governing equations in spatial and time domains. In the second method, which is semi-analytical, the Galerkin approach is first used to discretize nonlinear partial differential governing equations in the spatial domain, and the obtained set of fractional-order ordinary differential equations are then solved by the predictor–corrector method. The accuracy of the results for the linear and nonlinear vibrations of fractional viscoelastic nanobeams with different boundary conditions is shown. Also, by comparing obtained results for different values of some parameters such as viscoelasticity coefficient, order of fractional derivative and parameters of surface stress model, their effects on the frequency and damping of vibrations of the fractional viscoelastic nanobeams are investigated.     

Nonlinear modeling and analysis of electrically actuated viscoelastic microbeams based on the modified couple stress theory

A nonclassical nonlinear continuum model of electrically actuated viscoelastic microbeams is presented based on the modified couple stress theory to consider the microstructure effect in the framework of viscoelasticity. The nonlinear integral-differential governing equation and related boundary conditions of are derived based on the extended Hamilton's principle and Euler–Bernoulli hypothesis for viscoelastic microbeams with clamped-free, clamped-clamped, simply-supported boundary conditions. The proposed model accounts for system nonlinearities including the axial residual stress, geometric nonlinearity due to midplane stretching, electrical forcing with fringing effect. The behavior of the microbeam is simulated using generalized Maxwell viscoelastic model. A new generalized differential/integral quadrature method is developed to solve the resulting governing equation. The developed model is verified against elastic behavior and a favorable agreement is obtained. Efficiency of the developed model is demonstrated by analyzing the quasistatic pull-in phenomena of electrically actuated viscoelastic microbeams with different boundaries at various material length scale parameters and axial residual stresses in the framework of linear viscoelasticity.     

Constitutive Modelling of Resins in the Stiffness Domain

An analytic method for inverting the constitutive compliance equations of viscoelasticity for resins is developed. These equations describe the HWKK/H rheological model, which makes it possible to simulate, with a good accuracy, short-, medium- and long-term viscoelastic processes in epoxy and polyester resins. These processes are of first-rank reversible isothermal type. The time histories of deviatoric stresses are simulated with three independent strain history functions of fractional and normal exponential types. The stiffness equations are described by two elastic and six viscoelastic constants having a clear physic meaning (three long-term relaxation coefficients and three relaxation times). The time histories of axiatoric stresses are simulated as perfectly elastic.The inversion method utilizes approximate constitutive stiffness equations of viscoelasticity for the HWKK/H model. The constitutive compliance equations for the model are a basis for determining the exact complex shear stiffness, whereas the approximate constitutive stiffness equations are used for determining the approximate complex shear stiffness. The viscoelastic constants in the stiffness domain are derived by equating the exact and approximate complex shear stiffnesses. The viscoelastic constants are obtained for Epidian 53 epoxy and Polimal 109 polyester resins. The accuracy of the approximate constitutive stiffness equations are assessed by comparing the approximate and exact complex shear stiffnesses. The constitutive stiffness equations for the HWKK/H model are presented in uncoupled (shear/bulk) and coupled forms. Formulae for converting the constants of shear viscoelasticity into the constants of coupled viscoelasticity are given as well.     

A boundary element methodology for viscoelastic analysis: Part II without cells

In this study Kelvin and Boltzmann viscoelastic models are implemented in a two-dimensional boundary element atmosphere. This general methodology is based on differential constitutive relations for viscoelasticity, avoiding the use of relaxation functions. In this part of the study, important algebraic operations are introduced into the formulation allowing analysing viscoelastic problems without using internal cells. This improvement is very important to model infinite and semi-infinite regions. The formulation is verified comparing the numerical results with analytical solutions. An extension of the formulation to consider soil–structure interaction is presented in order to demonstrate the vast applicability of the technique.     

DECAY RATES AND CONVERGENCE OF SOLUTIONS TO SYSTEM OF ONE-DIMENSIONAL VISCOELASTIC MODEL WITH DAMPING

This paper considers the asymptotic behavior of solutions to the system of onedimensional viscoelastic model with damping and prove that the corresponding solutions time-asymptotically behave like nonlinear diffusion wave as in [4,11]. In addition, It is also shown that the system of one-dimensional viscoelastic model with damping is a viscosity approximation of a hyperbolic conservation laws with damping.     

DECAY RATES AND CONVERGENCE OF SOLUTIONS TO SYSTEM OP ONE-DIMENSIONAL VISCOELASTIC MODEL WITH DAMPING

This paper considers the asymptotic behavior of solutions to the system of one-dimensional viscoelastic model with damping and prove that the corresponding solutions time-asymptotically behave like nonlinear diffusion wave as in [4,11]. In addition, It is also shown that the system of one-dimensional viscoelastic model with damping is a viscosity approximation of a hyperbolic conservation laws with damping.     

Damping sources interactions in impact of viscoelastic composite plates with damping treated SMA wires,using a hyperbolic plate theory

In the present paper, a comprehensive study is made on effects of the viscoelastic and the phase-transformation-based dissipations and their interactions on impact responses of viscoelastic composite plates with damping treated (structural hierarchy) shape memory alloy (SMA) wires, for the first time. In contrast to almost all of the available researches, a high-order hyperbolic plate theory that includes not only odd but also even functions of the transverse coordinate, is proposed and employed here. While a hierarchical viscoelastic constitutive law is employed for both the orthotropic and SMA materials, Brinson's constitutive law is refined to include the loading fluctuations and structural hierarchy of the SMA wire, simultaneously. The traditional Hertz and Yang-Sun contact laws are modified accordingly. The resulting highly nonlinear piecewise-defined integro-differential finite element governing equations are solved by an iterative algorithm within each time step. The presented discussions show that in contrast to the common belief, the zero-shear traction condition on the top and bottom surfaces of the viscoelastic orthotropic plate cannot be satisfied by the available plate theories, even for the symmetric lamination schemes. Results show that the viscoelasticity and phase-transformation effects on the resulting dynamic responses are more pronounced for the low and high energy impacts, respectively.     

关于分数阶Maxwell淤泥模型及其在水波衰减问题中的应用

分数阶Maxwell模型可用来模拟粘弹性的海床淤泥．与传统的有理分式模型相比,分数阶Maxwell模型可以用更少的自由参数,较好地描述某些真实淤泥的流变特性．将该分数阶Maxwell模型用于研究淤泥与自由表面水波的相互作用,并得到了线形单色波的衰减率．从水波衰减率曲线中可观测到淤泥层的共振现象,共振时衰减率将达到峰值．对于线形单色波,其衰减率还可表示为各模态衰减率之和,从而可研究某一模态的运动对水波衰减的影响．模态分析表明,当某一模态运动引发共振时,总衰减率由该模态的模态衰减率决定．     

Unsteady flow of viscoelastic fluid with the fractional K-BKZ model between two parallel plates

In the present paper, we investigate the unsteady flow of a viscoelastic fluid between two parallel plates which is generated by the impulsively accelerated motion of the bottom plate. Based on the result of (Jaishankar and McKinley, 2014), the fractional K-BKZ constitutive equation is obtained from the fractional Maxwell model. Using respectively the fractional Maxwell model and fractional K-BKZ model, the unidirectional flows between two plates are simulated and compared. The velocity field and shear stress of the flows are calculated by developing efficient finite difference schemes. The results show that the fluid with the fractional Maxwell model gradually loses the viscoelasticity, but the fluid with the fractional K-BKZ model continues to preserve the viscoelasticity. The dependence of the flow velocity on various parameters of the fractional K-BKZ model is analyzed graphically.