曲面物理和力学

回顾了平坦空间中的经典物理和力学, 综述了卷曲空间中物理和力学的三个典型案例, 即生物膜力学、扁壳力学和曲面扩散动力学, 揭示了平坦空间与卷曲空间物理和力学的根本差异. 论文阐明了, 平坦空间的物理和力学只有一个基本微分算子------经典梯度, 而曲面物理和力学必然存在两个基本微分算子------经典梯度和形状梯度.     

?????????????е???

简要综述了生物膜力学与几何的新进展. 在生物膜力学中,着重介绍了基于微分算子的平衡 理论和几何约束理论;在生物膜几何中,重点评述了源于生物膜力学的新梯度算子及其积分 性质. 指出：新梯度算子可能在生物膜曲面上诱发新的驱动力;生物膜力学与几何是一个有 机整体,其背后存在着一个对称的几何体系,包括对称的微分算子以及对称的积分定理系统.     

基于物理信息的神经网络求解曲面上对流扩散方程

引入基于物理信息的神经网络PINNs(Physics-informed Neural Networks)并将其用于求解曲面对流扩散方程。区别于传统的神经网络模型,PINNs在建立模型过程中引入了自动微分技术,并将物理信息即偏微分方程信息编译其中,通过定义损失函数得到关于该模型中神经网络参数即权重和偏置的优化目标,随后利用已有的优化算法进行求解。显而易见,PINNs通过添加额外的物理信息约束放宽了对于数据量的要求,对于一个确定性模型显示出更好的鲁棒性。本文基于曲面微分算子与欧氏空间下标准微分算子的解析关系,引入两种曲面微分算子处理技术,即非本征技术和嵌入技术,并结合PINNs针对定义在高维复杂曲面上的对流扩散方程进行求解,多个数值算例证明了该方法的有效性、鲁棒性以及其在求解此类问题的潜力。     

基于物理信息神经网络的功能梯度材料稳态/瞬态热传导分析

采用物理信息神经网络PINN(Physics-informed Neural Networks)求解稳态和瞬态功能梯度材料(FGMs)热传导问题。该方法利用控制方程、边界及初始条件的残差构造损失函数,在无任何响应数据的情况下得到了更具泛化能力的神经网络模型,同时避免了传统数值方法在求解计算力学问题时所需的微分、积分公式推导以及繁重的建模和划分网格等前处理工作。本文探究了PINN及其域分解的扩展物理信息神经网络XPINN(eXtended Physics-informed Neural Networks)在求解稳态和瞬态FGMs热传导问题时的适用性,讨论了网络结构对预测结果的影响。研究结果表明,PINN/XPINN在解决几何复杂的稳态和瞬态FGMs热传导问题时仍具有较高的可靠性和简洁的求解流程,同时,为极端环境下求解复杂多场耦合和夹杂等问题提供了新思路。     

细分曲面边界元法的黏附吸声材料结构拓扑优化分析

等几何分析采用样条基函数构造几何模型和实施变量近似,实现了计算机辅助设计和辅助工程的无缝连接,并已广泛应用于弹性力学、电磁场和位势问题等领域.然而直接采用等几何方法难以构造复杂模型,限制了该方法在大规模实际工程问题上的应用.细分曲面法可用于克服这一问题,该方法对传统模型的离散网格进行细分和拟合操作,构造出极限光滑曲面,连续性更高,对复杂结构的适用性更强.该方法主要有以下优点:(1)适用于任意拓扑结构;(2)数值计算稳定;(3)实施简单;(4)局部细化与连续性控制.由于该方法在复杂结构模型构造方面具有较强的灵活性和便利性,已被广泛应用于航空航天、汽车、动画、游戏制作等建模领域.将细分曲面法与边界元法相结合进行结构声学分析,几何场与物理场均采用箱样条基函数进行插值近似.以黏附吸声材料结构的声散射问题为例,建立吸声材料分布拓扑优化数学模型,并采用移动渐进线算法进行设计变量更新,最终获得最优材料分布.      

