空间钢框架结构的改进双重非线性分析

为了探讨三维结构的高等分析方法,给出基于UL法的严格三维单元虚功增量方程,详细推导了考虑剪切变形影响的三维梁柱单元的几何非线性切线刚度矩阵和基于三维单元简化塑性区模型的双重非线性刚度方程,并编制空间钢框架结构的双重非线性分析程序。使用包括几何与材料双重非线性的数值算例来验证方法和计算机程序的可行性、精确度与有效性。利用该程序,只需一个单元/构件即可准确预测空间钢框架结构的极限承载力与失稳模态,提高了结构的非线性分析效率。     

橡胶隔振器高频动态特性的计算方法

该文采用一种由Mooney-Rivlin模型和多个Maxwell模型叠加组成的非线性粘弹性本构模型,用于计算橡胶隔振器的高频动态特性。该文给出了在时域和频域范围内拟合本构模型中粘弹性参数的方法,利用拟合得到的本构模型参数,对某款橡胶悬置跨点动态特性进行计算,并与实验结果进行对比。该文还建立了橡胶隔振器等效力学模型,分析了原点动刚度和跨点动刚度的区别,分析表明:使用跨点动态特性测试法可消除测试中附加惯性力的影响,适用于橡胶隔振器高频动特性的测试;同时,该文搭建了橡胶隔振器有限元模型,分别用于分析其跨点动刚度与原点动刚度,并将分析结果与实验结果进行对比,分析结果验证了有限元模型和力学模型的正确性。除此之外,该文还分析对比了时域（松弛、蠕变）和频域（简谐动态试验）拟合粘弹性参数方法的优缺点。     

改进共旋坐标法的Timoshenko梁单元非线性分析

为提高空间Timoshenko梁单元非线性问题的计算精度,在共旋坐标法的基础上,提出了一种改进的Timoshenko梁单元几何非线性分析方法。利用虚功原理得到改进空间梁单元的刚度矩阵;使用有限质点法中的逆向运动思路计算单元局部坐标系下的刚体旋转矩阵;根据整体坐标系与局部坐标系之间旋转角度的转化以及微分关系,求得空间梁单元的切线刚度矩阵;编制了相应的有限元程序,对多个经典的大变形结构进行几何非线性分析。计算结果印证了该文所提出改进方法的正确性,同时与传统共旋坐标法相比,具有更高的精度。     

基于隔离非线性的实体单元模型与计算效率分析

实体有限元模型计算中往往需要较多的计算单元与结点数量,且这些单元状态判定以及大规模的刚度矩阵分解将消耗大量的计算资源,计算效率低。该文基于隔离非线性法理论建立了线性四面体与六面体等参单元分析模型,采用直接积分格式的6积分点替代六面体等参单元的8高斯点作为非线性应变插值点,能够在保证计算精度的同时提高单元状态判定效率。控制方程采用Woodbury公式与组合近似法联合求解,使得整个求解过程只有矩阵回代以及矩阵与向量的乘积,进一步提高了求解效率。基于时间复杂度的计算效率分析表明:随着结点自由度数目的增加,该文方法的计算效率相对传统变刚度法显著提高,数值算例验证了实体单元模型的正确性以及算法的高效性。     

基于PARAFAC分析和SVM离心泵故障诊断方法

在多尺度平行因子分析理论的基础上,将原始信号经过多尺度小波分解得到三维时频信号,再经平行因子分析得到通道加载因子、时间加载因子和频率加载因子,通过实验分析,后二者可以明显地表征设备正常或故障状态,利用这一特征建立不同状态的离心泵与其对应的时间加载因子和频率加载因子的映射关系,并以此作为改进粒子群算法优化后的支持向量机分类器的特征向量进行故障分类。与小波包能量特征相比,所提的这种诊断方法用于离心泵故障诊断时提取特征更为简便,所提分类器的分类准确率有显著提高,而其复杂度却没有明显增加。     

充气薄膜褶皱分析的高效互补有限元列式

