Thermal-induced snap-through buckling of simply-supported functionally graded beams

The instability of functionally graded material (FGM) structures is one of the major threats to their service safety in engineering applications. This paper aims to clarify a long-standing controversy on the thermal instability type of simply-supported FGM beams. First, based on the Euler-Bernoulli beam theory and von Kármán geometric nonlinearity, a nonlinear governing equation of simply-supported FGM beams under uniform thermal loads by Zhang's two-variable method is formulated. Second, an approximate analytic solution to the nonlinear integro-differential boundary value problem under a thermal-induced inhomogeneous force boundary condition is obtained by using a semiinverse method when the coordinate axis is relocated to the bending axis (physical neutral plane), and then the analytical predictions are verified by the differential quadrature method (DQM). Finally, based on the free energy theorem, it is revealed that the symmetry breaking caused by the material inhomogeneity can make the simply-supported FGM beam under uniform thermal loads occur snap-through postbuckling only in odd modes; furthermore, the nonlinear critical load of thermal buckling varies non-monotonically with the functional gradient index due to the stretching-bending coupling effect. These results are expected to provide new ideas and references for the design and regulation of FGM structures.     

Geometrical nonlinear analysis of thin-walled composite beams using finite element method based on first order shear deformation theory

Based on a seven-degree-of-freedom shear deformable beam model, a geometrical nonlinear analysis of thin-walled composite beams with arbitrary lay-ups under various types of loads is presented. This model accounts for all the structural coupling coming from both material anisotropy and geometric nonlinearity. The general nonlinear governing equations are derived and solved by means of an incremental Newton–Raphson method. A displacement-based one-dimensional finite element model that accounts for the geometric nonlinearity in the von Kármán sense is developed to solve the problem. Numerical results are obtained for thin-walled composite beam under vertical load to investigate the effects of fiber orientation, geometric nonlinearity, and shear deformation on the axial–flexural–torsional response.     

Quasi-static and dynamical analyses of a thermoviscoelastic Timoshenko beam using the differential quadrature method

The quasi-static and dynamic responses of a thermoviscoelastic Timoshenko beam subject to thermal loads are analyzed. First, based on the small geometric deformation assumption and Boltzmann constitutive relation, the governing equations for the beam are presented. Second, an extended differential quadrature method(DQM)in the spatial domain and a differential method in the temporal domain are combined to transform the integro-partial-differential governing equations into the ordinary differential equations. Third, the accuracy of the present discrete method is verified by elastic/viscoelastic examples, and the effects of thermal load parameters, material and geometrical parameters on the quasi-static and dynamic responses of the beam are discussed. Numerical results show that the thermal function parameter has a great effect on quasi-static and dynamic responses of the beam. Compared with the thermal relaxation time, the initial vibrational responses of the beam are more sensitive to the mechanical relaxation time of the thermoviscoelastic material.     

Analytical solution for a strained reinforcement layer bonded to lip-shaped crack under remote mode III uniform load and concentrated load

In this paper, the analytical solution of stress field for a strained reinforcement layer bonded to a lip-shaped crack under a remote mode III uniform load and a concentrated load is obtained explicitly in the series form by using the technical of conformal mapping and the method of analytic continuation. The effects of material combinations, bond of interface and geometric configurations on interfaciai stresses generated by eigenstrain, remote load and concentrated load are studied. The results show that the stress concentration and interfaciai stresses can be reduced by rational material combinations and geometric configurations designs for different load forms.     

缝合式C/SiC复合材料非线性本构关系及断裂行为研究

C/SiC复合材料具有高比强度、高比模量和优良的热稳定性能等一系列优点, 广泛应用于航空航天领域中. 裂纹扩展进而引起的脆性断裂是其主要失效形式之一, 因而材料的断裂性能分析对材料的结构设计和应用有重要的指导意义. 本文开展了缝合式C/SiC复合材料简单力学试验和断裂试验, 研究了材料在不同载荷下的力学响应及断裂特征. 基于缝合式C/SiC复合材料简单力学试验, 建立了材料宏观非线性损伤本构方程, 并模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 本构方程采用简单函数描述了材料在复杂应力状态下的非线性应力-应变曲线, 并考虑了反向加载过程中造成的裂纹闭合. 基于商业有限元软件ABAQUS, 通过编写UMAT子程序实现非线性损伤本构方程, 采用单个单元验证了建立的本构方程的有效性. 在此基础上, 采用线弹性损伤本构和非线性损伤本构分别模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 采用非线性损伤本构方程模拟的力-位移曲线结果与试验结果更为吻合, 非线性损伤本构预测的失效载荷与试验失效载荷更为接近, 验证了所建立的非线性损伤本构方程的准确性, 为C/SiC复合材料断裂行为的研究提供了借鉴, 为缝合式C/SiC复合材料结构的设计和应用提供了理论基础.      

