An aeroelastic model for composite rotor blades with straight and swept tips. Part I: Aeroelastic stability in hover

An analytical model for predicting the aeroelastic behavior of composite rotor blades with straight and swept tips is presented. The blade is modeled by beam type finite elements along the elastic axis. A single finite element is used to model the swept tip. The non-linear equations of motion for the finite element model are derived using Hamilton's principle and based on a moderate deflection theory and accounts for: arbitrary cross-sectional shape, pretwist, generally anisotropic material behavior, transverse shears and out-of-plane warping. Numerical results illustrating the effects of tip sweep, anhedral and composite ply orientation on blade aeroelastic behavior are presented. It is shown that composite ply orientation has a substantial effect on blade stability. At low thrust conditions, certain ply orientations can cause instability in the lag mode. The flap-torsion coupling associated with tip sweep can also induce aeroelastic instability in the blade. This instability can be removed by appropriate ply orientation in the composite construction.     

The cut-and-try method in determining the saturate spacing of transverse cracks of composite laminates

In this paper, the saturate spacing of transverse cracks of the 90° ply is originally calculated by the 3-D finite element method. Thus, a new approach is put forward for predicting the saturate spacing of transverse cracks.     

765THE CUT-AND-TRY METHOD IN DETERMINING THE SATURATE SPACING OF TRANSVERSE CRACKS OF COMPOSITE LAMINATES

I.IntroductionTransversecracksfirstoccurinthe9o"plyandspreadalongtheentirewidth.Withstaticorcyclicloadadding,thenumberoftransversecracksinthepliescontinuous1yincreases,lastinguntilthedistributionofcracksreachesasaturatestate.Theoccurrenceoftransversecrack…     

Transverse Stiffness of a Sinusoidally Corrugated Plate∗

Abstract

An analytical model for the initial transverse stiffness of a sinusoidally corrugated plate is derived, incorporating deformations due to extension, shear, and bending. A nondimensional plot is developed for determining transverse stiffness based on thickness and corrugation for a range of plate geometries. This model shows that for most corrugated plates the transverse stiffness is dramatically decreased from that of an uncomigated plate of the same thickness. For thin plates, a simple approximate polynomial expression for initial transverse stiffness is obtained. For thick plates with a small degree of corrugation, transverse stiffness is not negligible relative to longitudinal stiffness. The exact model is verified using a linear-elastic two-dimensional finite element model.     

Stress state characterization of delamination cracks in [0/90] symmetric laminates by BEM

The presence in a laminate of plies oriented at 90° with respect to the preferred direction of load generates almost immediately the appearance in these plies of cracks transverse to the load (parallel to the fibres in the lamina). These cracks reach the interface with the neighbouring ply, which in this paper will be considered oriented 0° with respect to the direction of the load. This may cause the bifurcation of the crack, which now appears propagating as a delamination crack between the two plies. The objective of this study is to characterize the stress state at the tip of both, the transverse crack in the 90° ply reaching the interface with the 0° ply, and the delamination crack for different lengths of the debonding. The analysis is performed by means of the Boundary Element Method allowing contact, without or with friction, to take place between the faces of the crack. The plies are considered as equivalent homogeneous bodies under a generalized plane strain state. The results are compared with those predicted by the open and contact models of Interfacial Fracture Mechanics. Accurate knowledge of the stress state at the neighbourhood of the tips of the cracks studied is necessary to generate failure criteria based on Fracture Mechanics parameters to predict the appearance and growth of the type of damage described.     

Computational analysis of multi-stepped beams and beams with linearly-varying heights implementing closed-form finite element formulation for multi-cracked beam elements

The model where the cracks are represented by means of internal hinges endowed with rotational springs has been shown to enable simple and effective representation of transversely-cracked slender Euler–Bernoulli beams subjected to small deflections. It, namely, provides reliable results when compared to detailed 2D and 3D models even if the basic linear moment–rotation constitutive law is adopted.This paper extends the utilisation of this model as it presents the derivation of a closed-form stiffness matrix and a load vector for slender multi-stepped beams and beams with linearly-varying heights. The principle of virtual work allows for the simple inclusion of an arbitrary number of transverse cracks. The derived at matrix and vector define an ‘exact’ finite element for the utilised simplified computational model. The presented element can be implemented for analysing multi-cracked beams by using just one finite element per structural beam member. The presented expressions for a stepped-beam are not exclusively limited to this kind of height variation, as by proper discretisation an arbitrary variation of a cross-section’s height can be adequately modelled.The accurate displacement functions presented for both types of considered beams complete the derivations. All the presented expressions can be easily utilised for achieving computationally-efficient and truthful analyses.     

