Thermo-damage-viscoelastic constitutive model of HTPB composite propellant

The main goal of this work is to establish a thermo-damage-viscoelastic model for Hydroxyl-Terminated Polybutadiene (HTPB) composite propellant based on the thermodynamic theory and elastic–viscoelastic correspondence principles. The model will also consider the influence of temperatures. The parameter α which represent the damage evolution rate and the material constants a, b are defined as exponential functions of temperature T, i.e. α(T), a(T) and b(T). Relaxation tests and uniaxial constant rate tensile tests are used to acquire the model parameters, and C(S) curves of different rates under the same temperature states are considered to be overlapped in this paper, while noncoincidence under the different temperature states. Then, uniaxial constant rate tensile tests and multi-step tensile-relaxation tests are used to verify the accuracy of the model. The results show that, the model is highly accurate in describing the mechanical property of HTPB under various loading conditions, but some drawback in describing the relaxation property inside of the HTPB’s nonlinear viscoelastic segments.     

Reduced order nonlinear aeroelasticity of swept composite wings using compressible indicial unsteady aerodynamics

Nonlinear dynamic aeroelasticity of composite wings in compressible flows is investigated. To provide a reasonable model for the problem, the composite wing is modeled as a thin walled beam (TWB) with circumferentially asymmetric stiffness layup configuration. The structural model considers nonlinear strain displacement relations and a number of non-classical effects, such as transverse shear and warping inhibition. Geometrically nonlinear terms of up to third order are retained in the formulation. Unsteady aerodynamic loads are calculated according to a compressible model, described by indicial function approximations in the time domain. The aeroelastic system of equations is augmented by the differential equations governing the aerodynamics lag states to derive the final explicit form of the coupled fluid-structure equations of motion. The final nonlinear governing aeroelastic system of equations is solved using the eigenvectors of the linear structural equations of motion to approximate the spatial variation of the corresponding degrees of freedom in the Ritz solution method. Direct time integrations of the nonlinear equations of motion representing the full aeroelastic system are conducted using the well-known Runge–Kutta method. A comprehensive insight is provided over the effect of parameters such as the lamination fiber angle and the sweep angle on the stability margins and the limit cycle oscillation behavior of the system. Integration of the interpolation method employed for the evaluation of compressible indicial functions at any Mach number in the subsonic compressible range to the derivation process of the third order nonlinear aeroelastic system of equations based on TWB theory is done for the first time. Results show that flutter speeds obtained by the incompressible unsteady aerodynamics are not conservative and as the backward sweep angle of the wing is increased, post-flutter aeroelastic response of the wing becomes more well-behaved.     

Three-dimensional instability problems for viscoelastic composite materials and structural members

Investigations of viscoelastic composite materials and structural members carried out using the TDLTSDB and initial-imperfection method are reviewed. The investigations address the internal and surface loss of stability of layered and fibrous composites, the loss of stability of plates, and the delamination (buckling) of plates with cracks. Each of these problems is reviewed separately. New areas of further research are proposed. The review focuses on the investigations carried out by the author and his students Published in Prikladnaya Mekhanika, Vol. 43, No. 10, pp. 3–27, October 2007.     

Equivalent nonlinear beam model for the 3-D analysis of shear-type buildings: Application to aeroelastic instability

Tower buildings can be very sensitive to dynamic actions and their dynamic analysis is usually carried out numerically through sophisticated finite element models. In this paper, an equivalent nonlinear one-dimensional shear–shear–torsional beam model immersed in a three-dimensional space is introduced to reproduce, in an approximate way, the dynamic behavior of tower buildings. It represents an extension of a linear beam model recently introduced by the authors, accounting for nonlinearities generated by the stretching of the columns. The constitutive law of the beam is identified from a discrete model of a 3D-frame, via a homogenization process, which accounts for the rotation of the floors around the tower axis. The macroscopic shear strain in the equivalent beam is produced by the bending of columns, accompanied by negligible rotation of the floors. A coupled nonlinear shear–torsional mechanical model is thus obtained. The coupling between shear and torsion is related to a non-symmetric layout of the columns, while mechanical nonlinearities are proportional to the slenderness of the columns. The model can be used for the analysis of the response of tower buildings to any kind of dynamic and static excitation. A first application is here presented to investigate the effect of mechanical and aerodynamic coupling on the critical galloping conditions and on the postcritical behavior of tower buildings, based on a quasi-steady model of aerodynamic forces.     

