粘弹性力学的对应原理及其数值反演方法

积分变换是处理粘弹性混合边值问题的重要数学工具,积分变换的应用使粘弹性混合边值问题在象空间与相应弹性混合边值问题对应起来,从而使粘弹性混合边值问题的求解可以继承和借鉴弹性问题的求解方法,再利用积分反演方法就可求得时间域粘弹性边值问题的解．本文结合国内外的研究成果,就粘弹性力学中存在的各种对应原理及数值反演方法进行了归类和总结．结合在求解粘弹性边值问题中的应用,对各类方法的特点进行了评述,并指出存在的问题及发展新的数值方法的研究重点．      

分数阶导数粘弹性模型的有限元法

在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数值解。并以二维沥青路面结构为例进行了路面动态粘弹性响应分析。算例分析表明,该方法能够正确有效地进行路面动态粘弹性分析。     

无网格方法在平面粘弹性力学问题中的应用

本文综合应用无网格方法(EFGM)、线性粘弹性与弹性力学之间的对应原理,Laplace变换和逆变换等方法求解了拟静态平面弹性和粘弹性力学问题。首先,利用Laplace变换和逆变换推导了平面问题的粘弹性本构关系,建立了拟静态粘弹性平面问题的边值问题;其次,利用粘弹性与弹性力学之间的对应原理得到了Laplace变换域中平面问题的基本方程,在Laplace变换域中建立了相应的泛函,并得到了用无网格方法离散的控制方程;同时,求解了几个拟静态弹性和粘弹性平面问题,给出了它们的表达式和数值结果;最后,采用Laplace逆变换和数值逆变换,得到了粘弹性力学平面问题在物理空间中的解,并比较了由解析解和无网格数值方法所得到的数值结果,可以看到它们是非常吻合的。说明本文方法的正确性和有效性。     

粘弹性固体的精细积分有限元算法

粘弹性固体本构方程的数学表达式分为微分型和积分型两种，其数值求解主要是时域上离散计算。文中从微分型表达式出发导出其状态空间方程的数学表达式，通过严格推导论证了它与微、积分型表达式的等价性；引入状态空间方程，从而利用精细积分格式来求解粘弹性固体本构方程；给出了粘弹性固体本构方程的精细积分有限元算法，为求解粘弹性固体本构方程的数值解提供了一个新的途径，具有计算简便，求解精度高等优点。     

结构粘弹性分析的线法

本文将以现代常微分方程求解器为支撑软件的线法推广应用于含时间变量的粘弹性力学问题。基基本思想是先用拦普拉斯积分变换把问题的控制偏微分方程组和定解条件进行交换,得到拉氏变换空间中的定常问题,然后用传统的线法借助于常微分方程求解器得到该变换空间中的数值解,再应用拉氏数值逆变换求得原时空域中的数值解。文中还给出方法的基本原理和实施过程。算例表明,采用拉普拉斯积分变换于线法分析结构粘弹性力学等含时间变量的     

模拟微可压粘弹性流体的WCCBS_SU方法

针对微可压缩粘弹性流动问题,发展了微可压缩流的WCCBS方法,详细推导了基于Oldroyd-B本构模型的WCCBS_SU方法的求解过程。在流场微可压的条件下,分别对平面Poiseuille流和4：1粘弹性收缩流进行了数值模拟。Poiseuille流在不同We数下数值结果与解析解的比较,验证了本文方法具有较高的精度和较好...     

考虑频率依赖性的涂层复合结构固有特性求解

粘弹性阻尼材料的力学特性参数会随着频率的变化而改变,即具有频率依赖性,因而传统的动力学建模及分析方法不能满足实际涂层结构优化设计的需要。在简要介绍粘弹性阻尼材料频率依赖性的基础上,本文提出用特征向量增值法来求解涂层复合结构的固有特性,并详细推导了特征向量增值法的求解原理。由此,提出了特征向量增值法的计算流程,包括计算无阻尼系统的固有特性,用Fox and Kapoor或者Nelson方法计算复特征向量增量;用Rayleigh熵法求解复特征值。最后,以涂敷粘弹性阻尼材料的钛基薄板为例,求解了该复合结构的固有特性,并与经典的模态应变能法进行了比较,证明了所提方法的正确性。     

