Stress relaxation in broken fibers in unidirectional composites: modeling and application to creep rupture analysis

This paper describes a model of stress relaxation in broken fibers in unidirectional metal matrix composites reinforced with long brittle fibers. A cylindrical cell with a broken fiber embedded in a power law creeping matrix is employed, and the broken fiber is assumed to have a bilinear distribution of axial stress. Then, on the basis of energy balance in the cell under constant overall strain, a relaxation equation of interfacial shear stress acting on stress recovery segments is derived in a simple form. The relaxation equation is approximated rationally and integrated to obtain an analytical solution, which is shown to agree fairly well with the numerical analysis of Du and McMeeking. (Du, Z.-Z., McMeeking, R.M., 1995. Creep models for metal matrix composites with long brittle fibres. J. Mech. Phys. Solids 43, 701–726.) Moreover, the relaxation equation is combined with Curtin's model (Curtin, W.A., 1991. Theory of mechanical properties of ceramic-matrix composites. J. Am. Ceram. Soc. 74, 2837–2845.), so that rupture time in long term creep is evaluated analytically and explicitly on the assumption of global load sharing. It is shown that the resulting relation represents well the dependence of creep rupture time on applied stress observed experimentally on a unidirectional metal matrix composite.     

Generalization of the Krook kinetic relaxation equation

One of the most significant achievements in rarefied gas theory in the last 20 years is the Krook model for the Boltzmann equation [1]. The Krook model relaxation equation retains all the features of the Boltzmann equation which are associated with free molecular motion and describes approximately, in a mean-statistical fashion, the molecular collisions. The structure of the collisional term in the Krook formula is the simplest of all possible structures which reflect the nature of the phenomenon. Careful and thorough study of the model relaxation equation [2–4], and also solution of several problems for this equation, have aided in providing a deeper understanding of the processes in a rarefied gas. However, the quantitative results obtained from the Krook model equation, with the exception of certain rare cases, differ from the corresponding results based on the exact solution of the Boltzmann equation. At least one of the sources of error is obvious. It is that, in going over to a continuum, the relaxation equation yields a Prandtl number equal to unity, while the exact value for a monatomic gas is 2/3.In a comparatively recent study [5] Holway proposed the use of the maximal probability principle to obtain a model kinetic equation which would yield in going over to a continuum the expressions for the stress tensor and the thermal flux vector with the proper viscosity and thermal conductivity.In the following we propose a technique for constructing a sequence of model equations which provide the correct Prandtl number. The technique is based on an approximation of the Boltzmann equation for pseudo-Maxwellian molecules using the method suggested by the author previously in [6], For arbitrary molecules each approximating equation may be considered a model equation. A comparison is made of our results with those of [5].     

螺旋压缩弹簧应力松弛特性试验分析及其应用

本文通过应力松弛试验、理论推导及数值模拟研究了高温下螺旋压缩弹簧的应力松弛规律,并利用加速模型对工况下弹簧应力松弛服役寿命做出预测。首先,根据螺旋压缩弹簧的结构特点搭建了弹簧应力松弛连续动态测试装置,该装置不仅避免了传统测试方法存在的缺陷,而且能够保证试验过程中位移载荷恒定,并实时监测载荷变化。本文以某飞机舱门单锁机构中的螺旋压缩弹簧为试验对象进行了不同温度条件下的应力松弛试验,得到其松弛动力学曲线,并基于Arrhenius模型建立了弹簧在工况下的应力松弛服役寿命预测模型;其次,基于应力松弛和蠕变在本质上的一致性,结合金属材料蠕变规律并根据试验弹簧的受力特点,推导出用于描述试验材料松弛行为的蠕变本构方程,由试验结果获得该本构方程的材料常数;最后,通过该本构方程及材料常数,在ANSYS软件中对试验弹簧的松弛过程进行模拟,结果表明,模拟结果与试验结果误差不大于4%。因此,通过本文方法所建立的蠕变方程对弹簧在不同载荷条件下的应力松弛规律进行仿真分析具有一定的可行性与准确性。     

Relaxation and retardation functions of the Maxwell model with fractional derivatives

