中厚板统计分析的随机边界元法

本文建立了分析含随机材料参数并具厚度不均匀性的中厚板问题的随机边界元法,基于Taylor级数展开技术,分析和到广义位移的均值和一阶偏差的积分方程,其中将材料参数的随机性和厚度的不均匀性作为等效荷载处理,从而得到广义边界位移或面力的均值和协方差,并进一步求出部点广义位移和内力的均值和协方差,最后用本文方法计算了两个数例,并对所得结果进行了分析,探讨。     

The method of transformation of the boundaries for structure the admissible displacements

In this paper the method of transformation of the boundaries for structure the admissible displacements with various boundary conditions is presented. What is called the method of transformation of the boundaries is that, first we transform the actual system into the basic system and additional boundary forces and displacements on the basis of the superposition principle, then apply variational principles to the basic system, finally find the permissible displacement of the actual system by means of the method of transformation of the series.In this paper, we also give mixed energy principles under variation of boundary conditions. The mixed energy principles as the potential and complementary energy principles under variation of boundary conditions are all the chief theoretical fundamental of the method of transformation of the boundaries.Applying the method of transformation of the boundaries, we form the permissible displacements of rectangular plates of plane stress and bending problems with various edge conditions.Because the method of transformation of the boundaries is in progress to follow the variational principles and definite program to form permissible displacements, the difficulty in supposing and piecing together permissible displacements in the Rayleigh-Ritz method will be overcome.     

Load and Boundary Condition Calibration Using Full-field Strain Measurement

For critical load bearing structures, it is often necessary to experimentally determine the load distribution on the structure so that accurate finite element models can be developed for stress and fatigue life predictions. An inverse problem approach is presented here for computing or calibrating the loads and boundary conditions acting on a structure. This enables the creation of more accurate finite element models, especially for structures that have complicated load distribution and compliant boundary conditions. The method presented here involves minimizing the least square error between the strains computed using the finite element model and the strains and displacements obtained experimentally. The nodal loads and the compliance at fixed boundaries are treated as the variables in the optimization problem. The compliance is modeled as springs attached at the nodes that are on the boundary where the structure is restrained. The method is verified by computing the loads and boundary conditions when displacements, maximum shear strain or both are available at large number of points on the surface of the structure. The experimental data set was generated using the luminescent photoelastic coating (LPC) technique.     

Integral form of the skin friction coefficient suitable for experimental data

An integral method to evaluate skin friction coefficient for turbulent boundary layer flow is presented. The method replaces streamwise gradients with total stress gradients in the wall-normal direction and is therefore useful in cases when measurements at multiple streamwise locations are not available or feasible. It is also shown to be especially useful for experimental data with typical noisy shear stress profiles such as rough-wall boundary layer flows for which there are limited ways by which skin friction can be determined.     

Opening-Function Simulation of the Three-Dimensional Nonstationary Interaction of Cracks in an Elastic Body

Integral relations between three-dimensional dynamic displacements (stresses) in an infinite elastic body with arbitrarily located plane cracks and discontinuities in the displacements of the opposite crack faces are presented. The influence of opening cracks on each other is considered in the problem on crack faces loaded by pulse forces. This problem is reduced to a system of boundary integral equations of the wave-potential type in a time domain. The dynamic mode I stress intensity factors are determined for two coplanar elliptic cracks under forces in the form of the Heaviside function     

基于余能原理的有限变形问题有限元列式

利用基面力概念,推导了一种基于余能原理的有限变形问题显式有限元列式,可应用于结构的大位移、大转动问题。以基面力为状态变量来表达单元的余能,将有限变形情况下的单元余能分解为变形余能部分和转动余能部分,利用Lagrange乘子法推导出余能原理有限元的控制方程,编制了相应的非线性有限元程序。通过算例分析,说明该列式和程序的可靠性和精确性。     

任意外型薄板非线性后屈曲问题的边界元解

本文就薄板后屈曲问题建立一组新型的边界元计算公式，用这组公式求解能方便处理各种边界问题，另外文中将面内应力分解成基本部份和附加部份，并利用微分算子分解理论导得了挠度的一个不同形式的基本解，由于计算公式中，实现了面内位移和挠度的解耦，从而使迭代过程得到简化，文末还对圆板后屈曲路径进行了计算，得到了满意的结果。     

Fracture mechanics analysis of curved stiffened panels using BEM

A dual boundary element method is developed for a analysis of reinforced cracked shallow shells. Boundary integral equations are derived from the Betti’s reciprocal theorem for a cracked shallow shell with transverse frames and longitudinal stiffeners. The effect of frames and stiffeners are treated as a distribution of line body forces. The radial basis function is used to transform domain integrals to boundary integrals. Stress intensity factors are evaluated from crack opening displacements. The effect of curvature on the stress intensity factors is illustrated by numerical examples. Three examples are presented to demonstrate the accuracy of this method compared with solutions obtained using the finite element method.     

