非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程■,结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的v_(k+i/s)(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算e~(Ht),而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。     

非线性度对修正双步长显式法及常用逐步积分法的影响

为了掌握非线性度对逐步积分法的影响,研究了几种积分算法在不同非线性度振动系统中的响应。通过3个典型非线性算例,对修正双步长显式法、蛙跳式中心差分法、Newmark法、广义α法和精细积分法的计算精度和稳定性能等进行了比较。结果表明:非线性度对广义α法、精细积分法和Newmark法的稳定性有影响;高非线性度对Newmark法的计算稳定性影响最大;时间步长越小,算法精度和计算量越高;相同小步长情况下,精细积分法的精度最高,而修正双步长显式法的计算量最小;在时间步长较大时,低非线性度会引起精细积分法不稳定,修正双步长显式法的精度最高,修正双步长显式法在非线性系统中具有很强的鲁棒性。     

结构非线性动力分析混合时间积分并行算法

综合隐式和显式时间积分技术,对结构非线性动力反应分析提出一种并行混合时间积分算法.该算法采用区域分解技术.将并发性引入到算法中,即利用显式时间积分技术进行界面节点积分而利用隐式算法求解局部子区域.为实现并行混合时间积分算法,设计了灵活的并行数据信息流.编写了该算法的程序,在工作站机群实现了数值算例,验证了算法的精度和性能.计算结果表明该算法具有良好的并行性能,优于隐式算法.     

离散精细时程积分的自适应求积算法研究

基于Hamilton体系下的精细时程积分方法,通过对载荷项进行离散,应用中值法使载荷项在时间步长内为常值,从而将非齐次动力方程转化为齐次动力方程,避免了矩阵的求逆运算;基于积分区间逐次半分的思想实现了任意时间步长的自适应求积。数值算例结果表明:在同等时间步长的非齐次系统中,精细时程积分的最大误差为中心差分法的2.8%,为Newmark法的2.2%,最大求解误差仅为0.029%。这充分说明了本文的离散精细时程积分的自适应求积算法具有很好的收敛性。     

????????????????????????????

合隐式和显式时间积分技术,对结构非线性动力反应分析提出一种并行混合时间积分算 法. 该算法采用区域分解技术. 将并发性引入到算法中,即利用显式时间积分技术进行界面 节点积分而利用隐式算法求解局部子区域. 为实现并行混合时间积分算法,设计了灵活的并 行数据信息流. 编写了该算法的程序,在工作站机群实现了数值算例,验证了算法的精度和 性能. 计算结果表明该算法具有良好的并行性能,优于隐式算法.     

解对流方程的子域精细积分并行算法

基于子域精细积分的思想,针对对流方程初边值问题,首先提出了含参数a＞0的一族三层显格式和一族二层隐格式,它们的局部截断误差分别为O(a△t+△t2+△x2)和O(α△t+△t+△x2).当参数a≥(In△t-ln△x)/2△t时三层显格式是稳定的,而二层隐格式则对所有的参数α＞0都是无条件稳定的.然后,以二层隐格式为基础,设计了一种交替分组显武迭代(AGEI)方法,并证明了该迭代过程的收敛性.由于三层显格式和AGEI方法的整个计算过程都是显式的,所以非常适合于并行计算.文末的数值算例表明,上述方法具有很高的精确度和良好的实用性.     

RTM充模过程数值模拟的隐式有限元算法

建立了基于欧拉方法描述树脂传递模塑(RTM)工艺充模过程的基本数学方程,并采用有限元隐式时间积分方法对基本方程进行了数值求解．编制了基于隐式有限元算法及传统有限元控制体算法的程序,通过具体算例比较了这两种算法的优缺点．与传统的有限元控制体法相比,该文提出的隐式有限元算法能节省计算时间,特别适合于单元、节点数目多的情况．隐式有限元算法是一种纯有限元方法,不需要使用控制体积技术,采用该算法计算出的流动前沿与时间步长无关。     

瞬态热传导问题的精细积分-双重互易边界元法

采用双重互易边界元法结合精细积分法求解二维含热源的瞬态热传导问题。针对边界积分方程中热源项和温度关于时间导数项引起的域积分,采用双重互易法处理,将域积分转换为边界积分。采用边界元法将边界积分方程离散后,得到关于时间的微分方程组,并利用精细积分法处理其中的指数型矩阵;对于微分方程组中由边界条件和热源项引起的非齐次项,采用解析的方法计算。为了比较精细积分-双重互易边界元法的计算效果,同时使用有限差分法计算温度对时间的导数项。通过数值算例验证了本文方法的有效性和精确性。计算结果表明:时间步长对于精细积分-双重互易边界元法的结果影响较小,而有限差分法对时间步长比较敏感且只在时间步长选取较小时有效;当选取较大时间步长时,精细积分-双重互易边界元法依然具有良好的计算精度。     

高黏性流动有限元模拟的迭代稳定分步算法

分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中,已有的显式和半隐式分步算法由于 其显式本质,必须采用小时间步长计算,不但降低了计算效率,同时也常与为使分步算法稳 分步算法已被广泛应用于数值求解不可压缩N-S方程. Guermond等认为时间步长必须大于 某个临界值方能使算法稳定. 然而在高黏性流动模拟中,已有的显式和半隐式分步算法由于 其显式本质,必须采用小时间步长计算,不但降低了计算效率,同时也常与为使分步算法稳 定必须满足的最小时间步长要求冲突. 本文目的是构造一种含迭代格式的分步算法,它能在 保证精度的前提下大幅度地增大时间步长. 方腔流和平面Poisseuille流数值计算结果证实 了此特点,该方法被有效应用于充填流动过程的数值模拟.     

近似平衡多项式加速度动力显式算法

利用已知初始时刻的信息，建立一种可以取到任意阶高精度的多项式加速度单步隐式算法。在该隐式方法中，待采解方程纽系数矩阵中质量阵的系数远远大于阻尼阵和剐度阵的系数，略去非对角阻尼阵和非对角刚度阵对方程组的影响，得到一种近似平衡多项式加速度动力显式计算方法。此方法的精度主要由加速度多项式插值的项数、步长、质量阵的每件数、质量刚度比（质量阵和刚度阵的范数之比）决定。在此基础上给出了这种算法的通式，进行了精度分析，结果表明：如果时间步长h足够短，n次加速度近似平衡动力显式算法的精度可以达到O（hn＋1）。算例采用5次加速度近似平衡显式算法，计算结果的精确性证明了本算法的可行性。     

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.