Free vibration analysis of open thin deep shells with variable radii of curvature

In the present paper, free vibration of a thin open curved shell with parabolic curvature was studied. This shell has a curvature with variable radius in one direction. The equations of motion of this shell were inferred by first order shell theory. According to perpendicular nature of loading on shell of marine structures, the assumptions of Donnell–Mushtari–Vlasov can be used with an acceptable level of accuracy and the in-plane displacement along shell straight direction “x” can be neglected as compared to the displacement in two other directions. The natural frequencies and mode shapes related to the first five vibrational modes were extracted using semi-analytical methods including power series method, Galerkin method and beam function method. The results of the semi-analytical methods were validated against those obtained by using the finite element method. Out of the studied semi-analytical methods, Galerkin method was found to have an appropriate convergence in both natural frequency and mode shape. Adopting eight terms of the response series, Galerkin method has an appropriate convergence compared with the results of finite element.     

Semi-analytical solution to one-dimensional consolidation in unsaturated soils

In this paper, a series of semi-analytical solutions to one-dimensional consolidation in unsaturated soils are obtained. The air governing equation by Fredlund for unsaturated soils consolidation is simplified. By applying the Laplace transform and the Cayley-Hamilton theorem to the simplified governing equations of water and air, Darcy＇s law, and Fick＇s law, the transfer function between the state vectors at top and at any depth is then constructed. Finally, by the boundary conditions, the excess pore-water pressure, the excess pore-air pressure, and the soil settlement are obtained under several kinds of boundary conditions with the large-area uniform instantaneous loading. By the Crump method, the inverse Laplace transform is performed, and the semi-analytical solutions to the excess pore-water pressure, the excess pore-air pressure, and the soils settlement are obtained in the time domain. In the case of one surface which is permeable to air and water, comparisons between the semi-analytical solutions and the analytical solutions indicate that the semi-analytical solutions are correct. In the case of one surface which is permeable to air but impermeable to water, comparisons between the semi-analytical solutions and the results of the finite difference method are made, indicating that the semi-analytical solution is also correct.     

Accuracy Analysis of the Semi-Analytical Method for Shape Sensitivity Calculation∗

ABSTRACT

The semi-analytical method of design sensitivity analysis that is widely used for calculating derivatives of static response with respect to design variables for structures modeled by finite elements is studied in this paper. It is shown that the method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. Errors are shown to be associated with an incompatibility of the sensitivity field with the structure. An error index is developed to test the accuracy of the semi-analytical method. It characterizes the difference in errors between a general finite difference method and the semi-analytical method. A method for improving the accuracy of the semi-analytical method (when possible) is provided. Examples are presented to demonstrate the use of the error index.     

结构优化半解析灵敏度及误差修正改进算法

提出结构半解析灵敏度分析及其针对刚体位移的误差修正方法的改进算法, 构建灵敏度分析与误差修正项可分离形式. 该方法实现简便, 数值精度不受摄动步长与单元数目的影响. 首先从总体角度推得静力问题的误差修正半解析灵敏度分析方法, 提出了位移误差修正灵敏度列式, 并给出算法实施途径; 然后将此思路推广于自振频率、屈曲临界载荷问题, 提出了相应的计算步骤. 随后, 给出梁单元与壳单元误差修正项的具体推导方法, 并分别使用两种单元构建有限元模型进行算例测试. 结果表明, 该方法适用于多种分析类型, 数值精度不受单元数目与摄动步长的影响. 由于灵敏度分析与误差修正项可以分开计算, 该方法支持将误差修正项直接叠加于灵敏度求解结果进行误差修正, 使已有灵敏度分析程序得到充分利用. 尤其对于复杂工程结构的优化设计, 特别是形状优化设计以及尺寸、形状混合优化设计, 相比于原误差修正方法, 实现更为简便, 效率有所提升, 能为半解析灵敏度分析方法及其程序实现提供新的思路.      

Wave scattering by a floating porous elastic plate of arbitrary shape: A semi-analytical study

In this paper, a semi-analytical model based on linear potential flow theory and an eigenfunction expansion method is developed to study wave scattering by a porous elastic plate with arbitrary shape floating in water of finite depth. The water domain is divided into the interior and exterior regions, corresponding to the domain beneath the plate and the rest extending towards infinite distance horizontally, respectively. The unknown coefficients in the potential expressions are determined by satisfying the continuity conditions for pressure and velocity at the interface of the two regions, together with the conditions for the motion/force at the edge of the plate, where the Fourier series expansion method is employed to deal with the terms associated with the radius function. A plate with three shapes – circular, cosine and elliptical – and three edge conditions are considered. We find that the largest deflection of the plate with a simply supported edge and a clamped edge is more sensitive to the change in porosity when the porosity is small. The power dissipated by an elliptical plate with its major axis perpendicular to the incident wave direction is the largest among the case studies for the majority of the porosity values tested.     

