横向剪切对双模量复合材料叠层矩形板非线性弯曲的影响

本文探讨了动力松弛(DR)法在双模量复合材料叠层矩形板非线性弯曲问题中的应用。在分析中分别采用叠层板大挠度经典理论和计及大转动(在Von Karman意义上)的复合材料叠层板剪切变形理论。我们发现,对于考虑横向剪切变形的非线性弯曲问题,如何计算虚拟密度以控制数值计算的稳定性,仍然需要进一步研究。本文提出了一种虚拟密度的计算方法,从而保证了本课题数值计算的稳定性。文中介绍了用DR法求解双模量复合材料叠层板非线性弯曲的主要步骤,给出了由轻度双模量材料(Born-Epxy(B-E))和高度双模量材料(Aramid-Rubber(A-R)和Polyester-Rubber(P-R))的两层正交叠层简支矩形板在正弦分布载荷及均布载荷作用下的非线性弯曲特性的数值结果。将所得结果和小挠度分析结果及普通复合材料的结果作了比较,并分析了横向剪切变形对无量纲中心挠度的影响。     

Natural vibrations of a strip made of a composite material with a locally curved structure

The possibility of using the finite-element method for investigating two-dimensional problems on natural vibrations in the mechanics of composite materials with curved structures is considered. With the example of a hinge-supported strip made of a composite material with a locally curved structure, the influence of geometrical and mechanical parameters of the strip on its eigenfrequencies is examined. It is established that the presence of local curving in the structure of strip material decreases the magnitude of eigenfrequencies.Translated from Mekhanika Kompozitnykh Materialov, Vol. 41, No. 1, pp. 71–78, January–February, 2005.     

碳纤维/树脂基复合材料曲壁蜂窝夹芯结构的三点弯曲性能北大核心CSCD

为研究碳纤维复合材料(CFRP)曲壁蜂窝结构在三点弯曲载荷作用下的承载特性与失效模式,对不同芯层高度、面板厚度的结构进行了理论预报、数值模拟及试验.首先,根据夹芯结构的主要失效模式,提出了相应的理论预报公式,并绘制了失效机制图;其次,建立了CFRP曲壁蜂窝夹芯结构的有限元仿真模型,对其在三点弯曲载荷作用下的典型失效行为进行模拟;最后,通过模压成型工艺制备了不同尺寸的CFRP曲壁蜂窝夹芯结构,并将试验结果与理论、模拟结果进行比较.结果表明,蜂窝夹芯结构承载能力与芯层高度、面板厚度密切相关,结构芯层及面板刚度随其尺寸的减小而下降,导致结构失效模式由面芯脱黏失效变为面板压溃失效.     

Vibrations of Rotors Supported on Nonlinear Bearings

Rotors in electrical machines are supported by various types of bearings. In general, the rotor bearings have nonlinear stiffness properties and they influence the rotor vibrations significantly. In this work, this influence of these nonlinearities is investigated. A simplified finite element model using Timoshenko beam elements is set up for the heterogeneous structure of the rotor. A transversally isotropic material model is adopted for the rotor core stack. Imposing the nonlinear bearing stiffnesses on the model, the Newton-Raphson procedure is used to carry out a run up simulation. The spectral content of these results shows nonlinear effects due to the bearings. The rotor vibrations are further investigated in detail for various constant speeds. These results show non-harmonic vibrations of the rotor in a section of the investigated speed range. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Sound absorption of a quadratic and cubic nonlinearly vibrating curved panel absorber

This study investigated the sound absorption of a nonlinearly vibrating curved panel backed by a cavity. Very few studies on similar nonlinear structural-acoustic problems have been conducted to date, although there have been many on nonlinear plate or linear structural-acoustic problems. A curved panel is considered because the overall absorption bandwidth can be designed by appropriately adjusting the panel curvature, which is a key factor in controlling the structural resonant frequencies and absorption peaks. The theoretical formulation is developed based on the assumptions of quadratic and cubic nonlinear structural vibrations, and the linear acoustic pressure induced within the cavity. An approach based on the numerical integration method is developed to solve the nonlinear governing equation of the structural-acoustic problem. In the parametric study, the panel displacement amplitude converges with an increasing number of modes. The effects of excitation level, cavity depth, and damping factor are also examined. The quadratic and cubic nonlinearities and their effects on the sound absorption are also investigated. An experiment was conducted. The theoretical and experimental observations correspond reasonably with each other.     

