基于赫林格?赖斯纳变分原理的一致高效无网格本质边界条件施加方法

无网格法具有高阶连续光滑的形函数, 在结构分析中呈现出显著的精度优势. 但无网格形函数在节点处一般没有插值性, 导致伽辽金无网格法难以直接施加本质边界条件. 采用变分一致尼兹法施加边界条件的数值解具有良好的收敛性和稳定性, 因而得到了非常广泛的应用, 然而该方法仍然需要引入人工参数来保证算法的稳定性. 本文以赫林格?赖斯纳变分原理为基础, 建立了一种变分一致的本质边界条件施加方法. 该方法采用混合离散近似赫林格?赖斯纳变分原理弱形式中的位移和应力, 其中位移采用传统无网格形函数进行离散, 而应力则在背景积分单元中近似为相应阶次的多项式. 此时的无网格离散方程可视为一种新型的尼兹法施加本质边界条件, 其中修正变分项采用再生光滑梯度和无网格形函数进行混合离散, 稳定项则内嵌于赫林格?赖斯纳变分原理弱形式中, 无需额外增加稳定项, 消除了对人工参数的依赖性. 该方法无需计算复杂耗时的形函数导数, 并满足积分约束条件, 保证了数值求解的精度. 数值结果表明, 所提方法能够保证伽辽金无网格法的计算精度最优误差收敛率, 与传统的尼兹法相比明显提高了计算效率.      

二阶一致无网格与有限元耦合离散研究

为准确方便地施加本质边界条件,在连续掺混法（Continuous Blending Method, CBM）的框架下,通过增加一个边中节点,发展了采用二阶基底的无网格与二阶有限元的耦合离散方法。Galerkin弱形式的数值积分采用具二阶一致性的3点积分方法（Quadratically Consistent 3-point integration method,QC3）。与原本在QC3中采用的Nitsche法相比,所发展的耦合离散方法可像有限元法一样简单高效地施加本质边界条件,不向弱形式中引入额外项,也不依赖于任何人工参数。而且,数值结果还表明,QC3的计算精度也得到进一步提高。     

A meshfree method for transverse vibrations of anisotropic plates

A meshfree approach, called displacement boundary method, for anisotropic Kirchhoff plate dynamic analysis is presented. This method is deduced from a variational principle, which uses a modified hybrid functional involving the generalized displacements and generalized tractions on the boundary and the lateral deflection in the domain as independent variables. The discretization process is based on the employment of the fundamental solutions of the static problem operator for the expression of the variables involved in the functional. The stiffness and mass matrices obtained for the dynamic model are frequency-independent, symmetric and positive definite and their computation involves boundary integrals of regular kernels only. Due to its features, the final resolving system can be solved with the classical approaches by using standard numerical procedures. To assess the formulation, the free vibrations of some anisotropic plates were calculated and the results compared with those obtained using other solution techniques. The present results are in good agreement with those found in the literature showing the accuracy and effectiveness of the proposed approach.     

Extended meshfree analysis of transverse and inplane loading of a laminated anisotropic plate of general planform geometry

An extended meshfree method is presented for the analysis of a laminated anisotropic plate under elastostatic loading. The plate may be of any planform shape with its thickness profile composed of perfectly bonded uniform thickness layers of distinct anisotropic materials. Both transverse and inplane loads are considered using a first order shear deformation theory for flexural behavior and generalized plane stress for the membrane behavior. In this extended meshfree method, a rectangular domain is initially considered with the plate of arbitrary geometry inscribed within it. A particular solution in the form of an analytic generalized Navier solution (a compound double Fourier series) is used to capture the response due to the loading within the rectangular domain. Then, a homogeneous solution by meshfree analysis is added to treat the augmented boundary conditions on the actual contour of the plate. These augmented conditions are composed of the prescribed values and that of the particular solution evaluated around the plate’s contour.Concentrated transverse and inplane loads in the form of uniform loads over a very small patch are considered with this generalized Navier solution representation. When a meshfree portion is added to account for the boundary conditions, such solutions constitute the Green’s functions for the plate. The viability of these double Fourier series representations is shown by the convergence rates for the kinematic and force/moment fields. An additional example of a two layer ±30° angleply circular plate is given to illustrate the capability of this extended meshfree method.     

等几何分析中采用Nitsche法施加位移边界条件

等几何分析使用NURBS基函数统一表示几何和分析模型,消除了传统有限元的网格离散误差,容易构造高阶连续的协调单元.对于结构分析,选择合适的几何参数可以得到光滑的应力解,避免了后置处理的应力磨平.但是由于NURBS基函数不具备插值性,难以直接施加位移边界条件.针对这一问题,提出一种基于Nitsche变分原理的边界位移条件"弱"处理方法,它具有一致稳定的弱形式,不增加自由度,方程组对称正定和不会产生病态矩阵等优点.同时给出方法的稳定性条件,并通过求解广义特征值问题计算稳定性系数.最后,数值算例表明Nitsche方法在h细化策略下能获得最优收敛率,其结果要明显优于在控制顶点处直接施加位移约束.     