A four-parameter Maxwell model is formulated with fractional derivatives of different orders of the stress and strain using the Riemann-Liouville definition. This model is used to determine the relaxation and retardation functions. The relaxation function was found in the time domain with the help of a power law series; a direct solution was used in the Laplace domain. The solution can be presented as a product of a power law term and the Mittag-Leffler function. The retardation function is determined via Laplace transformation and is solely a power law type.The investigation of the relaxation function shows that it is strongly monotonic. This explains why the model with fractional derivatives is consistent with thermodynamic principles.This type of rheological constitutive equation shows fluid behavior only in the case of a fractional derivative of the stress and a first order derivative of the strain. In all other cases the viscosity does not reach a stationary value.In a comparison with other relaxation functions like the exponential function or the Kohlrausch-Williams-Watts function, the investigated model has no terminal relaxation time. The time parameter of the fractional Maxwell model is determined by the intersection point of the short- and long-rime asymptotes of the relaxation function.     

Shear relaxation characteristics of rock joints under stepwise loadings

The stress relaxation characteristic of rock mass is an important aspect of rheology and has important practical significance for rock engineering. In order to investigate the relaxation characteristic of rock joints with different slope ratios and normal stresses, a series of shear stress relaxation tests were conducted on artifical rock joints poured by cement mortar. Test results show that the relaxation curves can be divided into three stages, i.e. instantaneous relaxation stage, attenuation relaxation stage, and stable relaxation stage. Furthermore, the nonlinear Maxwell relaxation equation was obtained by using the relation between the viscosity coefficient and time, and the theoretical curves based on the empirical equation agreed well with the test results. Moreover, the change law of the initial viscosity coefficient was investigated. Accordingly, a stress relaxation method, termed as relaxation stress peak method, was proposed to determine the long-term strength of rock joints.     

Interrelation of creep and relaxation for nonlinearly viscoelastic materials: application to ligament and metal

Creep and stress relaxation are known to be interrelated in linearly viscoelastic materials by an exact analytical expression. In this article, analytical interrelations are derived for nonlinearly viscoelastic materials which obey a single integral nonlinear superposition constitutive equation. The kernel is not assumed to be separable as a product of strain and time dependent parts. Superposition is fully taken into account within the single integral formulation used. Specific formulations based on power law time dependence and truncated expansions are developed. These are appropriate for weak stress and strain dependence. The interrelated constitutive formulation is applied to ligaments, in which stiffness increases with strain, stress relaxation proceeds faster than creep, and rate of creep is a function of stress and rate of relaxation is a function of strain. An interrelation was also constructed for a commercial die-cast aluminum alloy currently used in small engine applications.     

Analysis of rheological deformation of soil by means of finite element method

A linear time-dependent viscoelastic behavior of soil is analyzed by the finite element method (FEM), which has great advantage in obtaining an approximate numerical solution of deformation or stress of a continuous body under complex boundary conditions, as known nowadays worldwide. A rheological three-element model, which is easily handled and represents rationally the actual behavior of soil, is suggested to obtain the rheological constants and the constitutive equation of soil. As actual examples of soil behavior, a stress relaxation of a soil block and a time-dependent sinkage of a rigid wheel are calculated by FEM and are also compared with test results and theoretical values.     

A non-linear uniaxial integral constitutive equation incorporating rate effects,creep and relaxation

A previously proposed first order non-linear differential equation for uniaxial viscoplasticity, which is non-linear in stress and strain but linear in stress and strain rates, is transformed into an equivalent integral equation. The proposed equation employs total strain only and is symmetric with respect to the origin and applies for tension and compression. The limiting behavior for large strains and large times for monotonic, creep and relaxation loading is investigated and appropriate limits are obtained. When the equation is specialized to an overstress model it is qualitatively shown to reproduce key features of viscoplastic behavior. These include: initial linear elastic or linear viscoelastic response: immediate elastic slope for a large instantaneous change in strain rate normal strain rate sensitivity and non-linear spacing of the stress-strain curves obtained at various strain rates; and primary and secondary creep and relaxation such that the creep (relaxation) curves do not cross. Isochronous creep curves are also considered. Other specializations yield wavy stress-strain curves and inverse strain rate sensitivity. For cyclic loading the model must be modified to account for history dependence in the sense of plasticity.     

考虑应力松驰复合型裂纹顶端塑性区分析