A nonlocal higher-order model including thickness stretching effect for bending and buckling of curved nanobeams

An analytical approach for static bending and buckling analyses of curved nanobeams using the differential constitutive law, consequent to Eringen’s strain-driven integral model coupled with a higher-order shear deformation accounting for through thickness stretching is presented. The formulation is general in the sense that it can be deduced to examine the influence of different structural theories, for static and dynamic analyses of curved nanobeams. The governing equations derived using Hamiltons principle are solved in conjunction with Naviers solutions. The formulation is validated considering problems for which solutions are available. A comparative study is made here by various theories obtained through the formulation. The effects various structural parameters such as thickness ratio, beam length, rise of the curved beam, and nonlocal scale parameter are brought out on bending and stability characteristics of curved nanobeams.     

Nonlinear thermoelastic analysis of layered structures

The nonlinear thermoelastic behavior of orthotropic layered slabs and cylinders including radiation boundaries, temperature-dependent material properties, and stress-dependent layer interface conditions is investigated. A one-dimensional finite element formulation employing quadratic layer and linear interface elements is used to perform the analyses. The transient heat conduction portion of the program is temporally discretized via an implicit linear time interpolation algorithm which includes Crank-Nicolson, Galerkin, and Euler backward differencing. The nonlinear heat conduction equations are iteratively evaluated using a modified Newton-Raphson scheme. Direct iteration between heat conduction and stress analysis is employed when stress-dependent interface thermal resistance coefficients are utilized. Verification problems are presented to demonstrate the accuracy of the finite element code.     

Modeling of Structures in the Mechanics of Composites

One possible general statement of a quasi-static problem in the mechanics of composites is considered. It is assumed that a composite is characterized not only by the heterogeneity of a regular structure, but also by the presence of imperfections, impurities, cracks, and the roughness of surfaces, which are partly taken into account by introducing appropriate couple stresses. Two statements, in displacements and in stresses, are considered together with the statement of the same problems in the case where the constitutive relations are linear integral operators. The boundary-value problem remains nonlinear due to the nonlinearity of a scattering function which enters into the heat equation. The theory of effective moduli for a nonpolar medium is discussed in more detail. The equilibrium equations for a homogeneous medium with reduced characteristics and the equation of heat inflow, introduced in nonlinear (in an explicit form) and linear variants, are examined. For a simple laminated composite, all effective mechanical and thermophysical characteristics are found in an explicit form. The effective material functions for a transversely isotropic medium are constructed on the basis of a unique dimensionless relaxation kernel with the use of several Il'ushin kernels. Based on the known solution of the boundary-value problem for the reduced medium, the stress and strain concentration tensors, at any point of a simple laminated composite, are also constructed in an explicit form. In this case, the changes in the structure are taken into account.     

Buckling Instability of a Thick Rectangular Plate of a Composite Material with a Spatially Periodically Curved Structure

Within the framework of the continuum approach for composite materials with a spatially curved structure developed by Akbarov and Guz', with the use of the three-dimensional linearized theory of stability, the buckling instability of a rectangular plate made of a composite material is investigated. Various edge conditions are considered, and, for obtaining a numerical result, a three-dimensional FEM modelling is developed. Uniaxial and biaxial precritical compression of the plate is analyzed. The numerical results presented illustrate the influence of problem parameters on the critical relative shortening of the plate.     

Follower loads for high order finite elements

Finite element modelling of hydrostatic compaction where the applied pressure acts normal to the deformed surface requires a geometric nonlinear formulation and follower load terms [1, 5, 7]. These concepts are applied to high order [6] (p-FEM) elements with hierarchic shape functions. Applying the blending function method allows to precisely describe curved boundaries on coarse meshes. High order elements exhibit good performance even for high aspect ratios and strong distortion and therefore allow an efficient discretization of thin-walled structures. Since high order finite elements are less prone to locking effects a pure displacement-based formulation can be chosen. After introducing the basic concept of the p-version the application of follower loads to geometrically nonlinear high order elements is presented. For the numerical solution the displacement based formulation is linearized yielding the basis for a Newton-Raphson iteration. The accuracy and performance of the high order finite element scheme is demonstrated by a numerical example. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Automatic solution of optimal-control problems. I. simplest problem in the calculus of variations