粘弹性夹芯层合旋转圆板的气动弹性动力稳定性分析

基于哈密尔顿原理,考虑作用于圆板上的空气动力载荷以及粘弹性复模量本构模型,本文建立了含有粘弹性夹芯层的旋转圆板的空气动力学模型,并给出了相应的稳定性的控制方程及边界条件。运用Galerkin方法,文中数值求解并给出了夹芯层合旋转圆板的动力学基本特征,以及结构几何参数、材料参数对于其前、后行波振动频率、阻尼和动力稳定性的影响。结果表明：通过优化合理选取粘弹性夹芯层合旋转圆板的几何参数和材料参数,可提高旋转圆板的临界转速和颤振速度,从而增强其稳定性。     

基于时域自适应算法的单向粘弹性节理岩体的等效分析

在粘弹性节理岩体的数值仿真中,必需合理地计及节理的影响.若将节理与岩体作为独立的材料组分考虑,在节理密集的情况下,时空两方面的计算开销可能是很难承受的.本文通过时域展开技术与一个简单的等效假定,得到了递推格式的粘弹性单向节理岩体的等效本构关系,及等效场的递推求解格式,并由此数值模拟了一个粘弹性节理岩体中的洞室,计算结果与考虑独立组分的ANSYS粘弹性分析进行了比较,二者的计算平均误差约为11.03%,但前者的计算效率为后者16.81倍.本文工作有可能为粘弹性节理岩体的数值模拟提供一条新的途径.     

设置粘弹性阻尼耗能框架结构减震性能分析

研究粘弹性阻尼器对框架结构的减震性能,分析该阻尼器的阻尼参数对框架结构地震反应的影响,给出设置阻尼器后结构地震作用的计算方法。对设置粘弹性阻尼器框架结构振动方程进行傅里叶变换,求解该振动方程的频域解,再通过傅里叶逆变换得出该方程的时域解,由此研究阻尼器的阻尼系数C_d和松弛时间系数η等对结构地震响应的影响。在此基础上,提出粘弹性阻尼器的阻尼系数在各楼层经济分布的实用计算公式,并对设置粘弹性阻尼器框架结构地震作用给出计算方法,分析结构各楼层的地震作用及结构各振型的底部剪力。研究表明,设置于框架结构的粘弹性阻尼器宜控制其松弛时间系数η≤0.5～0.7。     

设置连接阻尼器的相邻结构随机地震动响应的简明封闭解

针对两相邻结构间设置连接阻尼器对结构的减震影响问题,研究了基于Kanai-Tajimi谱地震动激励下的Kelvin型粘弹性阻尼器与相邻结构形成的组合体系的随机地震动系列响应(绝对位移及层间位移)的简明封闭解。首先,利用Kelvin型粘弹性阻尼器本构关系及Kanai-Tajimi谱的滤波方程,将组合体系基于复杂地震动激励精确转化为基于简明白噪声激励的运动方程;其次,利用复模态法获得了组合结构相对于地面的绝对位移、层间位移等系列响应方差及0阶~2阶谱矩的统一简明封闭解。最后,通过算例及与虚拟激励法进行对比,证明本文方法的正确性和简明性;通过与未设置阻尼装置结构体系的动力响应对比,说明了阻尼装置对相邻结构具有良好的减震效果,但局部楼层的层间位移及层间剪力会有所增加。     

General dynamic equation and dynamical characteristics of viscoelastic Timoshenko beams