Ｓｈｉｈ［１］应用奇异单元，获得了不考虑应力松驰小范围屈服条件下复合型裂纹尖端塑性区形状。Ｚ．Ｚ．Ｚｕ等［２］采用Ｒｉｃｅ［５］给出的裂纹尖端应力关系式，利用有限元分析获得了不考虑应力松驰下复合型裂纹尖端塑性区，本文基于静力学中内力与外力平衡条件，用线弹性的全场解代替局部解，给出了考虑应力松驰下复合型裂纹尖端塑性区边界方程，获得了考虑应力松驰下的任意方向的塑性区尺寸及塑性区形状     

一种新的玻尔兹曼方程可计算模型构造与分析

临近空间飞行器因各部件尺寸差异较大, 在高空高马赫数条件下飞行会出现多流区共存的多尺度复杂非平衡流动现象, 流场中的气体分子速度分布函数与当地位置、流场分子速度、气体密度、流动速度、温度、热流矢量、应力张量等相关. 通过分析玻尔兹曼方程的一阶查普曼?恩斯科近似解, 构造了一种同时考虑热流矢量和应力张量影响、满足玻尔兹曼方程高阶碰撞矩的跨流域统一可计算模型方程, 并在数学上分析了其守恒性、H定理等基本属性, 证明了新模型方程与玻尔兹曼方程的相容性, 给出了新模型与现有模型如沙克霍夫(Shakhov)模型等的递进关系, 基于碰撞动力学确定了各流域统一气体分子碰撞松弛参数表达式. 在气体动理论统一算法中采用新建模型及现有模型模拟了一维激波结构、二维近空间飞行环境平板和多体圆柱干扰流动, 并与直接模拟蒙特卡洛方法对比分析, 结果表明在流场中粘性效应显著的区域新建模型能更好地捕捉激波位置, 尤其是在激波内部新模型模拟的宏观参数分布与直接模拟蒙特卡洛方法结果符合更好, 验证了新模型的有效性和可靠性, 同时说明在非平衡显著的流动区域碰撞松弛模型受多参数共同作用的影响.      

液晶高分子各向异性粘弹性流体本构方程理论

将液晶高分子－各向异性流体的本构方程,建立在Oldroyd随体导数观点基础上。推广上随机Oldroyd B流体模型,提出共转OldroydB流体模型,同时将微观结构的影响通过宏观参数表示出来,使在宏观理论中包含微观结构的贡献,即引入取向物质函数,非线性各向异性黏度函数和各向异性松弛时间及推迟时间等,表征取向运动对黏度和松弛及推迟现象的影响,在此基础上开展了一类新的液晶高分子－Oldroyd型本构方程理论,由该类型本构方程得出的物质函数,液晶高分子流体的第一、第二法向应力差与实验结果一致,解释了液晶高分子溶液的第一、第二法向应力差的特殊流变学行为。     

黏弹性材料界面裂纹应力场奇性分析

研究两半无限大黏弹性体间Griffith界面裂纹在简谐载荷作用下裂纹尖端动应力场的奇异特性.通过引入裂纹张开位移和裂纹位错密度函数,相应的混合边值问题归结为一组耦合的奇异积分方程.渐近分析表明裂尖动应力场的奇异特征完全包含在奇异积分方程的基本解中.通过对基本解的深入分析发现黏弹性材料界面裂纹裂尖动应力场具有与材料参数和外载荷频率相关的振荡奇异特性.以标准线性固体黏弹材料为例讨论了材料参数和载荷频率对奇性指数和振荡指数的影响.     

‘Generalized Burgers' equation’ for nonlinear viscoelastic waves

This paper deals with nonlinear longitudinal waves in a viscoelastic medium in which the viscoelastic relaxation function has the form K(t) = const. t^-v (0<v<1). This sort of slow relaxation may be more appropriate for polymers than the often used exponential relaxation. For a far field evolution of unindirectional waves, a “generalized Burgers' equation” is obtained, which is of a form with the second derivative in the usual Burgers' equation replaced by the derivative of real order 1 + v. The steady shock solution and self-similar pulse solution to this equation are discussed. In both cases numerical solutions are presented and analytic results are obtained for the asymptotic behaviors of the solutions. It is found that both shock and pulse solutions rise exponentially, but in their tails they have slow, algebraic decay.     