An automatic method for obtaining the numerical solution for the simplest problem in the calculus of variations is described. The nonlinear two-point boundary-value Euler-Lagrange equation is solved using the Newton-Raphson method for obtaining successive approximations of the solution. The derivatives required for the solution of the problem are computed automatically using the table method. The user of the program need only input the integrand of the objective function in the calculus-of-variations problem and specify the boundary conditions. None of the derivatives usually associated with the Euler-Lagrange equation and the Newton-Raphson method need be calculated by hand. An example is given with numerical results. The automatic solution of the simplest problem in the calculus of variations in this paper is considered to be the first step in the automatic solution of more general optimal-control problems.     

Stress Analysis for a Thick Rectangular Plate of a Composite Material with a Spatially Periodically Curved Structure under Forced Vibration

Natural and forced vibrations of a thick rectangular plate fabricated from a composite material with a spatially periodically curved structure are investigated with the use of exact three-dimensional equations of motion of the theory of elastic anisotropic bodies. The investigations are carried out within the framework of the continuum approach developed by Akbarov and Guz'. It is supposed that the plate is clamped at all its edges and is loaded on the upper face with uniformly distributed normal forces periodically changing with time. The influence of curving parameters on the fundamental frequency of the plate and on the distribution of the normal stress acting in the thickness direction under forced vibration is studied. The corresponding boundary-value problems are solved numerically by employing the three- dimensional FEM modeling.     

Variational statement of deformation problems for a composite latticed plate with various types of lattices

The variational statement of various boundary value problems for tangential displacements and forces in a latticed plate with an arbitrary piecewise smooth contour is investigated. The lattice consists of several families of bars made of a homogeneous composite material with a matrix of relatively low shear stiffness. The energy method reduces the problem to the variational problem of minimizing the energy functional. The conditions on the plate contour are established under which the functional is minimal and positive definite, which ensures that the problem is well posed.     

Successive approximation procedure of optimal tracking control for nonlinear similar composite systems

The optimal tracking control (OTC) problem for a class of affine nonlinear composite systems with similar structure is considered. By using a modeling technique, the nonlinear similar composite system is first transformed into some quasi-decoupled subsystems. Then the high-order, strongly coupled, nonlinear two-point boundary value (TPBV) problem is transformed into a sequence of linear decoupled TPBV problems through a successive approximation procedure. The obtained OTC law consists of an accurate linear term and a nonlinear compensation term which is the limit of the adjoint vector sequence. A suboptimal tracking control law is obtained by truncating a finite iterative result of the adjoint vector sequence as its nonlinear compensation term.     

A local iteration scheme for nonlinear two-dimensional steady-state heat conduction: a BEM approach

An efficient algorithm is proposed to solve the steady-state nonlinear heat conduction equation using the boundary element method (BEM). Nonlinearity of the heat conduction equation arises from nonlinear boundary conditions and temperature dependence of thermal conductivity. Using Kirchhoff's transformation, the case of temperature dependence of thermal conductivity can be transformed to the nonlinear boundary conditions case. Applying the BEM technique, the resulting matrix equation becomes nonlinear. The nonlinearity, however, only involves the boundary nodes that have nonlinearboundary conditions. The proposed local iterative scheme reduces the entire BEM matrix equation to a smaller matrix equation whose rank is the same as the number of boundary nodes with nonlinear boundary conditions. The Newton-Raphson iteration scheme is used to solve the reduced nonlinear matrix equation. The local iterative scheme is first applied to two one-dimensional problems (analytical solutions are possible) with different nonlinear boundary conditions. It is then applied to a two-region problem. Finally, the local iterative scheme is applied to two cavity problems in which radiation plays a role in the heat transfer.     

Galerkin's procedure,quasilinearization, and nonlinear boundary-value problems

