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

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

基于三维荧光光谱法和PARAFAC对多环芳烃定性定量分析

三维荧光光谱法在研究多环芳烃(PAHs)类物质的荧光信息时起到了重要作用。多环芳烃类物质具有致癌性,难降解性,多由尾气排放,垃圾焚烧产生,危害着人类健康及环境,因此人们不断探索对多环芳烃检测的方法。实验选取多环芳烃中的苊和萘作为检测物质,利用FLS920荧光光谱仪,为避免荧光光谱仪本身产生的瑞利散射影响,设置起始的发射波长滞后激发波长40 nm,设置扫描的激发波长(λex)范围为:200~370 nm,发射波长(λem)范围为:240~390 nm,对多环芳烃进行荧光扫描获取荧光数据,采用三维荧光光谱技术结合平行因子算法对混合溶液中的苊和萘进行定性定量分析。实验选用的苊和萘均购于阿拉丁试剂官网,配制浓度为10 mg·L-1的一级储备液,再将一级储备液稀释,得到苊和萘浓度为0.5,1,1.5,2,2.5,3,3.5,4和4.5 mg·L-1的二级储备液,并将苊和萘进行混合。在进行光谱分析前需要对苊和萘的光谱进行预处理,采用空白扣除法扣除拉曼散射的影响,并采用集合经验模态分解(EEMD)消除干扰噪声。实验测得苊存在两个波峰,位于λex=298 nm,λem=324/338 nm处,萘存在一个波峰,位于λex=280 nm,λem=322 nm处。选用的PARAFAC算法对组分数的的选择很敏感,因此采用核一致诊断法预估组分数,估计值2和3的核一致值都在60%以上,分别对混合样品进行了2因子和3因子的PARAFAC分解,将分解后得到的激发发射光谱数据和各组分浓度数据进行归一化处理,并绘制光谱图,与归一化处理后的真实的激发发射光谱图和各组分浓度图进行对比。同时将PARAFAC得到的混合样本的预测浓度,通过计算回收率(R)和均方根误差(RMSEP)来判定定量分析的准确度。选择2因子时,各混合样品中苊和萘拟合度为95.7%和96.7%,平均回收率分别为101.8%和98.9%,均方根误差分别为0.0187和0.0316;选择3因子时,各混合样品中苊和萘拟合度为95.3%和95.8%,平均回收率分别为97%和102.5%,均方根误差分别为0.033和0.116,由三项指标可得选用2因子进行定性定量分析的效果明显好于选用3因子。分析实验结果表明,基于三维荧光光谱法和PARAFAC算法对混合样品进行定性定量分析,能够有效的判定混合样品的类别,同时能够成功的预测出混合样品的浓度。     

核工程钢板混凝土墙防撞击贯穿实用计算方法

针对带对拉钢筋的双钢板混凝土墙，研究了双钢板混凝土墙防贯穿计算方法，建立了基于能量法的防贯穿计算公式。在已知弹体和钢板混凝土墙体材料与几何相关参数时，利用防贯穿实用计算公式可对带对拉钢筋的核工程双钢板混凝土墙体贯穿速度以及弹体剩余速度进行计算，避免了复杂的结构抗撞击反应动力时程数值分析。为验证公式的可靠性，将公式计算结果与已有刚性弹体撞击双钢板混凝土墙实验结果及其动力有限元计算结果进行对比，结果表明:防贯穿实用计算公式可以正确判断双钢板混凝土墙的贯穿状态，实用计算公式给出的弹体剩余速度与实验结果符合良好。为进一步验证公式的适用范围，将公式计算结果与共10个工况的飞机发动机撞击双钢板混凝土墙的有限元计算结果进行了对比分析，结果表明:除1个工况计算结果偏差略超10%外，其余工况的偏差均在10%以内，说明该计算方法合理可行。     

基于CEEMD的重大事件对香港住宅价格影响的实证分析

