双柱状颗粒在线性剪切流场中的动力学分析

为研究柱状颗粒在线性剪切流场中的运动状态和受力情况,本文以颗粒长径比为2,颗粒之间的初始距离ΔS_Py=4D为例,基于直接力浸入边界法数值模拟了双柱状颗粒在三维线性剪切流场中的运动过程。根据模拟结果分析了柱状颗粒周围流场参数分布,在考虑壁面对颗粒的影响和颗粒之间相互影响的条件下,研究了颗粒的受力和运动的变化,探索了流体曳力导致柱状颗粒迁移和转动的规律。研究结果表明,双柱状颗粒在线性剪切流场中易向速度大的流体区域运动;前后两颗粒运动状态和轨迹不同,颗粒之间距离较近时,曳力会产生较大的波动;只有当颗粒在壁面附近时,滞后颗粒才能追上领先颗粒,两颗粒发生牵引、翻滚和分离过程。     

近源堆沉积土石体细观结构与抗剪强度指标关联分析

细观结构是认知土石体力学行为本质的关键科学问题。本文通过引入物理学和数学等方法,借助颗粒物质力学理论,从几何排列与接触力的空间分布来定量刻画水平固结与山前坡地堆积两种典型环境下土石体的细观结构特征,并建立其与抗剪强度指标的关联。研究表明,（1）两种环境的土石体在细观结构上存在较大差异。在几何排列上,水平固结环境下的土石体具有长程无序和短程有序的特点,坡地堆沉积环境下的土石体表现出了无序的无定形结构;在接触力与单位接触向量的空间分布上,两者较为相似,绝大多数接触力以小于均值接触力的形式存在,其概率密度曲线P（f）呈幂函数衰减;90%以上接触方位角集中在40°~160°和220°~340°范围内。（2）基于径向分布函数、接触力概率密度和单位接触向量分别定义细观结构的特征量K_a,K_s和K_o,发现三个特征量的增大对内摩擦角呈线性促进作用,对黏聚力呈非线性削弱作用。     

海南岛清澜-八门湾潟湖潮汐水动力的模拟研究

基于DELFT3D模型研究了清澜-八门湾潟湖的水位分布和潮汐波内部结构,以及由于人类开垦引起的海岸线变化（以1962年、1985年和2008年为例）对潟湖水动力特性的影响。结果表明,清澜-八门湾潟湖潮汐是由多种分潮耦合而成的复杂驻波,其中K₁,O₁,M₂,S₂和M₄分潮的影响最大。由湾外向湾内传输,由于湾内红树林和浅滩引起的底部损耗增加,M₂,S₂,K₁和O₁分潮幅值减弱;M₄分潮幅值增强,表现出明显的浅水增幅效应;M₂和S₂分潮相位在文教河和文昌江领域表现出明显的干湿效应。不同年代海岸线的研究表明,1985年和2008年间,人类复垦导致潟湖及其潮汐汊道附近的红树林和滩涂区域严重破坏,海岸线缩减,引发了水位降低、纳潮量减少和潟湖潮汐汊道底摩擦弱化,从而削弱了干湿、潮呛和浅水效应。后果提示持续的人类复垦活动将会引发清澜-八门湾潮汐水动力环境的进一步恶化,可能导致未来发生更大的自然灾害。     

压电压磁材料中正n边形孔边裂纹分析

基于电磁复合材料力学,运用Stroh型公式和复变函数方法,针对压电压磁材料中含正n边形孔边裂纹反平面问题进行了研究。利用Schwarz-Christoffel变换技术,结合Cauchy积分公式和留数定理,导出了磁电全非渗透型边界条件下任意正n边形裂纹尖端场强度因子和能量释放率的解析解。当缺失磁场时,所得解退化为已有结果,以此验证方法的有效性。通过数值算例,对比分析了n=3,n=4和n=5三种特殊情形对应的孔口边长、裂纹长度和受到的力、电和磁载荷对等效场强度因子和无量纲能量释放率的影响规律。研究结果发现,正n边形孔洞的尺寸和裂纹长度均会促进裂纹扩展,且前者的影响更显著一些;正n边形边的数量增加会阻止裂纹的扩展;在磁电全非渗透型边界条件下,机械载荷始终促进裂纹的扩展,电位移载荷可以促进或抑制裂纹的扩展,磁载荷对裂纹的扩展贡献较少。本研究结果适用于任意正n边形孔边裂纹求解问题,为压电压磁材料元器件的优化设计和断裂特性分析提供了新思路。     