Buckling analysis of a laminated cylindrical shellunder torsion subjected to mixed boundary conditions

A new efficient method is developed in this paper for buckling analysis of a cross-plylaminated cylindrical shell under torsion subjected to mixed boundary conditions. The transverseshear is taken into account by a first-order theory with a shear correction factor of 5/6. The mixedboundary conditions include conditions in forces as well as conditions in displacements, and theseforces and displacements are selected as basic unknowns. The other displacements and forces areexpressed in terms of the basic unknowns by taking inverse of a matrix composed of operators.The equations of buckled equilibrium in terms of the basic unknowns are solved with doubletrigonometric series which satisfy the mixed boundary conditions. Comparison of the obtainednumerical results with those given in the literature based on completely clamped boundaryconditions checks with the fact that the mixed boundary conditions yield appreciably lowerbuckling load and less circumferential wave number than the completely clamped boundaryconditions. The curves in the figures show how the difference in buckling loads between the twokinds of boundary conditions varies when the length and thickness of the shell vary.     

Bending tests to estimate the axial force in slender beams with unknown boundary conditions

This paper presents a static method for the axial load identification of slender prismatic beams with uncertain length and unknown boundary conditions as is typical of struts and ties of truss structures or tie-rods of arches and vaults. The proposed method requires the knowledge of the beam flexural rigidity only. Flexural displacements or curvatures are measured at five cross sections of the beam subjected to an additional concentred lateral load. Unlike analogous dynamic methods, any set of experimental data may be used in the identification algorithm. The proposed algorithm is verified by means of many numerical and experimental tests on beams having different boundary conditions. Excellent estimates of the axial forces are obtained when the greatest possible distance between sensors is adopted, even if in the presence of high values of the axial forces very accurate displacement measurements are required.     

Non-linear viscoelastodynamic equations of three-dimensional rotating structures in finite displacement and finite element discretization

This paper deals with the non-linear viscoelastodynamics of three-dimensional rotating structure undergoing finite displacement. In addition, the non-linear dynamics is studied with respect to geometrical and mechanical perturbations. On part of the boundary of the structure, a rigid body displacement field is applied which moves the structure in a rotation motion. A time-dependent Dirichlet condition is applied to another part of the boundary. For instance, this corresponds to the cycle step of a helicopter rotor blade. A surface force field is applied to the third part of the boundary and depends on the time history of the structural displacement field. For example, this might corresponds to general unsteady aerodynamics forces applied to the structure. The objective of this paper is to model the non-linear dynamic behavior of such a rotating viscoelastic structure undergoing finite displacements, and to allow small geometrical and mechanical (mass, constitutive equations) perturbations analysis to be performed. The model is constructed by the introduction of a reference configuration which is deduced from the non-linear steady boundary value problem. A constitutive equation deduced from the Coleman and Noll theory concerning the viscoelasticity in finite displacement is used. Thereafter, the weak formulation of the boundary value problem is constructed and discretized using the finite element method. In order to simplify the mathematical study of the equations, multilinear forms are introduced in the algebraic calculation and their mathematical properties are presented.     

Generalized expressions for boundary conditions of shells of revolution

For the boundary conditions of shells of revolution, traditionally, four out of the eight quantities which are the four displacements on the middle surface u, v, w and if together with the four corresponding forces, are given. when the generalized displacements on the nodal circles are used as basic unknowns, the number of unknowns on a nodal circle is more than four^[1][2][3][4]. In this case, how to deal with the boundary conditions is still a problem that has not been solved satisfactorily yet. In this paper,the relations between the generalized and nongeneralized quantities of a shell’s edge are derived according to the principle of virtual work. Seven types of common edges are studied and their expressions of boundary conditions in the form of generalized displacements or forces are qiven. The number of expressions for each type of edge may correspond with the number of unknowns used on a nodal circle. Kith these expressions, boundary conditions can be put directly into equations of motion of generalized displacement method so as to solve the generalized displacements. By so doing, the process of transformation and inverse transformation about unknowns in [2] is avoided. Not only is the argument simple and clear, but the calculation work is reduced.Having the set of generalized expressions of boundary conditions, the generalized displacement method of the shell of revolution may be more perfect in theory.     

THE STOCHASTIC BOUNDARY ELEMENT METHOD IN STATISTICAL ANALYSIS OF MODERATELY THICK PLATES