褶皱变形是柔性薄膜结构的一种常见的失稳模式,其数值模拟具有挑战性。基于连续体和张力场理论,提出了一种适用于充气薄膜结构褶皱分析的互补共旋有限元方法。采用共旋坐标法,将物体的大变形分解为结构整体坐标系下的刚体运动和单元局部坐标系下小应变变形,推导了一个空间三节点三角形膜单元的切线刚度矩阵。该刚度矩阵包含材料刚度、旋转刚度和平衡投影刚度矩阵三个部分,涵盖了随动载荷对单元刚度的影响。在单元局部坐标系下,依据双模量材料本构关系构造了一个褶皱模型,能够判断单元处于“张紧”“褶皱”或“松弛”状态。进一步通过建立等价的线性互补问题,消除了迭代求解过程中的内力振荡,改善了算法的稳定性。数值算例表明:该文方法能够准确地预测充气薄膜结构的位移、应力以及褶皱区域。较之已有的“拟动态”和“惩罚”方法,该方法不需要引入额外的求解技术来保证收敛,具有良好的稳定性。     

An enhanced layerwise finite element using a robust solid-shell formulation approach

The application of layerwise theories to correctly model the displacement field of sandwich structures or laminates with high modulus ratios, usually employs plate or facet shell finite element formulations to compute the element stiffness and mass matrices for each layer. In this work, a different approach is proposed, using a high performance hexahedral finite element to represent the individual layer mass and stiffness. This 8-node hexahedral finite element is formulated based on the application of the enhanced assumed strain method (EAS) to resolve several locking pathologies coming from the high aspect ratios of the finite element and the usual incompressibility condition of the core materials. The solid-shell finite element formulation is introduced in the layerwise theory through the definition of a projection operator, which is based on the finite element variables transformation matrix. The new finite element is tested and the implemented numerical remedies are verified. The results for a soft core sandwich plate are hereby presented to demonstrate the proposed finite element applicability and robustness.     

基于等几何分析的参数化曲梁结构非线性动力学降阶模型研究

模型降阶方法通过构造全阶模型的低阶近似模型有效提升了求解效率，同时也保留了原阶模型的主要信息从而保证了较高的计算精度。对于结构非线性以及参数化的模型降阶问题，常需要重复计算刚度矩阵等非线性以及参数依赖项，求解效率较低。此外，当参数化模型的几何形状改变时，往往需要重复进行CAD与有限元(FEA)模型的转换，这对于复杂结构较为耗时。等几何分析采用描述几何形状的非均匀有理B样条(NURBS)插值物理场，实现了CAD与FEA模型的统一，消除了两者之间繁琐的模型转换过程，其具有几何精确、高阶连续等优点，并且几何形状在细化过程中保持不变，非常适合于薄壁类结构的分析以及参数化表达。该研究结合等几何分析、特征正交分解(POD)以及离散经验插值方法(DEIM)研究参数化的平面曲梁结构的非线性动力学模型降阶问题。数值结果表明，基于等几何分析的POD-DEIM降阶模型能够显著提升平面曲梁结构的非线性动力学计算效率，并且该模型对于参数化以及变载荷等情形显示出了良好的适应性。     

A simplified co‐rotational method for quadrilateral shell elements in geometrically nonlinear analysis

This paper presents a simplified co‐rotational formulation for quadrilateral shell elements inheriting the merit of element‐independence from the traditional co‐rotational approach in literature. With the objective of application to nonlinear analysis of civil engineering structures, the authors further simplify the formulation of the geometrical stiffness using the small strain assumption, which is valid in the co‐rotational approach, with the warping effects considered as eccentricities. Compared with the traditional element‐independent co‐rotational method, the projector is neglected both in the tangent stiffness matrix and in the internal force vector for simplicity in formulation. Meanwhile, a quadrilateral flat shell element allowing for drilling rotations is adopted and incorporated into this simplified co‐rotational algorithm for geometrically nonlinear analysis involved with large displacements and large rotations. Several benchmark problems are presented to confirm the efficiency and accuracy of the proposed method for practical applications.     

