The viscoelastic behavior of short-fiber composite materials

This article examines the viscoelastic behavior of unidirectional short-fiber composites. Each fiber is approximated by an ellipsoidal elastic inclusion embedded in a viscoelastic matrix. The static correspondence principle which relates the viscoelastic problem to the elastic problem is applied to the modified expression of the effective elastic moduli of short-fiber composites originally derived by the authors. Numerical examples for various fiber aspect ratios are shown based upon the collocation approximation method. Comparisons of the result with the self-consistent method are also presented.     

求解一类非线性粘弹性问题的弹性回复对应原理

本文提出了一个全新的粘弹性理论体系,它把传统的线粘弹性理论作为特殊情况包括在内。其主要结果是两个用于求解一类物理非线性粘弹性问题的弹性回复对应原理。利用它,只要知道相应的非线性弹性问题的解,就可以求出非线性粘弹性问题的解答。对应原理不是基于本构关系的相似性,而是基于从粘弹性现时响应到瞬时弹性响应的可回复性^[1]。首先找到非线性粘弹性与非线性弹性本构关系之间的对应,然后导出了两个弹性回复对应原理。通过对改性聚丙烯材料的实验验证了该对应原理的正确性和对此类材料的实用性。     

高分子材料非线性粘弹性问题的解法

尽管实验限制在小变形范围之内，高分子材料仍表现出明显的非线性粘弹性性能。目前由于分析上的困难，通常只好采用线粘弹性假设。基于前文，提出了一个求解高分子材料物理非线性（含线性）粘弹性问题的新原理，通过对改性聚丙烯材料的实验表明了该原理的正确性和对于此类材料的实用性。     

粘弹结构动力学分析

本文研究由单一粘弹材料构成的粘弹结构的振动.在时间域和拉氏域内建立了它的动力学方程.给出了粘弹结构的"对偶弹性结构"概念,并用来讨论了粘弹结构的模态特性,证实位移形态可由实模态向量张成的向量空间内加以描述,它就是"对偶弹性结构"的模态向量,无需增广坐标.每个模态自由度含有数个反映蠕变的实特征值和一对反映衰减振动的复共轭特征值.分析了在模态向量空问内的动响应,给出了传递函数(频率响应)和脉冲响应.从而揭示了粘弹结构的动力学特性.     

The normal aging viscoelastic contact problem

The problem of a mechanically smooth rigid indentor of arbitrary asymmetric curved profile being pressed against the surface of a linear viscoelastic half-space is reconsidered. A mechanical argument is used to remove restrictions and extend the range of validity of the solution given in [2] which is simplified and given for an aging material. Still, the shape of the contact area is always the same as that of some member of the one parameter family of contact areas associated with the corresponding elastic problem and the viscoelastic solution is expressed in terms of the allied elastic solutions.     

The band gap of 1D viscoelastic phononic crystal

The band gap of one dimensional (1D) phononic crystal with viscoelastic host material is studied in this paper. The standard solid model is used to simulate the viscoelastic behavior of the host material and the fillers embedded in the host material are still assumed elastic material. The band gap problem in 1D phononic crystal leads to an eigenvalue problem by using the plane wave expansion method and the Bloch–Floquet wave theory in a periodic structure. An iterative algorithm is designed to obtain the band gap structure due to the dependence of elastic constants on frequency for the viscoelastic host material. A numerical example is given for steel/epoxy phononic crystal. The band gap of 1D phononic crystal is evaluated for different viscoelastic constants, namely, relaxation time, initial and final state elastic modulus. It is found that the viscoelastic constants of host material affect not only the location but also the width of band gaps.     

Dynamic stress concentration around a hole in a viscoelastic plate