In this paper a stochastic boundary element method(SBEM)isdeveloped to analyze moderately thick plates with random material parameters andrandom thickness.Based on the Taylor series expansion,the boundary integrationequations concerning the mean and deviation of the generalized displacements arederived,respectively.It is found that the randomness of material parameters isequivalent to a random load,so the mean and covariance matrices of unknowngeneralized boundary displacements and tractions can be obtained.Furthermore,themean and covariance of generalized displacements and forces at internal points can alsobe obtained.A numerical example has been worked out with the method proposed andnecessary analysis is made for the results.     

局部有限元模型边界条件的直接传递方法

提出一种将整体分析得到的节点力或节点位移直接传递到精细化局部有限元模型的方法,即部分混合单元法。沿精细化局部有限元模型周边建立一组过渡单元,该组过渡单元采用与整体模型一致的单元类型和模拟方式,其外侧边界上的节点与整体模型节点的相对坐标对应,内侧边界与精细化局部有限元模型采用基于面约束的方式连接。在外侧边界上根据节点坐标对应施加整体分析获得的节点力或节点位移,过渡单元就可直接将边界条件传递到精细化局部有限元模型。通过贵州红水河特大桥钢-混结合段的精细化有限元分析,验证了本文方法的实用性和有效性。     

A continuum model for size-dependent deformation of elastic films of nano-scale thickness

Ultra-thin elastic films of nano-scale thickness with an arbitrary geometry and edge boundary conditions are analyzed. An analytical model is proposed to study the size-dependent mechanical response of the film based on continuum surface elasticity. By using the transfer-matrix method along with an asymptotic expansion technique of small parameter, closed-form solutions for the mechanical field in the film is presented in terms of the displacements on the mid-plane. The asymptotic expansion terminates after a few terms and exact solutions are obtained. The mid-plane displacements are governed by three two-dimensional equations, and the associated edge boundary conditions can be prescribed on average. Solving the two-dimensional boundary value problem yields the three-dimensional response of the film. The solution is exact throughout the interior of the film with the exception of a thin boundary layer having an order of thickness as the film in accordance with the Saint-Venant’s principle.     

Elastostatics with non absolutely continuous data

The main problem of classical elastostatics is considered when given body forces and surface forces act on the body. The problem is studied in the case when the given forces are represented by general vector measures, i.e. when concentrated loads may be given on subsets either of the body or of its boundary. An existence and uniqueness (up to rigid displacements) theorem is proved. Conditions under which the energy integral is finite are given.     

Factorization method in the geometric inverse problem of static elasticity

The factorization method, which has previously been used to solve inverse scattering problems, is generalized to geometric inverse problems of static elasticity. We prove that finitely many defects (cavities, cracks, and inclusions) in an isotropic linearly elastic body can be determined uniquely if the operator that takes the forces applied to the body outer boundary to the outer boundary displacements due to these forces is known.     

Correct Solutions of Plane Elastic Problems for a Half-Plane

Continuous and integrable solutions and one-to-one relationships between boundary forces and displacements are found through the direct integration of the differential equations of the plane elastic problem for a half-plane with boundary conditions for either forces or displacements or with mixed boundary conditions. The necessary equilibrium conditions for forces and the compatibility conditions for displacements that ensure the correctness of the solutions are formulated     

On quadruple integral equations related to a certain crack problem

An exact solution of a four part mixed boundary value problem representing a three colinear crack system connected with specified crack opening displacements between the cracks is obtained. The three cracks thus become one with pressure and/or opening displacement prescribed on the crack face. From considerations of dual symmetry and a formulation based on Papkovich-Neuber harmonic functions, the boundary value problem is reduced to solving a quadruple set of integral equations. An exact solution of these equations is derived using a modified finite Hilbert transform technique. The closed form results for the stress distributions and the crack-tip stress intensity factors are presented. Limiting cases of the solution yield results which agree with well known solutions.     

Elasticity solutions for functionally graded annular plates subject to biharmonic loads

Based on England’s expansion formula for displacements, the elastic field in a transversely isotropic functionally graded annular plate subjected to biharmonic transverse forces on its top surface is investigated using the complex variables method. The material parameters are assumed to vary along the thickness direction in an arbitrary fashion. The problem is converted to determine the expressions of four analytic functions α (ζ), β (ζ), ? (ζ) and ψ (ζ) under certain boundary conditions. A series of simple and practical biharmonic loads are presented. The four analytic functions are constructed carefully in a biconnected annular region corresponding to the presented loads, which guarantee the single-valuedness of the mid-plane displacements of the plate. The unknown constants contained in the analytic functions can be determined from the boundary conditions that are similar to those in the plane elasticity as well as those in the classical plate theory. Numerical examples show that the material gradient index and boundary conditions have a significant influence on the elastic field.