Influence of rotation and generalized magneto-thermoelastic on Rayleigh waves in a granular medium under effect of initial stress and gravity field

Influence of rotation, relaxation times, magnetic field, initial stress and gravity field on attenuation coefficient (Imaginary part of frequency equation root) and Rayleigh waves velocity (the real part of frequency equation root) in an elastic half-space of granular medium is studied. The analytical solution is obtained by using Lame’s potential techniques. The numerical calculations are carried out for the frequency equation of Rayleigh waves velocity. The results are displayed graphically. Some results of previous investigations are deduced as special cases from this study.     

A state space formalism and perturbation method for generalized thermoelasticity of anisotropic bodies

A state space formalism for generalized anisotropic thermoelasticity accounting for thermomechanical coupling and thermal relaxation is developed, which includes the classical thermoelasticity as a special case. By properly grouping the field variables using matrix notations, the basic equations of thermoelasticity are formulated into a state equation and an output equation in terms of the state vector. To obtain the solution for a specific problem it suffices to solve the state equation under the prescribed conditions. For weak thermomechanical coupling an asymptotic solution can be obtained by using the method of perturbation with multiple scales. Propagation of plane harmonic thermoelastic waves in an anisotropic medium is studied within the context.     

Relaxation behavior and modeling

Load relaxation tests deliver several orders of magnitude of inelastic strain rate data while elastic strains are converted into inelastic strains [see Lemaitre and Chaboche, 1994. (Mechanics of Solid Materials, Oxford University Press, Cambridge p. 264)]. Hart used this test for providing information on the inelastic deformation behavior for modeling purposes. Characteristic relaxation curves were obtained with ductile metals and alloys at room and high temperature showing a scaling relation derived from Hart's theory. Subsequent testing with servo-controlled testing machines and strain measurement on the gage length showed that an increase of prior strain rate also increased the average relaxation rate. For relaxation tests starting in the flow stress region, the relaxation curves can be independent of the stress and strain at the start of the relaxation tests. For the modeling of these newly found relaxation behaviors and other phenomena the viscoplasticity theory based on overstress (VBO) has been introduced. It is shown that VBO admits a long-term (asymptotic) solution that can be used with sufficient accuracy for the flow stress region of the stress–strain diagram. The long-term solution predicts the observed relaxation behaviors and that the relaxation curves coincide when shifted along the stress axis. This behavior is observed for the recently obtained data and is confirmed by two sets of the Hart-type data when they are plotted according to the new method.     

A GENERALIZED PIEZOELECTRIC-THERMOELASTIC THEORY WITH STRAIN RATE AND ITS APPLICATION 1)

工程中大量材料的形变介于弹性与黏性之间, 既具有弹性固体特性, 又具有黏性流体特点, 即为黏弹性. 黏弹性使得材料出现很多力学松弛现象, 如应变松弛、滞后损耗等行为. 在研究受热载荷作用的多场耦合问题的瞬态响应时, 考虑此类问题中的热松弛和应变松弛现象, 对准确描述其瞬态响应尤为重要. 针对广义压电热弹问题的瞬态响应, 尽管已有学者建立了考虑热松弛的广义压电热弹模型, 但迄今, 尚未计入应变松弛. 本文中, 考虑到材料变形时的应变松弛, 通过引入应变率, 在Chandrasekharaiah广义压电热弹理论的基础之上, 经拓展, 建立了考虑应变率的广义压电热弹理论. 借助热力学定律, 给出了理论的建立过程并得到了相应的状态方程及控制方程. 在本构方程中, 引入了应变松弛时间与应变率的乘积项, 同时, 分别在本构方程和能量方程中引入了热松弛时间因子. 其后, 该理论被用于研究受移动热源作用的压电热弹一维问题的动态响应问题. 采用拉普拉斯变换及其数值反变换, 对问题进行了求解, 得到了不同应变松弛时间和热源移动速度下的瞬态响应, 即无量纲温度、位移、应力和电势的分布规律, 并重点考察了应变率对各物理量的影响效应, 将结果以图形形式进行了表示. 结果表明: 应变率对温度、位移、应力和电势的分布规律有显著影响.