自适应无网格热弹塑性接触模型研究

提出一种自适应无网格热弹塑性接触求解模型,求解接触问题的线性规划-增量初应力法与基于应变能梯度的自适应无网格法相结合,给出了模型计算理论和算法实现.通过圆柱体与弹塑性平面热弹塑性接触算例对模型进行验证.对是否考虑材料应变硬化,是否考虑摩擦力和热输入,是否考虑材料屈服强度温度相关等情况的两种算例进行了讨论.结果表明,该模型能有效地求解考虑不同情况下的热弹塑性接触问题,在较真实地模拟接触状况的同时,具有较高的计算精度和计算效率.     

无网格局部Petrov-Galerkin方法在弹塑性断裂力学问题中的应用

采用无网格局部Petroy-Galerkin方法来分析弹塑性断裂力学问题.这种无网格方法采用移动最小二乘法(MLS)来构造近似试函数和采用Heaviside函数作为加权残值法中的权函数,由于近似函数不满足KroneckerDelta条件,因此采用直接插值法来施加本质边界条件.如果不考虑体力,所形成的整体刚度矩阵只包含局部边界积分,而不包含局部域积分和奇异积分.采用增量Newton-Raphson迭代法来求解弹塑性增量形式的局部Petrov-Galerkin方程.数值算例结果表明,该文方法对于弹塑性断裂力学问题的求解是可行的和有效的,并且所得到的结果具有较好的精度.     

用无网格局部径向点插值法分析非均质中厚板的弯曲问题

用无网格局部径向点插值法分析了非均质中厚板的弯曲问题.利用虚位移原理推导了中厚板的离散系统方程.采用径向基函数耦合多项式基函数来近似试函数,用四次样条函数作为加权残值公式中的权函数.所构造成的形函数具有Kronecker delta性质,可以很方便地施加本质边界条件.此方法不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行,是一种真正的无网格方法.在计算过程中,取积分中的高斯点的材料参数来模拟问题域材料特性的变化.算例结果表明这种无网格方法具有效率高、精度高和易于实现等优点.     

基于混合网格多维梯度重构的热流预测方法研究

高超声速气动热环境的数值计算对算法和网格的敏感度极高. 随着高超声速飞行器外形日益复杂, 生成高质量的结构网格时间成本呈指数增加, 难以满足工程应用的需求. 非结构/混合网格因具有很强的复杂外形适应能力, 为了缩短任务周期, 有必要在非结构/混合网格上开展高精度的气动热环境数值计算方法研究. 梯度重构方法是影响非结构/混合网格热流计算精度的重要因素之一. 本文通过引入多维梯度重构方法, 发展了基于常规的非结构/混合网格的高精度热流计算方法, 对典型的高超声速Benchmark算例(二维圆柱)进行了模拟, 并与气动力计算广泛采用的Green-Gauss类方法和最小二乘类方法进行了对比. 计算结果表明, 多维梯度重构方法能有效提高非结构/混合网格热流预测精度, 其鲁棒性和收敛性更好. 最后将多维梯度重构方法应用于常规混合网格的三维圆柱和三维双椭球绕流问题, 得到了与实验值吻合较好的热流计算结果, 展现了良好的应用前景.      

非结构网格二阶有限体积法中黏性通量离散格式精度分析与改进

非结构网格二阶有限体积离散方法广泛应用于计算流体力学工程实践中,研究非结构网格二阶精度有限体积离散方法的计算精度具有现实意义. 计算精度主要受到网格和计算方法的影响,本文从单元梯度重构方法、黏性通量中的界面梯度计算方法两个方面考察黏性流动模拟精度的影响因素. 首先从理论上分析了黏性通量离散中的“奇偶失联”问题,并通过基于标量扩散方程的制造解方法验证了“奇偶失联”导致的精度下降现象,进一步通过引入差分修正项消除了“奇偶失联”并提高了扩散方程计算精度;其次,在不同类型、不同质量的网格上进行基于扩散方程的制造解精度测试,考察单元梯度重构方法、界面梯度计算方法对扩散方程计算精度的影响,结果显示,单元梯度重构精度和界面梯度计算方法均对扩散方程计算精度起重要作用;最后对三个黏性流动算例(二维层流平板、二维湍流平板和二维翼型近尾迹流动)进行网格收敛性研究,初步验证了本文的结论,得到了计算精度和网格收敛性均较好的黏性通量计算格式.      

