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

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

基于局部加密纯无网格法非线性Cahn-Hilliard方程的模拟

为数值求解描述不同物质间相位分离现象的高阶非线性Cahn-Hilliard(C-H)方程,发展了一种基于局部加密纯无网格有限点集法(local refinement finite pointset method,LR-FPM).其构造过程为:1)将C-H方程中四阶导数降阶为两个二阶导数,连续应用基于Taylor展开和加权最小二乘法的FPM离散空间导数;2)对区域进行局部加密和采用五次样条核函数以提高数值精度;3)局部线性方程组求解中准确施加含高阶导数Neumann边值条件.随后,运用LR-FPM求解有解析解的一维/二维C-H方程,分析粒子均匀分布/非均匀分布以及局部粒子加密情况的误差和收敛阶,展示了LR-FPM较网格类算法在非均匀布点情况下的优点.最后,采用LR-FPM对无解析解的一维/二维C-H方程进行了数值预测,并与有限差分结果相比较.数值结果表明,LR-FPM方法具有较高的数值精度和收敛阶,比有限差分法更易数值实现,能够准确展现不同类型材料间相位分离非线性扩散现象随时间的演化过程.     

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

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

悬链线竖向集中力作用下的构形分析和张力计算

目前,悬链线在竖向集中力和均布荷载共同作用下的构形分析和受力计算的理论仍不完善。针对这一问题,通过引入悬链线的几何约束方程、力平衡方程和超越方程,建立了竖向集中力与均布荷载共同作用下的非线性方程组。采用牛顿迭代法求解方程组,得到了悬链线的构形和受力情况。为了验证理论计算的正确性,进行了算例和试验验证。结果表明,算例的计算结果与文献结论保持一致,试验测得的构形和水平张力大小与理论计算的构形和水平张力大小吻合较好。本文的理论计算可以更加简单精确地计算出悬链线在竖向集中力和均布荷载共同作用下的构形和受力情况,为实际工程提供重要的理论指导。     

基于特征正交分解的一类瞬态非线性热传导问题的新型快速分析方法

提出了一种基于特征正交分解(POD)和有限元法的瞬态非线性热传导问题的模型降阶快速分析方法, 建立了导热系数随温度变化的一类瞬态非线性热传导问题有限元格式的POD降阶模型. 在隐式时间推进方法的基础上有效结合单元预转换方法和多级线性化方法发展了一种加速求解瞬态非线性热传导降阶模型的新型计算方法,并通过二维和三维算例验证了该方法的准确性和高效性. 研究结果表明: (1)降阶模型解的均方根误差在经过初始时段轻微的脉动后稳定于0.01%以下, 而其计算效率比有限元全阶模型提高2$\sim $3个数量级, 并且自由度数量(DOFs)愈大提高的幅度也愈加显著; (2)新型算法解决了常规算法在计算非线性降阶模型时加速性能差的问题, 即使是在DOFs比较小的时候也能够明显提高计算效率; (3)常数边界条件下得到的POD模态可以用来建立相同求解域在各种复杂时变边界条件下的瞬态非线性热传导降阶模型, 并对其传热过程和温度场进行快速准确的分析与预测, 具有很好的工程应用价值.      

一种新型介电弹性体仿生可调焦透镜的变焦分析

介电弹性体 (dielectric elastomer) 是电活性聚合物智能材料的一种,在外加电场作用下,可产生多种形式的响应.在驱动柔性透镜的变焦方面,相对于传统的机械操控变焦方法 显示出独特的优势.针对一款在电压激励下可高效调节焦距的介电弹性体仿人眼变焦透镜,该透镜由上下两层介电弹性薄膜和固定框架构成,并在封闭腔内充入盐水,上层薄膜涂覆环形柔性电极.在电压激励下,上层膜发生变形,由于盐水的体积保持恒定,引 起下层膜随之变形,使得透镜的焦距发生改变.采用 neo-Hookean 模型,利用变分原理导出了该透镜的控制方程、边界条件和连 续条件.利用打靶法求解了该非线性问题并高效地处理了非线性问题的界面连续条件. 理论分析结果与实验结果相吻合. 利用此模型开展了广泛的参数分析,研究表明,透镜的几何形状、初始焦距、介电弹性体薄膜的预拉伸率、涂覆的电极面积、材料的剪切模量等对透镜焦距的调节性能都有重要的影响.所建立的理论分析模型可为柔性仿生透镜的设计和参数优化提供有效的分析方法.      