Among the popular and successful techniques for solving boundary-value problems for nonlinear, ordinary differential equations (ODE) are quasilinearization and the Galerkin procedure. In this note, it is demonstrated that utilizing the Galerkin criterion followed by the Newton-Raphson scheme results in the same iteration process as that obtained by applying quasilinearization to the nonlinear ODE and then the Galerkin criterion to each linear ODE in the resulting sequence. This equivalence holds for only the Galerkin procedure in the broad class of weighted-residual methods.This work was supported in part by the National Science Foundation, Grant No. GJ-1075.     

Formulation and evaluation of a finite element model for the linear analysis of stiffened composite cylindrical panels

A finite element model for linear static and free vibration analysis of composite cylindrical panels with composite stiffeners is presented. The proposed model is based on a cylindrical shell finite element, which uses a first-roder shear deformation theory. The stiffeners are curved beam elements based on Timoshenko and Saint-Venant assumptions for bending and torsion respectively. The two elements are developed in a cylindrical coordinate system and their stiffness matrices result from a hybrid-mixed formulation where the element assumed stress field is such that exact equilibrium equations are satisfied. The elements are free of membrane and shear locking with correct satisfaction of rigid body motions. Several examples dealing with stiffened isotropic and laminated plates and shells with eccentric as well as concentric stiffeners are analyzed showing the validity of the models.     

Formation of Saddle-Shaped Composite Sheets

Based on the deformation model of an unbalanced multilayer composite, changes in bending curvatures of sheet-type composites with nonsymmetric structure relative to the midplane of the sheet, depending on the moisture of layers, are predicted. The bending curvatures of saddle-shaped sheets of wood-based composites are calculated with regard to the physical and mechanical properties, geometrical dimensions, orientation, and distribution of layers. The analytical results are compared with the bending curvatures found experimentally for a four-layered unbalanced composite made of birch veneer. The applied calculation model enables us to determine the values of bending curvatures of saddle-shaped wood composite sheets, which can be used in elaborating the technological recommendations.     

Dynamic modeling and analysis of rotating FG beams for capturing steady bending deformation

A dynamic analysis of rotating functionally gradient (FG) beams is presented for capturing the steady bending deformation by using a novel floating frame reference (FFR) formulation. Usually, the cross section of bending beams should rotate round the point at the neutral axis while centrifugal inertial forces are supposed to act on centroid axis. Due to material inhomogeneity of FG beams, centroid and neutral axes may be in different positions, which leads to the eccentricity of centrifugal forces. Thus, centrifugal forces can be divided into three componets: transverse component, axial component and force moment acting on the points of the neutral axis, in which transverse component and force moment can make the beam produce the steady bending deformation. However, this speculation has not been presented and discussed in existing literatures. To this end, a novel FFR formulation of rotating FG beams is especially developed considering centroid and neutral axes. The FFR and its nodal coordinates are used to determine the displacement field, in which kinetic and elastic energies can be accurately formulated according to centroid and neutral axes, respectively. By using the Lagrange's equations of the second kind, the nonlinear dynamic equations are derived for transient dynamics problems of rotating FG beams. Simplifying the nonlinear dynamic equations obtains the equilibrium equations about inertial and elastic forces. The equilibrium equations can be solved to capture the steady bending deformation. Based on the steady bending state, the nonlinear dynamic equations are linearized to obtain eigen-frequency equations. Transient responses obtained from the nonlinear dynamic equations and frequencies obtained from the eigen-frequency equations are compared with available results in existing literatures. Finally, effects of material gradient index and angular speed on the steady bending deformation and vibration characteristics are investigated in detail.     

Identification of parameters in nonlinear boundary-value problems

This paper presents a method for identification of parameters in nonlinear boundary-value problems. The successive approximations technique proposed uses the theory of Lagrange multipliers and the Newton-Raphson method. This method does not require storage of functions and is quadratically convergent. Numerical results are presented.This research was sponsored by the National Institutes of Health, Grant No. GM-16197-01. Computing assistance was obtained from the Health Sciences Computing Facility, University of California at Los Angeles, NIH Grant No. FR-3.