A fixed‐grid b‐spline finite element technique for fluid–structure interaction

We present a fixed‐grid finite element technique for fluid–structure interaction problems involving incompressible viscous flows and thin structures. The flow equations are discretised with isoparametric b‐spline basis functions defined on a logically Cartesian grid. In addition, the previously proposed subdivision‐stabilisation technique is used to ensure inf–sup stability. The beam equations are discretised with b‐splines and the shell equations with subdivision basis functions, both leading to a rotation‐free formulation. The interface conditions between the fluid and the structure are enforced with the Nitsche technique. The resulting coupled system of equations is solved with a Dirichlet–Robin partitioning scheme, and the fluid equations are solved with a pressure–correction method. Auxiliary techniques employed for improving numerical robustness include the level‐set based implicit representation of the structure interface on the fluid grid, a cut‐cell integration algorithm based on marching tetrahedra and the conservative data transfer between the fluid and structure discretisations. A number of verification and validation examples, primarily motivated by animal locomotion in air or water, demonstrate the robustness and efficiency of our approach. Copyright © 2013 John Wiley & Sons, Ltd.     

Nonlinear free vibrations of porous composite microplates incorporating various microstructural-dependent strain gradient tensors

The main objective of the present numerical analysis is to predict the nonlinear frequency ratios associated with the nonlinear free vibration response of porous composite plates at microscale in the presence of different microstructural gradient tensors. To achieve this end, by taking cubic-type elements into account, isogeometric models of porous composite microplates are obtained with and without a central cutout and relevant to various porosity patterns of distribution along the plate thickness. The established unconventional models have the capability to capture the effects of various unconventional gradient tensors continuity on the basis of a refined shear deformable plate formulation. For the simply supported microsized uniform porous functionally graded material (U-PFGM) plate having the oscillation amplitude equal to the plate thickness, it is revealed that the rotation gradient tensor causes to reduce the frequency ratio about 0.73%, the dilatation gradient tensor causes to reduce it about 1.93%, and the deviatoric stretch gradient tensor leads to a decrease of it about 5.19%. On the other hand, for the clamped microsized U-PFGM plate having the oscillation amplitude equal to the plate thickness, these percentages are equal to 0.62%, 1.64%, and 4.40%, respectively. Accordingly, it is found that by changing the boundary conditions from clamped to simply supported, the effect of microsize on the reduction of frequency ratio decreases a bit.

     

等几何壳体分析与形状优化

首先基于Reissner-Mindlin理论进行了三维壳体等几何分析,而后基于此对三维壳体进行形状优化,提出了形状优化中灵敏度的全解析计算方法,包括位移应变阵、雅克比阵和刚度阵等相对控制顶点位置的灵敏度解析计算公式;通过实例验证了壳体等几何分析和灵敏度全解析计算方法的有效性。与传统的基于网格的灵敏度半解析计算方法相比,基于NURBS的灵敏度全解析计算具有精确、计算效率高的特点,且可以避免优化迭代中的网格畸变。     

面内功能梯度三角形板等几何面内振动分析

基于平面应变理论,利用等几何有限元方法分析了弹性边界条件下面内功能梯度三角形板的面内振动特性.板的材料属性沿厚度方向呈均匀分布,而在面内方向呈任意指数梯度变化.采用非均匀有理B样条(NURBS)基函数对三角形结构进行等几何建模和位移描述,实现了三角形板几何设计和振动分析的无缝衔接.在三角形板边界上引入虚拟弹簧约束并通过调节虚拟弹簧刚度,实现任意边界条件的施加.通过不同的单元细化方案和对比算例,验证了等几何方法的灵活性、准确性和快速收敛性.系统研究了边界条件、材料属性和几何参数对三角形板振动特性的影响.同时给出了弹性边界条件下面内功能梯度三角形板的振动特性解,具有重要参考价值.     

摩擦接触问题的比例边界等几何B可微方程组方法

摩擦接触问题是计算力学领域最具挑战性的问题之一,接触系统的泛函具有非线性、非光滑的特点,导致接触算法的收敛性与精确性难以保证.因此将比例边界等几何分析(scaled boundary isogeometric analysis,SBIGA)与B可微方程组(B dierential equation,BDE)相结合,提出了求解二维摩擦接触问题的比例边界等几何B可微方程组方法.在比例边界等几何坐标变换的基础上,通过虚功原理推导了关于边界控制点变量的接触平衡方程,表示成B可微方程组形式的接触条件可被严格满足,求解B可微方程组的算法的收敛性有理论保证.此比例边界等几何B可微方程组方法(SBIGA-BDE)只需在接触体边界进行等几何离散,使问题降低一维,能精确描述接触边界,并可通过节点插入算法进行真实接触区域的识别.此外,由于几何建模和数值分析使用相同的基函数,节约了划分网格的时间.以赫兹接触问题和悬臂梁摩擦接触问题为例,通过与解析解及数值计算软件ANSYS计算结果进行对比,验证了该方法求解二维摩擦接触问题的有效性及高精度等特点.      

Navier-Stokes方程的最小二乘等几何方法

基于Hermite多项式的C¹型单元构造复杂,限制了最小二乘有限元法的应用.引入高阶光滑的非均匀有理B样条作为基函数简化C¹型单元构造,提出求解黏性不可压流动Navier-Stokes方程的最小二乘等几何方法.用Newton法或Picard法对Navier-Stokes方程线性化,用线性化偏微分方程的余量定义最小二乘泛函,导出最小二乘变分方程,用NURBS构造高阶光滑的有限维空间来近似速度场和压力场.计算表明:本文方法计算的二维顶盖驱动流数值解能准确描述流动状况,计算的二维通道内圆柱绕流全局质量损失由最小二乘有限元法的6%降为0.018%,该方法可用于Navier-Stokes方程的求解,并且具有较好的质量守恒性.     

Isogeometric analysis for the dynamic problem of curved structures including warping effects

Curved beam structures have been used in many civil, mechanical, aircraft, and aerospace constructions. The analysis is mainly based on solid and plate models due to the fact that traditional curved beam elements do not include nonuniform warping effects, especially in the dynamic analysis. In this article, independent warping parameters have been taken into account and the initial curvature effect is considered. Curved beam’s behavior becomes more complex, even for dead loading, due to the coupling between axial force, bending moments, and torque that curvature produces. In addition to these, the Isogeometric tools (b-splines or NURBS), either integrated in the Finite Element Method or in a Boundary Element–based Method called Analog Equation Method, have been employed in this contribution for the dynamic analysis of horizontally curved beams of open or closed (box-shaped) cross sections. Free vibration characteristics and responses of the stress resultants and displacements to moving loading have been studied.     

Improved quadratic isogeometric element simulation of one-dimensional elastic wave propagation with central difference method

Two improved isogeometric quadratic elements and the central difference scheme are used to formulate the solution procedures of transient wave propagation problems. In the proposed procedures, the lumped matrices corresponding to the isogeometric elements are obtained. The stability conditions of the solution procedures are also acquired. The dispersion analysis is conducted to obtain the optimal Courant-Friedrichs-Lewy (CFL) number or time-step sizes corresponding to the spatial isogeometric elements. The dispersion analysis shows that the isogeometric quadratic element of the fourth-order dispersion error (called the isogeometric analysis (IGA)-f quadratic element) provides far more desirable numerical dissipation/dispersion than the element of the second-order dispersion error (called the IGA-s quadratic element) when appropriate time-step sizes are selected. The numerical simulations of one-dimensional (1D) transient wave propagation problems demonstrate the effectiveness of the proposed solution procedures.