双重点最小二乘配点无网格法

配点类无网格法需要计算近似函数的二阶导数,因而在移动最小二乘(MLS)近似中至少要采用二次基函数。本文利用Voronoi图对双重点移动最小二乘近似法进行了改进,建立了基于Voronoi图的双重点移动最小二乘近似(VDG),并利用加权最小二乘法离散微分方程,导出了双重点最小二乘配点无网格法(MD GLS)。该方法将求解域用节点离散,并以节点为生成点建立Voronoi图,取Voronoi多边形的顶点为辅助点。近似函数及其二阶导数的计算过程可分解为两个步骤:首先用场函数节点值拟合辅助点处近似函数的一阶导数,再以辅助点处近似函数的一阶导数值拟合节点处近似函数的二阶导数。由于在每一步中只需计算MLS形函数及其一阶导数,这种近似方法需要较少的影响点和较小的影响域。同时借助于Voronoi结构的优良几何性质,可以快速地搜索影响点。研究表明,与基于MLS的加权最小二乘无网格法(MWLS)相比,这种方法可以显著提高计算效率,并且在精度和收敛性方面也有所改善。     

连续体结构屈曲约束的ICM方法拓扑优化

基于ICM(独立、连续、映射)方法解决具有屈曲约束的连续体拓扑优化问题。建立以结构重量为目标,以屈曲临界力为约束的拓扑优化模型;采用独立的连续拓扑变量,借助泰勒展式将目标函数作二阶近似展开;借助瑞利商、泰勒展式、过滤函数将约束化为近似显函数,避免了灵敏度的计算;将优化模型转化为对偶规划,并利用序列二次规划求解,减少了设计变量的数目,缩小了模型的求解规模。给出三个算例,结果表明:该方法可有效地解决屈曲约束的连续体拓扑优化问题,能够得到合理的拓扑结构,并有较高的计算效率。     

基于响应面方法的桁架截面敏度分析和优化

把响应面方法引入桁架截面优化中,将应力和位移约束近似表达为桁架截面倒变量的线性函数。为拟合响应面,基于中心复合和单纯形试验设计方法开发了中心对称和拟单纯形试验设计两种方法,既可保证约束近似精度,又降低了结构分析计算量。对于桁架结构重量目标函数,直接推出倒变量的二阶形式,以桁架重量最小为目标的优化问题构造为标准的二次规划模型。将响应面方法计算的位移对设计变量的敏度与莫尔积分方法的近似显式进行了对比。以MSC.Patran为平台的数值算例表明:结合响应面方法,应用序列二次规划对问题进行寻优,其收敛精度及稳定性都可获得保障。     

基于快速Fourier变换的可靠性分析方法

提出了一种用于非线性隐式功能函数可靠性分析的快速Fourier变换方法.首先利用高效多维数值积分方法或基于响应面法的卷积积分求解功能函数的分布函数的特征函数,也就是进行Fourier正变换;然后利用快速Fourier逆变换数值求解功能函数的概率密度函数,进而求得失效概率.给出的算例表明此法具有精度和计算效率较高的优点.     

改进的基于雅可比椭圆函数的随机平均法

改进了基于雅可比椭圆函数的随机平均法,用于预测高斯白噪声激励下硬弹簧及软弹簧系统的随机响应.引入包含雅可比椭圆正弦函数、余弦函数及delta函数的雅可比椭圆函数变换,导出关于响应幅值和相位的随机微分方程.应用随机平均原理,将响应幅值近似为Markov扩散过程,建立其平均的It随机微分方程.响应幅值的稳态概率密度由相应的简化Fokker-Planck-Kolmogorov方程解出,进而得到系统位移和速度的稳态概率密度.以Duffing-Van der Pol振子为例,研究了硬弹簧及软弹簧情形下的随机响应,通过与Monte-Carlo数值模拟结果比较证实了本文方法的可行性及精度.     

杆系结构几何优化中的广义变量近似方法

提出一种杆系结构几何优化的广义中间变量近似方法，首先引入一套广义中间变量，包括各杆件局部柔度特性、方向余弦及内力与位移关系矩阵各元素，结构响应与这些广义中间变量的关系比与设计变量本身成更好的线性关系。因此以广义中间变量做一介泰勒展开近似位移和内力。应力约束和屈曲约束由近似内力计算，近似问题由优化器在设计变量空间内求解。最后给出了几个算例，结果表明本文的方法。是十分有效的。     

杆系结构几何优化中的广义变量近似方法

提出一种杆系结构几何优化的广义中间变量近似方法。首先引入一套广义中间变量，包括各杆件局部柔度特性、方向余弦及内力与位移关系矩阵各元素，结构响应与这些广义中间变量的关系比与设计变量本身（即节点坐标和截面积）成更好的线性关系。因此，以广义中间变量做一阶泰勒展开近似位移和内力。应力约束和屈曲约束由近似内力计算。近似问题由优化器在设计变量空间内求解。最后给出了几个算例，结果表明了本文的方法是十分有效的。     

桁架拓扑优化的多点逼近遗传算法

提出一种基于多点逼近函数和遗传算法的桁架拓扑优化方法。该方法建立了包含连续尺寸和离散拓扑两类变量的优化模型，并通过构造多点逼近函数建立了结构优化问题的第一级序列显式近似，然后采用分层优化方法：在外层对拓扑变量采用遗传算法进行优化；在内层对尺寸变量通过可由对偶法求解的第二级序列近似问题进行优化。几个经典的桁架拓扑优化算例表明该方法能以较少的结构分析次数获得比较理想的概率意义上的最优解。     

含模糊随机变量的多模式系统广义失效概率计算的内插迭代线抽样方法

通过正则化基本变量的度量空间,定义了单个基本变量同时具有模糊和随机双重不确定性时的广义失效概率.在广义失效概率的计算中,模糊随机变量被等价变换为随机变量.从而使得广义失效概率的计算变换为随机失效概率的计算.当模糊随机变量的密度函数和隶属函数均为正态型时,推导了其等价概率密度函数的形式和参数.采用自适应线抽样方法对基本变量同时具有模糊和随机不确定性时的多模式广义失效概率进行了计算,并采用数值算例对自适应线抽样广义失效概率计算方法的效率和精度进行了验证.算例分析表明该方法的计算结果是合理的,并且由于自适应线抽样法具有较高的效率和精度,因而所提方法具有一定的工程意义.     

结构可靠性分析高效自适应重要抽样方法

提出了一种基于自适应Metropolis算法和快速高斯变换技术的结构可靠性分析高效自适应重要抽样方法. 该方法首先利用自适应Metropolis算法高效生成结构失效域样本, 然后运用自适应宽核密度估计方法构造重要抽样密度函数, 最后采用快速高斯变换加速重要抽样过程中核函数的计算. 与传统方法相比, 自适应Metropolis算法能够在相同计算量下提供更多结构失效域信息从而改善计算精度, 即为求得给定精度问题的解, 可有效减少样本生成过程中的结构分析次数, 提高方法的计算效率; 快速高斯(Gauss)变换大幅降低核密度估计的计算复杂度从而大幅缩减重要抽样的计算耗时. 通过数值算例可以看出该方法具有较高的计算精度和效率.      

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.