Multi-scale constitutive modeling of the small deformations of semi-crystalline polymers

A multi-scale constitutive model for the small deformations of semi-crystalline polymers such as high density Polyethylene is presented. Each macroscopic material point is supposed to be the center of a representative volume element which is an aggregate of randomly oriented composite inclusions. Each inclusion consists of a stack of parallel crystalline lamellae with their adjacent amorphous layers.Micro-mechanically based constitutive equations are developed for each phase. A viscoplastic model is used for the crystalline lamellae. A new nonlinear viscoelastic model for the amorphous phase behavior is proposed. The model takes into account the fact that the presence of crystallites confines the amorphous phase in extremely thin layers where the concentration of chain entanglements is very high. This gives rise to a stress contribution due to elastic distortion of the chains. It is shown that the introduction of chains’ elastic distortion can explain the viscoelastic behavior of crystalline polymers. The stress contribution from elastic stretching of the tie molecules linking the neighboring lamellae is also taken into account.Next, a constitutive model for a single inclusion considered as a laminated composite is proposed. The macroscopic stress-strain behavior for the whole RVE is found via a Sachs homogenization scheme (uniform stress throughout the material is assumed).Computational algorithms are developed based on fully implicit time-discretization schemes.     

On the predictive/fitting capabilities of the advanced differential constitutive equations for linear polyethylene melts

This contribution examines the capabilities of three differential constitutive models (XPP, PTT-XPP, and modified Leonov) in predicting rheological properties of two virtually linear polyethylene materials (HDPE Tipelin FS 450-26, mLLDPE Exact 0201) with specific attention to both steady as well as transient shear and uniaxial elongational flow situations. For each situation the (dis)advantages of the individual models are discussed and both, qualitative and quantitative model efficiency evaluation has been carried out.     

Multiscale micromechanical modeling of the constitutive response of carbon nanotube-reinforced structural adhesives

Appropriate formulations are developed to allow for the atomistic-based continuum modeling of nano-reinforced structural adhesives on the basis of a nanoscale representative volume element that accounts for the nonlinear behavior of its constituents; namely, the reinforcing carbon nanotube, the surrounding adhesive and their interface. The newly developed representative volume element is then used with analytical and computational micromechanical modeling techniques to investigate the homogeneous dispersion of the reinforcing element into the adhesive upon both the linear and nonlinear properties. Unlike our earlier work where the focus was on developing linear micromechanical models for the effective elastic properties of nanocomposites, the present approach extends these models by describing the development of a nonlinear hybrid Monte Carlo Finite Element model that allows for the prediction of the full constitutive response of the bulk composite under large deformations. The results indicate a substantial improvement in both the Young’s modulus and tensile strength of the nano-reinforced adhesives for the range of CNT concentrations considered.     

Drilling couples and refined constitutive equations in the resultant geometrically non-linear theory of elastic shells

It is well known that distribution of displacements through the shell thickness is non-linear, in general. We introduce a modified polar decomposition of shell deformation gradient and a vector of deviation from the linear displacement distribution. When strains are assumed to be small, this allows one to propose an explicit definition of the drilling couples which is proportional to tangential components of the deviation vector. The consistent second approximation to the complementary energy density of the geometrically non-linear theory of isotropic elastic shells is constructed. From differentiation of the density we obtain the consistently refined constitutive equations for 2D surface stretch and bending measures. These equations are then inverted for 2D stress resultants and stress couples. The second-order terms in these constitutive equations take consistent account of influence of undeformed midsurface curvatures. The drilling couples are explicitly expressed by the stress couples, undeformed midsurface curvatures, and amplitudes of quadratic part of displacement distribution through the thickness. The drilling couples are shown to be much smaller than the stress couples, and their influence on the stress and strain state of the shell is negligible. However, such very small drilling couples have to be admitted in non-linear analyses of irregular multi-shell structures, e.g. shells with branches, intersections, or technological junctions. In such shell problems six 2D couple resultants are required to preserve the structure of the resultant shell theory at the junctions during entire deformation process.     