细分曲面边界元法的黏附吸声材料结构拓扑优化分析

等几何分析采用样条基函数构造几何模型和实施变量近似,实现了计算机辅助设计和辅助工程的无缝连接,并已广泛应用于弹性力学、电磁场和位势问题等领域.然而直接采用等几何方法难以构造复杂模型,限制了该方法在大规模实际工程问题上的应用.细分曲面法可用于克服这一问题,该方法对传统模型的离散网格进行细分和拟合操作,构造出极限光滑曲面,连续性更高,对复杂结构的适用性更强.该方法主要有以下优点:(1)适用于任意拓扑结构;(2)数值计算稳定;(3)实施简单;(4)局部细化与连续性控制.由于该方法在复杂结构模型构造方面具有较强的灵活性和便利性,已被广泛应用于航空航天、汽车、动画、游戏制作等建模领域.将细分曲面法与边界元法相结合进行结构声学分析,几何场与物理场均采用箱样条基函数进行插值近似.以黏附吸声材料结构的声散射问题为例,建立吸声材料分布拓扑优化数学模型,并采用移动渐进线算法进行设计变量更新,最终获得最优材料分布.     

DCM 相位梯度提取技术在光力学图像信息识别中的应用

为解决相位去包裹对光力学图像处理的影响问题,本文将一种直接从数字编码傅立叶变换或相移条纹图中逐相素提取相位梯度的方法应用于光力学图像信息识别中。该方法是一个广泛应用于光纤传感器中的微分正交乘法器(DCM)的正、余弦相位检波技术。这一技术算法简单,不需计算反正切,不受相位去包裹方法标准的限制,非常适于应用在以相位梯度作为最终所需识别的信息的光力学图像处理中。本文以载波等达因条纹图应力场计算及四幅图相移等达因法求板的接触应力为例,对该技术在光力学图像信息识别中的应用进行了论证。     

纵横载荷作用下功能梯度梁的弯曲和过屈曲

基于一阶非线性梁理论和物理中面概念,导出了纵横向载荷作用下功能梯度材料(FGM)梁非线性弯曲和过屈曲问题的控制方程,并获得了该问题的精确解;据此解研究了梯度材料性质、外载荷、横向剪切变形以及边界条件等因素对功能梯度材料梁非线性力学行为的影响,分析中假设功能梯度材料性质只沿梁厚度方向,并按成分含量的幂指数函数形式变化。结果表明:纵横载荷共同作用下,功能梯度梁的弯曲构形将有无限多个;随着梯度指数的增大,梁的变形减小,临界载荷升高;随着长高比的增大,横向剪切变形的影响减小。     

具有切应力梯度的平行平板流动腔的构造

在近年来关于流体切应力与细胞力学行为之间关系的研究中,流体切应力梯度被认为是诱发动脉粥样硬化和内膜增生的重要因素之一.本文探讨如何利用常用的离体细胞力学实验工具--平行平板流动腔模拟具有梯度的定常流切应力环境.结果表明,根据Hele-Shaw流的原理和常用复势W(Z)=AZn(n＞1)的特性,可构造出具有各种切应力梯度分布的流动腔.与其它模拟切应力梯度的方法比较,本文的方法更加简洁、可行.     

功能梯度材料结构热力耦合概率分析方法

基于自洽平均微观力学W-T(Wakashima-Tsukamoto)模型和拉丁超立方抽样,建立一种功能梯度平板热力耦合概率分析方法.以材料物理属性和空间分布的不确定性随机参数作为概率分析的输入,以热应力作为输出,基于本征正交分解POD(Proper Orthogonal Decomposition)对热应力的随机时间历...     

Shape gradient and classical gradient of curvatures: driving forces on micro/nano curved surfaces