Delamination identification and response prediction in composites using bivariate Gamma function-based microgenetic algorithms

This study deals with an identification of stiffness reduction and response predictions occurred by the delamination damage in laminated composite plates under impact loads. Combined bivariate Gamma function and microgenetic algorithms are developed to determine the crack region due to the delamination and predict future responses. The validity of the proposed method was verified using impact-induced data obtained from a two-dimensional delamination finite element model. Examples indicated that the proposed approach is a feasible and advantageous method through which future dynamic responses can be predicted and the distribution of the degraded stiffness of laminated composite structures can be inspected for different measuring locations and fibre angles.     

复合材料副簧刚度的匹配设计方法

为了对复合材料副簧的刚度进行匹配设计,设计了包含有复合材料副簧的主副簧总成结构。采用集中载荷法计算复合材料副簧的等效载荷,根据原钢板弹簧的挠度变化来估算复合材料副簧的等效刚度。根据钢板弹簧设计理论,对复合材料副簧的等效刚度进行匹配设计。采用ABAQUS软件对设计的主副簧总成的总刚度进行有限元模拟,通过调整复合材料副簧的铺层数量来修正复合材料副簧的等效刚度。提出的匹配设计方法对复合材料板簧的推广应用具有重要意义。      

Coarse/fine mesh preconditioners for the iterative solution of finite element problems

A class of preconditioners built around a coarse/fine mesh framework is presented. The proposed method involves the reconstruction of the stiffness equations using a coarse/fine mesh idealization with relative degrees-of-freedom derived from the element shape functions. This approach leads naturally to effective preconditioners for iterative solvers which only require a factorization involving coarse mesh variables. A further extension is the application of the proposed method to super-elements in conjunction with substructuring (domain decomposition) techniques. The derivation of the coarse/fine mesh discretization via the use of transformation matrices, allows a straightforward implementation of the proposed techniques (as well as multigrid type procedures) within an existing finite element system.     

A new finite element scheme for instability analysis of thin shells

This paper describes a new finite element scheme for the analysis of instability phenomena of arbitrary thin shells. A computationally efficient procedure is proposed for calculating the non-linear stiffness and tangential stiffness matrices for a doubly-curved quadrilateral element defined by co-ordinate lines. The essential feature is the explicit addition of the non-linear terms into the rigid-body motion of the element. Thus the non-linear and tangential element stiffness matrices can easily be generated by transforming the generalized element stiffness matrix for linear analysis, and the non-linear terms of these matrices are separated into a number of component terms multiplied by the rigid-body rotations. These component terms can be stored permanently and used to calculate efficiently the non-linear and tangential stiffness matrices at each iteration. Illustrative examples are presented which confirm the validity of the present approach in the analysis of instability phenomena of thin plates and shells.     

基于薄层单元的螺栓连接结构接触面不确定性参数识别

提出了一种螺栓连接接触面不确定性参数识别方法,首先采用薄层单元对接触面进行参数化,然后根据不确定性识别方法识别薄层单元参数。以四螺栓搭接结构试验模型为研究对象,开展接触面不确定性参数识别方法仿真研究。采用Monte-Carlo方法构造待识别参数真实值样本,代入基准有限元模型中计算获得具有统计意义的仿真试验数据;采用不确定性参数识别方法预测薄层单元参数均值与标准差,仿真结果表明:该方法能够较为准确的模拟接触面法向和切向接触刚度,并显著提高连接结构的建模效率,建立反映真实结构动态性能统计特征的有限元模型。     

Topology optimization of plate/shell structures with respect to eigenfrequencies using a biologically inspired algorithm

