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

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

Kirchhoff-Love壳的等几何分析

论文利用等几何分析研究了基于Kirchhoff-Love理论的薄壳的静态问题.等几何分析采用等参思想,将精确描述几何形状的NURBS基函数同时作为场变量的插值函数,保证了在分析和网格优化过程中模型的几何精确性,并可以轻易地构造任意高阶连续的单元.该方法具有很高的数值精度.计算结果表明,在等几何分析中,NURBS单元的阶次越高,网格数越多,计算结果越精确.     

基于等几何分析的比例边界有限元方法

提出了一种具有比例边界有限元的半解析特性和等几何分析的几何特性的新方法。该新方法是在比例边界有限元框架中用NURBS曲线或曲面精确描述域边界几何形状,同时域边界位移场采用描述几何形状的NURBS形函数等参构造。这种新方法具有比例边界有限元固有的径向解析特性和NURBS的高阶连续性的优点。数值算例显示,与传统的比例边界有限元相比,基于等几何分析的比例边界有限元方法提高了域边界单元和域内应力场的连续性,减少了计算自由度。应用此方法可以用较少的计算自由度获得更高连续阶和更高精度的位移、应力和应变场。     

节点梯度光滑有限元配点法

配点法构造简单、计算高效, 但需要用到数值离散形函数的高阶梯度,而传统有限元形函数的梯度在单元边界处通常仅具有C$^{0}$连续性,因此无法直接用于配点法分析. 本文通过引入有限元形函数的光滑梯度,提出了节点梯度光滑有限元配点法. 首先基于广义梯度光滑方法,定义了有限元形函数在节点处的一阶光滑梯度值,然后以有限元形函数为核函数构造了有限元形函数的一阶光滑梯度,进而对一阶光滑梯度直接求导并用一阶光滑梯度替换有限元形函数的标准梯度,即完成了有限元形函数二阶光滑梯度的构造.文中以线性有限元形函数为基础的理论分析表明,其光滑梯度不仅满足传统线性有限元形函数梯度对应的一阶一致性条件,而且在均布网格假定下满足更高一阶的二阶一致性条件.因此与传统线性有限元法相比,基于线性形函数的节点梯度光滑有限元法的$L_{2}$和$H_{1}$误差均具有二次精度,即其$H_{1}$误差收敛阶次比传统有限元法高一阶, 呈现超收敛特性.文中通过典型算例验证了节点梯度光滑有限元配点法的精度和收敛性,特别是其$H_{1}$或能量误差的精度和收敛率都明显高于传统有限元法.      

基于IGA-SIMP法的连续体结构应力约束拓扑优化

建立了一种IGA-SIMP框架下的连续体结构应力约束拓扑优化方法。基于常用的SIMP模型,将非均匀有理B样条（NURBS）函数用于几何建模、结构分析和设计参数化,实现了结构分析和优化设计的集成统一。利用高阶连续的NURBS基函数,等几何分析（IGA）提高了结构应力及其灵敏度的计算精度,增加了拓扑优化结果的可信性。为处理大量局部应力约束,提出了基于稳定转换法修正的P-norm应力约束策略,以克服拓扑优化中的迭代振荡和收敛困难。通过几个典型平面应力问题的拓扑优化算例表明了本文方法的有效性和精确性。应力约束下的体积最小化设计以及体积和应力约束下的柔顺度最小化设计的算例表明,基于稳定转换法修正的约束策略可以抑制应力约束体积最小化设计中的迭代振荡现象,获得稳定收敛的优化解;比较而言,体积和应力约束下的柔顺度最小化设计的迭代过程更加稳健,适合采用精确修正的应力约束策略。     

基于改进位移模式的一维C1有限元超收敛算法

提出了基于改进位移模式的一维C1有限元超收敛算法。利用单元内部需满足平衡方程的条件,推导了超收敛计算解析公式的显式,即将高阶有限元解的位移模式用常规有限元解的位移模式表示。用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式。采用积分形式推导了单元刚度矩阵。该算法在前处理阶段使用了超收敛计算公式,在常规试函数的基础上,增加了高阶试函数,使得单元内平衡方程的残差减少,从而达到提高精度的目标。对于Hermite单元,本文的结点和单元的位移、导数都达到了h4阶的超收敛精度。     

基于IGA-SIMP法的连续体结构应力约束拓扑优化

建立了一种IGA-SIMP框架下的连续体结构应力约束拓扑优化方法。基于常用的SIMP模型,将非均匀有理B样条(NURBS)函数用于几何建模、结构分析和设计参数化,实现了结构分析和优化设计的集成统一。利用高阶连续的NURBS基函数,等几何分析(IGA)提高了结构应力及其灵敏度的计算精度,增加了拓扑优化结果的可信性。为处理大量局部应力约束,提出了基于稳定转换法修正的P-norm应力约束策略,以克服拓扑优化中的迭代振荡和收敛困难。通过几个典型平面应力问题的拓扑优化算例表明了本文方法的有效性和精确性。应力约束下的体积最小化设计以及体积和应力约束下的柔顺度最小化设计的算例表明,基于稳定转换法修正的约束策略可以抑制应力约束体积最小化设计中的迭代振荡现象,获得稳定收敛的优化解;比较而言,体积和应力约束下的柔顺度最小化设计的迭代过程更加稳健,适合采用精确修正的应力约束策略。     