Recent experiments and molecule dynamics simulations have shown that adhesion droplets on conical surfaces may move spontaneously and directionally. Besides, this spontaneous and directional motion is independent of the hydrophilicity and hydrophobicity of the conical surfaces. Aimed at this important phenomenon, a general theoretical explanation is provided from the viewpoint of the geometrization of micro/nano mechanics on curved surfaces. In the extrinsic mechanics on micro/nano soft curved surfaces, we disclose that the curvatures and their extrinsic gradients form the driving forces on the curved spaces. This paper focuses on the intrinsic mechanics on micro/nano hard curved surfaces and the experiment on the spontaneous and directional motion. Based on the pair potentials of particles, the interactions between an isolated particle and a micro/nano hard curved surface are studied, and the geometric foundation for the interactions between the particle and the hard curved surface is analyzed. The following results are derived: (a) Whatever the exponents in the pair potentials may be, the potential of the particle/hard curved surface is always of the unified curvature form, i.e., the potential is always a unified function of the mean curvature and the Gaussian curvature of the curved surface. (b) On the basis of the curvature-based potential, the geometrization of the micro/nano mechanics on hard curved surfaces may be realized. (c) Similar to the extrinsic mechanics on micro/nano soft curved surfaces, in the intrinsic mechanics on micro/nano hard curved surfaces, the curvatures and their intrinsic gradients form the driving forces on the curved spaces. In other words, either on soft curved surfaces or hard curved surfaces and either in the extrinsic mechanics or the intrinsic mechanics, the curvatures and their gradients are all essential factors for the driving forces on the curved spaces. (d) The direction of the driving force induced by the hard curved surface is independent of the hydrophilicity and hydrophobicity of the curved surface, explaining the experimental phenomenon of the spontaneous and directional motion.     

Operator-splitting methods for the incompressible Navier-Stokes equations on non-staggered grids. Part 1: First-order schemes

New implicit finite difference schemes for solving the time-dependent incompressible Navier-Stokes equations using primitive variables and non-staggered grids are presented in this paper. A priori estimates for the discrete solution of the methods are obtained. Employing the operator approach, some requirements on the difference operators of the scheme are formulated in order to derive a scheme which is essentially consistent with the initial differential equations. The operators of the scheme inherit the fundamental properties of the corresponding differential operators and this allows a priori estimates for the discrete solution to be obtained. The estimate is similar to the corresponding one for the solution of the differential problem and guarantees boundedness of the solution. To derive the consistent scheme, special approximations for convective terms and div and grad operators are employed. Two variants of time discretization by the operator-splitting technique are considered and compared. It is shown that the derived scheme has a very weak restriction on the time step size. A lid-driven cavity flow has been predicted to examine the stability and accuracy of the schemes for Reynolds number up to 3200 on the sequence of grids with 21 × 21, 41 × 41, 81 × 81 and 161 × 161 grid points.     

Variation transforming analysis (I)

This paper is concerned with operator variation from the transforming point of view, and presents some new concepts and new relations. Related problems and concepts include: convex operator, reciprocity set, reciprocity principles, unilateral variation principles, (H, H ₁,H ₂) — generalized solution and operator-differential equation, etc..Variation principle and variation method is a sort of basic concept and method in the analysis of many problems concerning mathematics, mechanics, physics and control theory. The change and transformation of variation or stationary value relations is the most fundamental form and process of these princeples and methods. This paper analyzes some properties of operator variation and stationary value from the transforming point of view and presents some new concepts, unifiedly treats and extends the reciprocity theorem of classical variation, variation theorem of quadratic functional and unilateral variation principles. Furthermore, we investigate the calculus of normed ring and the solution of a sort of operator differential equation.     

PERTURBATION SOLUTION OF THE WEAK-NONLINEAR DIFFERENTIAL EQUATION WITH δ-FUNCTION