IMPOSING DISPLACEMENT BOUNDARY CONDITIONS WITH NITSCHE＇S METHOD IN ISOGEOMETRIC ANALYSIS

等几何分析使用 NURBS 基函数统一表示几何和分析模型, 消除了传统有限元的网格离散误差, 容易构造高阶连续的协调单元. 对于结构分析, 选择合适的几何参数可以得到光滑的应力解, 避免了后置处理的应力磨平. 但是由于 NURBS 基函数不具备插值性, 难以直接施加位移边界条件. 针对这一问题, 提出一种基于 Nitsche 变分原理的边界位移条件“弱”处理方法, 它具有一致稳定的弱形式, 不增加自由度, 方程组对称正定和不会产生病态矩阵等优点. 同时给出方法的稳定性条件, 并通过求解广义特征值问题计算稳定性系数. 最后, 数值算例表明 Nitsche 方法在h细化策略下能获得最优收敛率, 其结果要明显优于在控制顶点处直接施加位移约束.}     

Strong-form framework for solving boundary value problems with geometric nonlinearity

In this paper, we present a strong-form framework for solving the boundary value problems with geometric nonlinearity, in which an incremental theory is developed for the problem based on the Newton-Raphson scheme. Conventionally, the finite element methods (FEMs) or weak-form based meshfree methods have often been adopted to solve geometric nonlinear problems. However, issues, such as the mesh dependency, the numerical integration, and the boundary imposition, make these approaches computationally inefficient. Recently, strong-form collocation methods have been called on to solve the boundary value problems. The feasibility of the collocation method with the nodal discretization such as the radial basis collocation method (RBCM) motivates the present study. Due to the limited application to the nonlinear analysis in a strong form, we formulate the equation of equilibrium, along with the boundary conditions, in an incremental-iterative sense using the RBCM. The efficacy of the proposed framework is numerically demonstrated with the solution of two benchmark problems involving the geometric nonlinearity. Compared with the conventional weak-form formulation, the proposed framework is advantageous as no quadrature rule is needed in constructing the governing equation, and no mesh limitation exists with the deformed geometry in the incremental-iterative process.     

Effect of magnetic field and non-homogeneity on the radial vibrations in hollow rotating elastic cylinder

In this paper, the natural frequencies of the radial vibrations of a hollow cylinder with different boundary conditions under influences of magnetic field, rotation and non-homogeneity have been studied. The solution of the problem is obtained by using technique of variables separation. In the present paper three different boundary conditions are considered, namely the free, fixed and mixed boundary conditions. The displacement and stresses components have been obtained in analytical form involving Bessel function of first and second kind and of order n. The determination is concerned with the eigenvalues of the natural frequencies of the radial vibrations for different boundary conditions in the case of harmonic vibrations. Numerical results are given and illustrated graphically for each case considered. Comparisons are made with the results in the absence of magnetic field, rotation and non-homogeneity. The results indicate that the effect of magnetic field, rotation and non-homogeneity are very pronounced.     

Fast approximations of dynamic stability boundaries of slender curved structures

Curved beams and panels can often be found as structural components in aerospace, mechanical and civil engineering systems. When curved structures are subjected to dynamic loads, they are susceptible to dynamic instabilities especially dynamic snap-through buckling. The identification of the dynamic stability boundary that separate the non-snap and post-snap responses is hence necessary for the safe design of such structures, but typically requires extensive transient simulations that may lead to high computation cost. This paper proposes a scaling approach that reveals the similarities between dynamic snap-through boundaries of different structures. Such identified features can be directly used for fast approximations of dynamic stability boundaries of slender curved structures when their geometric parameters or boundary conditions are varied. The scaled dynamic stability boundaries of half-sine arches, parabolic arches and cylindrical panels are studied.     

NEW TREATMENT OF ESSENTIAL BOUNDARY CONDITIONS IN EFG METHOD BY COUPLING WITH RPIM

One of major difficulties in the implementation of meshfree methods using the moving least square (MLS) approximation, such as element-free Galerkin method (EFG), is the imposition of essential boundary conditions as the approximations do not pass through the nodal parameter values. Another class of meshfree methods based on the radial basis point interpolation can satisfy the essential boundary conditions exactly since its approximation function passes through each node in an influence domain and thus its shape functions possess the properties of delta function. In this paper, a coupled element-free Galerkin(EFG)-radial point interpolation method (RPIM) is proposed to enhance their advantages and avoid their disadvantages. Discretized equations of equilibrium are obtained in the RPIM region and the EFG region, respectively. Then a collocation approach is introduced to couple the RPIM and the EFG method. This method satisfies the linear consistency exactly and can maintain the stiffness matrix symmetric. Numerical tests show that this method gives reasonably accurate results consistent with the theory.     