A semi-analytical elastic stress–displacement solution for notched circular openings in rocks

A semi-analytical plane elasticity solution of the circular hole with diametrically opposite notches in a homogeneous and isotropic geomaterial is presented. This solution is based on: (i) the evaluation of the conformal mapping function of a hole of prescribed shape by an appropriate numerical scheme and (ii) the closed-form solutions of the Kolosov–Muskhelishvili complex potentials. For the particular case of circular notches––which resemble to the circular cavity breakout in rocks––it is demonstrated that numerical results pertaining to boundary stresses and displacements predicted by the finite differences model FLAC^2D, as well as previous analytical results referring to the stress-concentration-factor, are in agreement with analytical results. It is also illustrated that the solution may be easily applied to non-rounded diametrically opposite notch geometries, such as “dog-eared” breakouts by properly selecting the respective conformal mapping function via the methodology presented herein. By employing a stress-mean-value brittle failure criterion that takes into account the stress-gradient effect in the vicinity of the curved surfaces in rock as well as the present semi-analytical solution, it is found that a notched hole, e.g. borehole or tunnel breakout, may exhibit stable propagation. The practical significance of the proposed solution lies in the fact that it can be used as a quick-solver for back-analysis of borehole breakout images obtained in situ via a televiewer for the estimation of the orientation and magnitude of in situ stresses and of strain–stress measurements in laboratory tests.     

功能梯度材料板件三维分析的半解析梯度有限元法

将半解析有限元与梯度有限元相结合,形成一种半解析梯度有限元来求解功能梯度材料板件问题。该方法兼有有限元法的适应性强、程序统一,半解析有限元法的节省单元与计算工作量,梯度有限元法的适应构件内部材料性能任意梯度分布等特点,并实现用一维数值计算给出构件三维分析结果。算例分析表明了方法的精度、功能与上述特点,充分揭示了功能梯度材料板件力学响应的三维形态。半解析梯度有限元法可推广应用到其他功能梯度材料面结构的各类分析中。     

Eigenproblems of two-dimensional acoustic cavities with smoothly varying boundaries by using the generalized multipole method

This paper presents a semi-analytical approach to solve the eigenproblem of a two-dimensional acoustic cavity with smoothly varying boundaries. The multipole expansion for the acoustic pressure is formulated in terms of Bessel and Hankel functions to satisfy the Helmholtz equation in the polar coordinate system. Rather than using the addition theorem, the multipole method and directional derivative are both combined to propose a generalized multipole method in which the acoustic pressure and its normal derivative with respect to non-local polar coordinates can be calculated. The boundary conditions are satisfied by uniformly collocating points on the boundaries. By truncating the multipole expansion, a finite linear algebraic system is acquired. The direct searching approach is applied to identify the natural frequencies using the singular value decomposition technique. Several numerical examples are presented, including those of an annulus cavity, a confocal elliptical annulus cavity and an arbitrarily shaped cavity with an inner elliptical boundary. The accuracy and numerical convergence of the proposed method are validated by comparison with results of the available analytical method and the commercial finite-element code ABAQUS. No spurious eigensolutions are found in the proposed formulation. Due to its semi-analytical character, excellent accuracy and fast rate of convergence are the main features of the proposed method.     

复合材料正交加筋板在横向载荷下挠度计算的半解析法

本文采用一种新的半解析法,即独特利用Heaviside函数建立与加筋板等效的变刚度模型来开展复合材料双向正交加筋板在横向载荷下的弯曲挠度分析．此模型可以准确地描述筋条在板面上的分布,以及由于筋条的存在而导致的板面刚度不均匀分布．使用Galerkin加权残值法求解该模型的控制方程,得到不同边界条件和载荷情况下的级数解．对于双向正交加筋板,将此半解析法的结果与传统均匀化方法和使用商业有限元软件ABAQUS建立的有限元模型所得到的弯曲挠度结果比较,验证了此方法的准确性和优越性．不同于传统均匀化方法,本双向正交加筋板的弯曲挠度半解析法可精确、有效地获取加筋间的局部弯曲挠度,可以促进复合材料结构的设计分析与优化的研究进展．     

A method of analysis for end and transitional effects in anisotropic cylinders