STUDY ON NONLINEAR CONSTITUTIVE RELATIONSHIP AND FRACTURE BEHAVIOR OF STITCHED C/SiC COMPOSITES$^{\bf 1)}$

C/SiC复合材料具有高比强度、高比模量和优良的热稳定性能等一系列优点, 广泛应用于航空航天领域中. 裂纹扩展进而引起的脆性断裂是其主要失效形式之一, 因而材料的断裂性能分析对材料的结构设计和应用有重要的指导意义. 本文开展了缝合式C/SiC复合材料简单力学试验和断裂试验, 研究了材料在不同载荷下的力学响应及断裂特征. 基于缝合式C/SiC复合材料简单力学试验, 建立了材料宏观非线性损伤本构方程, 并模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 本构方程采用简单函数描述了材料在复杂应力状态下的非线性应力-应变曲线, 并考虑了反向加载过程中造成的裂纹闭合. 基于商业有限元软件ABAQUS, 通过编写UMAT子程序实现非线性损伤本构方程, 采用单个单元验证了建立的本构方程的有效性. 在此基础上, 采用线弹性损伤本构和非线性损伤本构分别模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 采用非线性损伤本构方程模拟的力-位移曲线结果与试验结果更为吻合, 非线性损伤本构预测的失效载荷与试验失效载荷更为接近, 验证了所建立的非线性损伤本构方程的准确性, 为C/SiC复合材料断裂行为的研究提供了借鉴, 为缝合式C/SiC复合材料结构的设计和应用提供了理论基础.     

热载荷作用下梯度多孔材料梁弯曲及过屈曲力学行为的研究

基于Bernoulli-Euler梁理论,引入物理中面解耦了复合材料结构的面内变形与横向弯曲特性,研究了梯度多孔材料矩形截面梁在热载荷作用下的弯曲及过屈曲力学行为．假设沿梁厚度方向材料的性质是连续变化的,利用能量法推导了矩形截面梁的控制微分方程和边界条件,并用打靶法对无量纲化的控制方程进行数值求解．利用计算得到的结果分析了材料的性质、热载荷、边界条件对矩形截面梁非线性力学行为的影响．结果表明,对称材料模型下,固支梁与简支梁均显示出了典型的分支屈曲行为特征,而其临界屈曲热载荷值均会随着孔隙率系数的增加而单调增加．非对称材料模型下,固支梁仍显示出分支屈曲行为特征,但其临界屈曲热载荷不再随着孔隙率系数的变化而单调变化;而对于两端简支梁,发生了弯曲变形,弯曲挠度随载荷的增大而增大．     

Nonlinear stability and vibration of pre/post-buckled microstructure-dependent FGPM actuators

The present study deals with the study of the nonlinear stability and small free vibration of microstructure-dependent functionally graded piezoelectric material (FGPM) beams in pre/post-buckling regimes. The Timoshenko beam theory with various inplane and out-of-plane boundary conditions are considered under different types of mechanical and thermal loads. The beam is assumed to be under inplane mechanical, thermal, and electrical excitations. Each thermo-electro-mechanical property of the beam is graded across the thickness (i.e., height) of the beam, based on a power law model. The von Kármán type geometric nonlinearity is included to account for the large deflection behavior of the beam under inplane loads. The modified couple stress theory is included to account for the size effects. A weak-form, displacement-based, finite element formulation is developed to discretize the equations of motion. The resulting system of nonlinear algebraic equations is solved using Newton’s iterative method. The numerical results of frequencies and lateral deflections as a function of load parameters reveal the existence of bifurcation or critical states in some cases. The effects of load type, microstructural dependency, boundary conditions, beam geometry, composition rule of the constituents, and actuator voltage are discussed through the various parametric studies.     

Effective elastic properties of flexible chiral honeycomb cores including geometrically nonlinear effects