投影方向极值原理的滑裂面应力矢量极限平衡抗滑稳定计算理论

根据有限元成果的滑动面切向应力达到可用极限值的抗滑稳定极限平衡概念,建立投影方向极值原理单元滑面应力矢量极限平衡抗滑稳定计算理论模型,揭示不平衡虚拟耗散力、耗散能和极值条件方程的力学本质联系。根据投影极值方向平衡状态时的垂直方向还存在不平衡耗散力,提出复杂滑裂面极限平衡状态接近程度系数λ及矢量安全度系数K_λ评价抗滑稳定性能的最佳路径。通过边坡考题和坝基深层抗滑稳定算例的势能极值规律存在性和极值点的吻合度,证明了理论模型的合理可靠性。有限元应力矢量理论模型的建立及其安全度系数K_λ是极限平衡抗滑稳定计算理论发展的重要基础理论突破和全新理论研究成果。     

基于COMSOL的脆性材料相场断裂模型

相场法通过一系列微分方程描述材料断裂过程,避免了繁琐的裂纹面追踪,在模拟裂纹的萌生、扩展和分叉等方面具有优势。介绍了基于相场法的脆性材料断裂模型,给出了脆性材料断裂问题相场法控制方程的推导过程,提出了基于分步迭代法在COMSOL中实现脆性材料相场断裂模型的方法。再现了脆性材料单元模型和单边缺口平板受拉及受剪作用下的开裂过程,模拟的裂纹扩展路径与已有文献的结果相近,验证了程序的合理性。针对脆性材料相场断裂模型包含的诸多参数,采用Morris法对影响荷载-位移关系的脆性材料断裂模型参数进行了全局敏感性分析,结果表明,杨氏模量(E)、临界能量释放率(G_c)和位移增量(Δu_x)是影响模型荷载-位移关系输出结果的主要参数。基于COMSOL实现的相场断裂模型能够有效模拟脆性材料的裂纹萌生和扩展断裂过程,模型参数E,G_c和Δu_x对材料断裂性能的提升或模型参数反演效率的提高具有重要影响。     

凹形多面体离散单元的水平集函数接触算法

多面体模型理论上可构造任意颗粒形态,然而受单元接触算法的限制,仅用于凸形颗粒材料的离散元模拟。对于具有凹形特征的多面体单元,单个接触点的搜索算法难以精确计算单元间的作用力。考虑多面体单元间存在单个或多个接触点的计算特性,本文发展了适用于凸形和凹形多面体颗粒材料的水平集函数接触算法。该方法通过点-三角形单元距离计算方法和奇-偶数判定方法建立多面体单元的零水平集函数和空间水平集函数,并对水平集函数进行三线性插值,可得到多面体单元间的单个或多个接触点。为检验水平集函数接触算法的可靠性,对球形和凹形多面体颗粒材料的堆积和倒塌过程进行离散元模拟,并分析颗粒形状对堆积密度和休止角的影响规律。     

稀土钇对Ni-W-ZrO2复合镀层性能的影响

