The effects of age on viscoelasticity and deformation recovery of articular chondrocytes

本文旨在表征年龄对软骨细胞在自然生长过程中的黏弹性和恢复变形能力的影响. 结果表明：年龄对软骨细胞黏弹性及其恢复变形能力产生显著影响，老年组软骨细胞各 项黏弹性参数值均明显高于幼年和中年组软骨细胞(p<0.0001), 而后两组无显著差异 (p>0.05); 老年组软骨细胞蠕变达到平衡态所需时间t_E显著小于幼年和中年组 (p<0.05), 而幼年和中年组无显著差异(p>0.05). 老年组软骨细胞最大蠕变位移 L_M显著大于幼年和中年组(p<0.005), 而幼年和中年组无显著差异(p>0.05). 老年组软骨细胞恢复变形时间t_R显著大于幼年和中年组(p<0.005), 而幼年和中年组 无显著差异(p>0.05). 恢复变形前8s的分析发现，幼年组软骨细胞恢复变形率K_y 显著高于中年和老年组(p<0.005), 而中年组(K_a)和老年组(K_o)软骨细胞 的恢复变形率无显著差异(p>0.05); 此外, 实验发现老年组软骨细胞的残余变形L_R 比幼年和中年组显著增大(p<0.005), 而后两组无显著差异(p>0.05). 研究工作对 于软骨组织工程、软骨细胞与支架材料相互作用以及探讨OA发生过程中的力学生物学机制具 有理论意义.     

软骨细胞黏弹性及恢复变形与年龄相关性研究

本文旨在表征年龄对软骨细胞在自然生长过程中的黏弹性和恢复变形能力的影响.结果表明:年龄对软骨细胞黏弹性及其恢复变形能力产生显著影响,老年组软骨细胞各项黏弹性参数值均明显高于幼年和中年组软骨细胞(p<0.000 1),而后两组无显著差异(p>0.05);老年组软骨细胞蠕变达到平衡态所需时间t_E显著小于幼年和中年组(p<0.05),而幼年和中年组无显著差异(p>0.05).老年组软骨细胞最大蠕变位移L_M显著大于幼年和中年组(p<0.005),而幼年和中年组无显著差异(p>0.05).老年组软骨细胞恢复变形时间t_R显著大于幼年和中年组(p<0.005),而幼年和中年组无显著差异(p>0.05).恢复变形前8s的分析发现,幼年组软骨细胞恢复变形率K_y显著高于中年和老年组(P<0.005),而中年组(K_α)和老年组(K_o)软骨细胞的恢复变形率无显著差异(p>0.05);此外,实验发现老年组软骨细胞的残余变形L_R比幼年和中年组显著增大(p<0.005),而后两组无显著差异(p>0.05).研究工作对于软骨组织工程、软骨细胞与支架材料相互作用以及探讨OA发生过程中的力学生物学机制具有理论意义.     

心肌细胞大变形黏弹性模型及在实验中的应用

在衡量单个细胞力学行为的研究中,越来越多地采用结合实验的数值模拟方法. 在连续介质力学框架下,发展了一种新的心肌细胞本构模型,并与微管吮吸实验结合,探讨了心肌细胞的力学特性. 本构模型是对普遍使用的仅能用于小变形分析的标准线性固体模型的一种扩展,它将超弹性性能引入到黏弹性模型中,用以描述细胞的大变形黏弹性效应. 基于改进的本构模型,对心肌细胞微管吮吸实验过程进行了有限元模拟,并将计算结果与实验结果以及经典理论解进行了对比. 结果显示发展的本构模型适合细胞大变形问题的有限元数值模拟.      

深埋圆柱形隧洞黏弹性衬砌-土系统稳态振动的解析解

