A Simple Semi-Implicit Scheme for Partial Differential Equations with Obstacle Constraints

We develop a simple and efficient numerical scheme to solve a class of obstacle problems encountered in various applications. Mathematically, obstacle problems are usually formulated using nonlinear partial differential equations (PDE). To construct a computationally efficient scheme, we introduce a time derivative term and convert the PDE into a time-dependent problem. But due to its nonlinearity, the time step is in general chosen to satisfy a very restrictive stability condition. To relax such a time step constraint when solving a time dependent evolution equation, we decompose the nonlinear obstacle constraint in the PDE into a linear part and a nonlinear part and apply the semi-implicit technique. We take the linear part implicitly while treating the nonlinear part explicitly. Our method can be easily applied to solve the fractional obstacle problem and min curvature flow problem. The article will analyze the convergence of our proposed algorithm. Numerical experiments are given to demonstrate the efficiency of our algorithm.     

带有摩擦耗能元件的框架结构动力分析

本文提出粘性屈服模型来模拟摩擦耗能元件的力-速度关系,该模型是连续变化的,克服了库仑摩擦模型不连续导致数值计算复杂的缺点,在进行摩擦耗能体系的动力分析中,采用缩减自由度技术,并作适当的变换,则带有摩擦耗能元件体系的动力分析归结为求解微分代数方程,本文采用增量型Rosenbrock二级三阶半隐式Runge-Kutta法求解该方程,以考虑框架和支撑的材料和几何非线性。对带有摩擦耗能元件的钢框架进行了弹塑性动力分析,研究了支撑刚度与结构层刚度的比值、摩擦力的大小以及地震波类型等参数对体系的影响。     

基于MPS-FEM耦合方法对比研究刚性与弹性挡板对液舱晃荡的抑制作用

由流体冲击载荷引起的流固耦合问题广泛存在于船舶与海洋工程领域.例如:在特定激励频率下载液货舱内流体的非线性运动引起对舱壁的砰击作用,进而可能影响液舱围护系统的安全性.由于此类流固耦合问题通常涉及多学科知识,且流体自由面的变化具有强非线性特征,对研究人员带来较大挑战.考虑到Lagrange类方法在处理结构和流体自由面大变形问题上的优势,基于MPS-FEM耦合方法开发了流固耦合求解器.其中,采用MPS方法来数值模拟流体场瞬态变化,FEM方法来分析结构场的变形问题.此外,该求解器采用了弱耦合的方式来实现流体场和结构场之间的数据传递.为了验证该方法在处理流固耦合问题上的可靠性,首先数值研究了溃坝泄洪流与弹性挡板之间的流固耦合标准算例,数值结果与实验标准结果能够较好地吻合.此后,采用该求解器数值研究了带刚性挡板和弹性挡板的液舱晃荡问题,对比分析了多种激励频率下两种挡板对液舱内流体运动及舱壁上冲击压力的抑制效果.