Summary The transient hoop stresses which are generated at the circumference of a circular hole in a large viscoelastic plate when a radial pressure pulse acts on the hole boundary, are determined. An analytical/numerical approach is employed which is based on the use of Laplace transform and Hankel functions, and the Dubner-Abate-Crump technique for inverting the Laplace-transformed solution. In our formulation of the problem, a general linear-viscoelastic material for the plate is considered, but numerical results are extracted only for a three-parameter model (standard linear solid). The present work accompanies recent efforts by Georgiadis on analogous transient viscoelastodynamic problems involving finite-length cracks in a stress-wave environment. The mathematically simpler problem considered here demonstrates the influence of viscoelastic effects on the dynamic-stress-concentration in marked similarity with the dynamic-stress-intensity behavior of crack-tip viscoelastic stress fields. Besides that, the present problem generalizes the classical Kromm-Selberg-Miklowitz problem in the sense that it considers viscoelastic response. The latter solutions were not utilized here through the usual elastic/viscoelastic correspondence principle but, instead, an independent solution was derived by simple means.     

Topology optimization of unconstrained damping treatments for plates

A finite element model for composite plates consisting of an elastic isotropic base layer covered with viscoelastic treatment is presented. The composite plate undergoes bending vibrations in the lateral direction. In contrast to many previous publications, the viscoelastic treatment is not constrained from the top in order to better simulate real cases. The objective is to find the optimum distribution of viscoelastic treatment which maximizes the modal damping ratio (MDR) for a certain volume of treatment. Topology optimization is performed with two strategies: optimizing the whole domain of viscoelastic treatment and optimizing a unit cell of the periodic treatment. Numerical examples show that the presented model is able to increase the MDR by an order of magnitude compared to plain treatments.     

Direct boundary integral procedure for a Boltzmann viscoelastic plane with circular holes and elastic inclusions

A direct boundary integral method in the time domain is presented to solve the problem of an infinite, isotropic Boltzmann viscoelastic plane containing a large number of randomly distributed, non-overlapping circular holes and perfectly bonded elastic inclusions. The holes and inclusions are of arbitrary size and the elastic properties of all of the inclusions can, in general, be different. The method is based on a direct boundary integral approach for the problem of an infinite elastic plane containing multiple circular holes and elastic inclusions described by Crouch and Mogilevskaya [1], and a time marching strategy for viscoelastic analysis described by Mesquita and Coda [2–8]. Benchmark problems and numerical examples are included to demonstrate the accuracy and efficiency of the method.     

下卧饱和半空间粘弹性土层上基础阻抗函数的分析

由于地下水的影响,实际土层一般由含水土层和不含水土层组成。由于求解的复杂性,目前基础阻抗函数的求解中很少考虑这种地基情形。根据干土是饱和土的特殊情形这一事实,通过将液相压缩模量、孔隙比以及惯性耦合质量密度取为零,可以使干土的动力反应统一到饱和土动力反应求解方案中,从而可通过饱和地基上基础阻抗函数的一种求解方法,解决饱和土和干土并存时基础阻抗函数的求解问题。地基无限域的影响通过局部透射人工边界考虑。通过算例验证了该方法的可行性,并考察了干土层厚度以及基础埋深对动力阻抗的影响。     

A fast multipole boundary element method for 2D viscoelastic problems

A fast multipole formulation for 2D linear viscoelastic problems is presented in this paper by incorporating the elastic–viscoelastic correspondence principle. Systems of multipole expansion equations are formed and solved analytically in Laplace transform domain. Three commonly used viscoelastic models are introduced to characterize the time-dependent behavior of the materials. Since the transformed multipole formulations are identical to those for the 2D elastic problems, it is quite easy to implement the 2D viscoelastic fast multipole boundary element method. Besides, all the integrals are evaluated analytically, leading to highly accurate results and fast convergence of the numerical scheme. Several numerical examples, including planar viscoelastic composites with single inclusion or randomly distributed multi-inclusions, as well as the problem of a crack in a pressured viscoelastic plane, are presented. The results are verified by comparison with the developed analytical solutions to illustrate the accuracy and efficiency of the approach.     

