椭球粒子声辐射力计算及分析

针对水下椭球粒子,以声散射理论为基础,采用分波序列的方法,建立了椭球粒子声辐射力的理论计算模型。进而根据声辐射力计算公式,以刚性椭球粒子和液体椭球粒子为例,计算并分析不同Bessel波束作用下椭球粒子的轴向声辐射力函数特征。数值仿真计算结果表明,对于刚性椭球粒子,扁平椭球粒子相对于细长椭球粒子更有助于激发负声辐射力；对于液体椭球粒子,细长椭球粒子相对于扁平椭球粒子更加容易产生负声辐射力；对于不同介质的椭球粒子,不同的入射波束激发的负声辐射力的效果也存在明显的差异。该结果为复杂的尺寸和介质粒子声操控技术提供了理论的可行性。     

Duffing系统的主--亚谐联合共振

以Duffing系统为研究对象,研究在多频激励下同时发生主共振和1/3次亚谐共振的动力学行为与稳定性.首先,通过多尺度法得到系统的近似解析解,利用数值方法检验近似程度,结果吻合良好,证明了求解过程和解析解的正确性.然后,从解析解中导出稳态响应的幅频方程和相频方程,从幅频曲线以及相频曲线中发现系统最多存在7个不同的周期解,这种多解现象可用于对系统状态进行切换.基于Lyapunov稳定性理论,得到联合共振定常解的稳定条件,利用该条件分析了系统的稳定性,并与Duffing系统的主共振和1/3次亚谐共振单独存在时比较.最后,通过数值方法分析了非线性项和外激励对系统动力学行为与稳定性的影响,发现了联合共振特有的现象:刚度软化时,非线性项不仅影响系统的响应幅值,同时还影响系统的多值性和稳定性;刚度硬化时,非线性项对系统的影响与单一频率下主共振和1/3次亚谐共振类似,仅影响系统的响应幅值.这些结果对Duffing系统动力学特性的研究具有重要意义.     

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

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

正电子冲击下氦电离三重微分截面的理论计算

我们发展了一种正电子碰撞原子电离的畸变波Born近似方法， 在这个方法中，正负电子偶素通道通过一个ab initio的光学势附加到入射粒子和靶的相互作用势上，且通道对电离作用被第一次被考虑在正电子碰撞原子电离的过程中. 应用这个方法计算了在50 eV入射能量范围氦的电离的三重微分截面，计算结果和实验数据很好的符合.     

2:1内共振条件下变转速预变形叶片的非线性动力学响应

在工程实际中,涡轮机叶片的转速在很多应用场景下不是一个定常值,比如发动机在启动、变速、停机等工况下,转子输入与输出功率失衡,伴随产生扭振,产生速度脉冲. 另外,由于服役环境、安装误差等因素会引起叶片在所难免的预变形. 本文主要研究预变形叶片,在变转速条件下的非线性动力学行为. 考虑叶片转速由一定常转速和一简谐变化的微小扰动叠加而成. 应用拉格朗日原理得到变转速叶片的动力学控制方程,并采用假设模态法将偏微分方程转为常微分方程,通过引入无量纲,使方程更具有一般性. 运用多尺度方法求解了该参激振动系统,得到了在 2:1 内共振情形下的平均方程,进而获得系统的稳态响应. 详细研究温度梯度、阻尼以及转速扰动幅值等系统参数对叶片动力学响应的影响规律,同时考察了立方项在 2:1 内共振下对方程的影响. 对原动力方程进行正向、反向扫频积分来观察其跳跃现象,并对解析解进行验证. 结果发现参数的变化对叶片均有不同程度影响,在 2:1 内共振下立方项对系统响应的影响很小,解析解与数值解吻合很好.      

基于局部加密纯无网格法非线性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波垂直入射下, 进行了海域场地和海洋工程结构的动力响应分析. 首先,对比分析了线性自由场和非线性自由场输入情形的海域场地非线性反应,结果表明线性自由场输入时反应不合理,自由场分析和场地分析应该采用相一致的本构模型. 然后,对比分析了海床分别为线性和非线性情形时,海域场地以及海水-海床-结构体系的反应特征. 与线性海床情形相比,非线性对海床反应的影响主要由如下两方面因素控制: 一方面,非线性导致饱和海床模量减小, 饱和海床与基岩间的波阻抗比减小,由基岩到饱和海床间的反射系数和透射系数增加, 导致反应增大; 另一方面,非线性导致阻尼加大, 使海床反应减小. 对于本文算例而言,阻尼对非线性海床结果的影响占主导作用.      

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

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

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

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