考虑周长约束的圆柱薄壳初始几何缺陷随机场建模方法

利用随机场对圆柱薄壳结构的初始几何缺陷进行建模,并据此建立了一种用于含初始几何缺陷轴压圆柱薄壳屈曲分析的随机分析方法。首先,指出已有将圆柱薄壳初始几何缺陷表征为二维高斯随机场的方法会导致与实际不相符的初始几何缺陷,如圆柱周长显著增大或缩小的几何缺陷。其次,提出一种考虑周长不变约束的随机场建模方法,以剔除与实际不相符的随机几何缺陷。最后,基于所建立的初始几何缺陷随机场模型,利用非干涉多项式混沌展开法进行圆柱薄壳的随机屈曲分析,给出临界屈曲载荷的概率分布。数值试验结果表明,基于随机场理论的初始几何缺陷建模方法可有效刻画几何缺陷对结构承载能力的影响,而提出的约束随机场建模方法又能有效减小结果的分散性。     

几何缺陷浅拱的动力稳定性分析

研究了几何缺陷对粘弹性铰支浅拱动力稳定性能的影响。从达朗贝尔原理和欧拉-贝努利假定出发推导了粘弹性铰支浅拱在正弦分布突加荷载作用下的动力学控制方程,并采用Galerkin截断法得到了可用龙格-库塔法求解的无量纲化非线性微分方程组。同时引入能有效追踪结构动力后屈曲路径的广义位移控制法,对含几何缺陷浅拱的响应曲线进行几何、材料双重非线性有限元分析。用这两种方法分析了前三阶谐波缺陷对浅拱动力稳定性能的影响,其中动力临界荷载由B-R准则判定。主要结论有:材料粘弹性使浅拱动力临界荷载增大且结构响应曲线与弹性情况差别很大;二阶谐波缺陷影响显著,它使动力临界荷载明显下降且使得浅拱粘弹性动力临界荷载可能低于弹性动力临界荷载。     

Buckling and post-buckling of long pressurized elastic thin-walled tubes under in-plane bending

The present paper focuses on the structural stability of long uniformly pressurized thin elastic tubular shells subjected to in-plane bending. Using a special-purpose non-linear finite element technique, bifurcation on the pre-buckling ovalization equilibrium path is detected, and the post-buckling path is traced. Furthermore, the influence of pressure (internal and/or external) as well as the effects of radius-to-thickness ratio, initial curvature and initial ovality on the bifurcation moment, curvature and the corresponding wavelength, are examined. The local character of buckling in the circumferential direction is also demonstrated, especially for thin-walled tubes. This observation motivates the development of a simplified analytical formulation for tube bifurcation, which considers the presence of pressure, initial curvature and ovality, and results in closed-form expressions of very good accuracy, for tubes with relatively small initial curvature. Finally, aspects of tube bifurcation are illustrated using a simple mechanical model, which considers the ovalized pre-buckling state and the effects of pressure.     

Vibrations of plates with initial geometrical periodical imperfections interacting with a periodic elastic foundation

Summary  The aim of this contribution is to formulate an engineering theory describing the effect of periodic structure of the elastic foundation on the dynamic behaviour of a plate with geometrical periodical imperfections, Fig. 1. It is shown that the obtained solutions depend on the structure length parameter l of the elastic foundation and the wave length of imperfections. The general results are illustrated by an example and compared with the known theories of orthotropic plates. Received 7 February 1999; accepted for publication 6 October 1999     

薄壁圆筒壳初始几何缺陷不确定性量化的极大熵方法

针对薄壁圆筒壳结构轴压屈曲载荷的缺陷敏感性以及真实几何缺陷的不确定性,提出一种基于实测缺陷数据和极大熵原理的初始缺陷建模与屈曲载荷预测方法。首先,将初始几何缺陷视为二维随机场,并利用实测缺陷数据和Karhunen-Loève展开法将初始缺陷的随机场建模转化为随机向量的建模;其次,利用极大熵方法确定随机向量的概率分布;最后,基于所构建的初始缺陷随机模型,利用MCMC抽样方法和确定性屈曲分析方法,进行随机屈曲分析并给出基于可靠度的屈曲载荷折减因子。数值算例表明,与直接假设随机场相关结构的方法相比,本文方法的结果是对薄壁圆筒壳屈曲载荷的一个更无偏估计。     

Dynamic, multimode buckling of thin-walled columns subjected to in-plane pulse loading

The present paper deals with dynamic, coupled buckling of long, prismatic columns simply supported at the ends. This investigation concerns thin-walled structures of a square cross-section with or without intermediate stiffeners under in-plane pulse loading. The dynamic load of a rectangular shape has been assumed in the analysis. The structures are composed of rectangular plates interconnected along longitudinal edges. A plate model is adopted in the analysis. The material of the structure is isotropic. The problem has been investigated on the basis of the disturbance theory. The dynamic critical load factor DLF has been determined using the Budiansky and Hutchinson criterion. The results obtained with the analytical-numerical method (ANM), which employs the asymptotic perturbation theory, have been compared with the finite element method (FEM).     