In this paper, a governing differential equation of viscoelastic Timoshenko beam including both extension and shear viscosity is developed in the time domain by direct method. To measure the complex moduli and three parameters of standard linear solid, the forced vibration technique of beam is successfully used for PCL and PMMA specimens. The dynamical characteristics of viscoelastic Timoshenko beams, especially the damping properties, are derived from a considerable number of numerical computations. The analyses show that the viscosity of materials has great influence on dynamical characteristics of structures, especially on damping, and the standard linear solid model is the better one for describing the dynamic behavior of high viscous materials.     

Fully resolved simulations of viscoelastic suspensions by an efficient immersed boundary-lattice Boltzmann method

An efficient immersed boundary-lattice Boltzmann method (IB-LBM) is proposed for fully resolved simulations of suspended solid particles in viscoelastic flows. Stress LBM based on Giesekus and Oldroyd-B constitutive equation are used to model the viscoelastic stress tensor. A boundary thickening-based direct forcing IB method is adopted to solve the particle–fluid interactions with high accuracy for non-slip boundary conditions. A universal law is proposed to determine the diffusivity constant in a viscoelastic LBM model to balance the numerical accuracy and stability over a wide range of computational parameters. An asynchronous calculation strategy is adopted to further improve the computing efficiency. The method was firstly applicated to the simulation of sedimentation of a single particle and a pair of particles after good validations in cases of the flow past a fixed cylinder and particle migration in a Couette flow against FEM and FVM methods. The determination of the asynchronous calculation strategy and the effect of viscoelastic stress distribution on the settling behaviors of one and two particles are revealed. Subsequently, 504 particles settling in a closed cavity was simulated and the phenomenon that the viscoelastic stress stabilizing the Rayleigh–Taylor instabilities was observed. At last, simulations of a dense flow involving 11001 particles, the largest number of particles to date, were performed to investigate the instability behavior induced by elastic effect under hydrodynamic interactions in a viscoelastic fluid. The elasticity-induced ordering of the particle structures and fluid bubble structures in this dense flow is revealed for the first time. These simulations demonstrate the capability and prospects of the present method for aid in understanding the complex behaviors of viscoelastic particle suspensions.     

带支撑广义Maxwell粘弹性阻尼耗能结构风振响应分析

工程中粘弹性阻尼器的安装通过支撑与结构进行连接,但在安装粘弹性阻尼器的耗能结构随机响应分析中,为了简化计算过程,将支撑的刚度看作无穷大或者忽略支撑刚度的影响。实际上,对支撑刚度的影响加以考虑更能符合工程实际。针对考虑支撑刚度影响的粘弹性阻尼耗能结构风振响应分析过程复杂的问题,提出了一种求解考虑支撑影响的广义Maxwell粘弹性阻尼耗能结构基于Davenport谱风振响应的简明解析解法。在广义Maxwell粘弹性阻尼器微分型本构模型基础上,给出了考虑支撑刚度的粘弹性阻尼器等效本构关系。将粘弹性阻尼器等效本构关系与结构运动方程联立,采用复模态法将其解耦,获得结构风振响应的统一表达式。将耗能结构在Davenport风速谱下的系列响应功率谱密度函数分解为频域响应函数与Davenport功率谱密度函数的乘积形式,基于随机振动理论中谱矩的定义,对响应功率谱密度函数积分后获得无积分项的系列响应谱矩表达式。在算例中通过与虚拟激励法计算结果对比,验证了所提解法的准确性,并分析了支撑刚度在耗能系统中的影响。     

Approximate and numerical analysis of nonlinear forced vibration of axially moving viscoelastic beams

Steady-state periodical response is investigated for an axially moving viscoelastic beam with hybrid supports via approximate analysis with numerical confirmation. It is assumed that the excitation is spatially uniform and temporally harmonic. The transverse motion of axially moving beams is governed by a nonlinear partial-differential equation and a nonlinear integro-partial-differential equation. The material time derivative is used in the viscoelastic constitutive relation. The method of multiple scales is applied to the governing equations to investigate primary resonances under general boundary conditions. It is demonstrated that the mode uninvolved in the resonance has no effect on the steady-state response. Numerical examples are presented to demonstrate the effects of the boundary constraint stiffness on the amplitude and the stability of the steady-state response. The results derived for two governing equations are qualitatively the same,but quantitatively different. The differential quadrature schemes are developed to verify those results via the method of multiple scales.     