通过脉冲电沉积技术在铝合金基体表面制备了Ni-W-ZrO₂复合镀层,利用光学显微镜、扫描电镜、X射线衍射仪、能谱仪、维氏硬度仪、电化学工作站和马弗炉等设备分别对复合镀层的截面形貌、元素分布、表面硬度、耐蚀性和抗高温氧化性分别进行了表征和测试,研究了Y（NO₃）₃添加量对复合镀层性能的影响.研究结果表明：所制备复合镀层截面平整、成分分布均匀,基体与镀层之间结合良好,ZrO₂颗粒弥散分布于Ni-W基质金属之中,镀层厚度约为130 μm;Y（NO₃）₃的加入能显著提高复合镀层的硬度、抗高温氧化性和耐蚀性;当硝酸Y添加量为1.5 g/L时,Ni-W-ZrO₂复合镀层的综合力学性能达到最佳.     

柱塞泵多目标优化设计及CFD仿真分析

液压泵噪声是液压系统的主要噪声源,针对轴向柱塞泵的流致振动噪声,提出一种改善泵配流特性的设计方案。首先,根据柱塞泵的工作原理对柱塞腔压力特性和泵出口流量特性准确建模并求解。通过分析压力冲击和流量脉动对错配角（φ₀）的响应,得φ₀=4°为佳。利用一种多目标遗传算法（NSGA-Ⅱ）,以减小压力超调量和流量脉动率为目标,对三角槽结构进行了优化;并获得该多目标优化问题的Pareto最优解集,通过对最优解集的分析知,深度角θ₁=16°且宽度角θ₂=85°时较为理想。最后,为了验证模型的正确性,建立流体域计算流体动力学（CFD）模型,对比两种模型计算结果发现吻合较好,能够相互验证。利用CFD分析结果可视化的特点,从柱塞泵流场的角度,进一步分析了泵压力冲击以及流量脉动产生的原因。     

基于转捩/尺度适应模型与FW-H声学方程的气动噪声数值模拟研究

采用计算流体力学(CFD)同"声类比"相结合的方法进行噪声模拟,利用CFD数值模拟 MD30P30N多段翼型失速攻角附近流场以及其气动特性来校核近场精度,进一步通过结合可穿透数据面的FW-H声学方程进行气动噪声分析。为了准确捕捉近场流场信息,为噪声预测提供可靠的声源精度,本文基于k-ω剪应力输运(SST)湍流模型,建立了尺度适应分离流(SAS)模型,并采用γ-Re_θt转捩模型耦合k-ω SST湍流模型建立了边界层转捩数值模拟技术;充分利用尺度适应模型在边界层表现为雷诺平均(RANS)方法这一特点,将γ-Re_θt转捩模型与SAS模型结合,建立针对包含转捩、分离现象复杂流场的数值模拟技术。文中以RANS为控制方程,分别采用全湍流k-ω SST模型、γ-Re_θtSAS转捩/分离流模型对多段翼进行近场数值模拟计算,结合可穿透数据面的FW-H声学方程进行气动噪声预测,在分析其对流场及气动噪声影响的基础上,得出了几点结论。     

大学物理《力学》中的自由度分析

对《力学》中的物体自由度进行多方面分析,以深化教学、提高学生正 确分析物理问题的能力.使用实际教学分析的研究方法,在《力学》范围内讨论自由度与坐标、 自由与约束的关系并得以下结论: (1) 同一物体的自由度随其所在的``空间'不同而不同, 不因坐标系的选取不同而 异, 在同类参考系中不因参考系的动静而有别;(2)自由度遵循叠加原理. 讨论了质点系的总自由度及相关计算问题,并指出研究《力学》中自由度的意义.     

Application of the theory of branching of solutions of nonlinear equations in the estimation of pole perturbance of linear closed-loop systems

When the state and input matrices of a multivariable linear time-invariant system are perturbed, the problem of the estimation of pole perturbance of closed-loop system is considered by making use of the theory of branching of solutions for nonlinear equations.     

两自由度振动系统的斜碰撞分析

