水力机组制动过程的分段递推数值解法

水力机组制动过程的数学模型为非线性一阶常微分方程,目前尚不法得到该方程的精确解析解.本文提出的分段递推法通过对该方程在小分段内进行合理的线性化处理,用拉普拉斯数学变换求出小分段内的近似解析解,得到递推算式,利用递推算式得到全部数值解.使用分段递推法对XAJ水力机组在各种制动方式下进行数值求解,得到了制动过程,并分析了各种制动方式的特点.同时,在自由制动和飞逸过程时与龙格—库塔法进行了比较,并在电气制动时用实验数据进行了验证.结果显示分段递推法是精确的,能满足科学研究和工程应用需要,是求解类似方程的一种简便通用方法.     

含多余广义坐标的拉格朗日方程及实例应用

 推导了含有多余广义坐标的离散动力系统的拉格朗日方程,利用该方程给出了一类双线摆在平衡位置小幅振动时近似频率的计算公式,与传统的拉格朗日乘子法相比,更适合对动力系统动力学性能的定性分析.     

流固耦合的周期加强板的振动及声辐射研究

研究了流体负载下的无穷大双周期加强板, 在周期谐振力作用下的振动响应和声辐射,并提出了一种基于有限元和空间波数法的半解析半数值方法. 首先利用有限元的方法对周期结构进行单元离散, 并将结构对薄板的作用力等效为节点力的作用. 然后通过周期结构的振动方程, 结合薄板与结构的位移边界条件, 建立了节点力与薄板节点位移的函数方程. 最后应用空间波数法和傅里叶变换, 并采用数值计算的方法求解出薄板的节点位移, 得到了周期加强板关于离散节点位移的振动和辐射声压方程. 在数值算例中, 对该方法的正确性进行了验证, 并且分析了周期结构对薄板的振动和声辐射的影响.     

一种O(n)算法复杂度的递推绝对节点坐标法研究

相比于传统的浮动坐标法,绝对节点坐标法(absolute nodal coordinate formulation,ANCF)在处理柔性体非线性大变形问题上具有显著优势,但是对于ANCF的求解目前主要依据拉格朗日方程等分析力学原理建立微分代数方程(differential algebraic equation,DAE)进行,其算法复杂度为O(n2)或O(n3)(n为系统自由度),且求解过程存在位置或速度的违约问题.据此,研究了一种O(n)算法复杂度的递推绝对节点坐标法(recursive absolute nodal coordinate formulation,RANCF).该方法采用ANCF描述大变形柔性体,借鉴铰接体递推算法(articulatedbody algorithm,ABA)思路建立多柔体系统逐单元的运动学和动力学递推关系,得到微分形式的系统动力学方程(ordinary differential equation,ODE).在ODE方程中,系统广义质量阵为三对角块矩阵,通过恰当的矩阵处理,可以得到逐单元求解该方程的递推算法.在此基础上,给出了RANCF算法的详细流程,并对流程中每个步骤进行了细致的算法效率分析,证明了RANCF是算法复杂度为O(n)的高效算法.RANCF方法保留了ANCF对大转动、大变形多柔体系统精确计算的优点,同时极大地提升了算法效率,特别在处理高自由度复杂多柔体系统中具有显著优势.并且该方法采用ODE求解,无DAE的违约问题,因此具有更高的算法精度.最后,在算例部分,通过MSC.ADAMS仿真软件、能量守恒测试、算法复杂度曲线对RANCF的正确性、计算精度和计算效率进行了验证.     

基于黎曼几何的非线性振子动力学分析递推解析算法

基于黎曼几何和变分原理,推导了黎曼流形上非线性耗散动力系统的二阶微分动力学方程,并运用流形收缩的概念将动力学方程离散化,进而建立了相应的递推求解格式。选取3个自治非线性阻尼振子系统,分别采用递推解析算法和龙格库塔法求解微分动力学方程,并比较分析了不同的时间步长下两种算法的计算耗时。结果表明,与龙格库塔法相比,基于黎曼几何的递推算法不仅能得到每一时步的解析表达式,而且计算耗时短,计算效率高。基于黎曼流形的动力学方程递推算法为非线性动力学系统的解析求解提供了新思路。     