基于混合网格多维梯度重构的热流预测方法研究

高超声速气动热环境的数值计算对算法和网格的敏感度极高.随着高超声速飞行器外形日益复杂,生成高质量的结构网格时间成本呈指数增加,难以满足工程应用的需求.非结构/混合网格因具有很强的复杂外形适应能力,为了缩短任务周期,有必要在非结构/混合网格上开展高精度的气动热环境数值计算方法研究.梯度重构方法是影响非结构/混合网格热流计算精度的重要因素之一.本文通过引入多维梯度重构方法,发展了基于常规的非结构/混合网格的高精度热流计算方法,对典型的高超声速Benchmark算例(二维圆柱)进行了模拟,并与气动力计算广泛采用的Green-Gauss类方法和最小二乘类方法进行了对比.计算结果表明,多维梯度重构方法能有效提高非结构/混合网格热流预测精度,其鲁棒性和收敛性更好.最后将多维梯度重构方法应用于常规混合网格的三维圆柱和三维双椭球绕流问题,得到了与实验值吻合较好的热流计算结果,展现了良好的应用前景.     

非结构网格二阶有限体积法中黏性通量离散格式精度分析与改进

非结构网格二阶有限体积离散方法广泛应用于计算流体力学工程实践中,研究非结构网格二阶精度有限体积离散方法的计算精度具有现实意义.计算精度主要受到网格和计算方法的影响,本文从单元梯度重构方法、黏性通量中的界面梯度计算方法两个方面考察黏性流动模拟精度的影响因素.首先从理论上分析了黏性通量离散中的"奇偶失联"问题,并通过基于标量扩散方程的制造解方法验证了"奇偶失联"导致的精度下降现象,进一步通过引入差分修正项消除了"奇偶失联"并提高了扩散方程计算精度;其次,在不同类型、不同质量的网格上进行基于扩散方程的制造解精度测试,考察单元梯度重构方法、界面梯度计算方法对扩散方程计算精度的影响,结果显示,单元梯度重构精度和界面梯度计算方法均对扩散方程计算精度起重要作用;最后对三个黏性流动算例(二维层流平板、二维湍流平板和二维翼型近尾迹流动)进行网格收敛性研究,初步验证了本文的结论,得到了计算精度和网格收敛性均较好的黏性通量计算格式.     

大变形问题分析的局部Petrov-Galerkin法

在微机电系统(MEMS)的建模和模拟研究中,大变形或大移动要充分予以考虑.用有限元法分析这类问题,由于难以避免的网格畸变,使模拟效率精度降低甚至失效,无网格方法(Meshless Method)则能在分析这类问题时显示出明显的优势,无网格局部Petrov-Galerkin(MLPG)法被誉为是一种有发展前景的真正无网格法.本文进一步发展了MLPG法,通过对任意的离散分布节点采用局部径向基函数构造插值形函数和Heaviside权函数,分析方程采用局部加权弱形式离散,建立了变量仅依赖于初始构型的完全Lagrange分析格式,最后用Newton-Raphson法迭代求解.文中分析了悬臂梁典型算例和微机电开关非线性大变形问题,通过与有限元结果的比较,表明本文提出的大变形问题无网格局部Petrov-Galerkin法具有稳定性好及收敛性快等优点.     

基于修正权函数的多裂纹无网格模拟

本文采用了一种基于不连续场修正权函数的无网格方法来处理二维平面多裂纹问题。相较于传统的无网格断裂不连续场和奇异场模拟方法,修正权函数法算法简便易实现。采用修正权函数处理多裂纹时,只需要对每一段裂纹周围节点的权函数进行修正,就能同时模拟多裂纹不连续位移场和多裂尖奇异场。本文采用基于不连续场修正权函数的无单元Galerkin方法（EFGM）,对Y型裂纹板、十字型裂纹板和孔边双裂纹板进行了分析。数值结果表明,在不引入扩展基函数情况下,通过修正权函数法能够得到精度较高的应力强度因子解,能较好地拟合多裂纹的裂尖奇异场。     

应变能中耦合项对旋转悬臂梁振动频率的影响

大范围运动悬臂梁的动力学建模问题对动力学特性分析及控制系统设计具有极其重要的作用. 当前研究多采用一次近似模型,其忽略了由轴向和横向变形所产生的应变能中的耦合项,然而这些项对动力学特性会产生影响. 通过讨论应变能的选取方式,计入了应变能中的耦合项;利用哈密尔顿原理建立结构的耦合振动模型;再借助瑞利—里兹法,以无大范围运动时的振型函数作为基本解组,得到了结构振动广义特征方程并求解. 通过数值算例对比分析,指出考虑应变能耦合项得到的频率与不考虑应变能耦合项得到的频率存在明显差别.      

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.     