In this article, a novel approach is presented to perform topology optimization in a simple and explicit way. The method capitalizes on the use of a bio-inspired algorithm to represent topology, leading to more flexible optimization solutions along with explicit structure representation. To avoid remeshing upon design changes, a special treatment called the enhanced stiffness transformation approach (ESTA) is introduced to transform the stiffness and mass matrices of the growing stiffener into a set of equivalent stiffness and mass matrices. In this way, stiffeners are separated from the finite element mesh and can grow in an arbitrary direction to form an optimized layout solution. Notably, this approach incorporates more geometric information into topology optimization, which improves the clarity of stiffener layouts. Finally, the effectiveness of the proposed method is illustrated with two examples of maximum eigenfrequency design of plate/shell structures.     

A solid-shell layerwise finite element for non-linear geometric and material analysis

The application of layerwise theories to correctly model the displacement field of sandwich structures or laminates with high modulus ratios usually employs plate or facet-shell finite element formulations to compute the element stiffness and mass matrices for each layer. In this work an alternative approach is proposed, using a high performance hexahedral finite element to represent the individual layer mass and stiffness. This eight-node hexahedral finite element is formulated based on the application of the enhanced assumed strain method (EAS) to solve several locking pathologies coming from the high aspect ratio of the finite element and the usual incompressibility condition of the core materials. The solid-shell finite element formulation is introduced in the layerwise theory through the definition of a projection operator, based on the finite element variables transformation matrix. The non-linear geometric and material capabilities are introduced into the finite element formulation, allowing for the representation of large displacements, large deformation and material non-linear behaviors. The developed formulation is numerically tested and benchmarked, being validated by using published experimental results obtained from sandwich specimens.     

An analytical stiffness derivative extended finite element technique for extraction of crack tip Strain Energy Release Rates

An analytical approach combined with the extended finite element method (XFEM) is proposed to extract the Strain Energy Release Rates within the classical stiffness derivative technique. The proposed idea hinges on the following two XFEM properties: (i) the crack is mesh independent, i.e. there is no need for mesh perturbations in the vicinity of the crack and (ii) the asymptotic crack tip field is embedded in the mathematical formulation of the stiffness matrix. By employing these properties we show that the derivative of the stiffness matrix with respect to the crack extension can be computed in a closed form and on the fly during the analysis. Thus the virtual crack extension, and the error inherent in the finite difference scheme of the classical stiffness derivative technique is completely avoided. Numerical results on few benchmark problems show that this method is comparable to the J-integral method.     

计及二阶效应的梁杆系统动力响应分析方法研究

基于动力刚度法和有限元理论提出了一种考虑二阶效应计算梁杆动力响应的新方法。通过求解轴向力作用下Bernoulli-Euler 梁横向和轴向挠度自由振动微分方程,利用位移边界条件反解出待定系数,得到了动态精确形函数;使用经典有限元方法推导了考虑截面自身旋转惯量的质量阵和考虑二阶效应的刚度阵,该质量阵和刚度阵各元素均为轴力和圆频率的超越函数;建立了杆系结构瞬态动力学分析的动力平衡方程,给出了稳定和高效的求解方案。对几个典型的算例进行了计算分析,并与通用软件ANSYS 的计算结果进行了比较。计算结果表明:该分析梁杆系统动力响应的新方法具有较高的计算精度和效率,特别是能够准确地计入轴力对于梁杆动力响应的影响。     

Bidirectional evolutionary structural optimization for nonlinear structures under dynamic loads

This study aims to develop efficient numerical optimization methods for finding the optimal topology of nonlinear structures under dynamic loads. The numerical models are developed using the bidirectional evolutionary structural optimization method for stiffness maximization problems with mass constraints. The mathematical formulation of topology optimization approach is developed based on the element virtual strain energy as the design variable and minimization of compliance as the objective function. The suitability of the proposed method for topology optimization of nonlinear structures is demonstrated through a series of two- and three-dimensional benchmark designs. Several issues relating to the nonlinear structures subjected to dynamic loads such as material, geometric, and contact nonlinearities are addressed in the examples. It is shown that the proposed approach generates more reliable designs for nonlinear structures.