研究斜碰撞振动系统动力学的一个关键问题是对系统在碰撞前后的状态进行合理描述和正确计算．针对两弹性体斜碰撞问题，基于瞬间碰撞假设，提出了采用步进冲量来分析和求解斜碰撞前后的状态关系；并以弹簧摆和振子组成的两自由度斜碰撞振动系统为例，具体介绍了该算法如何实现．用解析方法讨论了该系统在斜碰撞过程中可能出现的各种力学现象，将冲量步进算法得到的数值解与解析结果进行对比，取得了完全一致的结果．该数值方法能适应多种斜碰撞问题的计算．     

One type of solutions of conduction of nervous pulses and the stability of action potential

This paper reports a type of laws which governs action poten-tial of nervous impulses,and it is discussed by general form-nonlinear dispersive process.We find that the nervous wave is aslowly varying amplitude solitary wave in the small dispersivecase.we prove that the solitary wave is not generated in the or-dinary dispersion,but a travelling wave with varying amplitudesmay be obtained.The stability of various possible action poten-tials and bifurcation in overdamped case are also discussed inthis paper.     

The distribution of zeroes of solutions of neutral equations

1.IntroductionThestudyofoscillator}'behaviorofneutraldifferentialequationsisarelativealnewfieldandisveryinterestinginapplications.Inrecentyears.theoscillatorytheoryfortheseequationshasbeenextensivelydeveloped.WerefertotherecentbookbyGyoriandLadaSI'I,Baino…     

沙漠地貌演化过程的湍流大涡模拟研究

采用大涡模拟和混合粒径群体沙粒运动算法结合内置边界方法来模拟沙波纹与沙丘的演化和发展,在平坦沙区形成沙波纹的同时,在沙丘上同样也出现了沙波纹的结构,即沙波纹的发展出现了多尺度现象.此外,还分析了粒径对沙波纹与沙丘自身形态演化和发展的影响,以及沙丘对平坦沙区沙波纹形成发展的作用.     

CHARACTERISTICS OF SUBDIFFERENTIALS OF FUNCTIONS

ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ（郭兴明）ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ￥ＧｕｏＸｉｎｇｍｉｎｇ（ＲｅｃｅｉｖｅｄＪｕｎｅ１６，１９９５...     

Motion equations of the dynamic theory of consolidation from the point of view of the theory of transient pipe flows

An elastic fluid-saturated porous medium is modeled as a bundle of parallel cylindrical tubes aligned in a direction parallel to the fluid movement. The pore space is filled with viscous compressible liquid. A cell model and the theory of transient pipe flow are used to derive one-dimensional governing equations in such media. All macroscopic constants in these equations are defined by the individual material constants of the fluid and solid. The interaction force includes an additional term not found in Biot's theory.     

Identification of viscoelastic characteristics of composite materials on the basis of results of an experimental and theoretical analysis of the dynamic behavior of hemispherical shells

A method is proposed for determining the stiffness and rheological characteristics of composite materials, which is based on minimizing the disagreement between experimental data and results of numerical simulations of deformation of hemispherical shells under explosive loading. The damping characteristics of randomly reinforced polymer materials are analyzed with the use of this method. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 126–133, May–June, 2006.     

On some problems of stability of systems with an infinite number of degrees of freedom

The problem of stability of equilibria of a physical pendulum with a nonstretchable thread attached to it is considered from the standpoint of the Lagrange theorem on stability and its inversion. Specific difficulties which one faces when studying an infinite dimensional mechanical system are discussed. A new approach to the study of stability with respect to two metrics is suggested. The influence of resonant phenomena on the motion of the shortened (linearized) system is considered.

---

Sommario Si considera il problema di stabilità di punti di equilibrio di un pendolo fisico con un filo inestensibile appeso ad esso dal punto di vista del teorema di Lagrange sulla stabilità e sulla sua inversione. Difficoltà specifiche relative allo studio di un sistema meccanico di dimensione infinita sono discusse. Si suggerisce un nuovo metodo per lo studio della stabilità rispetto a due metriche. L'influenza di fenomeni di risonanza sul moto del sistema ridotto (linearizzato) è considerata.