Viscoelastic contact problems with friction

A mechanically rough rigid indentor of abitrary asymmetric curved profile is pressed against the surface of a linear viscoelastic half-space whose behaviour in shear and dilatation is similar. In the first place complete adhesion is assumed and the viscoelastic solution is found to be expressible in terms of corresponding elastic solutions. The shape of the contact area and the tangential component of displacement of a point of the surface of the viscoelastic body inside the contact area are unaffected by the viscoelasticity of the body. Fully adhesive contact between two viscoelastic bodies may be treated similarly. Also, certain frictional elastic contact problem solutions, where the contact area contains regions of slip as well as regions of adhesion, are extended to viscoelasticity.     

Estimation of piezoelastic and viscoelastic properties in laminated structures

An inverse method for material parameter estimation of elastic, piezoelectric and viscoelastic laminated plate structures is presented. The method uses a gradient based optimization technique in order to solve the inverse problem, through minimization of an error functional which expresses the difference between experimental free vibration data and corresponding numerical data produced by a finite element model. The complex modulus approach is used to model the viscoelastic material behavior, assuming hysteretic type damping. Applications that illustrate the influence of adhesive material interfaces and viscoelastic parameter identification are presented and a few simulated test cases aid the interpretation of results.     

Study of viscoelastic effect on the frequency shift of microcantilever chemical sensors

Microcantilevers coated with a chemically sensitive layer are increasingly being used in chemical detection systems. The sensitive coating, often a polymer, absorbs specific molecules, which can be detected by monitoring the shift in the mechanical resonant frequency. Usually, the frequency shift resulting from molecular absorption is interpreted as a mass loading effect. However, mass loading is not the only effect that has an impact on the frequency shift; the viscoelastic properties of the sensitive coating also are affected by the sorption process. Sorption-induced modulus changes are typically difficult to characterize. However, it is known that the sorption of analyte molecules in a polymer coating results in the plasticization of the coating. In most cases, the polymer becomes more rubbery with increasing concentration of analyte molecules, i.e., the coating becomes softer with increasing loss modulus, and the storage modulus decreases. Using a new analytical model developed for the resonant frequency expression of a hybrid microcantilever (elastic base and viscoelastic layer), the effects of the modification of the storage and loss moduli of the sensitive layer on the resonant frequency are examined. The main conclusion of this analytical study is that, even if the sensitive coating moduli are small compared to the base cantilever's Young's modulus, the effect of the change in the viscoelastic coating properties could contribute significantly to the overall frequency shift (8-23% in the simulations depending on the coating thickness, with even higher contributions for other sets of problem parameters).     

Efficient techniques for predicting viscoelastic behavior of sublaminates

Even for linear elastic behavior, stress analysis of thick laminated composites can be very computation intensive if every lamina is modeled discretely. In such cases, modeling of individual lamina is impractical and the homogenization method for sublaminates becomes essential. In the current work, 3D homogenization formulas for an elastic sublaminate, which were derived by the authors in previous work, were utilized to determine the 3D effective properties for a viscoelastic sublaminate. The properties were determined by three methods that exploited the 3D elastic homogenization formulas: (i) quasi-elastic method, (ii) correspondence principle, and (iii) direct time integration of the incremental viscoelastic equations. The finite element method with discrete modeling of the plies was used to obtain reference solutions. The effective viscoelastic properties obtained using the three methods based on the elastic homogenization formulas were in very good agreement with the reference solution. Among these methods, the quasi-elastic method was found to be both accurate and the simplest method in determining the effective properties. The methods were also used to predict the stress response of a sublaminate to different strain histories. The direct time integration method using the 3D elastic homogenization formulas performs accurately and efficiently for this problem.     

点弹性支承粘弹性矩形薄板的基频分析计算