衬砌和土体具有黏弹性性质.将土骨架和衬砌结构视为具有分数阶导数本构的黏弹性体,在频率域内研究了深埋圆柱形隧洞衬砌和土体系统的动力特性.基于黏弹性理论,根据界面连续性条件,分别得到了黏弹性土体和衬砌结构的径向位移、应力等的解析表达式.在此基础上,对比分析了经典弹性土和弹性衬砌系统、分数导数黏弹性衬砌和土体系统的动力特性.考察了土体和衬砌的模量比、衬砌厚度、分数导数阶数、材料参数比对系统动力响应的影响.结果表明:经典弹性土和弹性衬砌系统与分数导数黏弹性衬砌和土体系统的动力特性存在较大差异.随着分数导数阶数的增加,衬砌的径向位移和环向应力幅值明显减小;土体的黏性对系统动力特性的影响大于衬砌黏性的影响.     

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

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

黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.      

热环境中黏弹性功能梯度材料及其结构的蠕变

基于经典的对应原理, 将 Mori-Tanaka 方法等细观力学结果推广于定常温度环境下的黏弹性情形. 根据泊松比与时间呈弱相关的特点, 给出 Laplace 象空间中功能梯度材料的松弛模量和热膨胀系数, 并直接建立耦合热应变的多维黏弹性本构关系. 在此基础上, 求解黏弹性功能梯度圆柱薄壳在热环境中的轴对称弯曲蠕变变形问题. 考虑材料热物参数的温度相关性, 首先确定稳态温度场, 导出相空间中轴对称弯曲变形的解析解, 采用数值反演得到蠕变变形. 算例表明, 蠕变初期, 热环境的影响明显, 随着时间增加, 热应力松弛, 影响逐渐消失. 当圆柱薄壳受轴压时, 相比于两端固支, 两端简支的端部变形更加明显. 通过圆柱薄壳的轴对称弯曲求解, 给出体积含量呈任意分布的黏弹性功能梯度结构在热机载荷下的蠕变分析途径.      

黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.     

黏弹性土层中单桩纵向振动的对比分析

在频率域内研究了黏弹性土层中端承桩纵向振动的动力特性.将土骨架视为具有分数阶导数本构关系的黏弹性体,基于黏弹性理论,采用平面应变模型给出了分数阶导数黏弹性土层的动力阻抗.考虑桩纵向振动时的横向惯性效应,将桩等效为Rayleigh-Love杆,得到了桩头动力复刚度和导纳的解析表达式.通过数值计算,分析了不同模型土条件下桩头动刚度因子和阻尼随激励频率的动力响应.同时,研究了Rayleigh-Love和Euler-Bernoulli两种模型桩动力特性的差异.分析了桩-土界面连续性模型和相对滑移模型对黏弹性土层中桩纵向振动的影响.结果表明:1随着阶数和材料参数比的增加,桩头刚度因子和阻尼明显减小;2对于大直径桩,随着外荷载激励频率的增加,桩横向效应对刚度因子和阻尼有显著影响.3连续性模型条件下桩头的刚度因子和阻尼在共振时的振幅小于相对滑移模型条件.     

层状分数阶黏弹性饱和地基与梁共同作用的时效研究

饱和地基与梁共同作用问题的研究在力学领域及工程界都具有重要意义. 采用分数阶Merchant模型研究饱和地基的流变固结, 该模型比常用整数阶黏弹性模型更能精确反映地基的时变特征. 基于层状正交各向异性黏弹性饱和地基的固结解答, 采用有限元法与边界元法耦合的方法, 研究梁与分数阶黏弹性饱和地基的共同作用问题. 依据Timoshenko梁理论将梁离散为若干单元, 进而得到梁的总刚度矩阵方程; 将黏弹性地基固结问题的精细积分解答作为边界积分的核函数, 采用边界元法建立地基柔度矩阵方程; 结合梁与地基接触面的位移协调条件以及力的平衡条件, 通过有限元法与边界元法的耦合, 最终求得层状分数阶黏弹性饱和地基与Timoshenko梁共同作用的解答. 将本文地基退化为Kelvin地基进行计算, 并与已有文献中的算例进行对比, 二者具有很好的一致性. 在此基础上, 探讨分数阶次和地基成层性对梁与黏弹性饱和地基共同作用的影响. 结果表明: 分数阶次高的黏弹性饱和地基的固结速率明显更快; 对于层状地基, 加固表层土体能有效控制地基整体沉降, 并减小差异沉降. 实际工程中, 应充分考虑饱和地基流变及土体分层性的影响, 以准确分析梁与地基的共同作用过程.      

