9/11: WHAT DID WE KNOW AND WHEN DID WE KNOW IT?

In February 1998, Osama Bin Laden published a signed statement calling for a fatwa against the United States for its having ‘declared war against God’. As we now know, the fatwa resulted in the unprecedented attack of 9/11. The issue of whether or not 9/11 was in any way predictable culminated in the public debate between Richard Clarke, former CIA Director George Tenet and the White House. This paper examines whether there was any evidence of a structural change in the terrorism data at or after February 1998 but prior to June 2001, controlling for the possibility of other breaks in earlier periods. In doing so, we use the standard Bai–Perron procedure and our sequential importance sampling (SIS) Markov Chain Monte Carlo (MCMC) method for identifying an unknown number of breaks at unknown dates. We conclude that sophisticated statistical time‐series analysis would not have predicted 9/11.     

伪谱法模拟各向异性/非均匀介质中弹性波的传播

数值模拟一直是研究弹性波在复杂介质中传播问题的有效手段之一。各向异性问题和非均匀问题日益引起关注。伪谱法是近年来发展较快、建立在弹性波动方程基础上的一种数值模拟方法。它的主要特征是将波动方程中的空间微分变换成频域中的乘法运算 ,从而只需要较少的空间格点就可以得到较高的计算精度。介绍了伪谱法及其稳定性条件 ;编制了一个适合于小型工作站或者PC机的模拟地震波在二维各向异性均匀 /非均匀介质中传播的数值模拟程序 ,模拟了各向异性 /非均匀介质中弹性波传播理论 ,分析了裂纹以及各向异性界面对波传播的影响     

弯管稀相气力输送CFD-DEM法数值模拟

利用CFD-DEM方法对稀相气固两相流在带有弯管的气力输送管道内的流动特性数值模拟。通过CFD求解连续气相流场,使用DEM求解离散颗粒相运动及受力,为提高仿真速度,模型忽略了空隙度对气相的影响。仿真结果显示了颗粒在弯管内形成颗粒绳及在垂直管道内颗粒绳分散的过程,并获取了颗粒-颗粒、颗粒-壁面的碰撞信息。对比弯管上下两部分的碰撞情况,颗粒及壁面在弯管下半部分受到的撞击和磨损更严重。通过对气力输送各类工作参数的研究发现,其对管道内气固两相流的流动特性和碰撞情况有着不同程度的影响。气流速度对颗粒绳分散影响甚微,但对弯管内颗粒碰撞强度有明显影响。随着颗粒质量流量的增加,形成的颗粒绳更紧凑、分散速度更慢,并在弯管中形成了阻碍颗粒-壁面碰撞的防护层。弯径比的增加也能加强颗粒绳的紧凑度,减缓颗粒绳的分散速度。     

非高斯噪声下信号盲检测算法*

针对实际甚低频和超低频接收机不仅受非高斯噪声的影响,同时受到接收机内部和外部环境中高斯噪声影响的问题,对噪声采用高斯尺度混合分布和高斯分布的混合模型建模,根据混合模型的性质,设计了一种基于马尔可夫链蒙特卡罗方法的信号盲检测算法。盲检测算法在贝叶斯层次模型下,采用Gibbs抽样和M-H抽样更新参数,同步检测信道衰落系数、噪声模型参数和信号。算法迭代效率快、精度高。通过与最优检测性能比较,盲检测算法性能优异,对甚低频和超低频信号接收具有重要的现实意义。     

动能拦截器姿控固体小火箭点火算法设计

使用固体姿控小火箭是实现动能拦截器快响应和高精度姿态控制的最佳方案之一。针对一种新型动能拦截器姿控小火箭布局,提出了点火组合混合搜索算法。描述了动能拦截器姿控小火箭的配置方案,分析了弹体自旋需求。设计了一种结合目标排序法和区间搜索法的点火组合混合搜索算法:当可用小火箭个数较少时,采用目标排序法;当可用小火箭个数较多时,采用区间搜索法。指令力矩近似仿真结果及姿态控制数值仿真结果表明:该算法能够有效地近似指令力矩,实现快速高精度的姿态跟踪。     

