《数值流形法》序

得知郑宏教授的《数值流形法》即将脱稿,激动之情难以言表.我完全同意郑教授将流形法的特点归结为“升阶”和“切割”.受郑教授的委托,我愿在此将建立流形法的历程做一简单的回忆.其实早在1985年———也就是在我完成不连续变形分析法(DDA)之前,“升阶”在我心中已完全成熟.“升阶”似乎能取得一致收敛,弥补有限元只能按能量收敛的缺憾.我没有急于发表是因为一直没能想到令我满意的解决“不连续”的办法.起初曾想过在DDA块体上引入高阶多项式,或在块体上再划网格,也就是用分片多项式来提高DDA块体的变形精度.但我对这些方案都不甚满意,因为实际计算时多项式的阶数会受到限制,引入网格又会招致有限元既有的毛病——解的精度受控于网格质量.     

三维位势问题的梯度边界积分方程的新解法

三维位势问题的边界元分析中,关于坐标变量的边界位势梯度的计算是一个困难的问题. 已有一些方法着手解决这个问题,然而,这些方法需要复杂的理论推导和大量的数值计算. 本文提出求解一般边界位势梯度边界积分方程的辅助边值问题法. 该方法构造了与原边界值问题具有相同解域的辅助边值问题,该辅助边值问题具有已知解,因此通过求解此辅助边值问题,可获得梯度边界积分方程对应的系统矩阵,然后将此系统矩阵应用于求解原边值问题,求解过程非常简单,只需求解一个线性系统即可获得原边值问题的解. 值得注意的是,在求解原边值问题时,不再需要重新计算系统矩阵,因此辅助边值问题法的效率并不很差. 辅助边值问题法避免了强奇异积分的计算,具有数学理论简单、程序设计容易、计算精度高等优点,为坐标变量梯度边界积分方程的求解提供了一个新的途径. 3个标准的数值算例验证了方法的有效性.      

基于波长调制-直接吸收光谱方法的CO分子1567 nm处谱线参数高精度标定

直接吸收光谱(DAS)可直接测量分子吸收率函数,并通过拟合吸收率函数确定待测气体参数.波长调制-直接吸收光谱(WM-DAS)在DAS基础上,结合了波长调制光谱(WMS)中谐波分析思想,利用傅里叶变换复现吸收率函数,可有效提高吸收率函数的测量精度.本文利用WM-DAS方法结合长光程气体吸收池,在室温低压条件下,对CO分子1567 nm处R5-R11近红外弱吸收谱线吸收率函数进行了精确复现,其拟合残差标准差低至3×10^-5,随后根据测得的吸收率函数对谱线的碰撞展宽、Dicke收敛以及速度依赖的碰撞展宽系数等光谱参数进行了高精度标定,并将其与高灵敏度的连续波腔衰荡光谱(CW-CRDS)测量结果进行了比较,实验结果表明该方法与CW-CRDS测量结果具有高度一致性,更具有系统简单、测量速度快、对环境要求低等优点.     

Mass and isotopic concentrations of water-insoluble refractory carbon in total suspended particulates at Mt.Waliguan Observatory(China)

Mass concentration and isotopic values δ¹³C and ¹⁴C are presented for the water-insoluble refractory carbon（WIRC） component of total suspended particulates（TSP）,collected weekly during 2003,as well as from October 2005 to May 2006 at the WMO-GAW Mt.Waliguan（WLG） site.The overall average WIRC mass concentration was（1183±120）ng/m³（n = 79）,while seasonal averages were 2081 ±1707（spring）,454±205（summer）,650±411（autumn）,and 1019±703（winter） ng/m3.Seasonal variations in WIRC mass concentrations were consistent with black carbon measurements from an aethalometer,although WIRC concentrations were typically higher,especially in winter and spring.The δ¹³C PDB value（-25.3 ± 0.8）%.determined for WIRC suggests that its sources are C₃ biomass or fossil fuel combustion.No seasonal change in δ¹³C PDB was evident.The average percent Modern Carbon（pMC） for ¹⁴C in WIRC for winter and spring was（67.2 ± 7.7）%（n = 29）.Lower pMC values were associated with air masses transported from the area east of WLG,while higher pMC values were associated with air masses from the Tibetan Plateau,southwest of WLG.Elevated pMC values with abnormally high mass concentrations of TSP and WIRC were measured during a dust storm event.     