Determination of the first ply failure load for a cross ply laminate subjected to uniaxial tension through computational micromechanics

A method for determining the in situ strength of fiber-reinforced laminas for three types of transverse loading including compression, tension and shear is presented. In the framework of this method, an analysis of local stresses that are responsible for the coalescence of matrix cracks is carried out by using a multi-fiber unit cell model and finite element method. The random distribution of fibers, fiber–matrix decohesion and matrix plastic deformations are taken into account in the micromechanical simulations. The present study also shows that the nonlinear hardening behavior of matrix reflects more realistically the influence of plastic deformations on the in situ transverse strength of lamina than the perfectly plastic behavior of matrix. The prediction of the in situ transverse strength is verified against the experimental data for a cross ply laminate subjected to uniaxial tension.     

任意斜裂纹转子的耦合振动研究

对包含不同类型裂纹(横裂纹、横-斜裂纹以及任意斜裂纹)的转子的耦合振动进行研究,以揭示裂纹转子在不同方向上刚度参数的变化规律及其交叉耦合机理,特别是由此引发的振动特征. 对于包含不同类型裂纹的转子轴段,采用六自由度Timoshenko梁单元模型对其进行单元建模,并基于应变能理论推导计算柔度参数和刚度矩阵. 在此基础上, 采用纽马克-$\beta$数值算法求解裂纹转子的运动方程,获得裂纹转子在单故障或多故障激励(不平衡激励、扭转激励或不平衡激励加扭转激励)作用下的耦合振动响应,进而分析耦合振动谱特征. 与横裂纹和横-斜裂纹相比,任意斜裂纹使转子刚度矩阵的交叉耦合效应更显著,导致转子发生更强烈的弯-扭耦合甚至是纵-弯-扭耦合振动.无论是在不平衡激励还是扭转激励作用下, 弯曲振动与扭转振动幅度都更大. 而且,包含不同类型裂纹的转子的耦合振动特征频率,例如旋转基频与二倍频、扭转激励频率及其边带成分的幅值,对裂纹面方向角具有不同的敏感性. 所得的这些研究结果,可以为转子裂纹的特征参数辨识与诊断提供理论依据.      

复合材料刚度性能表征的协同损伤力学模型

针对复合材料层合板的弥散型损伤,提出一个刚度性能表征的协同损伤力学模型. 该模型兼顾了微观物理损伤响应和宏观材料刚度性能表征. 从微观角度,建立细观RVE 模型求解裂纹表面张开位移和滑开位移,以此定义损伤张量,并在宏观上通过对材料应变和损伤表面位移进行均匀化处理,建立单向板或层合板的损伤刚度矩阵和损伤张量之间的联系. 以基体裂纹为例,详细分析并建立了横向裂纹和纵向裂纹的损伤本构. 计算了[±θ/90₄]_S 铺层层合板中基体横向裂纹对刚度性能的影响,结果表明该方法能够准确地预测复合材料层合板由损伤导致的刚度性能衰减.      

Numerical simulation of interlaminar damage propagation in CFRP cross-ply laminates under transverse loading

This paper proposes a numerical simulation of interlaminar damage propagation in FRP laminates under transverse loading, using the finite element method. First, we conducted drop-weight impact tests on CFRP cross-ply laminates. A ply crack was generated at the center of the lowermost ply, and then a butterfly-shaped interlaminar delamination was propagated at the 90/0 ply interface. Based on these experimental observations, we present a numerical simulation of interlaminar damage propagation, using a cohesive zone model to address the energy-based criterion for damage propagation. This simulation can address the interlaminar delamination with high accuracy by locating a fine mesh near the damage process zone, while maintaining computational efficiency with the use of automatic mesh generation. The simulated results of interlaminar delamination agreed well with the experiment results. Moreover, we demonstrated that the proposed method reduces the computational cost of the simulation.     

复合材料层板分层伴以横向裂纹扩展的三维有限元分析