基于黎曼几何的非线性振子动力学分析递推解析算法

基于黎曼几何和变分原理，推导了黎曼流形上非线性耗散动力系统的二阶微分动力学方程，并运用流形收缩的概念将动力学方程离散化，进而建立了相应的递推求解格式。选取3个自治非线性阻尼振子系统，分别采用递推解析算法和龙格库塔法求解微分动力学方程，并比较分析了不同的时间步长下两种算法的计算耗时。结果表明，与龙格库塔法相比，基于黎曼几何的递推算法不仅能得到每一时步的解析表达式，而且计算耗时短，计算效率高。基于黎曼流形的动力学方程递推算法为非线性动力学系统的解析求解提供了新思路。     

考虑局部非线性结构振动的模型预测控制

针对具有局部非线性的结构振动主动控制问题,首先,利用Newmark-β法推导了局部非线性结构动力响应的时域显式递推表达式。然后,基于该表达式,一方面建立了局部非线性结构振动预测模型,并借鉴模型的线性化思想对非线性环节进行处理,提出考虑局部非线性结构振动的模型预测控制算法;另一方面,在获取各仿真时刻的预测初始状态时,为避免直接对非线性动力学方程所有维度的完全迭代求解,推导了时域显式降维迭代方法,从而使在线计算效率得到显著提高。最后,分别以局部非线性弹簧振子系统在简谐激励作用下的无控动力学分析和含粘滞阻尼器毗邻框架建筑结构在地震激励作用下的受控振动为例,进行数值仿真,结果表明,本文算法计算准确且控制效果较好。     

热环境下壁板非线性颤振分析

基于一阶活塞气动力理论,采用Von Karman大变形应变-位移关系建立了无限展长壁板热环境下颤振方程,采用伽辽金方法对方程进行离散处理.取温度为分叉参数,研究壁板颤振时的分叉及混沌等复杂动力学特性.结果表明:温度载荷降低了系统的颤振临界动压,改变了颤振特性.在整个分岔参数范围内,系统呈现出较为复杂的变化,包括衰减振动、极限环振动、拟周期振动和混沌型振动.当考虑材料热效应时,系统的颤振动压将进一步降低,其响应也表现出更为丰富的非线性动态力学行为.     

索-梁耦合系统解的稳定性分析

研究了在惯性参考系中弹性斜拉索与悬臂梁耦合结构的非线性动力学问题,利用Hamilton原理建立了索-梁耦合系统的非线性动力学方程,利用Galerkin方法将索-梁耦合系统的非线性运动偏微分方程离散为一组常微分方程,然后利用多尺度法分析研究索-梁耦合动力学系统的非线性振动,对耦合系统解的稳定性进行了分析,用Runge-Kutta法对数学模型进行数值计算,并提出对工程有实际意义的结论.     

离散系统通用动力学方程求解算法的研究进展

离散系统中的颗粒物在凝并、破碎、冷凝/蒸发、成核、沉积等事件作用下颗粒尺度分布的时间演变由通用动力学方程所描述.该方程为一典型的部分积分微分方程,普通数值方法难以求解.本文详细介绍了求解通用动力学方程的矩方法、分区法、离散法、离散-分区法、MonteCarlo方法等几种算法的原理、优缺点和最新的研究进展,并着重介绍了MonteCarlo算法,包括基于时间驱动Monte Carlo方法、基于事件驱动MonteCarlo方法、常数目法、常体积法以及多重Monte Carlo算法.      

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

The 7th International Conference on Fluid Mechanics May 24-27,2015,Qingdao,China

正http://www.icfm7.org First Announcement and Call for PapersThe objective of International Conference on Fluid Mechanics(ICFM)is to provide a forum for researchers to exchange new ideas and recent advances in the fields of theoretical,experimental,computational Fluid Mechanics as well as interdisciplinary subjects.It was successfully convened by the Chinese Society of Theoretical and Applied Mechanics(CSTAM)in Beijing(1987,     

INFORMATION FOR CONTRIBUTORS

Contributions： The Journal, Acta Mechanica Solida Sinica, is pleased to receive papers from engineers and scientists working in various aspects of solid mechanics. All contributions are subject to critical review prior to acceptance and publication.