Buckling analysis of cylindrical shells with random geometric imperfections

In this paper the effect of random geometric imperfections on the limit loads of isotropic, thin-walled, cylindrical shells under deterministic axial compression is presented. Therefore, a concept for the numerical prediction of the large scatter in the limit load observed in experiments using direct Monte Carlo simulation technique in context with the Finite Element method is introduced. Geometric imperfections are modeled as a two dimensional, Gaussian stochastic process with prescribed second moment characteristics based on a data bank of measured imperfections. (The initial imperfection data bank at the Delft University of Technology, Part 1. Technical Report LR-290, Department of Aerospace Engineering, Delft University of Technology). In order to generate realizations of geometric imperfections, the estimated covariance kernel is decomposed into an orthogonal series in terms of eigenfunctions with corresponding uncorrelated Gaussian random variables, known as the Karhunen-Loéve expansion. For the determination of the limit load a geometrically non-linear static analysis is carried out using the general purpose code STAGS (STructural Analysis of General Shells, user manual, LMSC P032594, version 3.0, Lockheed Martin Missiles and Space Co., Inc., Palo Alto, CA, USA). As a result of the direct Monte Carlo simulation, second moment characteristics of the limit load are presented. The numerically predicted statistics of the limit load coincide reasonably well with the actual observations, particularly in view of the limited data available, which is reflected in the statistical estimators.     

Closed form solutions for free vibrations of rectangular Mindlin plates

A new two-eigenfunctions theory, using the amplitude deflection and the generalized curvature as two fundamental eigenfunctions, is proposed for the free vibration solutions of a rectangular Mindlin plate. The three classical eigenvalue differential equations of a Mindlin plate are reformulated to arrive at two new eigenvalue differential equations for the proposed theory. The closed form eigensolutions, which are solved from the two differential equations by means of the method of separation of variables are identical with those via Kirchhoff plate theory for thin plate, and can be employed to predict frequencies for any combinations of simply supported and clamped edge conditions. The free edges can also be dealt with if the other pair of opposite edges are simply supported. Some of the solutions were not available before. The frequency parameters agree closely with the available ones through pb-2 Rayleigh-Ritz method for different aspect ratios and relative thickness of plate.     

Perturbation initial parameter method for solving the geometrical nonlinear problem of axisymmetrical shells

In the previous paper[7], the author presented a System of First-Order Differential Equations for the problem of axisynrm’trically loaded shells of revolution with small elastic. strains and arbitrarily large axial deflections, and a Method of Variable-Characteristic Nondimensionization with a Scale of Load Parameter. On this basis, by taking the weighted root-mean-square deviation of angular deflection from linearity as perturbation parameter, this paper pressents a perturbation system of nondimensional differential equations for the problem, thus transforms the geometrical nonlinear problem into several linear problems. This paper calculates these linear problems by means of the initial parameter method of numerical integration. The numerical results agree quite well with the experiments[4].     

Localized analysis of thin-walled structure's buckling/initial post-buckling and its accuracy

A method of localization is proposed to lower the high order of equations in FEM calculation for the stability of a complex thin-walled structure. The localized analysis enables us to obtain both the upper and lower limits for the bifurcating point in a whole linear-elastic structural system, as well as an approximate solution to asymptotic post-buckling problem. Some numerical examples are included. Project supported by National Natural Science Foundation of China     

PERTURBATION ANALYSIS FOR MAGNETO-PLASTIC INSTABILITY OF FERROMAGNETIC BEAM-PLATES WITH GEOMETRIC IMPERFECTION

The magneto-plastic instability of a ferromagnetic beam-type plate with simple supports and small initial imperfection is analytically investigated in this paper for that the plastic deformation of the plate with a linear-strain hardening relation is co…     

A note on the determination of critical loads for flexible noncircular cylindrical shells with clamped edges

An approach is proposed to determine the upper and lower bounds of the critical load for flexible noncircular long cylindrical shells with clamped longitudinal edges under nonuniform loading. An exact analytical solution is designed. Its critical points are analyzed __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 11, pp. 78–87, November 2005.     

非均匀介质中流动和传热问题的有限分析法

数值求解非均匀介质中的输运问题广泛应用于科学计算和工程领域．介质的强非均匀性给相关问题的准确求解带来极大的困难．近年来,本课题组将有限分析法拓展到该领域,建立了非均匀介质中输运问题的有限分析法．该算法基于网格奇点邻域内类拉普拉斯方程局部解析解构建,算法具有很高的精度,且不依赖于介质的非均匀性强度．不管相邻网格传导率差异如何,仅需对原始网格进行很少地细分就可以获得非常准确的计算结果,因此与其他传统数值算法相比,可以大幅提高计算精度和效率．该算法可广泛应用于求解非均匀多孔介质中的渗流、复合材料中的热传导及电场分布等问题．     

On computations of stress intensity factors for stiffened half planes with imperfections and cracks