Influence of Anisotropy on the Response Characteristics of Finite Cylinders Under Free Vibrations

A numerical–analytical approach is proposed to solve a problem on the free vibrations of cylindrical bodies. The approach is based on three-dimensional elastic theory and the semianalytic finite-element method. The free vibrations of isotropic and anisotropic solid cylinders of finite length are examined. It is studied how boundary conditions and mechanical and geometrical parameters affect the distribution of dynamic properties. The efficiency of the approach proposed is tested by comparing results produced by different approaches     

Nonlinear modeling of tuned liquid dampers (TLDs) in rotating wind turbine blades for damping edgewise vibrations

Tuned liquid dampers (TLDs) utilize the sloshing motion of the fluid to suppress structural vibrations and become a natural candidate for damping vibrations in rotating wind turbine blades. The centrifugal acceleration at the tip of a wind turbine blade can reach a magnitude of 7–8g. This facilitates the use of a TLD with a relatively small fluid mass and with feasible geometric dimensions to mitigate the lightly-damped edgewise vibrations effectively. In the present paper, modal expansions are carried out directly on the velocity field and the free surface of the sloshing liquid in the rotating coordinate system. A formulation has been proposed leading to coupled nonlinear ordinary differential equations, which have been obtained through the Galerkin variational approach together with the modal expansion technique. Two models, with one sloshing mode and three sloshing modes, have been studied in the numerical simulation. It is shown that the one-mode model is able to predict the sloshing force and the damped structural response accurately, since the primary damping effect on the structure is achieved by the first sloshing mode of the fluid. Although it is unable to predict the fluid free-surface elevation equally well, the one-mode model can still be utilized for the design of TLD. Parametric optimization of the TLD is carried out based on the one-mode model, and the optimized damper effectively improves the dynamic response of wind turbine blades.     

Non-linear vibrations of shallow circular cylindrical panels with complex geometry. Meshless discretization with the R-functions method

Geometrically non-linear forced vibrations of a shallow circular cylindrical panel with a complex shape, clamped at the edges and subjected to a radial harmonic excitation in the spectral neighborhood of the fundamental mode, are investigated. Both Donnell and the Sanders–Koiter non-linear shell theories retaining in-plane inertia are used to calculate the elastic strain energy. The discrete model of the non-linear vibrations is build using the meshfree technique based on classic approximate functions and the R-function theory, which allows for constructing the sequences of admissible functions that satisfy given boundary conditions in domains with complex geometries; Chebyshev orthogonal polynomials are used to expand shell displacements. A two-step approach is implemented in order to solve the problem: first a linear analysis is conducted to identify natural frequencies and corresponding natural modes to be used in the second step as a basis for expanding the non-linear displacements. Lagrange approach is applied to obtain a system of ordinary differential equations on both steps. Different multimodal expansions, having from 15 up to 35 generalized coordinates associated with natural modes, are used to study the convergence of the solution. The pseudo-arclength continuation method and bifurcation analysis are applied to study non-linear equations of motion. Numerical responses are obtained in the spectral neighborhood of the lowest natural frequency; results are compared to those available in the literature. Internal resonances are also detected and discussed.     

Forced vibrations of orthotropic shells: nonclassical boundary-value problems

The forced vibrations of a cylindrical orthotropic shell are studied. Two types of boundary conditions on the outer surface are examined considering that the displacement vector prescribed on the inner surface varies harmonically with time. Asymptotic solutions of associated dynamic equations of three-dimensional elasticity are found. Amplitudes of forced vibrations are determined and conditions under which resonance occurs are established. Boundary-layer functions are defined. The rate of their decrease with distance from the ends inside the shell is determined. A procedure of joining solutions for the internal boundary-layer problem is outlined in the case for the, if clamping boundary conditions are prescribed at the ends     

几何精确NURBS有限元中边界条件施加方式对精度影响的三维计算分析

非均匀有理B样条(NURBS)有限元法把计算机辅助几何设计(CAGD)中的NURBS几何构形方法与有限元方法有机结合起来,有效消除了有限元离散模型的几何误差,提高了计算精度。但是由于NURBS基函数不是插值函数,直接在控制节点上施加位移边界条件会引起较大误差。本文详细讨论了NURBS基函数的插值特性,在NURBS有限元分析中采用罚函数法施加位移边界条件,提高了收敛率和计算精度。结合典型三维弹性力学问题,对两种施加位移边界条件的方法进行了对比和分析。计算结果表明,直接施加位移边界条件会导致收敛率和精度的明显降低,而基于罚函数法的NURBS有限元分析则能达到最优收敛率,并具有更高的精度。     

