Application of general variational methods with discontinuous fields to bending,buckling, and vibration of beams

The effectiveness of consistent variational statements for discontinuous fields is illustrated for bending, buckling, and vibration of elastic beams. General variational statements are developed, from which finite-element formulations with piece-wise constant and piece-wise linear trial functions, are obtained. Results are illustrated by means of numerical examples.     

Application of a refined plate bending element to buckling problems

The buckling stiffness matrix of a refined plate bending element is derived for various continuous and non-uniform distributions of in-plane forces. Elements of this matrix are expressed in an explicit form and can, therefore, be readily used in a finite element computer program for solving plate buckling problems with various loading and edge conditions. A number of problems are solved and the results obtained are compared with analytical and/or numerical results obtained by using other procedures. It is shown that the refined plate bending element is superior to simpler elements when used for solving plate buckling problems.     

Application of a trigonometric finite difference procedure to numerical analysis of compressive and shear buckling of orthotropic panels

A numerical analysis developed for the buckling of rectangular orthotropic layered panels under combined shear and compression is described. This analysis uses a central finite difference procedure based on trigonometric functions instead of using the conventional finite differences which are based on polynomial functions. Inasmuch as the buckle mode shape is usually trigonometric in nature, the analysis using trigonometric finite differences can be made to exhibit a much faster convergence rate than that using conventional differences. Also, the trigonometric finite difference procedure leads to difference equations having the same form as conventional finite differences; thereby allowing available conventional finite difference formulations to be converted readily to trigonometric form. For two-dimensional problems, the procedure introduces two numerical parameters into the analysis. Engineering approaches for the selection of these parameters are presented and the analysis procedure is demonstrated by application to several isotropic and orthotropic panel buckling problems. Among these problems is the shear buckling of stiffened isotropic and filamentary composite panels in which the stiffener is broken. Results indicate that a break may degrade the effect of the stiffener to the extent that the panel will not carry much more load than if the stiffener were absent.     

基于PCL的折叠式夹层板结构参数化建模

折叠式夹层板具有优越的力学性能,在有限元分析中,由于夹层板结构复杂,建模工作量大.针对这些问题,基于PCL对Patran进行二次开发,通过输入夹层板尺寸和网格数量等参数实现V形折叠式夹层板结构的参数化建模.编制折叠式夹层板建模界面,实现夹层板结构有限元分析参数化,提高建模效率,为夹层板结构设计提供支撑.     

Natural frequencies of vibration of cantilever sandwich beams

Theoretical and experimental studies were made in obtaining the natural frequencies of cantilever sandwich beams subjected to only gravity forces. The method of minimizing the total energy of the system was used for determining the frequencies. A vibration system made by Unholtz-Dickie was utilized to set the beam in vibration. Resonance occurred when the frequency of the shaker coincided with the natural frequency of the beam. The resonance frequencies were measured by transducers mounted at various locations on the beam. A total of sixteen beams of various lengths, thickness and core density were tested.

It was found that the natural frequency of a cantilever sandwich beam depends largely upon the thickness, length, core density and stiffness of the beam. In addition, the natural frequency has a nonlinear variation with the mode and for any particular mode, the value of the frequency increases as the length of the beam decreases.

Design factors were developed based upon the ratios of the theoretical frequencies of homogeneous beams having the same thicknesses and stiffnesses of that of sandwich beams and of the frequencies experimentally determined for similar sandwich beams.     

Numerical analysis of sandwich beams

This paper describes the development of a numerical model for the physical nonlinear analysis of simply supported sandwich beams, specifically with foamed-concrete cores and concrete faces. The long-term behaviour is included in view of creep and shrinkage of both faces and core. The structural behaviour of sandwich beams is described by a fourth-order differential equation in the deformation w and a second-order differential equation in the shear deformation of the core γ_k. The flexural stiffness of the core is taken into account. The general-solution procedure is based on the finite-difference method, together with a successive-substitution algorithm using the secant flexural moduli of the core and faces and the secant shear modulus of the core. The option of tension stiffening is incorporated to represent the nonlinear behaviour of reinforced concrete in tension. The tension stiffening is numerically calculated from a distributed tensile load instead of a load acting on both ends of the reinforced bar. Creep and shrinkage are calculated separately from the differential equations with an algorithm based on increments of time. With the presented model, the time-dependent deflection along the axis of the beam and the state of stress in every fibre can be calculated.     

Nonlinear vibration and stability of a shallow unsymmetrical orthotropic sandwich shell of double curvature with orthotropic core

In this paper, nonlinear behaviours for a shallow unsymmetrical, orthotropic sandwich shell of double curvature with orthotropic core having different elastic characteristics have been studied by a new set of uncoupled differential equations. The face sheet may be of unequal thickness of different materials. However, a restriction that the elements radii of curvature be large compared to the overall thickness of the sandwich has been imposed.