北京工业大学实验班大学物理课后沙龙

<正>为了把"教学做"融为一体,北京工业大学实验班大学物理课开展了"课后沙龙"活动。学生在第四节下课后因为食堂拥挤不急于离开教室,教师因势利导,利用课后10-20分钟时间,请两位同学给大家介绍他所感兴趣的科技知识或对生活的感悟等等。每周一次,一次两人,已经坚持开展了两届。对愿     

含硫氨基酸与本征及缺陷石墨烯表面的相互作用

半胱氨酸及蛋氨酸是人体的两种含硫氨基酸,在生物活性中发挥着巨大的作用.本研究采用密度泛函理论方法对以上两种氨基酸在本征及缺陷石墨烯表面的吸附机理进行了详细研究.主要考虑了两种吸附体系:半胱氨酸及蛋氨酸平躺在两种石墨烯表面;两种氨基酸垂直地放置于两种石墨烯表面,且含硫的基团靠近表面.研究结果表明,半胱氨酸及蛋氨酸初始构型对它们之间的相互作用有一定的影响.两种氨基酸平躺时有较大的吸附能.此外,吸附能的结果显示两种氨基酸可以更好的与缺陷石墨烯表面紧密结合.同时,蛋氨酸与本征及缺陷石墨烯相互作用均大于半胱氨酸与本征及缺陷石墨烯相互作用.模拟结果有望为含硫氨基酸的石墨烯传感器提供有用的指导.     

颗粒材料三维应力路径下的接触组构特性

颗粒材料的宏观应力变形特征与其微观接触力、组构等紧密相关.一般而言,强接触系统属于颗粒内部体系的传力结构,其对应的组构张量是影响宏观应力性质的重要因素.细观数值方法(如离散单元法)能够反映物理试验的基本规律,并且可以方便地提取宏微观数据来研究颗粒体系的应力变形机制.采用离散单元法(discrete element method,DEM)进行一系列等$p$等$b$应力路径下颗粒材料的真三轴试验,在此基础上研究了三维应力路径下颗粒材料的宏微观力学参数的演化过程、三维组构张量与应力张量多重联系以及强接触体系反映的宏观应力特征.研究表明:颗粒体系偏应力峰值状态和临界状态均存在与加载路径无关的宏微观特征;三维应力路径下组构张量与应力张量存在非共轴性,但其联合不变量演化过程表现出加载路径无关的特征;与弱接触系统的组构张量相比,强接触系统的组构张量更能反映宏观应力张量的特征;强弱接触体系的组构张量对颗粒体系宏观响应的贡献不同,其分界点存在一定取值范围,但采用平均接触力较为简单合理.      

基于双状态Χ~2检测和快速强跟踪AUKF的组合导航算法

针对GNSS/SINS/摄影测量组合导航中某个子系统发生故障时,整个导航系统易受到故障数据污染的问题,提出了一种基于快速强跟踪AUKF的双状态卡方(Χ~2)检测数据融合方法。首先,采用快速强跟踪AUKF算法进行滤波;然后,引入卡方检验通过检测UKF子滤波器输出的状态向量来定位故障参数;最后,采用强跟踪滤波准确跟踪状态矢量突变以增强系统鲁棒性,并根据自适应因子实时调整预测协方差阵以修正增益矩阵,使滤波结果不受异常信息的干扰。将提出的改进算法与常规算法分别应用于无人机着陆导航系统,结果显示:与传统UKF相比,提出的算法得到的位置误差减少了62.6%以上;与强跟踪UKF相比,导航误差也至少减小了32.6%。     

含参广义集值强向量平衡问题的稳定性

本文借助集合极限的性质和弱f-性假设证明了含参广义集值强向量平衡问题解集映射的下半连续性,其方法不同于最近文献(Zhao,2016和Meng,2018).此外,建立了含参广义集值强向量平衡问题解集连通性的充分条件,并举例验证了所得结果的正确性.本文得结果推广和改进了已有文献(Gong,2008,Xu,2009,Chen,2010,Xu,2013和Zhao,2013)中相应结果.