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

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

????????????????????????м???

面板堆石坝的三维非线性有限元分析,对并行计算提出了较高的要求. 在建模阶段挖掘出模 型内在的并行性,由相邻坝段构成子结构,既方便灵活又能满足区域分解法的要求. 编制了 适用于机群的并行多波前法计算程序,并应用于西龙池下库沥青混凝土面板堆石坝的有限元 数值分析,其计算结果和串行完全相同,通过对加速比和并行效率进行分析,结果表明采用 并行多波前法求解面板堆石坝的三维非线性有限元问题有效且算法简单,易于实现,无需辅 助的任务分配程序.     

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

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

基于集群的并行有限元分析研究

利用面向对象的方法实现了基于集群的并行有限元分析.首先介绍构建高性能计算集群的系统分析新方法,并针对建成集群系统的特点,对并行消息传递库MPI进行了并行语义分析,进而对其主要的实现函数进行面向对象的重构而建立起OO消息传递库,以此为基础分别利用两种不同的策略实现了面向对象的并行有限元分析(OOParaFEA: Object-Oriented Parallel FEA),一种是通过对传统的基于域分解方法的并行PCG算法进行改造而在集群平台实现,另一种是对已有的基于子结构方法的串行有限元分析程序加以扩展,加入系统方程组并行求解器而达到有限元分析并行化的目的.多个分析算例表明,基于集群进行面向对象的并行有限元分析可以有效提高计算效率,为进一步的网络化CAD/CAE研究奠定良好基础.     

有限元瞬态响应分析软件的一种并行开发方法

瞬态响应分析是有限元动力分析的重要内容之一，而串行计算机上运行的有限元软件在解题规模和速度上都受到很大限制。为此，基于系统集成思想对串行有限元软件进行并行开发。分析得出了瞬态响应分析并行开发的重点——位移方程组的瞬态响应，给出了线性方程组并行求解的思路和实现方法。用一个实例系统的实现验证了上述开发思路，从而也为并行应用软件开发探索了一条新的途径；最后对并行求解程序进行了算例验证。     

二维自适应非结构网格DSMC并行算法研究

研究了二维自适应非结构网格DSMC并行算法实现的过程.首先提出了一类非结构网格自适应策略,有效降低了网格尺度对计算结果的影响,提高了流场的分辨率;然后基于PC-CLUSTER群机并行体系结构与消息传递库MPI并行环境,利用分区并行思想,设计了非结构网格DSMC并行算法,节约了计算时间.利用For-tran90的动态分配内存技术编制了通用计算程序;最后对过渡流域高超声绕流进行了数值模拟,计算结果初步验证了算法的可行性与有效性.     

基于PANDA的并行显式有限元程序开发

冲击响应数值模拟在军用与民用领域有广泛的应用需求,但现有串行分析程序难以满足求解规模和效率的需求,因此需要开发并行显式有限元程序.本文在分析了显式有限元串行基本算法的基础上,设计了相应的并行算法,采用模块化程序开发模式基于PANDA程序框架开发了并行显式有限元分析程序DynPack,并进行了算例验证.算例表明,DynP...     

一种裂纹扩展的有限元仿真分析方法及其实现

为分析含多裂纹体在任意载荷下的断裂过程问题,给出了一种基于有限元的仿真方法及程序.基于商业有限元软件的应力分析结果和合适的断裂力学准则,判断裂纹的扩展以及扩展角;通过对每一个裂纹尖端扩展影响区域进行单元重划分,实现了多裂纹任意扩展的几何描述有限元模型更新.由于只对裂纹扩展区局部进行有限单元更新,网格划分和前期网格结果到新网格的映射的工作量相对较小,典型的裂纹扩展分析实例表明本文方法和程序的有效性.     

浅谈并行应用程序开发

在并行机上实现一个并行程序,要先谋而后动,全局考虑。从4个方面描述了对程序进行并行分解以及并行程序设计步骤。这些方法和步骤从实际工作中总结出来,对并行程序设计者有一定的参考价值。     

基于GPU并行计算的浅水波运动数值模拟

利用有限体积法求解描述水流运动的二维浅水方程组,模拟洪水波运动传播过程,并通过GPU并行计算技术对程序进行加速,建立了浅水运动高效模拟方法。数值模拟结果表明,基于本文提出的GPU并行策略以及通用并行计算架构(CUDA)支持,能够实现相比CPU单核心最高112倍的加速比,为利用单机实现快速洪水预测以及防灾减灾决策提供有效支撑。此外,对基于GPU并行计算的浅水模拟计算精度进行了论证,并对并行性能优化进行了分析。利用所建模型模拟了溃坝洪水在三维障碍物间的传播过程。     

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.