A simple approach used in the present analysis can be applied for stability as well as vibration. For the symmetrical case, where the face sheets are of equal thickness and of same materials, these equations can be shown to reduce to those given by Grigolyuk in 1957. Numerical results of a square rectangular simply supported curved plate, and of a rectangular sandwich cylindrical shell under mechanical and dynamic loading, have been computed and compared with other known results.     

功能梯度材料椭圆板的非线性热振动及屈曲

采用弹性理论建立了功能梯度材料板的静力平衡方程,利用静力平衡方程确定了功能梯度材料板的中性面位置,在此基础上推导出了功能梯度材料板在均匀温度场中的非线性振动及屈曲微分方程组,求得了功能梯度材料椭圆板的非线性振动及屈曲的近似解,讨论分析了中性面位置、梯度指数、温度等因素对功能梯度材料椭圆板非线性振动及屈曲的影响.把该方法计算结果与有限元计算结果进行了比较,验证了该方法的计算结果是可靠的.算例分析表明,中性面位置对均匀温度场中功能梯度材料椭圆板的非线性振动及屈曲有一定影响.     

功能梯度材料圆板的非线性热振动及屈曲

采用弹性理论建立了功能梯度材料板的静力平衡方程,利用静力平衡方程确定了功能梯度材料板的中性面位置,在此基础上推导出了功能梯度材料板在均匀温度场中的非线性振动及屈曲微分方程组,求得了功能梯度材料圆板的非线性振动及屈曲的近似解,讨论分析了中性面位置、梯度指数、温度等因素对功能梯度材料圆板非线性振动及屈曲的影响.把该方法计算结果与有限元计算结果进行了比较,验证了该方法的计算结果是可靠的.算例分析表明,中性面位置对均匀温度场中功能梯度材料圆板的非线性振动及屈曲有一定影响.     

非线性地基上正交异性矩形板的非线性热振动

研究了非线性地基上正交异性矩形板的非线性固有热振动.采用常规的L-P法分析非线性地基上正交异性矩形板的非线性热振动难以得到高精度的近似解,为此,先对该强非线性振动系统进行参数变换,将该强非线性振动系统转化为弱非线性振动系统.然后采用改进的L-P法进行求解,得到了强非线性振动系统的高精度近似解.此外,讨论了温度、地基特征参数、长宽比等因素对非线性地基上正交异性矩形板非线性热振动固有频率的影响,得到了非线性地基上正交异性矩形板热振动频率随温度下降、地基特征参数变大、长宽比变大而增大的结论.     

Vibration and buckling analysis of axisymmetric polar orthotropic circular plates

The vibration and stability analysis of polar orthotropic circular plates using the finite element method is discussed. In order to formulate the eigenvalue problems associated with the vibration and stability analyses, the clement stiffness, mass, and stability coefficient matrices are presented. By assuming the static displacement function, which is an exact solution of the polar orthotropic circular plate equation, approximates the vibration and buckling modes, the mass and stability coefficient matrices are readily derived from the given displacement function. Results showing the effects of orthotropy on natural frequencies and buckling loads are compared with their isotropic counterpart.     

Trefftz-unsymmetric finite element for bending analysis of orthotropic plates

This work develops a new four-node quadrilateral displacement-based Trefftz-type plate element for bending analysis of orthotropic plates within the framework of the unsymmetric finite element method (FEM). In the present formulation, the modified isoparametric interpolations are employed to formulate the element’s test functions in which the deflection is effectively enriched by the nodal rotation degrees of freedom (DOFs). Meanwhile, the element’s trial functions are determined based on the Trefftz functions that can a prior satisfy the governing equations of orthotropic Mindlin–Reissner plates. Numerical benchmark tests reveal that the new unsymmetric plate element is free of shear locking problem and can produce satisfactory results for both the displacement and stress resultant. In particular, it exhibits quite good tolerances to the gross mesh distortion.

     

Evolutionary natural frequency optimization of thin plate bending vibration problems

This paper extends the evolutionary structural optimization method to the solution for maximizing the natural frequencies of bending vibration thin plates. Two kinds of constraint conditions are considered in the evolutionary structural optimization method. If the weight of a target structure is set as a constraint condition during the natural frequency optimization, the optimal structural topology can be found by removing the most ineffectively used material gradually from the initial design domain of a structure until the weight requirement is met for the target structure. However, if the specific value of a particular natural frequency is set as a constraint condition for a target structure, the optimal structural topology can be found by using a design chart. This design chart describes the evolutionary process of the structure and can be generated by the information associated with removing the most inefficiently used material gradually from the initial design domain of a structure until the minimum weight is met for maintaining the integrity of a structure. The main advantage in using the evolutionary structural optimization method lies in the fact that it is simple in concept and easy to be included into existing finite element codes. Through applying the extended evolutionary structural optimization method to the solution for the natural frequency optimization of a thin plate bending vibration problem, it has been demonstrated that the extended evolutionary structural optimization method is very useful in dealing with structural topology optimization problems.     

