一种改进的R-L分数阶导数定义在固体推进剂粘弹性本构模型中的应用

针对Riemann-Liouville分数阶导数定义的不足之处进行了改进,利用改进过的分数阶导数定义,建立了类Kelvin体粘弹性本构模型,并应用于某固体推进剂上,对本构方程中的三个参数进行了求解,与经典的prony级数模型进行比较,采用分数阶导数的类Kelvin体粘弹性本构模型与实验结果能很好地吻合.     

松弛与蠕变力学特性转换关系的研究

本文首先论述了拉伸应力松弛模量Ｅ（ｔ）与拉伸蠕变柔量Ｄ（ｔ）之间的近似与精确转换关系式，提出一种诊断Ｄ（ｔ）的数值积分算法，通过应力松弛模量实验给出一种固体推进剂Ｅ（ｔ）和Ｄ（ｔ）的实测与转换计算结果，计算结果有较小的误差具有重要的实用价值。     

非定常微分型黏弹性本构模型

在应力作用下, 材料的力学参数随着微观结构的变化而变化, 需要考虑参数的时间效应. 利用黏滞系数随时间变化的黏性元件, 构造出非定常Maxwell模型、非定常Kelvin模型和非定常Zener模型. 求解非定常模型的微分型本构方程得到它们的松弛模量、蠕变柔量和卸载方程. 结果表明, 可以把常见的经验松弛函数和经验蠕变函数视为非定常微分型本构模型.      

边界元法在流—固相互作用分析中的应用

本文采用边界元法研究水与粘弹性介质在地震波入射条件下的相互作用问题。首先,采用半空间格林函数建立液体介质区域的边界元方程,然后考虑液—固界面的连续条件,获得水对固体介质的影响矩阵,最后用此影响矩阵修正固体介质的边界元方程。通过对简单问题的计算表明,采用本文方法求解流—固相互作用问题,其计算量与单独求解固体介质的波动问题时相当。     

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

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

有关高填充聚合物损伤粘弹介质的非线性本构方程与脱湿损伤模型

本文从连续损伤力学的角度,提出了复合固体推进剂损伤粘弹介质的非线性本构方程,并对单轴定应变速率拉伸情况提出了脱湿损伤的演变模型.理论计算与实验结果吻合较好.     

固体推进剂非线性拉伸行为的实验研究

固体推进剂非线性拉伸行为的实验研究王晓明，沈亚鹏（西安交通大学，７１００４９）关键词本构关系，粘弹性，非线性，推进剂１引言Ｈａｄｄａｄ（１）和杨挺青（２）对粘弹性非线性本构方程作了比较全面的综述，其中最简单，而又有利于实验确定材料常数的本构方程是其中...     

含RDX四组元HTPB固体推进剂的冲击起爆特性研究

为了研究含RDX四组元HTPB固体推进剂的冲击起爆行为和在低温条件下的适应性,在常温和低温条件下,对该固体推进剂进行了冲击加载拉氏分析实验。采用锰铜压力计测量了推进剂中不同位置处的压力变化历程,采用电离探针测量了固体推进剂的爆速。分析了固体推进剂的爆轰成长规律,获得了推进剂的临界起爆压力、爆速、爆压和爆轰成长距离等爆轰特征参量。通过对比不同条件下的特征参量发现,低温对固体推进剂的冲击起爆特性影响较小。此外,还对固体推进剂的冲击起爆过程进行了数值模拟,标定了固体推进剂点火增长模型的反应速率方程参数和推进剂的未反应JWL状态方程参数。     

纤维沥青混凝土FRA黏弹性预测模型

纤维沥青混凝土FRA(Fiber reinforced asphalt concrete)的力学行为取决于沥青混合料和纤维的物理属性及几何特征。通过三参数固体模型,考虑纤维的几何特性,构造了一种预测FRA黏弹性行为模型,推导出该模型的本构方程、蠕变柔量和松弛模量表达式。采用该模型研究了纤维体积分、比半径和弹性模量等对FRA松弛模量和蠕变柔量的影响。结果表明,纤维掺入量和比半径对FRA的松弛模量和蠕变柔量影响较大,而弹性模量影响较小。通过试验,验证了FRA预测模型具有较高的精度。在研究试验范围内,SMA-13沥青混凝土的聚丙烯腈纤维(PAN)最佳掺入量和长径比分别为0.3%～0.4%和600～900。     