Preface

This special issue of PARTICUOLOGY is devoted to the first UK-China Particle Technology Forum taking place in Leeds, UK, on 1-3 April 2007. The forum was initiated by a number of UK and Chinese leading academics and organised by the University of Leeds in collaboration with Chinese Society of Particuology, Particle Technology Subject Group （PTSG） of the Institution of Chemical Engineers （IChemE）, Particle Characterisation Interest Group （PCIG） of the Royal Society of Chemistry （RSC） and International Fine Particle Research Institute （IFPRI）. The forum was supported financially by the Engineering and Physics Sciences Research Council （EPSRC） of United Kingdom,     

基于鲁棒滤波的捷联导引头视线角速度估计方法

针对捷联导引头无法直接获取视线角速度等信息的问题,研究了鲁棒滤波在大气层外飞行器捷联导引头视线角速度估计中的应用。为了建立非线性滤波估计模型,考虑目标视线角速度的慢变特性,采用一阶马尔科夫模型建立了状态方程;推导了视线角速度的解耦模型,并建立了量测方程;考虑到实际应用中存在系统噪声统计特性失准的问题,基于Huber-Based鲁棒滤波方法,设计了视线角速度滤波器,并完成了基于Huber-Based滤波方法和扩展卡尔曼滤波方法的数学仿真。仿真结果表明Huber-Based滤波方法的视线角、视线角速度及视线角加速度估计精度分别达到0.1140'、0.1423'/s、0.0203'/s2,而扩展卡尔曼滤波方法的视线角、视线角速度及视线角加速度估计精度仅分别为0.6577'、0.6415'/s、0.0979'/s~2。仿真结果证明了该方法可以有效地估计出相对视线角速度等信息,并且在非高斯噪声的条件下,依然可获得较高的估计精度,具有一定的鲁棒性。     

Guide for authors

正Each of the sections below provides essential information for authors.We recommend that you take the time to read them before submitting a contribution to Acta Mechanica Sinica.We hope our guide to authors may help you navigate to the appropriate section.How to prepare a submission This document provides an outline of the editorial process involved in publishing a scientific paper in Acta Mechanica     

Use specification of multiscale materials for life spanned over macro-, micro-, nano-, and pico-scale

Multiscale material intends to enhance the strength and life of mechanical systems by matching the transmitted spatiotemporal energy distribution to the constituents at the different scale, say—macro, micro, nano, and pico,—, depending on the needs. Lower scale entities are, particularly, critical to small size systems. Large structures are less sensitive to microscopic effects. Scale shifting laws will be developed for relating test data from nano-, micro-, and macro-specimens. The benefit of reinforcement at the lower scale constituents needs to be justified at the macroscopic scale. Filling the void and space in regions of high energy density is considered.Material inhomogeneity interacts with specimen size. Their combined effect is non-equilibrium. Energy exchange between the environment and specimen becomes increasingly more significant as the specimen size is reduced. Perturbation of the operational conditions can further aggravate the situation. Scale transitional functions and/or f_j/j+1 are introduced to quantify these characteristics. They are represented, respectively, by , and (f_mi/ma,f_na/mi,f_pi/na). The abbreviations pi, na, mi, and ma refer to pico, nano, micro and macro.Local damage is assumed to initiate at a small scale, grows to a larger scale, and terminate at an even larger scale. The mechanism of energy absorption and dissipation will be introduced to develop a consistent book keeping system. Compaction of mass density for constituents of size 10⁻¹², 10⁻⁹, 10⁻⁶, 10⁻³ m, will be considered. Energy dissipation at all scales must be accounted for. Dissipations at the smaller scale must not only be included but they must abide by the same physical and mathematical interpretation, in order to avoid inconsistencies when making connections with those at the larger scale where dissipations are eminent.Three fundamental Problems I, II, and III are stated. They correspond to the commonly used service conditions. Reference is made to a Representative Tip (RT), the location where energy absorption and dissipation takes place. The RT can be a crack tip or a particle. At the larger size scales, RT can refer to a region. Scale shifting of results from the very small to the very large is needed to identify the benefit of using multiscale materials.