加强型Bi-2212超导股线有限元拉伸性能分析

为探究加强型Bi-2212超导股线的力学性能以及不锈钢包带对超导股线力学性能的影响，运用Ansys软件建立加强型Bi-2212超导股线的三维模型，采用有限元分析方法，对其机械结构进行仿真，模拟高温高压热处理后轴向受拉形变量和应力应变分布情况。仿真结果显示，在同样施加100 N拉力的情况下，与超导原线（Bi-2212超导线）相比，密绕Ni₈₀Cr₂₀包带后的加强型Bi-2212超导股线轴向受拉形变程度减小，约10%;应力应变均降低，约30%。采用Ni₈₀Cr₂₀作为外包带制成的加强型Bi-2212超导股线力学性能显著，是未来超导线发展方向之一。     

基于内聚力模型的快速荷载下CFRP-混凝土界面剪切剥离模拟

鉴于纤维增强复合材料(Fiber Reinforced Polymer, FRP)已在混凝土结构加固中普遍应用, 能有效提高结构的承载和变形能力. 为揭示动载作用下FRP-混凝土界面的剥离机理, 本文基于ABAQUS平台, 采用内聚力模型模拟CFRP-混凝土界面层, 实现了快速荷载下CFRP-混凝土界面剥离的高效模拟. 结果表明: CFRP表面应变在加载过程中由加载端向自由端方向传递, 随着加载速率的提高, 界面承载能力也随之提高; 界面峰值剪应力也存在显著的应变率效应; 模拟结果与试验结果基本吻合, 说明了新方法的有效性.     

梁柱简化模型在自立塔塔线体系动态响应中的适用性研究

对输电塔进行合理简化可以提高塔线体系动力学仿真的效率。本文给出自立塔梁柱简化模型的计算方法，并提出利用梁柱简化模型计算方法建立自立塔塔线体系整体模型，同时采用桁梁混合模型建立精细化塔线体系整体模型，对两种模型塔线体系静力特性及振型和固有频率等动力特性进行对比分析。以脱冰工况为例，采用生死单元技术将施加在输电线节点上的集中质量单元杀死来模拟脱冰，实现对塔线体系动力学响应的有限元模拟，研究塔线体系简化模型在动态响应中的适用性。结果表明，两种模型弯曲变形误差小，低阶的振型相同，固有频率值误差小，动力特性基本相同；脱冰工况下，自立塔节点位移和塔材内力时程曲线一致，在提高计算效率的情况下，能有效保证计算精度。     

基于银行下侧风险规避的第三方B2B平台激励策略研究

以电子仓单融资为例, 基于银行下侧风险规避角度, 研究联合授信和委托授信下当第三方B2B平台存在行为隐匿的道德风险时, 银行对B2B平台的激励策略设计问题。研究发现:B2B平台的最优努力水平随收益分配比例、回购比例的增大而减小, 随质押率、贷款利率、产品采购量、损失补偿比例的增大而增大;同时银行为规避违约风险, 需设置质押率、贷款利率和贷款额上限及回购比例下限, 并且银行最优收益分配比例与损失补偿比例、最优损失补偿比例与贷款损失率均成正相关关系。此外, 随着B2B平台工作效率的提高, 联合授信下最优收益分配比例将减小, 最优损失补偿比例将增大, 最终近似于委托授信下的最优损失补偿比例。最后给出数值分析。     

电子结构计算的数值方法与理论

第一原理电子结构计算已成为探索与研究物质机理、理解与预测材料性质的重要手段和工具.虽然第一原理电子结构计算取得了巨大的成功，但是如何利用高性能计算机又快又好地计算大规模体系，如何从数学角度理解电子结构模型的合理性与计算的可靠性和有效性，依然充满各种挑战.基于密度泛函理论的第一原理电子结构计算的核心数学模型为Kohn-Sham方程或相应的Kohn-Sham能量泛函极小问题.近年来，人们分别从非线性算子特征值问题的高效离散及Kohn-Sham能量泛函极小问题的最优化方法设计两个方面对电子结构计算的高效算法设计及分析展开了诸多研究.本文重点介绍我们小组在电子结构计算的方法与理论方面的一些进展，同时简单介绍该领域存在的困难与挑战.     

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

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

基于SWT方法的钢绞线索微动疲劳特性分析

为得到钢绞线索丝间接触区的应力场分布并预测微动疲劳裂纹萌生位置和微动疲劳寿命，本文利用参数化方法建立了精细化的钢绞线拉索有限元模型，包括整索模型和不同层丝间接触区域的局部精细化子模型．分析了钢绞线索在两种交变荷载工况下的应力场变化情况，并基于多轴疲劳SWT(Smith-Watson-Topper)临界平面法进行了疲劳特性分析和疲劳寿命预测．主要结论如下：钢绞线索内接触区边缘处的微动幅值较大，中心处几乎没有相对滑动，微动疲劳的初始裂纹萌生点位于接触区域边缘；经不同区域子模型分析比较，在轴向循环荷载作用下，外层钢丝的接触区域比内层钢丝更易发生微动疲劳损伤；在横向位移循环荷载作用下，同层钢丝因位置角度不同而产生了较大的疲劳特性差异，且相比轴向循环拉伸，该工况下最不利单丝的微动疲劳寿命更低；与非接触区域相比，接触区的疲劳寿命大幅降低，微动现象对钢绞线索的抗疲劳性能有明显降低作用．     

《数值流形法》序

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

平面波输入下海水-海床-结构动力相互作用分析

海洋工程结构的地震反应分析是保证海洋工程结构地震安全的重要环节.由于其所处的复杂环境, 该问题涉及到流固耦合和土-结相互作用.本文基于海水、饱和海床、基岩流固耦合统一计算框架,采用Davidenkov模型和修正的Masing准则考虑饱和海床的非线性,在脉冲SV波垂直入射下, 进行了海域场地和海洋工程结构的动力响应分析. 首先,对比分析了线性自由场和非线性自由场输入情形的海域场地非线性反应,结果表明线性自由场输入时反应不合理,自由场分析和场地分析应该采用相一致的本构模型. 然后,对比分析了海床分别为线性和非线性情形时,海域场地以及海水-海床-结构体系的反应特征. 与线性海床情形相比,非线性对海床反应的影响主要由如下两方面因素控制: 一方面,非线性导致饱和海床模量减小, 饱和海床与基岩间的波阻抗比减小,由基岩到饱和海床间的反射系数和透射系数增加, 导致反应增大; 另一方面,非线性导致阻尼加大, 使海床反应减小. 对于本文算例而言,阻尼对非线性海床结果的影响占主导作用.      

基于六面体网格的向量式有限元分析及应用

基于向量式有限元基本原理，给出了八节点六面体等参实体单元的基本公式，通过投影方式将空间曲面六面体转换为投影六面体，采用参考面的逆向运动求解节点纯变形，通过单元形函的虚功方程计算节点内力；针对坐标模式和内力积分模式等关键问题提出了有效的处理方案。编制了六面体实体单元的数值计算程序，并进行工程结构算例分析。结果表明，所编制程序可有效模拟实体结构的静力、动力及大变形大位移行为分析，验证了本文理论和程序的正确性和实用性。