New matrix method for analyzing vibration and damping effect of sandwich circular cylindrical shell with viscoelastic core

Based on the linear theories of thin cylindrical shells and viscoelastic materials, a governing equation describing vibration of a sandwich circular cylindrical shell with a viscoelastic core under harmonic excitation is derived. The equation can be written as a matrix differential equation of the first order, and is obtained by considering the energy dissipation due to the shear deformation of the viscoelastic core layer and the interaction between all layers. A new matrix method for solving the governing equation is then presented With an extended homogeneous capacity precision integration approach. Having obtained these, vibration characteristics and damping effect of the sandwich cylindrical shell can be studied. The method differs from a recently published work as the state vector in the governing equation is composed of displacements and internal forces of the sandwich shell rather than displacements and their derivatives. So the present method can be applied to solve dynamic problems of the kind of sandwich shells with various boundary conditions and partially constrained layer damping. Numerical examples show that the proposed approach is effective and reliable compared with the existing methods.     

Stability of axially accelerating viscoelastic beams: multi-scale analysis with numerical confirmations

Stability is investigated for an axially accelerating viscoelastic beam. The material time derivative is used in the viscoelastic constitutive relation, not simply the partial time derivative. The method of multiple scales is applied directly to the governing equation without discretization. When the axial speed is characterized as a simple harmonic variation about the constant mean speed, the instability conditions are presented for axially accelerating viscoelastic beams constrained by simple supports with rotational springs in parametric resonance. The finite difference schemes are developed to solve numerically the equation of axially accelerating viscoelastic beams with fixed supports for the instability regions in the principal parametric resonance. The numerical calculations confirm the analytical results. Numerical examples show the effects of the constraint stiffness, the mean axial speed, and the viscoelasticity.     

An incompressible SPH method for simulation of unsteady viscoelastic free-surface flows

In this paper, an incompressible smoothed particle hydrodynamics (SPH) method is presented to solve unsteady free-surface flows. Both Newtonian and viscoelastic fluids are considered. In the case of viscoelastic fluids, both the Maxwell and Oldroyd-B models are investigated. The proposed SPH method uses a Poisson pressure equation to satisfy the incompressibility constraints. The solution algorithm is an explicit predictor-corrector scheme and employs an adaptive smoothing length based on density variations. To alleviate the numerical difficulties encountered when fluid is highly stretched, an artificial stress term is incorporated into the momentum equation which reduces the risk of unrealistic fractures in the material. Two challenging test cases, the impacting drop and the jet buckling problems, are solved to demonstrate the capability of the proposed scheme in handling viscoelastic flows with complex free surfaces. The jet buckling test case was solved for a wide range of Weissenberg numbers. It was shown that in all cases the method is stable and fairly accurate and agrees well with the available data.     

A simplified viscoelastic constitutive equation applied to a dynamic finite deformation problem

A single integral type of constitutive equation for finite viscoelastic deformation is proposed. A special case, which is a viscoelastic generalization of the constitutive equation for a neo-Hookean elastic solid, is used to consider the finite deformation problem of shock wave propagation resulting from the sudden application of compressive force at the end of a semi-infinite viscoelastic bar. An approximate method is used to determine the shock front path and shock strength when the viscoelastic dissipative effect is small.     

Optimized finite difference iterative scheme based on POD technique for 2D viscoelastic wave equation

This study develops an optimized finite difference iterative(OFDI) scheme for the two-dimensional(2D) viscoelastic wave equation. The OFDI scheme is obtained using a proper orthogonal decomposition(POD) method. It has sufficiently high accuracy with very few unknowns for the 2D viscoelastic wave equation. Existence, stability, and convergence of the OFDI solutions are analyzed. Numerical simulations verify efficiency and feasibility of the proposed scheme.