In this paper, starting from some fundamental properties of Heaviside function and δ-function, making use of singular perturbation methods we provide a method of finding the asymptotic analytic solution of equation where M is an n -order linear differential operator, f(u) is a polynomial. By means of this method, we discuss some examples concretely. The results can be explained satisfac-torily in physics. If we deal with linear problem by this method, the result will agree with that drawn from theorem of impulse.     

Eigenfunction expansion method and its application to two-dimensional elasticity problems based on stress formulation

This paper proposes an eigenfunction expansion method to solve twodimensional （2D） elasticity problems based on stress formulation. By introducing appropriate state functions, the fundamental system of partial differential equations of the above 2D problems is rewritten as an upper triangular differential system. For the associated operator matrix, the existence and the completeness of two normed orthogonal eigenfunction systems in some space are obtained, which belong to the two block operators arising in the operator matrix. Moreover, the general solution to the above 2D problem is given by the eigenfunction expansion method.     

Interfacial PIV to resolve flows in the vicinity of curved surfaces

PIV measurements near a wall are generally difficult due to low seeding density, low velocity, high velocity gradient, and strong reflections. Such problems are often compounded by curved boundaries, which are commonly found in many industrial and medical applications. To systematically solve these problems, this paper presents two novel techniques for near-wall measurement, together named Interfacial PIV, which extracts both wall-shear gradient and near-wall tangential velocity profiles at one-pixel resolution. To deal with curved walls, image strips at a curved wall are stretched into rectangles by means of conformal transformation. To extract the maximal spatial information on the near-wall tangential velocity field, a novel 1D correlation function is performed on each horizontal pixel line of the transformed image template to form a “correlation stack”. This 1D correlation function requires that the wall-normal displacement component of the particles be smaller than the particle image diameter in order to produce a correlation signal. Within the image regions satisfying this condition, the correlation function yields peaks that form a tangential velocity profile. To determine this profile robustly, we propose to integrate gradients of tangential velocity outward from the wall, wherein the gradient at each wall-normal position is measured by fitting a straight line to the correlation peaks. The capability of Interfacial PIV was validated against Particle Image Distortion using synthetic image pairs generated from a DNS velocity field over a sinusoidal bed. Different velocity measurement schemes performed on the same correlation stacks were also demonstrated. The results suggest that Interfacial PIV using line fitting and gradient integration provides the best accuracy of all cases in the measurements of velocity gradient and velocity profile near wall surfaces.     

基于超导的全张量重力梯度敏感头设计

基于超导的迈斯纳效应与超导量子干涉技术,结合柔性并联机构理论,设计一种基于超导的全张量重力梯度敏感头。敏感头采用6个完全相同的同时具有移动与转动自由度的敏感结构,对称的布置在正六面体外,对称面的两个敏感结构相对旋转90°成垂直状态。轴向间距使两个敏感结构直接测量重力梯度轴向分量,相互垂直使两个敏感结构既可以测量重力梯度交叉分量,也可以测量共模角加速度。利用超导线圈的电感变化响应质量块位移,进一步通过超导回路将其转变为磁场变化,并由超导量子干涉器进行检测。敏感结构采用8分支的柔性并联机构支承,构成空间对称的形式,可以实现对称的力学特性,保证各处的柔性铰链产生均匀变形,减少非对称的偏移,避免单一铰链的应力集中,具有沿轴移动刚度与绕轴转动刚度小、非设计的寄生误差方向刚度大的优势。在惯性系下的全张量重力梯度值可由坐标变换得到,可以预期得到1E的测量精度。     

From the second gradient operator and second class of integral theorems to Gaussian or spherical mapping invariants

By combining of the second gradient operator, the second class of integral theorems, the Gaussian-curvature-based integral theorems and the Gaussian （or spherical） mapping, a series of invariants or geometric conservation quantities under Gaussian （or spherical） mapping are revealed. From these mapping invariants important transformations between original curved surface and the spherical surface are derived. The potential applications of these invariants and transformations to geometry are discussed