粘弹性夹层圆板自由振动的理论解

研究了粘弹性夹层圆板的自由振动特性.基于经典弹性薄板理论和Kelvin-Voigt粘弹性本构方程,建立了粘弹性夹层圆板振动控制方程.采用分离变量法导出了粘弹性夹层圆板的自然频率及振型解析表达式,计算了固支和简支粘弹性夹层圆板的自然频率,并与有限元计算结果进行比较;讨论了粘弹性夹层圆板的夹心层比率对自然频率及衰减系数的影响.研究表明：（1）随着夹心层厚度的增大,系统频率先增大后减小,高阶时该趋势表现更为明显;（2）随着夹心层厚度的增大,衰减系数一直增大,高阶时该趋势表现更为明显.     

Lateral buckling analysis of beams by the fem

Two new finite element formulations for the calculation of the lateral buckling load for elastic straight prismatic thin-walled open beams under conservative static loads, are presented. The stability criterion used is based on the positive definiteness of the second variation of the total potential energy. One formulation is suitable for sections where the initial bending is about a dominant major axis. The other finite element formulation takes account of initial bending curvature and essentially takes the form of a quadratic eigenvalue problem. Both formulations are tested with problems that have classical solutions or experimentally determined results and are shown to be accurate.     

Probabilistic free vibration analysis of beams subjected to axial loads

A stochastic finite-element-based algorithm for the probabilistic free vibration analysis of beams subjected to axial forces is proposed in this paper through combination of the advantages of the response surface method, finite element method and Monte Carlo simulation. Uncertainties in the structural parameters can be taken into account in this algorithm. Three response surface models are proposed. Model I: star experiment design using a quadratic polynomial without cross-terms; Model II: minimum experiment design using a quadratic polynomial with cross-terms; Model III: composite experiment design using a quadratic polynomial with cross-terms.A separate set of finite element data is generated to verify the models. The results show that the Model II is the most promising one in view of its accuracy and efficiency. Probabilistic free vibration analysis of a simply supported beam is performed to investigate the effects of various parameters on the statistical moments of the frequency response of beams. It is found that the geometric properties of beams have significant effects on the variation of frequency response.     

Static,vibration and buckling analysis of axisymmetric circular plates using finite elements

The static, vibration, and buckling analysis of axisymmetric circular plates using the finite element method is discussed. For the static analysis, the stiffness matrix of a typical annular plate element is derived from the given displacement function and the appropriate constitutive relations. By assuming that the static displacement function, which is an exact solution of the circular plate equation ?²?²W = 0, closely represents the vibration and buckling modes, the mass and stability coefficient matrices for an annular element are also constructed. In addition to the annular element, the stiffness, mass, and stability coefficient matrices for a closure element are also included for the analysis of complete circular plates (no center hole). As an extension of the analysis, the exact displacement function for the symmetrical bending of circular plates having polar orthotropy is also given.     

An element for static, vibration and buckling analysis of thick laminated plates

A conforming finite element formulation of the equations governing composite multilayered plates using Reddy's higher-order theory is presented. The element has eight degrees of freedom, u₀, v₀, w, ∂w/∂x, ∂w/∂y, ∂²w/∂x∂y, γ_x, γ_y, per node. The transverse displacement of the present element is described by a modified bicubic displacement function while the in-plane displacements and shear-rotations are interpolated quadraticly. The element is evaluated for its accuracy in the analysis of static, vibration, and buckling of anisotropic rectangular plates with different lamination schemes and boundary conditions. The conforming finite element described here for the higher-order theory gives fairly accurate results for displacements, stresses, buckling loads, and natural frequencies.     

变厚度智能硬夹心板振动分析

研究由形状记忆合金与普通钢材制成的硬夹心板的振动控制方法,求出硬夹心夹层板的平衡方程,并分析了变厚度智能夹心板的振动问题.算例表明该方法对于夹心板的振动能够有效地控制.     

Refined sandwich model for the vibration of beams with embedded shear piezoelectric actuators and sensors

This work extends a previously presented refined sandwich beam finite element (FE) model to vibration analysis, including dynamic piezoelectric actuation and sensing. The mechanical model is a refinement of the classical sandwich theory (CST), for which the core is modelled with a third-order shear deformation theory (TSDT). The FE model is developed considering, through the beam length, electrically: constant voltage for piezoelectric layers and quadratic third-order variable of the electric potential in the core, while mechanically: linear axial displacement, quadratic bending rotation of the core and cubic transverse displacement of the sandwich beam. Despite the refinement of mechanical and electric behaviours of the piezoelectric core, the model leads to the same number of degrees of freedom as the previous CST one due to a two-step static condensation of the internal dof (bending rotation and core electric potential third-order variable). The results obtained with the proposed FE model are compared to available numerical, analytical and experimental ones. Results confirm that the TSDT and the induced cubic electric potential yield an extra stiffness to the sandwich beam.