Study of static and dynamic stability of flexible rods in a geometrically nonlinear statement

We study static and dynamic stability problems for a thin flexible rod subjected to axial compression with the geometric nonlinearity explicitly taken into account. In the case of static action of a force, the critical load and the bending shapes of the rod were determined by Euler. Lavrent’ev and Ishlinsky discovered that, in the case of rod dynamic loading significantly greater than the Euler static critical load, there arise buckling modes with a large number of waves in the longitudinal direction. Lavrent’ev and Ishlinsky referred to the first loading threshold discovered by Euler as the static threshold, and the subsequent ones were called dynamic thresholds; they can be attained under impact loading if the pulse growth time is less than the system relaxation time. Later, the buckling mechanism in this case and the arising parametric resonance were studied in detail by Academician Morozov and his colleagues.In this paper, we complete and develop the approach to studying dynamic rod systems suggested by Morozov; in particular, we construct exact and approximate analytic solutions by using a system of special functions generalizing the Jacobi elliptic functions. We obtain approximate analytic solutions of the nonlinear dynamic problem of flexible rod deformation under longitudinal loading with regard to the boundary conditions and show that the analytic solution of static rod system stability problems in a geometrically nonlinear statement permits exactly determining all possible shapes of the bent rod and the complete system of buckling thresholds. The study of approximate analytic solutions of dynamic problems of nonlinear vibrations of rod systems loaded by lumped forces after buckling in the deformed state allows one to determine the vibration frequencies and then the parametric resonance thresholds.     

Natural vibrations of thin conical panels of variable thickness

A numerical analytical approach is proposed to study the natural vibrations of thin isotropic conical shells with varying thickness. The approach is based on the spline-approximation of unknown functions. Open shells (panels) with different boundary conditions are considered. The effect of variable thickness on the dynamic characteristics is analyzed     

On modified displacement version of the non-linear theory of thin shells

We discuss the non-linear theory of thin shells expressed in terms of displacements of the shell reference surface as the only independent field variables. The formulation is based on the principle of virtual work postulated for the reference surface. In our approach: (1) the vector equilibrium equations are represented through components in the deformed contravariant surface base, and using the compatibility conditions the resulting tangential equilibrium equations are additionally simplified, (2) at the shell boundary the new scalar function of displacement derivatives is defined and new sets of four work-conjugate static and geometric boundary conditions are derived, as well as (3) for prescribed shell geometry all non-linear shell relations are generated automatically by two packages set up in Mathematica. The displacement boundary value problem and the associated homogeneous shell buckling problem are generated exactly without using any additional approximations following from errors of the constitutive equations. Both problems are extremely complex and available only in the computer memory. Such an approach allows us to account also for those a few supposedly small terms which may be critical for finding the correct buckling load of shells sensitive to imperfections. This approach is used in the accompanying paper by Opoka and Pietraszkiewicz [Opoka, S., Pietraszkiewicz, W., 2009. On refined analysis of bifurcation buckling for the axially compressed circular cylinder. International Journal of Solids and Structures, 46, 3111–3123.] to perform the refined numerical analysis of bifurcation buckling for the axially compressed circular cylinder.     

Effect of ventricle motion on the dynamic behaviour of chorded mitral valves

An Immersed Boundary (IB) model is employed to investigate the dynamic behaviour of a novel chorded mitral prosthesis, which is in the early stages of its development, under physiological flow conditions. In vivo magnetic resonance images (MRIs) of the left ventricle are analysed to determine the relative motion of the mitral annulus and the papillary muscle regions of the ventricle. The dynamic boundary conditions are incorporated into IB simulations to test the valve in a more realistic dynamic geometric environment. The IB model has successfully identified the effect of the dynamic boundary conditions on the mechanical behaviour of the valve and revealed the strengths and weaknesses of the current mitral design. The mechanical performance of the prosthesis is compared with recent studies of native porcine valves; differences in mechanical behaviour are observed. Potential improvements for the design of the prosthesis are proposed.     

无网格动力分析的循环卷积神经网络代理模型

在无网格动力分析中,除了无网格形函数本身构造复杂引入的计算成本,还需要逐步递推求解每个时间步的动力响应,因而计算效率较为低下.本文通过研究无网格离散数据与机器学习训练样本、无网格动力分析递推计算过程与循环卷积神经网络序列信息传递模式之间的本征联系,构建了与无网格法相匹配的循环卷积神经网络设计方法,进而提出了一种无网格动...