本文采用三维有限元模型,对典型铺设的[0_2/±45_2/90_2]_s碳/环氧复合材料层板中分层伴以横向裂纹的产生和扩展导致的层间应力分布进行了分析。计算结果表明层间裂纹首先在90°层中部出现并开裂至相邻界面处而产生横向裂纹,横向裂纹的出现引起局部分层按三角形状扩展;并指出分层损伤过程是一个主导性的稳定扩展过程,这是导致刚度下降的主要因素。最后,数值计算结果与实验结果比较,两者是吻合的。     

NUMERICAL SIMULATION OF CRACK GROWTH OF REINFORCED CONCRETE FLEXURAL MEMBERS BASED ON XFEM1)

基于最大主应力准则,利用扩展有限元法对钢筋混凝土三点弯曲梁的裂纹扩展进行数值模拟。通过与已有模型试验结果的比对,验证了该方法的准确性与有效性。研究发现,带裂纹钢筋混凝土构件的裂纹扩展与力学性能不但受到单裂纹的尺寸、位置、倾角的影响,而且还受多裂纹间位置关系的影响;其中,裂纹横向分布位置和裂纹纵向分布位置是影响带裂纹钢筋混凝土梁裂缝扩展的主要因素。     

FRP-混凝土界面粘结行为的参数影响研究

FRP-混凝土界面的粘结性能对FRP加固混凝土结构力学行为和破坏模式有着重要影响。本文对表征FRP-混凝土界面粘结性能的三个重要参数(界面初始刚度、最大剪应力、界面破坏能)开展研究,通过13个单剪试件的试验考察了混凝土强度、胶层厚度和粘结长度等因素对界面粘结行为的影响,根据试验结果拟合了界面破坏能、最大剪切应力与胶层剪切刚度、混凝土强度之间的函数关系。在试验研究基础上,构建了外贴FRP-混凝土界面粘结的有限元模型。通过有限元分析考察了界面破坏能等三个参数不变的前提下,不同的局部粘结滑移本构关系对界面粘结行为的影响;进而研究了其中一个参数变化时引起的界面粘结性能改变。研究结果表明:界面粘结承载力随着胶层厚度增加而逐渐提高;胶层厚度与界面破坏能成正比,与峰值剪应力成反比;当界面破坏能等三个参数保持不变时,局部粘结滑移本构关系对FRP-混凝土界面粘结性能的影响较小;三个参数中的一个增大时将延缓界面破坏的过程。     

考虑粘结面厚度的局部粘结单元法

传统无厚度粘结单元法CFEM (Cohesive finite element method)在模拟脆性材料断裂方面具有很强的优势,但也存在很大问题.一是单元尺寸增大,收敛性变差;二是单元尺寸变小,模型刚度发生折减.为了克服这两个问题,发展了考虑厚度的局部粘结单元法,即在裂纹可能扩展区插入具有一定厚度的粘结面单元.粘结面单元采用拓展虚内键本构(Augmented virtual internal bond)描述.由于考虑了厚度,粘结面交叉处会形成多边形空缺.为了弥补这一空缺,将其看作多边形键元胞,采用离散虚内键模型(Discretized virtual internal bond)对其建模,保证了模型的几何完整性.模拟结果表明,本文方法有效,克服了传统CFEM方法的刚度折减问题,提高了计算稳定性和收敛性.     

A representative volume element based on translational symmetries for FE analysis of cracked laminates with two arrays of cracks

A methodology is proposed for the construction of a representative volume element (RVE) for analysis of laminated composites containing two arrays of ply cracks running in different directions. The only requirement is that the cracks in any ply are uniformly spaced, and if more than one ply of a given orientation is cracked, then the crack spacing of individual plies must only be in exact multiples of each other. The spacing of cracks in the two directions can be fully independent. The RVE is constructed through a systematic consideration of translational symmetries present in the cracked laminate. As a result, the boundary conditions on the RVE can be imposed without compromising accuracy. Examples of the application of the RVE methodology are given to illustrate its broad capability and a finite element (FE) stress analysis is performed for these cases to illustrate results such as the crack surface displacements, local stress fields and RVE-averaged elastic properties. For one case, the average properties are compared with experimental results, showing good agreement.     

The effect of transverse normal strain in contact of an orthotropic beam pressed against a circular surface