Stabilization in a two-dimensional chemotaxis-Navier–Stokes system

This paper deals with an initial-boundary value problem for the system $$\left\{ \begin{array}{llll} n_t + u\cdot\nabla n &=& \Delta n -\nabla \cdot (n\chi(c)\nabla c), \quad\quad & x\in\Omega, \, t > 0,\\ c_t + u\cdot\nabla c &=& \Delta c-nf(c), \quad\quad & x\in\Omega, \, t > 0,\\ u_t + \kappa (u\cdot \nabla) u &=& \Delta u + \nabla P + n \nabla\phi, \qquad & x\in\Omega, \, t > 0,\\ \nabla \cdot u &=& 0, \qquad & x\in\Omega, \, t > 0,\end{array} \right.$$ which has been proposed as a model for the spatio-temporal evolution of populations of swimming aerobic bacteria. It is known that in bounded convex domains ${\Omega \subset \mathbb{R}^2}$ and under appropriate assumptions on the parameter functions χ, f and ?, for each ${\kappa\in\mathbb{R}}$ and all sufficiently smooth initial data this problem possesses a unique global-in-time classical solution. The present work asserts that this solution stabilizes to the spatially uniform equilibrium ${(\overline{n_0},0,0)}$ , where ${\overline{n_0}:=\frac{1}{|\Omega|} \int_\Omega n(x,0)\,{\rm d}x}$ , in the sense that as t→∞, $$n(\cdot,t) \to \overline{n_0}, \qquad c(\cdot,t) \to 0 \qquad \text{and}\qquad u(\cdot,t) \to 0$$ hold with respect to the norm in ${L^\infty(\Omega)}$ .     

Incompressible 3D Navier–Stokes Equations as a Limit of a Nonlinear Parabolic System

In this paper, we consider the Cauchy problem for a nonlinear parabolic system ${u^\epsilon_t - \Delta u^\epsilon + u^\epsilon \cdot \nabla u^\epsilon + \frac{1}{2}u^\epsilon\, {\rm div}\, u^\epsilon - \frac{1}{\epsilon}\nabla\, {\rm div}\, u^\epsilon = 0}$ in ${\mathbb {R}^3 \times (0,\infty)}$ with initial data in Lebesgue spaces ${L^2(\mathbb {R}^3)}$ or ${L^3(\mathbb {R}^3)}$ . We analyze the convergence of its solutions to a solution of the incompressible Navier?CStokes system as ${\epsilon \to 0}$ .     

Global Compensated Compactness Theorem for General Differential Operators of First Order

Let A ₁(x, D) and A ₂(x, D) be differential operators of the first order acting on l-vector functions ${u= (u_1, \ldots, u_l)}$ in a bounded domain ${\Omega \subset \mathbb{R}^{n}}$ with the smooth boundary ${\partial\Omega}$ . We assume that the H ¹-norm ${\|u\|_{H^{1}(\Omega)}}$ is equivalent to ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_1u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ and ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_2u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ , where B _i  = B _i (x, ν) is the trace operator onto ${\partial\Omega}$ associated with A _i (x, D) for i = 1, 2 which is determined by the Stokes integral formula (ν: unit outer normal to ${\partial\Omega}$ ). Furthermore, we impose on A ₁ and A ₂ a cancellation property such as ${A_1A_2^{\prime}=0}$ and ${A_2A_1^{\prime}=0}$ , where ${A^{\prime}_i}$ is the formal adjoint differential operator of A _i (i = 1, 2). Suppose that ${\{u_m\}_{m=1}^{\infty}}$ and ${\{v_m\}_{m=1}^{\infty}}$ converge to u and v weakly in ${L^2(\Omega)}$ , respectively. Assume also that ${\{A_{1}u_m\}_{m=1}^{\infty}}$ and ${\{A_{2}v_{m}\}_{m=1}^{\infty}}$ are bounded in ${L^{2}(\Omega)}$ . If either ${\{B_{1}u_m\}_{m=1}^{\infty}}$ or ${\{B_{2}v_m\}_{m=1}^{\infty}}$ is bounded in ${H^{\frac{1}{2}}(\partial\Omega)}$ , then it holds that ${\int_{\Omega}u_m\cdot v_m \,{\rm d}x \to \int_{\Omega}u\cdot v \,{\rm d}x}$ . We also discuss a corresponding result on compact Riemannian manifolds with boundary.     

Two-Phase Inertial Flow in Homogeneous Porous Media: A Theoretical Derivation of a Macroscopic Model

The purpose of this article is to derive a macroscopic model for a certain class of inertial two-phase, incompressible, Newtonian fluid flow through homogenous porous media. Starting from the continuity and Navier–Stokes equations in each phase β and γ, the method of volume averaging is employed subjected to constraints that are explicitly provided to obtain the macroscopic mass and momentum balance equations. These constraints are on the length- and time-scales, as well as, on some quantities involving capillary, Weber and Reynolds numbers that define the class of two-phase flow under consideration. The resulting macroscopic momentum equation relates the phase-averaged pressure gradient to the filtration or Darcy velocity in a coupled nonlinear form explicitly given by

On the viscoelastic characterization of the Jeffreys–Lomnitz law of creep

In 1958, Jeffreys (Geophys J?R Astron Soc 1:92–95) proposed a power law of creep, generalizing the logarithmic law earlier introduced by Lomnitz, to broaden the geophysical applications to fluid-like materials including igneous rocks. This generalized law, however, can be applied also to solid-like viscoelastic materials. We revisit the Jeffreys–Lomnitz law of creep by allowing its power law exponent α, usually limited to the range 0?≤?α?≤?1 to all negative values. This is consistent with the linear theory of viscoelasticity because the creep function still remains a Bernstein function, that is positive with a completely monotone derivative, with a related spectrum of retardation times. The complete range α?≤?1 yields a continuous transition from a Hooke elastic solid with no creep $\left(\alpha \,\to\, -\infty\right)$ to a Maxwell fluid with linear creep $\left(\alpha \,=\,1\right)$ passing through the Lomnitz viscoelastic body with logarithmic creep $\left(\alpha\, =0\right)$ , which separates solid-like from fluid-like behaviors. Furthermore, we numerically compute the relaxation modulus and provide the analytical expression of the spectrum of retardation times corresponding to the Jeffreys–Lomnitz creep law extended to all α?≤?1.     

Determination of axial forces during the capillary breakup of liquid filaments – the tilted CaBER method

The capillary breakup extensional rheometry (CaBER) is a versatile method to characterize the elongational behavior of low-viscosity fluids. Commonly, data evaluation is based on the assumption of zero normal stress in axial direction ( $\upsigma_{\rm zz}=0$ ). In this paper, we present a simple method to determine the axial force using a CaBER device rotated by 90° and analyzing the deflection of the filament due to gravity. Forces in the range of 0.1–1,000?μN could be assessed. Our study includes experimental investigations of Newtonian fructose solutions and silicon oil mixtures (viscosity range, 0.9–60?Pa s) and weakly viscoelastic polyethylene oxide (PEO, $M_{\rm w}=10^{6}$ ?g/mol) solutions covering a concentration range from c?≈?c* (critical overlap concentration) up to c?>?c _e (entanglement concentration). Papageorgiou’s solution for the stress ratio $\upsigma_{\rm zz}/\upsigma_{\rm rr}$ in Newtonian fluids during capillary thinning is experimentally confirmed, but the widely accepted assumption of vanishing axial stress in weakly viscoelastic fluids is not fulfilled for PEO solutions, if c _e is exceeded.     

含金属微粉泡沫材料中爆炸波衰减规律实验研究

设计并制备了一种新型的聚氨酯泡沫材料,研究了爆炸波在该材料中的衰减规律.材料的主要设计原理是在聚氨酯材料中均匀混合了10μm量级的金属微粉,以提高骨架的吸热能力,从而提高其抗爆性能.搭建了实验平台并分别测量了爆炸波在相同孔隙率下含金属微粉和不含金属微粉的聚氨酯材料,同时测量了不同金属微粉含量的聚氨酯材料中的传播特性.实验结果表明:含金属微粉的聚氨酯材料具有更好的抗爆性能.     

Regularity and Asymptotic Behavior of Nonlinear Stefan Problems

We study the following nonlinear Stefan problem $$\left\{\begin{aligned}\!\!&u_t\,-\,d\Delta u = g(u) & &\quad{\rm for}\,x\,\in\,\Omega(t), t > 0, \\ & u = 0 \, {\rm and} u_t = \mu|\nabla_{x} u|^{2} &&\quad {\rm for}\,x\,\in\,\Gamma(t), t > 0, \\ &u(0, x) = u_{0}(x) &&\quad {\rm for}\,x\,\in\,\Omega_0,\end{aligned} \right.$$ where ${\Omega(t) \subset \mathbb{R}^{n}}$ ( ${n \geqq 2}$ ) is bounded by the free boundary ${\Gamma(t)}$ , with ${\Omega(0) = \Omega_0}$ , μ and d are given positive constants. The initial function u ₀ is positive in ${\Omega_0}$ and vanishes on ${\partial \Omega_0}$ . The class of nonlinear functions g(u) includes the standard monostable, bistable and combustion type nonlinearities. We show that the free boundary ${\Gamma(t)}$ is smooth outside the closed convex hull of ${\Omega_0}$ , and as ${t \to \infty}$ , either ${\Omega(t)}$ expands to the entire ${\mathbb{R}^n}$ , or it stays bounded. Moreover, in the former case, ${\Gamma(t)}$ converges to the unit sphere when normalized, and in the latter case, ${u \to 0}$ uniformly. When ${g(u) = au - bu^2}$ , we further prove that in the case ${\Omega(t)}$ expands to ${{\mathbb R}^n}$ , ${u \to a/b}$ as ${t \to \infty}$ , and the spreading speed of the free boundary converges to a positive constant; moreover, there exists ${\mu^* \geqq 0}$ such that ${\Omega(t)}$ expands to ${{\mathbb{R}}^n}$ exactly when ${\mu > \mu^*}$ .     

Gradient Bounds for Elliptic Problems Singular at the Boundary

Let Ω be a bounded smooth domain in \({{R}^N, N \geqq 2}\), and let us denote by d(x) the distance function d(x, ?Ω). We study a class of singular Hamilton–Jacobi equations, arising from stochastic control problems, whose simplest model is

$ - \alpha \Delta u+ u + \frac{\nabla u \cdot B (x)}{d (x)}+ c(x) |\nabla u|^2=f (x) \quad {\rm in}\,\Omega, $

where f belongs to \({W^{1,\infty}_{\rm loc} (\Omega)}\) and is (possibly) singular at \({\partial \Omega, c\in W^{1,\infty} (\Omega)}\) (with no sign condition) and the field \({B\in W^{1,\infty} (\Omega)^N}\) has an outward direction and satisfies \({B\cdot \nu\geqq \alpha}\) at ?Ω (ν is the outward normal). Despite the singularity in the equation, we prove gradient bounds up to the boundary and the existence of a (globally) Lipschitz solution. We show that in some cases this is the unique bounded solution. We also discuss the stability of such estimates with respect to α, as α vanishes, obtaining Lipschitz solutions for first order problems with similar features. The main tool is a refined weighted version of the classical Bernstein method to get gradient bounds; the key role is played here by the orthogonal transport component of the Hamiltonian.