加权采样的序列蒙特卡罗滤波技术

在处理非线性滤波问题时,常用的几种通过模型近似或计算近似的滤波方法均在某种特定的情况下具有适用性.近年来随着计算机处理能力的快速发展,提出了序列蒙特卡罗滤波方法,因其在处理复杂的非线性和非高斯问题时表现出强大的潜力而引起广泛关注.详述了序列蒙特卡罗滤波算法的基本思想和原理,对其关键技术进行了归纳分析,并指出了该方法亟待解决的一些难点问题.     

基于蒙特卡洛和multi-Agents的舰艇编队维修保障能力优化求解

针对舰艇编队的维修保障能力评估,采用蒙特卡洛方法和智能计算技术进行了建模,并在分布交互仿真环境下进行了求解.通过将舰艇编队维修系统按维修单位的功能划分成多个子系统,每个子系统抽象为能进行局部优化的Agent,实现了基于Agent的建模.在此基础上,重点讨论了基于HLA的仿真环境,通过应用RTI的服务确保了求解过程的时空一致性.仿真试验结果表明:通过采用蒙特卡洛方法、智能技术及分布交互仿真技术相结合的方法,能有效地对舰艇编队的维修保障能力进行评估.     

加权随机模拟的时变可靠性分析方法

针对结构时变可靠性的随机模拟分析方法计算代价大的问题,在极值方法的基础上提出基于加权随机模拟的时变可靠性分析策略。时变可靠性分析需要计算在不同时间处的失效概率,通常需要进行多次可靠性分析,计算代价巨大。所提方法通过对常规静态可靠性的随机模拟方法进行改进拓展,运用加权策略分别发展了加权蒙特卡洛法和加权重要抽样法,使之能够高效分析计算时变可靠性问题。所提方法仅需一次常规可靠性分析模拟,即可得到时变失效概率函数估计。采用管状悬臂梁和十杆桁架两个算例进行验证。结果表明,基于加权思想的分析方法在能确保精确度的前提下能够大幅度减小计算量,提高计算效率。     

蒙特卡洛剖分投影法的爆震弹破片平均比动能评估

为了评估爆震弹封装壳体破片致死、致伤半径,并削弱其平均比动能,基于蒙特卡洛剖分投影法建立破片平均比动能计算模型。通过LS-DYNA及自编程联合仿真方法对爆震弹自然破片、半预制破片比动能计算模型进行仿真求解,获得了全破片全时域的质量分布、初始速度、垂直靶分布、平均比动能阈值及安全半径等指标。结果表明:在相同装药参数下,半预制破片相比于自然破片,质量分布、初始速度阈值更低且更集中;在小于等于2.5 m范围内半预制破片平均比动能阈值比自然破片更低,大于2.5 m后半预制破片平均比动能阈值比自然破片更高。相比于自然破片,采用半预制破片可显著减小爆震弹破片的致死半径,但并不能明显增大其安全半径。     

考虑黏弹性材料随机性的被动约束层阻尼梁动力学分析

采用传递函数方法研究了阻尼层黏弹性材料随机性对被动约束层阻尼(PCLD)梁动力学特性的影响。由Hamilton原理建立了PCLD梁六阶运动微分方程,通过引入状态向量,建立了系统的状态空间方程,利用传递函数方法得到了梁的固有频率和损耗因子。以黏弹性材料分数导数模型中的参数作为基本的随机变量,并假设其服从正态分布,使用Monte Carlo直接抽样法考察了材料模型参数的随机性对结构固有频率和模态损耗因子的影响。计算结果表明黏弹性材料参数的随机性对梁动力学特性的变异系数影响较大,模态损耗因子的变异系数最大值是材料参数变异系数的4.5倍。