缝合式C/SiC复合材料非线性本构关系及断裂行为研究

C/SiC复合材料具有高比强度、高比模量和优良的热稳定性能等一系列优点, 广泛应用于航空航天领域中. 裂纹扩展进而引起的脆性断裂是其主要失效形式之一, 因而材料的断裂性能分析对材料的结构设计和应用有重要的指导意义. 本文开展了缝合式C/SiC复合材料简单力学试验和断裂试验, 研究了材料在不同载荷下的力学响应及断裂特征. 基于缝合式C/SiC复合材料简单力学试验, 建立了材料宏观非线性损伤本构方程, 并模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 本构方程采用简单函数描述了材料在复杂应力状态下的非线性应力-应变曲线, 并考虑了反向加载过程中造成的裂纹闭合. 基于商业有限元软件ABAQUS, 通过编写UMAT子程序实现非线性损伤本构方程, 采用单个单元验证了建立的本构方程的有效性. 在此基础上, 采用线弹性损伤本构和非线性损伤本构分别模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 采用非线性损伤本构方程模拟的力-位移曲线结果与试验结果更为吻合, 非线性损伤本构预测的失效载荷与试验失效载荷更为接近, 验证了所建立的非线性损伤本构方程的准确性, 为C/SiC复合材料断裂行为的研究提供了借鉴, 为缝合式C/SiC复合材料结构的设计和应用提供了理论基础.      

Ba7AgGa5S15:一种新型混合金属硫化物的非线性光学性能研究

采用高温熔融-自发结晶法成功获得一种新型混合金属硫化物Ba7AgGa5S15.该化合物结晶于非中心对称的P31c空间群(No.159),晶胞参数为a=0.964 53(10) nm,c=1.805 9(4) nm,Z=2.其结构是由[Ga4 S10] T2超四面体与[AgS4]四面体共顶点连结形成的含有18元环孔道的三维网状框架,孤立的[Ga(2)S4]四面体填充在孔道中,Ba2+填充在该三维框架结构的空隙当中.第一性原理计算研究了该化合物的电子结构、态密度、双折射率、二阶非线性光学系数,以及倍频密度.结果 表明,该化合物具有大的光学带隙(3.76 eV),其带隙主要由S3p,Ba5d和Ga4s轨道决定;其d33方向上的倍频系数约为AgGaS2的0.4倍,主要倍频贡献来源于[AgS4]和[GaS4]四面体.该研究表明在Ag-Ga-S体系中引入Ba2+,形成的Ba7 AgGa5S15表现出比AgGaS2更宽的带隙,有利于产生高的激光损伤阈值(LDT).     

可变序结构下向量优化中的一个新非线性标量化函数及其应用

在具有可变序结构的一般拓扑向量空间中定义了一个新的非线性标量化函数,讨论了该函数的主要性质.同时作为应用,通过该函数构造出了一族半范数和一类赋范线性空间,并在最后建立了该非线性标量化函数和半范数的上、下半连续性结论.     

具有非线性传染率和时滞的随机SIQR计算机病毒模型的渐近行为

研究了一类具有非线性发生率和时滞的随机SIQR计算机病毒模型.首先证明了该系统具有唯一的全局正解,然后通过构造适当的Lyapunov函数并利用伊藤公式,分析了该模型的解在无病平衡点附近及地方病平衡点附近的渐近行为,最后通过数值模拟对随机系统解的渐近行为做了进一步的分析并给出了结论.