应用互补集成经验模态分解（CEEMD）方法对香港1997―2018年的住宅价格月度数据进行了分解，将经过重构后的数据分成高频序列、低频序列与残差项。将BP多断点检测应用于低频序列，并结合样本时段内的重大事件进行实证分析。结果表明：1997年亚洲金融风暴对房价的影响大于2008年金融危机；外部经济体的救市政策间接地影响香港房价；在经济不景气的大环境下“孙九招”政策没有立即见效；资本投资者入境计划、住房供给调整与按揭贷款调整对房价的影响较为显著；税收调整对房价影响不显著、对交易量影响显著；SARS爆发使住宅价格下降约1%。     

形状记忆合金椭圆孔口问题的有限单元法

基于Lagoudas形状记忆合金(SMA)三维本构模型,假设材料为各向同性,推导了SMA平面应力状态的增量型本构方程,继而编写了ABAQUS用户自定义材料(UMAT)子程序,研究了在双向拉伸情况下,外载荷、温度、椭圆孔口长短轴之比对超弹性SMA椭圆孔口板中应力诱发马氏体相变区的影响。数值结果表明:应力诱发马氏体相变首先发生在椭圆孔口长轴端点部位,在外加载荷作用下逐渐扩展到板内,并由内向外形成马氏体相区、相变混合区和奥氏体相区;SMA板内应力诱发马氏体完全相变区面积与施加外载荷成正相关,与温度成负相关;随着椭圆孔口长短轴之比增大,SMA板内应力诱发马氏体完全相变区面积呈现出先减小后增大的趋势;拉应力差值相同时,相较于拉应力沿椭圆孔口长轴方向较大的情况,当拉应力沿椭圆孔口短轴方向较大时,SMA板内完全相变区面积较大,椭圆孔口周边应力集中现象更明显。     

基于图像的二维剪纸自动生成方法

中国剪纸的设计极具挑战性, 要求画面简洁、直观, 还需要表达特定的文化内涵, 且整张剪纸须整体连通。提出了一种基于图像的二维剪纸自动生成方法, 能够将任意数码照片自动转化为剪纸图形。首先，利用图像分割方法建立区域连接图; 接着, 基于该连接图对颜色、边界对比度和区域连通性进行数学建模, 并获得优化目标函数; 最后, 通过模拟退火算法求解目标方程, 自动生成保持图像内容的剪纸图形。还开发了连通性后处理和区域指定等用户交互工具, 允许用户在自动生成的剪纸图形中方便地加入个人设计。实验表明, 所生成的剪纸图形画面简洁、整体连通。 本方法在降低剪纸设计难度的同时还可满足个性化的设计需求, 有助于传播和传承我国的民间剪纸艺术。     

虚拟震源测深法及其在INDEPTH-3台阵中的应用

基于中美合作项目INDEPTH第3期在青藏高原布设的台站，使用虚拟震源测深法研究青藏高原中部的地壳厚度。结果显示，拉萨地体和羌塘地体的地壳结构存在巨大差异。拉萨地体的地壳厚度大约为57 km，与艾里均衡说预测的地壳厚度基本一致，说明拉萨地体的地壳结构比较简单。羌塘地体的地壳厚度为60~75 km，向北有增厚趋势，明显较艾里均衡说预测的地壳厚，说明羌塘地体地壳结构比较复杂，原因有可能是羌塘地体下存在高温流体和低速带，或者与印度板块岩石圈在班公湖-怒江缝合带以北向下俯冲有关。     

紫外吸收光谱法和电化学法研究原儿茶醛与牛血清蛋白的相互作用及共存金属离子的影响

利用紫外可见吸收光谱法和电化学法研究了中药小分子原儿茶醛与牛血清蛋白的相互作用,考察了温度和共存金属离子Cu~(2+),Ni~(2+),Zn~(2+)对二者结合常数的影响情况,计算了原儿茶醛与牛血清蛋白作用过程的ΔG,ΔH,ΔS。结果显示,原儿茶醛与牛血清蛋白的结合常数随着温度的升高和金属离子的存在都有不同程度的降低,二者的结合是一个焓和熵共同驱动的自发过程,之间的作用力主要表现为静电引力。运用循环伏安法测试了原儿茶醛与牛血清蛋白相互作用的电化学行为,结果表明,原儿茶醛与牛血清蛋白结合位点数为1,两者之间形成一种电活性的超分子化合物。