Flexible chiral honeycomb cores generally exhibit nonlinear elastic properties in response to large geometric deformation, which are suited for the design of morphing aerospace structures. However, owing to their complex structure, it is standard to replace the actual core structure with a homogenized core material presenting reasonably equivalent elastic properties in an effort to increase the speed and efficiency of analyzing the mechanical properties of chiral honeycomb sandwich structures. As such, a convenient and efficient method is required to evaluate the effective elastic properties of flexible chiral honeycomb cores under conditions of large deformation. The present work develops an analytical expression for the effective elastic modulus based on a deformable cantilever beam under large deformation. Firstly, Euler–Bernoulli beam theory and micropolar theory are used to analyze the deformation characteristics of chiral honeycombs, and to calculate the effective elastic modulus under small deformation. On that basis, the expression for the effective elastic modulus is improved by including the stretching deformation of the chiral honeycomb structure for a unit cell under conditions of large deformation. The effective elastic moduli calculated by the respective analytical expressions are compared with the results of finite element analysis. The results indicate that the analytical expression obtained under consideration of the geometric nonlinearity is more suitable than the linear expressions for flexible chiral honeycomb cores under conditions of high strain and low elastic modulus.     

平面梁杆结构几何非线性分析的一种简便方法

本文提出了一种新的几何非线性分析方法，适用于结点位移任意大，单元刚体转角任意大、单元局部弯曲比较小的平面梁杆结构。文中的刚度矩阵和附加荷载列阵都是以显式形式给出的，可直接应用。     

Analytical solution for a strained reinforcement layer bonded to lip-shaped crack under remote mode Ⅲ uniform load and concentrated load

In this paper,the analytical solution of stress field for a strained reinforcement layer bonded to a lip-shaped crack under a remote mode Ⅲ uniform load and a concentrated load is obtained explicitly in the series form by using the technical of conformal mapping and the method of analytic continuation.The effects of material combinations,bond of interface and geometric configurations on interfacial stresses generated by eigenstrain,remote load and concentrated load are studied.The results show that the stress concentration and interfacial stresses can be reduced by rational material combinations and geometric configurations designs for different load forms.     

Hamilton体系下功能梯度梁的热冲击动力屈曲分析

在Hamilton体系下,基于Euler梁理论研究了功能梯度材料梁受热冲击载荷作用时的动力屈曲问题;将非均匀功能梯度复合材料的物性参数假设为厚度坐标的幂函数形式,采用Laplace变换法和幂级数法解析求得热冲击下功能梯度梁内的动态温度场:首先将功能梯度梁的屈曲问题归结为辛空间中系统的零本征值问题,梁的屈曲载荷与屈曲模态分别对应于Hamilton体系下的辛本征值和本征解问题,由分叉条件求得屈曲模态和屈曲热轴力,根据屈曲热轴力求解临界屈曲升温载荷。给出了热冲击载荷作用下一类非均匀梯度材料梁屈曲特性的辛方法研究过程,讨论了材料的梯度特性、结构几何参数和热冲击载荷参数对临界温度的影响。     

Mathematical programming approaches for the safety assessment of semirigid elastoplastic frames

This paper presents two complementary mathematical programming based approaches for the accurate safety assessment of semirigid elastoplastic frames under quasistatic loads. The inelastic behavior of the flexible connections and material plasticity are accommodated through piecewise linearized nonlinear yield surfaces. As is necessary for this class of structures, geometric nonlinearity is taken into account. Moreover, only a 2nd-order geometric approximation is included as this is sufficiently accurate for practical structures. The work described has a twofold contribution. First, we develop an algorithm that can robustly and efficiently process the complete (path-dependent) nonholonomic response of the structure in a stepwise (path-independent) holonomic fashion. The governing formulation is cast in mixed static-kinematic variables and leads naturally to what is known in the mathematical programming literature as a mixed complementarity problem (MCP). The novelty of the proposed algorithm is that it processes the MCP directly without using some iterative (and often cumbersome) predictor–corrector procedure. Second, in the spirit of simplified analyses, the classical limit analysis approach is extended to compute the limit load multiplier under the simultaneous influence of joint flexibility, material and geometric nonlinearities, and limited ductility. Our formulation is an instance of the challenging class of optimization problems known as a mathematical program with equilibrium constraints (MPEC). Various nonlinear programming based algorithms are proposed to solve the MPEC. Finally, four numerical examples, concerning practical structures and benchmark cases, are provided to illustrate application of the analyses as well as to validate the accuracy and robustness of the proposed schemes.     

Shock wave method for monitoring crack repair processes in reinforced concrete structures