基于光学诊断的某新型含铝推进剂燃烧特性分析

为探究某新型含铝固体推进剂燃烧特性和规律，在模拟固体发动机的高压条件下，采用可调功率激光器结合高速摄影、发射光谱等光学诊断技术对该新型含铝固体推进剂开展了系统的点火及燃烧过程研究。通过对该推进剂的点火延迟、退移速率、燃烧温度以及团聚物颗粒尺寸的定量测量和分析，明确了该推进剂的点火延迟量级；证实此推进剂的退移速率严格遵循Summerfield燃速公式；判断出其最高燃烧温度高于3 300 K，且随压力增大而升高；通过对燃烧过程中发光凝聚相产物面积的量化分析得出推进剂产物中团聚物粒径尺寸受环境参数的影响规律。     

Effect of the silicon-carbide micro- and nanoparticle size on the thermo-elastic and time-dependent creep response of a rotating Al–SiC composite cylinder

The history of stresses and creep strains of a rotating composite cylinder made of an aluminum matrix reinforced by silicon carbide particles is investigated. The effect of uniformly distributed SiC micro- and nanoparticles on the initial thermo-elastic and time-dependent creep deformation is studied. The material creep behavior is described by Sherby’s constitutive model where the creep parameters are functions of temperature and the particle sizes vary from 50 nm to 45.9 μm. Loading is composed of a temperature field due to outward steady-state heat conduction and an inertia body force due to cylinder rotation. Based on the equilibrium equation and also stress-strain and strain-displacement relations, a constitutive second-order differential equation for displacements with variable and time-dependent coefficients is obtained. By solving this differential equation together with the Prandtl–Reuss relation and the material creep constitutive model, the history of stresses and creep strains is obtained. It is found that the minimum effective stresses are reached in a material reinforced by uniformly distributed SiC particles with the volume fraction of 20% and particle size of 50 nm. It is also found that the effective and tangential stresses increase with time at the inner surface of the composite cylinder; however, their variation at the outer surface is insignificant.     

考虑收缩徐变影响的新旧混凝土叠合梁单元刚度矩阵改进

混凝土自身的收缩徐变会在新旧混凝土叠合梁中使应力重分布.为了计算重分布应力,首先推导以挠度表达的叠合梁非线性微分方程,然后通过求解该微分方程,引入位移形函数、刚度形函数和等效节点载荷形函数,最后得出新混凝土梁、旧混凝土梁和Goodman弹性夹层三合一的叠合梁改进型单元刚度矩阵和等效节点载荷,从而为收缩徐变影响下的混凝土的内力计算提供了一种有效的新方法.文中还进行了实例验证分析,并从中得出了一些有益的结论.     

Dynamic penny-shaped cracks in multilayer sandwich composites

A point force method is proposed for obtaining the dynamic elastic response of a multilayer sandwich composite in the presence of a penny-shaped crack under a harmonic loading. The sandwich composite is a multilayered solid whose lower half is the mirror image of the upper half with the center plane as the mirror. The crack is lying on the mirror plane of the composite. The solution of the mode I dynamic crack problem is formulated by integrating the Green’s function of a time-harmonic surface normal point force over the crack surface with an unknown point force distribution. The dual integral equations of the unknown point force distribution are established by considering the boundary conditions, which can be reduced to a Fredholm integral equation of the second kind. A complete solution of the crack problem under consideration can be obtained by solving this Fredholm integral equation. It will be shown that the results obtained by this approach are the same as some existing solutions.     

Fracture of a burning solid propellant attenuated by a crack-type cavity

A plane steady problem of fracture mechanics for a burning deformable solid propellant attenuated by a crack-type cavity with a burning surface is considered. The crack-type cavity is assumed to have tip zones with bonds between the faces, and mixed boundary conditions are imposed on the propellant charge boundary. The problem of equilibrium of a propellant charge containing a crack-type cavity reduces to solving a system of nonlinear singular integrodifferential equations with a Cauchy-type kernel. Based on the solution obtained, the normal and tangential forces in the tip zones of the crack-type cavity are found. Local conditions determining stability (safety) of the solid propellant burning regime are found.     

