非光滑稀疏约束优化问题的最优性条件及算法

针对目标函数非光滑的稀疏约束优化问题,给出基本可行性和λ-平稳性两个必要最优性条件,利用所给出的必要最优性条件构造出稀疏次梯度投影算法.在理论上分析了算法的收敛性,证明了由该算法所产生序列的任意聚点都是λ-平稳点.最后,通过两个数值实例验证了算法的收敛性、有效性和优化能力.     

基于卷积神经网络的HL-2A装置H模辨识研究

基于HL-2A装置的放电实验数据，利用卷积神经网络和时间窗口算法开发了高约束(H)模时段的识别算法，得到了可靠的高成功率的高约束模时段识别结果。算法中，选取206次放电实验数据中等离子体储能及氘α通道信号作为双通道原始数据进行学习，得到一个深度为21层的二分类卷积神经网络。该网络模型经过其他474次放电数据的测试集检验，高约束模识别的正确率达到了98.17%。     

聚变堆相关的液态金属膜流初步实验研究

在磁约束核聚变堆的面对等离子部件设计中,液态金属锂膜流因具有带走杂质、保护面对等离子固壁等优点而被认为是优选方案之一. 然而,如何克服聚变堆中强磁场环境下产生的磁流体力学效应并形成大面积均匀铺展锂膜流动是目前亟需解决的问题.本文通过搭建室温液 态镓铟锡回路和高温液态锂回路,开展了两种不同特性的液态金属膜流实验, 并采用传统可视化方法获得了展向磁场存在时镓铟锡和锂在导电底板形成的液膜流动表面特征.实验结果 表明: 无磁场时,两种液态金属膜流流动表面波动特性与常规流体膜流均一致, 即随着流动雷诺数的增加表面波动变得更为混乱; 而展向磁场存在时,镓铟锡膜流表面波动变得更为规则, 且沿着磁场方向平行排列,表现为拟二维波动的特征; 而锂膜流却产生了明显的磁流体 力学阻力效应,表现为在流动方向局部产生锂滞留现象, 且滞留点随雷诺数增大向下游移动. 最后通过膜流受力分析,进一步阐述了锂膜流受到比镓铟锡膜流更为严重磁流体力学效应影响的原因.      

稀疏张量约束的低剂量CT图像重建

提出了一种稀疏张量约束重建算法,该方法利用非局部相似的先验信息,将CT图像分割成一系列图像块组;采用张量的多维低秩分解方法,将这一先验信息引入低剂量CT重建中,构造目标函数;通过重建图像更新和图像块组张量稀疏编码两个步骤,交替迭代求解目标函数。基于仿真数据和临床数据的实验结果验证了该算法的有效性,实验结果表明:与经典重建算法相比,所提算法在抑制噪声的同时,能更好地保持重建图像的细节,获得更高质量的图像。     

具有指数型保护障碍的带约束的动态保护基金的定价

动态基金保护可以确保投资者在投资期间内基金价格永不会低于某个障碍水平.为了减少发行人所承担的资产的下行风险,Han等(2016)提出了一个带约束的动态基金保护,产品约定只有当标的基金价值达到某个事先给定的上保护水平时才启动保护.受此启发,考虑了具有指数型上下保护水平的带约束的动态保护基金的定价.利用布朗运动的反射原理,给出了具有指数型保护水平的带约束的动态保护基金的定价公式.     

广角任意反射面速度干涉仪的光学性质研究

研究了广角任意反射面速度干涉仪(VISAR)的光学性质。阐述了广角VISAR的原理，指出广角诊断靶中椭球镜的作用是将靶丸内表面成虚像在靶丸中心附近。使用Zemax模拟了成像弯曲对动态干涉条纹形成的影响，提出使用异形光纤面板进行像面矫正。研究了工程误差对干涉仪成像的影响，若要取得良好的成像效果，椭球镜的位置偏差不得多于30 μm，倾角不得超过4°，长轴方向加工误差需小于0.1 μm，短轴方向误差需小于4 μm，镜面反射率应高于70%。讨论了广角VISAR光学研究的进一步发展方向如影响动态条纹的更多可能因素、像面矫正的其他方法、物与像面的光学对应等。     