A nondestructive method for monitoring the crack state in reinforced concrete structures based on the recoding of wave processes in these structures under shock actions is proposed. The essence of the method and its possibilities are demonstrated by an example of the study of the behavior of a reinforced concrete beam with a crack at various stages of crack development and repair. Numerical simulation was used to study variations in the wave front characteristics in the crack area. A quantitative criterion was formulated, which permits estimating the concrete integrity or the existence of crack in it and monitoring the variations in the crack state in the process of loading the structure and the crack repair. The criterion is determined as the ratio of the amplitudes of the first half-waves of the acceleration wave front registered in regions on the opposite shores of the crack. The criterion value is independent of the amplitude of the shock action and the beam fixation conditions and is solely determined by the mechanical state of the material used to repair the crack. The criterion adequacy was demonstrated by comparing the results of numerical simulation with experimental data. A cycle of numerical experiments were carried out, which, for each duration of the shock action, permits determining the optimal values of the distance between the pulse application point and the acceleration recording points at which the criterion is most sensitive to the crack state.     

On a probabilistic stability theory for imperfection sensitive structures

The concept of almost sure sample stability and sample stability in probability are formulated for elastic systems. Using a Koiter type approach these concepts are used in the analysis of imperfection sensitive structures. The applied load and the initial geometric imperfections are introduced into the analysis as random quantities. A compressed beam of finite length on a nonlinear elastic foundation is used in an example calculation.     

随机梁式结构静力损伤识别的一种改进方法

针对已有的损伤识别方法会出现损伤误识别的问题,本文在已有方法的基础上发展了一种随机梁式结构静力损伤识别的改进方法。假定静力荷载下梁式结构初始模型参数（如弹性模量和几何尺寸等）及测量误差为随机量,给出已有的基于随机有限元模型的梁式结构静力损伤识别方法,并进一步提出了一种改进方法。该方法通过设定损伤概率指标的阈值和反复迭代对结构损伤识别进行改进。数值算例和简支梁静力试验表明,考虑初始模型的不确定性以及静力响应测量误差,本文方法相较已有方法可以更有效地识别梁式结构的损伤。     

钢管混凝土柱-钢梁节点的力学性能分析

基于弹塑性有限元理论建立了钢管混凝土柱-钢梁节点荷载-位移全过程非线性有限元模型,在单元分析中采用改进的AUL表述推导得到梁柱单元刚度矩阵方程,同时考虑了材料的物理非线性和单元的几何非线性,并编制了非线性有限元程序NLFEACFST。采用该模型对相关研究者和作者进行的节点试验进行了分析,理论计算结果与试验结果比较表明,该模型具有很好的适用性和精度。在理论分析模型得到试验结果验证的基础上,对典型的中柱节点进行了荷载-位移全过程非线性特性分析,并对影响节点承载力和荷载-位移骨架曲线的因素进行了参数分析,为进一步从理论研究钢管混凝土框架结构的力学性能创造了条件。     

远距离爆炸荷载作用下钢框架几何相似律研究

根据Π定理推导了远距离爆炸荷载作用下钢框架原型结构与缩比模型的几何相似律表达式。基于已有的钢框架子结构爆炸实验,采用AUTODYN建立了钢框架子结构数值模型,验证了流固耦合方法在结构爆炸响应分析中的可靠性。在此基础上,对比了流固耦合方法和解析爆炸边界方法在钢框架远距离爆炸数值模拟中的准确性和计算效率,结果表明,解析爆炸边界方法可以合理地模拟远距离爆炸荷载作用下钢框架的动态响应,且计算效率较高。最后,采用该方法分析了具有不同相似比的两层三跨钢框架结构在远距离爆炸荷载作用下的动态响应及毁伤效应,结果表明:该结构的动态响应和毁伤效应符合几何相似规律。     

Uniaxial tension of a plate with a circular hole at finite strain

Some kinematic characteristics of a geometrically nonlinear deformation process obtained experimentally at uniaxial tension of a plate with a circular hole are given. The plate is made of a highly elastic material. The geometric linearization method is used to determine the numerical characteristics of the deformation process under study.     

Stress analysis of restrained bending for box beam with rectangular cross-section filled by honeycomb core

On the basis about studying free bending for box beam with rectangular crosssection filled by honeycomb core, supplementary displacements and stresses of restrained bending for such beam were analyzed. The hypothesis for separated variables was adopted to solve displacement. According to this, three aspect equations of geometrical, physical and balance were obtained. With Galerkin‘s method, it is summed up as two-order ordinary differential equations with the attenuation character. Analysis makes clear that attenuation speed of stress is concerned with a big load or a small one, geometric dimensions of crosssection of beam, and physical parameter of material.