从三维粘弹性本构关系出发,导出了具有多个点弹性支承的Kelvin型粘弹性矩形薄板的运动微分方程.针对方程中出现的二维广义δ函数,采用积分方程法导出了具有多个点弹性支承的四边简支Kelvin型粘弹性矩形薄板自由振动的复特征方程,分析了材料的无量纲延滞时间、点弹性支承的弹性系数和支承位置对矩形薄板的固有频率的影响.     

A boundary integral method for a dynamic, transient mode I crack problem with viscoelastic cohesive zone

We consider the problem of the dynamic, transient propagation of a semi-infinite, mode I crack in an infinite elastic body with a nonlinear, viscoelastic cohesize zone. Our problem formulation includes boundary conditions that preclude crack face interpenetration, in contrast to the usual mode I boundary conditions that assume all unloaded crack faces are stress-free. The nonlinear viscoelastic cohesive zone behavior is motivated by dynamic fracture in brittle polymers in which crack propagation is preceeded by significant crazing in a thin region surrounding the crack tip. We present a combined analytical/numerical solution method that involves reducing the problem to a Dirichlet-to-Neumann map along the crack face plane, resulting in a differo-integral equation relating the displacement and stress along the crack faces and within the cohesive zone.     

粘弹性复合材料中的渐近均匀化方法

主要研究了由各向同性线弹性加强体和各向同性线粘弹性基体组成的复合材料的问题。在已有的线弹性多层材料的渐近均匀化方法的基础上,应用弹性-粘弹性对应原理,在Carson域中求解粘弹性问题,通过两次运用均匀化方法,得到一类单向强化复合材料的有效模量的表达式。反演可得到单向强化复合材料的有效松弛模量在时间域中的表达式,并且与其它结果进行了比较。     

Modeling of the contact between a rigid indenter and a heterogeneous viscoelastic material

In this paper, the contact problem between a rigid indenter and a viscoelastic half space containing either isotropic or anisotropic elastic inhomogeneities is solved. The model presented here is 3D and based on semi-analytical methods. To take into account the viscoelastic properties of the matrix, contact and subsurface problem equations are discretized in the spatial and temporal dimensions. A conjugate gradient method and the fast Fourier transform are used to solve the normal problem, contact pressure, subsurface problem and real contact area simultaneously. The Eshelby’s formalism is applied at each step of the temporal discretization to account for the effect of the inhomogeneity on pressure distribution and subsurface stresses. This method can be seen as an enrichment technique where the enrichment fields from heterogeneous solutions are superimposed to the homogeneous viscoelastic problem solution. Note that both problems are fully coupled. The model is validated by comparison with a Finite Element Model. A parametric analysis of the effect of elastic properties and geometrical features of the inhomogeneity is proposed. The model allows to obtain the contact pressure distribution disturbed by the presence of inhomogeneities as well as subsurface and matrix/inhomogeneity interface stresses at every step of the temporal discretization.     

饱和横观各向同性分数导数黏弹性土中半封闭衬砌振动响应

在频率域内研究了饱和横观各向同性分数导数黏弹性土体中深埋圆形隧道半封闭衬砌振动响应问题。根据土体在长期沉积过程中存在各向异性的特点,将土骨架视为具有分数导数本构关系的横观各向同性黏弹性体,采用饱和多孔介质理论和弹性理论,利用衬砌内边界应力协调以及土体和衬砌界面处应力和位移连续,得到了简谐荷载作用下饱和横观各向同性黏弹性土和弹性衬砌的位移、应力和孔隙水压力解析表达式。考察了饱和经典弹性土、饱和分数导数性黏弹性土和饱和经典黏弹性土三种条件下饱和黏弹性土和衬砌各参数的影响,表明：横观各向同性面的弹性模量和衬砌厚度对系统动力响应的影响与分数导数阶数和土骨架的黏性有关;渗透系数较小时,系统存在明显的共振现象。另外,在三种条件下半封闭衬砌振动响应存在较大差异。