多元统计分析中一类矩阵迹函数最小化问题的有效算法

研究来源于多元统计分析中的一类矩阵迹函数最小化问题$\min c+ tr(AX)+\sum\limits_{j=1}^{m}tr(B_j X C_jX^{T}),\ \ {\rm s. t.} \ X^TX=I_p,$其中$c$为常数, $A\in R^{p\times n}\ (n\geq p)$, $B_j\in R^{n\times n}, C_j\in R^{p\times p}$为给定系数矩阵. 数值实验表明已有的Majorization算法虽可行, 但收敛速度缓慢且精度不高. 本文从黎曼流形的角度重新研究该问题, 基于Stiefel流形的几何性质, 构造一类黎曼非单调共轭梯度迭代求解算法, 并给出算法收敛性分析.数值实验和数值比较验证所提出的算法对于问题模型是高效可行的.     

核电结构土-结相互作用分析分区混合计算方法

土-结构相互作用分析是核电结构进行抗震设计和安全评估的重要环节.在核电结构的土-结相互作用分析中,阻尼和非线性是影响结构反应的重要因素. 若采用频域分析,可以方便考虑阻尼,但需通过等效线性化来考虑非线性,不适合于强震作用下的土体非线性.若采用时域分析的逐步积分方法,适合于考虑非线性,但材料阻尼一般采用瑞利阻尼模型,除了紧靠指定阻尼比的少数几个振型外,其他振型的反应将受到瑞利阻尼模型所确定的大阻尼所抑制,造成地震反应与真实情形有较大差异.若采用时域分析的模态叠加法,可合理计入阻尼效应,但模态叠加法不能考虑非线性.因此,如何合理考虑阻尼和非线性是核电结构土-结相互作用分析需要关注的问题.基于此,本文提出一种模态叠加和时步积分结合的土-结相互作用分区算法.其中,出于安全性考虑,地震作用下核电主体结构一般不允许进入非线性,因此结构可采用模态叠加方法,以便合理考虑结构阻尼;土体和基础采用显式时步积分法,可考虑土体非线性;通过人工边界条件考虑无限域的影响 (辐射阻尼).通过简单算例对该方法进行了验证,并用于CAP1400核电结构的土-结相互作用分析中,对比分析了采用模态阻尼和瑞利阻尼时核电结构和场地反应的差异,结果表明结构阻尼模型对场地的反应影响不大,但对结构反应影响明显,在实际工程中应合理选取阻尼模型.      

大型偏心结构吊装欠约束问题研究

吊装施工过程中被吊模块的水平度是作业要求的重要指标,通常需要增加配重调平。传统有限元方法需要补充约束以消除单元刚体位移,且需要重复计算平衡方程来求解调平载荷,效率不高。将模块的运动分解为随动坐标系的整体运动以及相对该坐标系的弹性变形,可将欠约束问题化为多体系统的静平衡问题。基于虚功率原理推导了吊装平顺时刻的节点力平衡方程以及相应的切线刚度矩阵,并将配重表示为基础配重与载荷系数相乘的形式。通过对节点力平衡方程求导,得到一组以载荷系数为自变量的微分方程,通过求解微分方程并结合水平度判据,可快速搜寻满足水平度要求的载荷系数。数值算例表明,该方法在解决偏心模块吊装欠约束问题方面具有明显的优势,在确定配重载荷方面具有较快的速度和合理的精度。     

三维无网格Galerkin结构分析方法及其应用研究

针对无网格Galerkin法在三维复杂几何形状的结构分析中存在的刚度矩阵的稀疏存储实现难、六面体背景网格适应性差等问题,本文采用逐节点对法组装刚度矩阵,利用CSR格式存储刚度矩阵,提出了一种基于四面体背景积分的改进的三维无网格Galerkin法。通过采用罚函数法施加对称约束和给定位移值约束,并推导出了施加这两种位移约束的统一格式。利用所提算法完成了三维悬臂梁的计算,所得结果与其理论解相吻合;完成了轴流式风机轮毂的结构分析,得到的位移与应力分布结果与其有限元解相吻合。这表明本文所提方法能满足工程应用中的计算要求,并适用于具有复杂形状的几何模型分析。