Vibration analysis of hard-coated composite beam considering the strain dependent characteristic of coating material

The strain dependent characteristics of hard coat-ings make the vibration analysis of hard-coated composite structure become a challenging task. In this study, the mod-eling and the analysis method of a hard-coated composite beam was developed considering the strain dependent char-acteristics of coating material. Firstly, based on analyzing the properties of hard-coating material, a high order polynomial was adopted to characterize the strain dependent character-istics of coating materials. Then, the analytical model of a hard-coated composite beam was created by the energy method. Next, using the numerical method to solve the vibra-tion response and the resonance frequencies of the composite beam, a specific calculation flow was also proposed. Finally, a cantilever beam coated with MgO+Al2 O3 hard coating was chosen as the study case;under different excitation levels, the resonance region responses and the resonance frequencies of the composite beam were calculated using the proposed method. The calculation results were compared with the experiment and the linear calculation, and the correctness of the created model was verified. The study shows that compared with the general linear calculation, the proposed method can still maintain an acceptable precision when the excitation level is larger.     

Fundamentals of two-phase flow by the method of irreversible thermodynamics

The paper explores thermodynamic aspect of modelling two-phase systems by the methods of irreversible thermodynamics in both classical (CIT) and extended (EIT) formulation. The conservation laws for two-phase model-continuum are derived. Then, the entropy production is analysed for two-fluid and homogeneous systems. Different equations of state are taken into consideration, namely that corresponding to the accompanying equilibrium state of physical element and more complex resulting from EIT. Obtained expressions for rate of entropy production per unit volume allow to identify the dissipative mechanisms in the two-phase system and suggest the forms of phenomenological relations to be adopted in the constitutive equations.     

利用结构基因法研究非线性周期性板结构的等效力学性能

非线性周期性板结构是一类在智能复合材料领域具有巨大应用潜力的结构,因其构成材料的非线性特性,以及结构中经常包含增强纤维、肋板和空洞等复杂微结构导致的材料几何非线性,利用常规的有限元方法进行建模和分析较为困难.本文提出了一种结构基因法,通过提取非线性周期性板结构的最小模型单元作为其结构基因,将异质周期性板结构等效为均质板结构,便捷地求解了非线性周期性板结构的微观力学性能和整体等效力学性能.算例表明,结构基因方法可用来分析复杂非线性复合材料结构问题,计算结果精度足够,为复合材料微观力学研究提供了有价值的参考.     

CONSTITUTIVE RELATION OF UNSATURATED SOIL BY USE OF THE MIXTURE THEORY(Ⅰ)-NONLINEAR CONSTITUTIVE EQUATIONS AND FIELD EQUATIONS

The nonlinear constitutive equations and field equations of unsaturated soils were constructed on the basis of mixture theory. The soils were treated as the mixture composed of three constituents. First, from the researches of soil mechanics, some basic assumptions about the unsaturated soil mixture were made, and the entropy inequality of unsaturated soil mixture was derived. Then, with the common method usually used to deal with the constitutive problems in mixture theory, the nonlinear constitutive equations were obtained. Finally, putting the constitutive equations of constituents into the balance equations of momentum, the nonlinear field equations of constituents were set up. The balance equation of energy of unsaturated soil was also given, and thus the complete equations for solving the thermodynamic process of unsaturated soil was formed. Foundation items: the National Natural Science Foundation of China (59678003); Special Research Plan of the Education Department of Shaanxi Province (01JK178) Biographies: HUANG Yi (1936-) ZHANG Yin-ke (1964-)     

Structure of the constitutive relations for hereditarily elastic materials reinforced by hard fibers

The problem of simplifying the nonlinear hereditary elasticity relations is considered for strongly anisotropic materials such as fiber-reinforced composites. This is done using their property that the material stiffness is high along the reinforcement and is low in the cross-sectional direction. The material is considered transversally isotropic. The simplification is performed by analyzing asymptotic representations of creep relations. Relations of various degrees of accuracy for various types of composites and stress states are obtained. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 120–127, May–June, 2005.     