In this paper we consider uniform extension problems for joined two half-planes with different thickness and material behavior and one of which contains an elliptical hole, the other contains a crack. Along the boundary of these half-planes there is a stiffening stringer. Computational formulas are given in power series form by complex variable—pertubation method. Results obtained here give extension to those of Handbook of stress intensity factors. Numerical results of special cases in this paper coincide with those of refs. [1], [3].     

Exact solution for in-plane static problems of circular beams made of functionally graded materials

Exact analytical solutions of in-plane static problems of circular beams with uniform cross-section made of functionally graded material (FGM) are obtained. Material properties are assumed to be varying arbitrarily through the thickness. The effects of axial extension and shear deformations are considered. The differential equation system is solved exactly using the initial values method. The circumferential stress distribution on the cross-section is also obtained. The results are compared with those of rather complex approaches in the literature, such as elasticity approach, and the comparison shows an excellent agreement. Effects of power law exponent and radius-to-height ratio of the beam on circumferential stress distribution and displacements are investigated.     

Near crack line elastic–plastic analysis for a infinite plate loaded by two pairs of point tensile forces

The near crack line analysis method is used to investigate a center crack loaded by two pairs of point tensile forces in an infinite plate in an elastic–perfectly plastic solid, and the analytical solutions are obtained in this paper. These solutions include: the elastic–plastic stress field near the crack line, the law that the length of the plastic zone along the crack line is varied with an external loads and the bearing capacity of an infinite plate with a center crack. The results of this paper are sufficiently precise near the crack line because the assumptions of the small scale yielding theory have not been used and no other assumptions have been taken.     

The application of graph theory to the limit analysis for torsion of thin-walled bars with complex cross-sections

Graph theory is employed in this paper as a means to establish the topological model of complex thin-walled cross-sections. On this basis, the upper and lower bound theorems of the plastic limit analysis are applied to the analysis of the plastic limit shear flows on the cross-section of thin-walled bars under St. Venant torsion. Corresponding mathematical programming problems are formulated and their duality is shown. After solving the linear programming problem corresponding to the lower bound theorem, the limit torsional moment of a thin-walled cross-section can be calculated according to the shear stress distribution in the limit state. The formula for calculating the limit torsional moment is given. Furthermore, the limit state of thin-walled cross-sections under St. Venant torsion is also discussed and the concept of the limit tree is introduced. A computer program has been developed by the author. Results calculated by the program for typical complex cross-sections are given.     

A comparative study on algorithms of boundary element solution for free torsional vibration of bodies of revolution

In this paper the functionu=rsin in cylindrical coordinates (r,,z) is introduced into the equation for free torsional vibration of bodies of revolution (where=v / r represents the angle of twist). With the static fundamental solution (–1 /R) a mixed BEM / FEM equation is derived. The domain integral term in the equation is discretized by Serendipity elements instead of commonly used constant value finite elements in the literature. The equation is an algebraic eigenvalue one. The dynamic fundamental solution (e ^1R /R) is also used for deriving the other mixed BEM / FEM equation. An appropriate iterative solution procedure is described. An algebraic eigenvalue equation can be obtained and its solution accuracy is almost interior meshing independent. A number of examples are studied. The results show the good economy and high accuracy of the algorithms proposed.The Project is Supported by National Natural Science Foundation of China.     

Two-fluid model for transient analysis of slug flow in oil wells

In this work it is presented a transient, one-dimensional, adiabatic model for slug flow simulation, which appears when liquid (mixture of oil and water) and gas flow simultaneously through pipes. The model is formed by space and time averaged conservation equations for mass, momentum and energy for each phase, the numerical solution is based on the finite difference technique in the implicit scheme. Velocity, pressure, volumetric fraction and temperature profiles for both phases were predicted for inclination angles from the horizontal to the vertical position (unified model) and ascendant flow. Predictions from the model were validated using field data and ten correlations commonly used in the oil industry. The effects of gas heating or cooling, due to compression and expansion processes, on the predictions and numerical stability, were studied. It was found that when these effects are taken into account, a good behavior of temperature predictions and numerical stability are obtained. The model presents deviations lower than 14% regarding field data and it presents better predictions than most of the correlations.     

Elastoplastic solution for an eccentric crack loaded by two pairs of point tensile forces

The strain energy density theory and the near crack line analysis method are applied to investigate an eccentric crack loaded by two pairs of tensile point forces in a finite plate. The minimum values of SED in the vicinity of the crack tip are determined, the initial growth orientation of crack are determined. Obtained is the elastic-plastic solution near the crack line of an eccentric crack loaded by two pairs of point tensile forces under large scale yielding condition. More specifically, the near field solution contains the unit normal vector of the elastic-plastic boundary and the elastic-plastic stress field. The length of the plastic zone along the crack line is found to vary with the external load and the bearing capacity of a finite plate with an eccentric crack loaded by two pairs of tensile point forces. Compared with small scale yielding condition, the normalized load obtained is higher than those under small scale yielding condition when the length of the plastic zone is the same.