We present a rigorous verification study and an extension to an existing semi-analytical finite element formulation for analysis of end and transition effects in prismatic cylinders. End and transition effects in stressed cylinders are phenomena associated with the difference between results that are predicted by the Saint–Venant solutions and the actual point-wise conditions. These differences manifest themselves as self-equilibrated stress states. Notwithstanding certain well-known exceptions (e.g., restrained torsion of open thin-walled sections), such effects in isotropic cylinders are usually confined to a very small neighborhood of a terminal boundary or transition zone, and are typically neglected. For anisotropy, as in the case of most smart/active and composite material systems, they can persist much further into the interior of the structure, and need to be quantified to design geometry transition zones and to fully understand the delamination effects. In the semi-analytical approach, we first discretize the governing equations within the cross-sectional plane of the cylinder. The end-solution fields satisfy the homogeneous form of the resulting semi-analytical system of ordinary differential equations. This leads to an algebraic eigenvalue problem, and an eigenfunction expansion of the stress and displacement fields due to end effects. Unique to the present study, we formulate a procedure to quantify the transitional effects for end-to-end connected cylinders for which the displacement and stress continuity along the transition interface need to be enforced. The semi-analytical approach has several distinct advantages: (i) It is computationally efficient, as only the cross-sectional geometry is discretized; (ii) it can be applied to arbitrary cross-sectional geometries and the most general form of anisotropy; and (iii) it yields direct measures for the decay lengths (or decay rates) of any end-or transition-solution field. Analytical solutions to end-effect problems are scarce. Those that exist are for simple geometry and material constitution. We use these analytical solutions, as well as solutions obtained using three-dimensional finite element models, to verify our approach and to assess its efficiency.     

Transient swelling-induced finite bending of hydrogel-based bilayers: analytical and FEM approaches

Hydrogels with their time-dependent intrinsic behaviors have recently been used widely in soft structures as sensors/actuators. One of the most interesting structures is the bilayer made up of hydrogels which may undergo swelling-induced bending. In this work, by proposing a semi-analytical method, the transient bending of hydrogel-based bilayers is investigated. Utilizing nonlinear solid mechanics, a robust semi-analytical solution is developed which captures the transient finite bending of hydrogel-based bilayers. Moreover, the multiphysics model of the hydrogels is implemented in the finite element method (FEM) framework to verify the developed semi-analytical procedure results. The effects of different material properties are investigated to illustrate the nonlinear behavior of these structures. The von-Mises stress contour extracted from FEM shows that the critical area of these soft structures is at the interface of the layers which experiences the maximum stress, and this area is most likely to rupture in large deformations.     

Parametric study of the distribution of the tensile stresses in pavilion structures constituted by four sectors of barrel shells

A semi-analytical approach is followed in order to obtain an approximate solution for an analytical model describing the static behaviour of pavilion shells. With the aim to parametrically investigate the distribution of the linear elastic tensile stresses, a linear elastic isotropic behaviour of the material is considered. Inside the pavilion shell, a family of arches and a family of straight beams can be recognized. This assumption justifies the separation of the variables which is at the basis of the proposed semi-analytical approach. The sensitivity of the behaviour of the shells to the values of the mechanical and geometrical parameters characterizing the system is investigated. Comparisons with results obtained by finite element models are performed to confirm both the validity of the semi-analytical approach and the provided results. Since only a linear elastic isotropic material has been taken into account, the results do not claim to describe the behaviour of masonry pavilion vaults.     

Buckling of a coating strip of finite width bonded to elastic half-space

A solution for buckling of a stiff strip of finite width bonded to a compliant elastic half-space and subjected to uniform axial compression is presented. Approximate semi-analytical and finite element solutions are obtained and compared with a two-dimensional case of a plate on elastic foundation. The comparison demonstrates that the two-dimensional solution can be applied to predict the buckling wavelength and critical compressive strain when the width of the strip is equal to or larger than the buckling wavelength. For narrow strips, the wavelength is smaller and critical strain is higher than that of a plate on foundation.     

On the circular footing plates on two-parameters foundation under arbitrary loads

Modelling soils by two-parameter foundation model, this paper calculates the distributions of displacements of circular footing plates on soils and reactions of soils under arbitrary loads using semi-analytical finite element method. And it improves F.Z. Vlazov’s solution in the case of axisymmetry. The results agree well in comparison with those by F E M. At the same time, the boundary conditions of circular plates on soils are discussed.     

Nonlinear vibration analysis of isotropic cantilever plate with viscoelastic laminate