多梁式波形钢腹板工字钢组合梁荷载横向分布系数研究

为探讨多梁式波形钢腹板工字钢组合梁横向分布系数计算方法,同时考虑波形钢腹板抗扭刚度、剪切变形及钢混滑移效应,对传统的偏心压力法、修正偏心压力法及刚接梁法进行修正,并结合有限元模型通过一座典型的4主梁波形钢腹板工字钢组合梁桥对比分析了上述方法的适用性,最后基于参数分析研究了横隔板数量及刚度对其横向荷载分布系数的影响。结果表明,考虑剪切变形及滑移效应的刚接梁法得到的荷载横向分布系数与有限元值符合最好;当桥梁宽跨比小于2时,宜采用刚接梁法计算各主梁荷载横向分布系数,当宽跨比大于2时,宜采用更为简洁的修正偏心压力法进行计算;横隔板的设置可改善各主梁荷载横向均匀分布,但跨间横隔板间距和刚度对其荷载横向分布系数影响较小,实桥设计时可仅在端部及跨中位置布置横隔板,横隔板刚度及其余部位横隔板数量可根据结构稳定性要求进行布置。     

复合材料层合板自由振动分析的无网格自然邻接点Petrov-Galerkin法

基于一阶剪切变形理论,提出了复合材料层合板自由振动分析的无网格自然邻接点Petrov-Galerkin法。计算时在复合材料层合板中面上仅需要布置一系列的离散节点,并利用这些节点构建插值函数。在板中面上的局部多边形子域上,采用加权余量法建立复合材料层合板自由振动分析的离散化控制方程,并且这些子域可由Delaunay三角形方便创建。自然邻接点插值形函数具有Kronecker delta函数性质,因而无需经过特别处理就能准确地施加本质边界条件。对不同边界条件、不同跨厚比、不同材料参数和不同铺设角度的复合材料层合板,由本文提出的无网格自然邻接点Petrov-Galerkin法进行自由振动分析时均可得到满意的结果。数值算例结果表明,本文方法求解复合材料层合板的自由振动问题是行之有效的。     

钢-混凝土组合梁桥徐变的三维虚拟层合元分析

采用三参数粘弹性模型,对钢-混凝土组合梁桥在长期荷载作用下的徐变性能进行了研究,并采用三维虚拟层合单元法研究这种效应以及截面由不同材料组成时结构的受力特点。在实验研究的基础上,拟合了粘弹性模型中的参数,并对钢-混凝土组合梁桥在长期荷载作用下的竖向变形和轴向应变进行了研究,通过对计算和试验结果进行对比分析验证了该法的准确性。     

CONSTITUTIVE RELATION OF UNSATURATED SOIL BY USE OF THE MIXTURE THEORY (Ⅱ)-LINEAR CONSTITUTIVE EQUATIONS AND FIELD EQUATIONS

The linear constitutive equations and field equations of unsaturated soils were obtained through linearizing the nonlinear equations given in the first part of this work. The linear equations were expressed in the forms similar to Biot’s equations for saturated porous media. The Darcy’s laws of unsaturated soil were proved. It is shown that Biot’s equations of saturated porous media are the simplification of the theory. All these illustrate that constructing constitutive relation of unsaturated soil on the base of mixture theory is rational.     

Application of the coordinatewise descent method on a unit interval for weight optimization of structures made of composite materials

Problems of weight optimization of anisotropic structures with restrictions on strength, stability, etc., are solved by the coordinatewise descent method on a unit interval, combined with effective methods of nonlinear programming (dichotomy and golden section methods). __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 210–218, July–August, 2009.     

钱令希教授发现的梁与拱结构响应函数族关系矩阵

阐发了钱令希教授《梁与拱结构函数分布图与感应图之间的连锁关系》一文的学术思想,指出从本质上看,结构响应函数族关系矩阵是材料力学的三函数关系的横向拓展和向梁与拱结构的提升。可见,钱令希教授的工作是力学领域独特的基础研究之一。此外,还介绍了(1)采用Green函数法求解了钱教授没有研究的在梁上分布载荷的情况。(2)作者及其团队首次把这一工作写进了《材料力学》教科书中。(3)给出了应用本文方法的两个例题。