Material point method applied to multiphase flows

The particle-in-cell method (PIC), especially the latest version of it, the material point method (MPM), has shown significant advantage over the pure Lagrangian method or the pure Eulerian method in numerical simulations of problems involving large deformations. It avoids the mesh distortion and tangling issues associated with Lagrangian methods and the advection errors associated with Eulerian methods. Its application to multiphase flows or multi-material deformations, however, encounters a numerical difficulty of satisfying continuity requirement due to the inconsistence of the interpolation schemes used for different phases. It is shown in Section 3 that current methods of enforcing this requirement either leads to erroneous results or can cause significant accumulation of errors. In the present paper, a different numerical method is introduced to ensure that the continuity requirement is satisfied with an error consistent with the discretization error and will not grow beyond that during the time advancement in the calculation. This method is independent of physical models. Its numerical implementation is quite similar to the common method used in Eulerian calculations of multiphase flows. Examples calculated using this method are presented.     

Probing the Contraction and Association of Polystyrene Chains in Semidilute Solution by Non‐Radiative Energy Transfer

The interchain contraction and association of polystyrene in semidilute solution near the Θ temperature are studies using NET. The results reveal that during a drop of temperature, polystyrene chains pass through different states: a smooth change of conformation at temperatures above the Θ temperature, an accelerated contraction and association, aggregation of polymer chains, and then further continuous contraction below the Θ temperature. A pronounced hysteresis is observed during repeated heating. A viscosity study of the same solution was compared with the NET results. It indicates that NET can effectively be used to study conformational changes like contraction and association of polymer chains in a semidilute solution.

     

Isotropic polynomial invariants of Hall tensor

The Hall tensor emerges from the study of the Hall effect, an important magnetic effect observed in electric conductors and semiconductors. The Hall tensor is third-order and three-dimensional, whose first two indices are skew-symmetric. This paper investigates the isotropic polynomial invariants of the Hall tensor by connecting it with a second-order tensor via the third-order Levi-Civita tensor. A minimal isotropic integrity basis with 10 invariants for the Hall tensor is proposed. Furthermore, it is proved that this minimal integrity basis is also an irreducible isotropic function basis of the Hall tensor.     

两端铰接的细长柔性圆柱体涡激振动响应特性数值研究

目前细长柔性圆柱体涡激振动响应的研究方法主要包括实验方法、计算流体动力学方法以及半经验模型方法. 鉴于实验方法研究成本较高、计算流体动力学方法计算时间较长,本文基于尾流振子模型对线性剪切来流下两端铰接的细长柔性圆柱体涡激振动响应特性进行了半经验模型方法研究. 先建立了柔性圆柱体结构振子以及尾流振子之间的耦合模型,紧接着基于二阶精度中心差分格式对耦合模型先离散后迭代进行求解. 对不同剪切参数下柔性圆柱体涡激振动响应的振动波长、振动频率、振动位移以及响应频率随时间的变化特性等参数进行了分析. 分析结果表明:圆柱体的涡激振动响应由驻波和行波混合组成. 当无量纲弯曲刚度较小时,在圆柱体两端附近,驻波占主导;而在圆柱体中间段附近,行波占主导. 当无量纲弯曲刚度较大时,在圆柱体整个长度区间上均为驻波占主导. 随着剪切参数的增大,振动位移以及振动波长均逐渐减小,而振动频率和频率带宽均逐渐增大.      

多自由度Maxwell阻尼器减震结构(Ⅱ)——等效阻尼

对多自由度带支撑Maxwell阻尼器减震结构的等效阻尼比进行了系统研究。构建了耗能结构一般微分和积分混合地震响应方程组;基于与多自由度随机平均法分析完全相同的等效准则,建立了耗能结构各振型等效阻尼比的一般解析计算式。与一些典型问题的模态应变能法的计算精度进行对比分析,结果显示:本文方法的计算结果在特定条件下与精确解完全相同;在一个二自由度系统中相对位移响应方差的误差分别为0.99%和0.45%,优于应变能法的8%和5.4%,表明了本文方法的有效性。从而为建立此类耗能结构等效阻尼分析的完备解析解法、直接应用反应谱法进行实际工程设计提供了依据。     

失谐弱耦合卫星天线结构振动分析及预测控制

为了研究弱耦合卫星天线结构的振动控制,建立了该结构的简化计算模型,并针对该模型研究了弱耦合卫星天线结构动力学性能的特殊性:结构失谐时的振动模态局部化现象;针对失谐前后的结构,采用预测控制方法进行了振动控制,并与二次线性最优控制(LQR)方法的振动抑制效果进行了对比. 仿真结果表明:弱耦合星载天线结构参数的微小失谐会导致结构振动产生明显的模态局部化;采用预测控制方法进行结构振动控制的效果明显优于LQR控制方法,且在失谐导致的模型失配时,预测控制方法对结构振动亦有较好的抑制;在进行此种结构的振动主动控制时必须考虑到结构失谐的影响.      

不同边界条件正交各向异性中厚板的振动与稳定分析

本文采用高阶剪切变形理论对正交各向异性中厚矩形板进行振动与稳定分析,数值计算采用样条有限点法,得出了六种不同边界条件矩形板的自振频率和屈曲载荷,并与相应的经典板理论的结果进行比较.结果说明横向剪切变形对复合材料层合板的影响与板的各向异性程度、板的宽厚比(b/h)、层合板的层数和板的支承条件有关,它随着层合板各向异性程度的增加而增加,随着层合板宽厚比的增加而逐渐消失.     

winkler地基上弹性圆薄板受刚体撞击的动力响应

本文针对工程实际中所遇到的撞击问题,在事先仅知撞击体初始速度的条件下,研究分析了半无限粘弹性Winkler地基上的弹性圆薄板受刚体撞击的动力响应问题,推导出了关于撞击力F(t)的非线性Volterra积分方程,给出了薄板位移响应W(r,θ,t)的一般表达式,并给出了相应的数值求解方法。作为实例,本文对周边固定的弹性圆薄板在圆心处受刚球撞击问题进行了分析计算,并对某些参数的影响进行了讨论。     

混合状态Hamiltonian元的半解析解和叠层板的计算

本文给出混合状态方程的一种强有力的半离散半解析解法,并给出了矩形域和叠层板的算例,讨论了可能的发展和应用。     

The peridynamic differential operator for solving time-fractional partial differential equations

Nonlinear Dynamics - In this paper, the numerical solution of time-fractional convection diffusion equations (TF-CDEs) is considered as a generalization of classical ones, nonexponential relaxation...