The contact problem of a straight orthotropic beam pressed onto a rigid circular surface is considered using beam theories that account for transverse shear and transverse normal deformations. The circular nature of the rigid surface emphasizes the difference between Euler Bernoulli theory behavior, where point loads develop at the edge of contact, and the higher order theories that predict non-singular pressure distributions. While Timoshenko beam theory is the simplest theory that addresses this behavior, the prediction of a maximum value of pressure at the edge of contact contradicts the elasticity theory result that contact pressure must drop to zero. Transverse normal strain is therefore introduced, both to study this fundamental discrepancy and to include an important effect in many contact problems. To investigate this effect, higher order beam theories that account for both constant and linear transverse normal strain through the beam thickness are derived using the principle of virtual work. The resulting orthotropic beam theories depend on the bending stiffness (EI), shear stiffness (GA), axial stiffness (EA₁) and transverse normal stiffness (EA₂), which are independent stiffness parameters that can differ by orders of magnitude. The above mentioned contact problem is then solved analytically for these theories, along with the Timoshenko beam model which assumes zero transverse normal strain. The results for different orthotropic materials show that inclusion of transverse normal deformation has a significant effect on the contact pressure solution. Furthermore, the solution using higher order beam theories encompasses the two extremes of a Hertz-like contact pressure when the half contact length is smaller than the thickness of the beam, and the Timoshenko beam theory case when the half contact length is much larger than the thickness. Concerning the behavior of the pressure at the edge of contact, adherence to the boundary conditions required by the principle of virtual work, shows that while the pressure does tend to zero, it does not become zero unless artificially enforced. In this regard the solution for the case of linear strain is better than that for constant strain. All beam solutions are validated with plane elasticity solutions obtained using the commercial finite element software ABAQUS.     

Mechanical Testing of Solid–Solid Interfaces at the Microscale

In order to determine the influence of internal interfaces on the material’s global mechanical behavior, the strength of single interfaces is of great interest. The experimental framework presented here enables quantitative measurements of the initiation and propagation of interfacial cracks on the microscale. Cantilever beams are fabricated by focused ion beam milling out of a bulk sample, with an interface of interest placed close to the fixed end of the cantilever. Additionally, a U-notch is fabricated at the location of the interface to serve as a stress concentrator for the initiation of the crack. The cantilevers are then mechanically deflected using a nanoindentation system for high resolution load-displacement measurements. In order to determine the onset and propagation of damage, the stiffness of the cantilevers is recorded by partial unloads during the test as well as by making use of a continuous stiffness technique. A finite element model is used to normalize the load and stiffness in order to establish the framework for comparisons between different interfaces.     

Effects of interface contacts on the magneto electro-elastic coupling for fiber reinforced composites

Imperfect bonding between constituents is studied where displacements, electric and magnetic static potentials are considered to have a jump proportional to the normal component of the mechanical traction, electric displacement and magnetic flux. This condition may model various interface damages or the thin glue layer between two adjacent phases. They are termed as the mechanically compliant, dielectrically weakly capacitance and magnetically weakly inductance at the interface. It is shown that while the more imperfect the interface is, the overall properties become weaker, such as longitudinal shear stiffness, out-of-plane piezoelectric coupling, and magnetoelectric coupling. Out-of-plane piezomagnetic coupling, transverse dielectric permittivity and transverse dielectric permeability exhibit no influence by imperfect bonding. The imperfect interface proposed is mimicked by the springs, capacitors and inductances for the mechanical, electric and magnetic interaction between the phases and are highly sensitive to the interphase properties. The results are compared mainly with the self consistent model reported in the literature and good agreements are shown.     

A new super convergent thin walled composite beam element for analysis of box beam structures

In this paper, a new composite thin wall beam element of arbitrary cross-section with open or closed contour is developed. The formulation incorporates the effect of elastic coupling, restrained warping, transverse shear deformation associated with thin walled composite structures. A first order shear deformation theory is considered with the beam deformation expressed in terms of axial, spanwise and chordwise bending, corresponding shears and twist. The formulated locking free element uses higher order interpolating polynomial obtained by solving static part of the coupled governing differential equations. The formulated element has super convergent properties as it gives the exact elemental stiffness matrix. Static and free vibration analyses are performed for various beam configuration and compared with experimental and numerical results available in current literature. Good correlation is observed in all cases with extremely small system size. The formulated element is used to study the wave propagation behavior in box beams subjected to high frequency loading such as impact. Simultaneous existence of various propagating modes are graphically captured. Here the effect of transverse shear on wave propagation characteristics in axial and transverse directions are investigated for different ply layup sequences.