Influence of the cross-sectional shape of tensile bars on their creep rupture strength in a corrosive medium

The creep rupture strength of tensile bars in a corrosive medium is studied. The time to failure of the bars is analyzed as a function of their cross-sectional shape. It is shown that the influence of the corrosive medium on the creep rupture strength of the bars is determined by the diffusion of its elements into the bar material, resulting in a decrease in its creep rupture strength. The diffusion of the corrosive medium into a bar is studied using an approximate method for solving the diffusion equation taking into account the motion of the diffusion front. Bars with different cross-sectional shapes are considered using Rabotnov kinetic theory. It is shown that at the same tensile stress, the bar with a circular cross-section has minimum time to failure.     

尼龙6/蒙脱土复合材料蠕变行为研究

针对经典黏弹性模型不能很好分析黏弹性材料的蠕变行为问题,运用分数阶导数的类标准线性体模型与Prony级数模型研究了黏弹性材料尼龙6/蒙脱土复合材料的蠕变行为.采用原位聚合法制备了尼龙6/蒙脱土复合材料,在室温环境下对其进行蠕变实验.然后运用分数阶导数的类标准线性体模型和Prony级数模型对复合材料的蠕变实验数据进行分析...     

黏弹性板的大挠度蠕变屈曲

研究了具有初始小挠度受轴向压载黏弹性板的蠕变屈曲问题，在建立控制方程时，利用了von Karman非线性应变-位移关系，并考虑了初始挠度，用标准线性固体模型描述材料的黏弹性特性，在求解非线性积分方程时，利用梯形公式计算记忆积分式，将非线性积分方程化为非线性代数方程进行数值求解，得到了结构的蠕变变形过程，又将问题退化到小挠度情况进行研究，得到了挠度随时间扩展的解析解，分析了瞬时失稳临界载荷、持久临界载荷的物理意义，讨论了考虑几何非线性对黏弹性板蠕变屈曲的影响。     

Creep performance of PVC aged at temperature relatively close to glass transition temperature

In order to predict the mechanical performance of the polyvinyl chloride(PVC) at a high operating temperature,a series of short-term tensile creep tests(onetenth of the physical aging time) of the PVC are carried out at 63 C with a small constant stress by a dynamic mechanical analyzer(DMA).The Struik-Kohlrausch(SK) formula and Struik shifting methods are used to describe these creep data for various physical aging time.A new phenomenological model based on the multiple relaxation mechanisms of an amorphous polymer is developed to quantitatively characterize the SK parameters(the initial creep compliance,the characteristic retardation time,and the shape factor) determined by the aging time.It is shown that the momentary creep compliance curve of the PVC at 63 C can be very well fitted by the SK formula for each aging time.However,the SK parameters for the creep curves are not constant during the aging process at the elevated temperatures,and the evolution of these parameters and the creep rate versus aging time curves at the double logarithmic coordinates have shown a nonlinear phenomenon.Moreover,the creep master curves obtained by the superposition with the Struik shifting methods are unsatisfactory in such a case.Finally,the predicted results calculated from the present model incorporating with the SK formula are in excellent agreement with the creep experimental data for the PVC isothermally aged at the temperature relatively close to the glass transition temperature.     

基于时序参数法的BFRP蠕变行为预测

本文研究了基于时序参数法的玄武岩纤维增强复合材料(BFRP)板弯曲蠕变行为．为了提高模型理论预测的精度,文中将整个实验时长的数据进行了合理的细化分组,且每组时长尽可能相同,然后对细化分组的数据组分别进行直接外推法拟合,得到各数据组对应的拟合方程参数,最后获得了BFRP板弯曲蠕变行为的理论预测模型．对比模型预测值和实验值可以发现：模型预测值随着时间的增长,越来越接近于实验值,最后二者的相对误差几乎为零．该结果表明：文中建立的蠕变预测模型具有足够高的精度,可以应用于BFRP的蠕变行为的预测．