基于Legendre多项式的双压电层合板三维压电谱元法模拟

谱单元作为一种高阶单元具有计算效率高和精度高的特点。本文在基于Legendre正交多项式的三维谱单元基础上提出了三维压电谱单元模型,用于压电层合板静力和动力性能模拟研究。压电层合板结构中的位移和电势自由度均离散到三维压电谱单元Gauss-Lobatto-Legendre(GLL)配置节点上,并且未对电势沿压电层厚度方向上的变化做任何假定。通过提高谱单元中沿压电层合板厚度方向上的形函数阶数的方法,来削弱三维谱单元在模拟薄板结构中出现的剪切闭锁现象。为验证单元的计算精度,取双压电层合板结构进行静力和动力行为模拟验证,并将计算结果与现有文献中的其它单元模型及有限元结果进行对比。结果表明,压电谱单元可有效模拟压电层合板的静力和动力行为,且提高谱单元形函数阶数可提高数值模拟精度。     

考虑结构自重的基于NURBS插值的3D拓扑描述函数法

在许多如大坝、桥梁等大型土木工程结构中,结构的自重是初始设计阶段必须考虑的重要载荷之一,因此研究自重载荷作用下的结构拓扑优化设计问题具有十分重要的意义.针对考虑自重载荷作用的拓扑优化问题所面临的主要困难,总结了现有处理考虑自重载荷的拓扑优化问题的三类主要方法;提出一种基于非均匀有理B样条(non-uniform rational B-splines,NURBS)基函数插值的拓扑描述函数方法,基于此方法研究了考虑设计依赖自重载荷作用的2D/3D结构优化设计问题.在列式下,高阶NURBS基函数被同时用于三维NURBS实体片中的几何场、位移场及设计变量场插值,实现了几何模型、分析模型和优化模型的有效统一,确保了位移场及设计变量场的高阶连续性;详细推导了基于NURBS基函数插值的考虑自重载荷作用的三维结构拓扑优化模型及其灵敏度列式,并采用移动渐进线方法(method of moving asymptotes,MMA)进行了优化求解;多个算例验证了方法的有效性和稳定性,结果表明,优化迭代过程稳健,收敛快,能够有效地克服自重载荷作用下连续体结构拓扑优化中经常遇到的低密度区域材料的寄生效应及目标函数的非单调性等问题.     

一类指数型间断函数的嵌入非连续等参有限元法

在嵌入非连续有限元的基本思想下,提出一类附加位移形函数———指数型间断函数,来模拟由于非连续结构,如裂纹和节理,所导致的位移不连续规律,该附加函数是以到间断处的垂直距离为自变量,且随距离的增大而呈指数衰减的函数.指数型间断函数具有在数学上的便于积分和求导的优点,且比阶梯间断函数更能反映实际破裂后的变形情况.本文用弱解形式推导了嵌入非连续有限元格式,编制了二维4节点和三维8节点的嵌入非连续等参有限元程序,并分别给出了算例.算例表明在模拟裂纹追踪时,指数型间断函数的嵌入非连续等参有限元法可行且有效.     

平面广义四节点等参元GQ4及其性能探讨

广义节点有限元是将传统有限元方法中的节点广义化,在不增加节点个数的前提下,仅通过提高广义节点的插值函数的阶次,从而达到提高有限元解精度的目的.与现有的p型和hp型有限元不同,在这种新的有限元中,节点自由度全部定义在节点处,在理论与程序实现上与传统有限元方法具有很好的相容性,传统有限元方法是这种新方法的广义节点退化为0阶时的特殊情形.文中主要讨论了这一新方法的四节点等参元（记为GQ4）的形式.对GQ4进行的各种数值试验表明,所发展的广义四节点等参单元具有精度高且无剪切自锁与体积自锁等的特点.     

二阶非自伴两点边值问题Garlerkin有限元的超收敛算法

提出了基于改进位移模式的二阶非自伴两点边值问题Garlerkin有限元的超收敛算法. 用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式,基于Garlerkin 方法,采用积分形式推导了单元平衡方程. 对于线性单元,本文给出了有代表性的算例,结点和单元的位移、导数都达到了h⁴阶的超收敛精度.     

Transient heat conduction analysis using the NURBS-enhanced scaled boundary finite element method and modified precise integration method

The Non-uniform rational B-spline(NURBS)enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper.The scaled boundary finite element method is a semi-analytical technique,which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction.In this method,only the boundary is discretized in the finite element sense leading to a reduction of the spatial dimension by one with no fundamental solution required.Nevertheless,in case of the complex geometry,a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often unavoidable in the conventional finite element approach,which leads to huge computational efforts and loss of accuracy.NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape.In the proposed methodology,the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions,while the straight part of the boundary is discretized by the conventional Lagrange shape functions.Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analysis and the solution is obtained using the modified precise integration method.The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion.Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method.The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.     