The nonlinear vibration of an isotropic cantilever plate with viscoelastic laminate is investigated in this article. Based on the Von Karman’s nonlinear geometry and using the methods of multiple scales and finite difference, the dimensionless nonlinear equations of motion are analyzed and solved. The solvability condition of nonlinear equations is obtained by eliminating secular terms and, finally, nonlinear natural frequencies and mode-shapes are obtained. Knowing that the linear vibration of this type of plate does not have exact solution, Ritz method is employed to obtain semi-analytical nonlinear mode-shapes of transverse vibration of this plate. Airy stress function and Galerkin method are employed to reduce nonlinear PDEs into an ODE of duffing type. Stability of plate and chaotic behavior are investigated by Runge–Kutta method. Poincare section diagrams are in good agreement with results of Lyapunov criteria.     

FINITE LAYER ANALYSIS FOR SEMI-INFINITE SOILS (Ⅰ)——GENERAL THEORY

In the present paper a finite layer method is studied for the elastodynamics of transverse isotropic bodies. With this method, semi-infinite soils can be considered as an transverse isotropic half-space, its material functions varying with depth. Dividing the half.space into a series of layers in the direction of depth the material fimetioms in each layer are simulated by exponential fumctions Consequently, the fundamental equations to be solved can be simplified if the fouricr transform with repsect to coordinates is used. We have obtained the relationship between the "layer forces" and "layer displacements". This finite layer method, in fact, can also be called a semi-analytical method. It possesses those advantages as the usual semi-analytical methods do, and can be used to analyse the problem of the interaction between soils and structures.     

ELASTO-PLASTIC COUPLED ANALYSIS OF BURIED STRUCTURE AND SOIL MEDIUM BY PERTURBATIONAL SEMI-ANALYTIC METHOD

I.IntroductionMoreandmoreelasto-plasticanalysesareusedinCivilEngineering.TheMohr-Coulombyieldconditionisoftenappliedtorockandsoilmediumtostudythesituationofdestruction.Theperturbationalmethodwasfirstusedinthestudyofcelestialmechanics,thenithasbeenusedinsolidmechanicstocarryoutthegeometricallyandphysicallynonlinearanalysis[']'12].Forlarge?scalcproblems,thegeneralnumericalmethodsmaynotbeapttothiscomplexcalculation,butsemi-analyticalmethodscandoitwell.Touseanalyticalfunctionsinoneortwodirectio…     

Axisymmetric deformations of circular rings made of linear and Neo-Hookean materials under internal and external pressure: A benchmark for finite element codes

The axisymmetric deformations of thick circular rings are investigated. Four materials are explored: linear material, incompressible Neo-Hookean material and Ogden's and Bower's forms of compressible Neo-Hookean material. Radial distributed forces and a displacement-dependent pressure are the external loads. This problem is relatively simple and allows analytical, or semi-analytical, solution; therefore it has been chosen as a benchmark to test commercial finite element software for various material laws at large strains. The solutions obtained with commercial finite element software are almost identical to the present semi-analytical ones, except for the linear material, for which commercial finite element programs give incorrect results.     

基于黏弹性的微注射成型流动模拟

应用有限元方法研究了微注射成型中瞬态、可压缩、非牛顿熔体流动的黏弹性对流动前沿及流动平衡的影响。基于Phan-Thien-Tanner模型建立了熔体流动的本构方程,利用Hele-Shaw假设和简化建立了瞬态、可压缩、非牛顿熔体流动的连续性方程、动量方程、能量方程;为了有效地描述微注射成型的尺寸效应,采用了边界滑移和表面张力边界条件。通过分部积分和待定系数法导出了带有边界信息的变分方程和求解应力分量的半解析公式,构造了有限元离散求解及超松驰迭代算法。模拟结果表明:熔体的黏弹性对浇口附近的压力和后续的熔体流动前沿有重要影响;与黏性模型相比,黏弹性模型可以控制模拟压力的快速增长,减少不同型腔之间的充填差异,与短射实验结果也更吻合。     

A semi-analytical approach for the non-linear large deflection analysis of laminated rectangular plates under general out-of-plane loading

A semi-analytical approach for the geometrically non-linear analysis of rectangular laminated plates with general inplane and out-of-plane boundary conditions under a general distribution of out-of-plane loads is developed. The analysis is based on the elastic thin plate theory with geometrically non-linear von Kármán strains. The solution of the non-linear partial differential equations is reduced to an iterative sequential solution of non-linear ordinary differential equations using the multi-term extended Kantorovich method. The efficiency, accuracy, and convergence of the proposed method are examined through a comparison with other semi-analytical methods and with finite element analyses. The capabilities of the approach and its applicability to the non-linear large deflection analysis of plate structures are demonstrated through various numerical examples. Emphasis is placed on combinations of lamination, boundary, and loading conditions that cannot be analyzed using alternative semi-analytical methods.