A piecewise beam element based on absolute nodal coordinate formulation

The element created in this investigation is based on the it absolute nodal coordinate formulation (ANCF) which has been successfully used in flexible multibody system dynamic and integration of computer aid design and analysis (ICADA). When modeling a B-spline curve with ANCF beam element, it is the common manner to convert this curve into a series of Bézier curves because the systematical conversion between ANCF beam element and a Bézier curve has already been built. In order to avoid the constrain equation produced in this method and to express a B-spline curve using only one element, an alternative approach is developed, leading to the piecewise ANCF (PANCF) beam element. It is demonstrated that when two ANCF beam elements are connected according to a particular continuity, they can constitute a PANCF element. Besides, a new PANCF element will be produced when an ANCF element is connected to an existing PANCF element. The continuity condition can be automatically ensured by the selection of nodal coordinates and the calculation of the piecewise continuous shape functions. The algorithm for converting a B-spline curve to a PANCF beam element is then given. There also are discussions on the features of PANCF element. When two neighboring segments of PANCF element have the same assumed length, the position vector at the interface cannot be expressed by the other coordinates so the position vector is preserved in the \(C^{2}\) continuous situation. Two examples are given to conclude the interpolation and continuity properties of the shape function and to demonstrate the feasibility of this PANCF in the ICADA.     

A layerwise C0-type higher order shear deformation theory for laminated composite and sandwich plates

A novel layerwise C⁰-type higher order shear deformation theory (layerwise C⁰-type HSDT) for the analysis of laminated composite and sandwich plates is proposed. A C⁰-type HSDT is used in each lamina layer and the continuity of in-plane displacements and transverse shear stresses at inner-laminar layer is consolidated. The present layerwise theory retains only seven variables without increasing the number of variables when the number of lamina layers are intensified. The shear stresses through the plate thickness derived from the constitutive equation of the present theory have the same shape as those calculated from the equilibrium equation. In addition, the artificial constraints are added in the principle of virtual displacements (PVD) and are certainly fulfilled through a penalty approach. In this paper, two C⁰-continuity numerical methods, such as the Finite Element Method (FEM) and Bézier isogeometric element (BIEM) are utilized to solve a discrete system of equations derived from the PVD. Several numerical examples with various geometries, aspect ratios, stiffness ratios, and boundary conditions are investigated and compared with the 3D elasticity solution, the analytical, as well as, numerical solutions based on various plate theories.     

基于改进位移模式的一维有限元超收敛算法

将多尺度方法的思想与超收敛计算的解析公式结合起来,提出了改进有限元位移模式的算法。利用超收敛计算的解析公式,将高阶有限元解的位移模式用常规有限元解的位移模式表示。用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式,采用积分形式推导了单元刚度矩阵。该算法在前处理和后处理两个阶段都使用超收敛计算公式,在常规试函数的基础上,增加了高阶试函数,使得单元内平衡方程的残差减少,从而达到提高精度的目标。对于线性单元,本文结点和单元的位移、导数都达到了h4阶的超收敛精度。     

A Direct Method for the Generation of Shape Functions of Any Desired Continuity for Rectangular Finite Elements

Abstract

A direct method is presented for the generation of shape functions for rectangular finite elements, with any desired continuity of the shape functions across the interelement boundaries, i.e., the shape functions can be of any prescribed class C^N, N = 0, 1, 2, 3 …. The method is illustrated with seven examples, four of which are previously published elements. The other three examples represent new elements, two of which are conforming C¹ elements, and the third one is of class C². One of the C1 elements and the C² element reproduce exactly all the terms up to the sextic polynomial.     

Multipatch discontinuous Galerkin isogeometric analysis of composite laminates

This paper is concerned with the isogeometric analysis (IGA) of composite laminates under cylindrical bending. Non-uniform rational B-splines (NURBS) are employed as basis functions for both geometric and computational implementations. In order to account for multiple domains, each lamina is modeled as a single NURBS patch. This multipatch representation corresponds to decomposition of the computational domain (composite laminate) into non-overlapping subdomains. As NURBS patches are discontinuous across their boundaries, a standard FEA-like procedure does not work for multipatch IGA; an additional numerical technique is required for coupling NURBS patches. Therefore, in this paper, one of the discontinuous Galerkin (DG) methods, namely symmetric interior penalty Galerkin formulation, is employed to allow for discontinuities. For numerical calculations, a composite laminate with stacking sequences $$0^{\circ }{/}90^{\circ }$$ and $$0^{\circ }{/}90^{\circ }{/}0^{\circ }$$, respectively, is adopted. The stresses are calculated along the thickness of the composite laminate, subjected to a sinusoidal load, and they are compared with the analytical solutions. It is shown that DG–IGA gives a better approximation in comparison with the standard IGA.