基于GPU并行的锥体导管架平台结构冰激振动DEM-FEM耦合分析

该文为分析海冰与锥体海洋平台的相互作用,采用离散元（DEM）-有限元（FEM）耦合方法建立冰激海洋平台结构的耦合模型。通过具有粘结-破碎性能的球体离散单元对海冰的漂移及破碎现象进行计算,海洋平台锥体部分采用平板型壳单元构造,其整体构架及锥体内部的加劲肋采用梁单元构造,即建立壳单元与梁单元组合的锥体海洋平台有限元模型。为提高DEM-FEM耦合算法的计算规模和效率,发展了离散单元与平板型壳单元接触算法及GPU并行环境下参数传递算法。基于此耦合模型分别讨论了平台结构的冰载荷、冰激振动以及锥体应力分布,并与相关实测数据进行对比,为寒区锥体海洋平台的结构设计提供有益的参考。     

复合材料板壳结构的杂交/混合有限元分析

本文建立了一个9节点Lagrange退化壳杂交/混合有限元模型,以用于复合材料板壳结构的有限元分析。该有限元模型基于修正的Hellinger-Reissner原理,位移和应力均采用分离法思想,使得这一单元不仅具有下列优点:1.位移和应力计算精度都比较高;2.消除了多余零能变形模式;3.具有厚薄通用性;4.具有几何不变性,并且较之一般杂交/混合单元计算工作量进一步降低。单元还考虑了横向剪切影响。计算实例表明,本单元关于复合材料浅壳和深壳的解都与参考解吻合很好,且收敛很快。      

改进的有限元混合法用于板壳分析

本文从广义变分原理出发,给出了考虑横剪切变形板弯曲混合变分原理的具体形式,在Herrmann 三角形混合单元的基础上,导出了一个挠度按二次多项式、内力矩按线性分布的12个参数的三角形混合板元(×12单元)以及组合平面应力单元的三角形壳元。通过典型结构的分析计算,结果表明×12单元收敛速度快,计算精度高,位移和内力收敛速度比较协调一致,计算剪力方便。     

扁壳广义协调曲面矩形元

本文从修正的扁壳胡海昌-鹫律原理泛函出发,引入两方面的广义协调条件(单元边界位移的积分型协调条件,膜应变与位移之间的积分型协调条件),使泛函退化为扁壳势能原理泛函,在此基础上导出一个具有二十个自由度的扁壳曲面矩形元。此单元对厚扁壳和薄扁壳都通用,不出现剪切闭锁和薄膜闭锁现象,具有良好的性能。     

Glass/Epoxy复合材料叠层板面内剪切连续损伤模型

为确定S2玻璃纤维/环氧树脂(S2-Glass/Epoxy) 叠层复合材料面内剪切应力-应变关系,对S2-Glass/Epoxy 叠层复合材料面内剪切拉伸载荷下的弹、塑性连续损伤本构模型及应用进行了研究。基于平面应力状态下的连续损伤力学模型,通过典型面内剪切拉伸实验,分别建立了忽略塑性应变和考虑塑性应变的两种连续损伤力学(CDM)模型,并确定相关参数。通过ABAQUS/Explicit 用户子程序VUMAT接口,分别采用两种CDM模型对S2-Glass/Epoxy 叠层复合材料面内剪切拉伸实验进行有限元数值计算,与实验结果对比,验证模型可靠性,并分析单元类型对有限元计算结果的影响。研究结果表明: 忽略塑性应变的CDM模型可以很好地预测复合材料面内剪切失效强度,但不能较好地预测其非线性力学响应; 考虑塑性应变,将塑性硬化与损伤耦合后的CDM模型则能较好的预测复合材料非线性力学响应和面内剪切失效强度; 该平面应力状态下建立的CDM模型可用于壳单元进行复合材料有限元数值计算,横向剪切作用导致传统壳单元数值计算的载荷位移曲线略低于平面应力单元计算结果; 减缩积分算法有利于提高有限元数值计算结果的准确性。     

复合材料大变形任意加筋壳单元

构造了用于复合材料偏心加筋壳形结构大变形分析的任意加筋壳单元。在此模型中,肋骨连同壳的整体被视为一个单元偏心加筋壳单元。肋骨可放在壳单元内的任意位置和任意方向。所构造单元的特点是在网格划分时,可不必考虑肋骨的位置,这就给网格划分带来了很大的灵活性。在壳和肋骨的方程中,引用Von-Karman大变形理论计及几何非线性的影响,按照Mindlin-Reissner一阶剪切变形理论考虑横向剪切变形。     

几何非线性复合材料层合固体壳单元

通过定义广义应力,提出了一个改进的刚度矩阵,以克服固体壳元的厚度自锁问题,并能保证沿复合材料层合结构厚度方向上的连续应力分布;将应力插值函数分为低阶和高阶两部分,建议了一个新的非线性变分泛函,推导了一个用于几何非线性分析的九节点固体壳单元,该单元的计算精度和效率基本上与九节点减缩积分单元相当,与同类型其他单元相比,该单元显著提高了计算效率。     

实体退化板单元及其在板的振动分析中的应用

经典板壳单元是由板壳理论构造出来的,而经典的板壳理论是在空间弹性理论的基础上考虑板壳的基本假定得来的。在空间等参数单元的基础上,直接引入板壳的基本假定,修改空间等参数单元的弹性矩阵,从而构造出适合于厚薄板壳分析的20结点实体退化板单元,并将其应用于开口圆柱薄壳的静力分析和厚薄板的固有振动分析。数值算例表明,该单元收敛快,稳定性好,具有较高的精度。此外,该单元还可以用于曲边变厚度板、壳体及层合板的振动分析。     

功能梯度压电压磁球对称问题的静力响应

研究径向载荷作用下功能梯度压电压磁空心球壳的空间球对称电磁弹耦合静力学问题。假设压电压磁空心球壳的材料参数沿球厚度方向呈幂函数分布, 在球坐标系下, 由材料的参数方程、本构方程、几何方程和平衡方程导出在外激励作用下空心球壳体的应力、电势、磁势等物理量的解析解。结构内表面材料设为常用BaTiO₃-CoFeO₄复合材料, 分别对球壳内表面受力以及内外表面存在电势差或磁势差的情况进行数值讨论, 分别给出不同梯度参数下结构内部径向应力、环向应力、电势和磁势的分布。结果表明梯度参数的选取对功能梯度压电压磁球壳的性能有很大的影响。     

具有无序尺寸肋骨配置的加筋圆柱壳振动局域化

为控制周期加筋圆柱壳在通带振动的远距离传播问题,拟采取肋骨尺寸无序配置的方式实现将振动局限于振源附近,为此开展了具有肋骨尺寸无序配置的加筋圆柱壳振动局域化研究。为定量预报局域化因子,将加筋圆柱壳各周向模式振动沿轴向的传播等效为耦合振子链的振动传播问题。以加筋圆柱壳有限元振动分析的结果为输入,使用结合波数分析的参数辨识技术给出了等效振子固有频率参数和耦合参数,利用耦合振子链的无序局域化因子公式预报肋骨尺寸无序配置加筋圆柱壳的局域化因子,用于判断振动局域化的发生。为预报局域化因子所需要的肋骨尺寸无序度参数使用了参考模型辨识技术。针对具有无序尺寸肋骨配置的加筋圆柱壳所开展的振动局域化分析结果表明,肋骨尺寸变化能导致单元对地刚度发生改变,从而导致单元固有频率发生改变。因此,肋骨尺寸无序配置可实现振动局域化,具体程度同振动模式、肋骨尺寸无序度参数和频率相关。     

An optimal versatile partial hybrid stress solid‐shell element for the analysis of multilayer composites

In the present contribution we propose an optimal low‐order versatile partial hybrid stress solid‐shell element that can be readily employed for a wide range of geometrically linear elastic structural analyses, that is, from shell‐like isotropic structures to multilayer anisotropic composites. This solid‐shell element has eight nodes with only displacement degrees of freedom and only a few internal parameters that provide the locking‐free behavior and accurate interlaminar shear stress resolution through the element thickness. These elements can be stacked on top of each other to model multilayer composite structures, fulfilling the interlaminar shear stress continuity at the interlayer surfaces and zero traction conditions on the top and bottom surfaces of composite laminates. The element formulation is based on the modified form of the well‐known Fraeijs de Veubeke–Hu–Washizu multifield variational principle with enhanced assumed strains formulation and assumed natural strains formulation to alleviate the different types of locking phenomena in solid‐shell elements. The distinct feature of the present formulation is its ability to accurately calculate the interlaminar shear stress field in multilayer structures, which is achieved by the introduction of the assumed interlaminar shear stress field in a standard enhanced assumed strains formulation based on the Fraeijs de Veubeke–Hu–Washizu principle. The numerical testing of the present formulation, employing a variety of popular numerical benchmark examples related to element patch test, convergence, mesh distortion, shell and laminated composite analyses, proves its accuracy for a wide range of structural analyses.Copyright © 2012 John Wiley & Sons, Ltd.     

A mixed solid‐shell element for the analysis of laminated composites

This paper presents a versatile low order locking‐free mixed solid‐shell element that can be readily employed for a wide range of linear elastic structural analyses, that is, from thick isotropic structures to multilayer anisotropic composites. This solid‐shell element has eight nodes with only displacement degrees of freedom and few assumed stress parameters that provide very accurate interlaminar stress calculations through the element thickness. These elements can be stacked on top of each other to model multilayer structures, fulfilling the interlaminar stress continuity at the interlayer surfaces and zero traction conditions on the top and bottom surfaces of the laminate. The element formulation is based on the well‐known Fraeijs de Veubeke–Hu–Washizu mixed variational principle with enhanced assumed strains formulation and assumed natural strains formulation to alleviate the different types of locking phenomena in solid‐shell elements. The distinct feature of the present formulation is its ability to accurately calculate the interlaminar stress field in multilayer structures, which is achieved by the introduction of a constraint equation on the interlaminar stresses in the Fraeijs de Veubeke–Hu–Washizu principle‐based enhanced assumed strains formulation. The intelligent computer coding of the present formulation makes the present element appropriate for a wide range of structural analyses. To assess the present formulation's accuracy, a variety of popular numerical benchmark examples related to element convergence, mesh distortion, and shell and laminated composite analyses are investigated and the results are compared with those available in the literature. These benchmark examples reveal that the proposed formulation provides very good results for the structural analysis of shells and multilayer composites. Copyright © 2011 John Wiley & Sons, Ltd.     

Two‐ and three‐dimensional transition element families for adaptive refinement analysis of elasticity problems

In this paper, new enhanced assumed strain (EAS) and hybrid stress transition element families are developed for 2D and 3D adaptive refinement analysis of elasticity problems. The EAS element families are based on some existing incompatible transition element families. By using the EAS method and the previous incompatible modes, the B ‐matrix columns associated with the EAS modes can be directly designed such that their domain integrals vanish automatically and they can be computed more efficiently. For 2D hybrid stress transition element families, it is possible to derive different stress fields that lead to rank‐sufficient transition elements. However, the task becomes intractable for 3D hybrid stress transition elements in which many combinations of mid‐side and mid‐face nodes are possible. This paper proposes to use hybrid stress transition element families in which the assumed stress fields are linearly complete. The new 2D element family is more accurate than the 2D rank‐sufficient element family. The new 3D element family is more accurate than the one with additional bilinear stress modes. Numerical examples reveal that the most accurate transition element families are the newly developed hybrid stress families followed by the EAS families, the incompatible families and then the compatible families. Copyright © 2008 John Wiley & Sons, Ltd.     

Stress resultant geometrically non-linear shell theory with drilling rotations. Part III: Linearized kinematics

A consistent formulation of the geometrically linear shell theory with drilling rotations is obtained by the consistent linearization of the geometrically non-linear shell theory considered in Parts I and II of this work. It was also shown that the same formulation can be recovered by linearizing the governing variational principle for the three-dimensional geometrically non-linear continuum with independent rotation field. In the finite element implementation of the presented shell theory, relying on the modified method of incompatible modes, we were able to construct a four-node shell element which delivers a very high-level performance. In order to simplify finite element implementation, a shallow reference configuration is assumed over each shell finite element. This approach does not impair the element performance for the present four-node element. The results obtained herein match those obtained with the state-of-the-art implementations based on the classical shell theory, over the complete set of standard benchmark problems.     

Mixed shell element for seven‐parameter formulation

A new mixed shell element is developed for a seven‐parameter formulation in this paper. The mixed shell element is constructed by assuming stress field and displacement field together. Assumed stress field and assumed displacement field can be combined by stress–strain relationship with Hu‐Washizu functional. The developed mixed shell element can provide more flexible stiffness than other commercial softwares. Additionally, seven‐parameter shell formulation is used instead of Reissner/Mindlin formulation, since it can provide the thickness change. Even though some commercial engineering software are not proper for very thick shell structure, the developed mixed shell element for seven‐parameter formulation can be used without distinction of thick shell and thin shell. An example of shell models with different thickness is provided with solid model. Static and modal analyses are also performed for verification. Copyright © 2005 John Wiley & Sons, Ltd.     

带有沙漏控制的相对自由度壳元

针对一点积分的八节点相对自由度壳单元存在的沙漏现象，提出采用拟应变法解决该问题的方法，并对锁死问题进行研究。给出了带有沙漏控制的八节点相对自由度壳元内的坐标、位移插值公式，推导了拟应变的表达式，通过Hu-Washizu变分原理，建立了有限元求解方程。利用Wilson非协调位移模式，单元的计算精度得到了明显改善。算例表明：基于八节点相对自由度壳单元，本文给出的沙漏控制算法能够有效的解决线性静力问题，并且具有较高的计算精度。     

An eight‐node hybrid‐stress solid‐shell element for geometric non‐linear analysis of elastic shells

This paper presents eight‐node solid‐shell elements for geometric non‐linear analysis of elastic shells. To subdue shear, trapezoidal and thickness locking, the assumed natural strain method and an ad hoc modified generalized laminate stiffness matrix are employed. A selectively reduced integrated element is formulated with its membrane and bending shear strain components taken to be constant and equal to the ones evaluated at the element centroid. With the generalized stresses arising from the modified generalized laminate stiffness matrix assumed to be independent from the ones obtained from the displacement, an extended Hellinger–Reissner functional can be derived. By choosing the assumed generalized stresses similar to the assumed stresses of a previous solid element, a hybrid‐stress solid‐shell element is formulated. Commonly employed geometric non‐linear homogeneous and laminated shell problems are attempted and our results are close to those of other state‐of‐the‐art elements. Moreover, the hybrid‐stress element converges more readily than the selectively reduced integrated element in all benchmark problems. Copyright © 2002 John Wiley & Sons, Ltd.     

内参型附加非协调位移基本项的推导和应用

在协调元位移模式基础上附加内参项是构造非协调元的一种常用方法。目前一般是先假设非协调位移模式(不能保证其通过小片试验)，然后按照一定的方法进行修改，从而形成能够保证收敛的非协调位移场，可是构造过程往往较复杂。本文从广义协调条件出发，首次推导了平面问题内参任意阶次附加非协调位移基本项通用公式，形式简单，便于工程人员直接应用于工程实践。根据通用公式，本文以Q8协调元为基础，发展了一个新的非协调元，数值试验表明它能够保证收敛，有较高精度，抗畸变能力强，从而证明了本文方法的可行性。     

A quadrilateral hybrid stress element for Mindlin plates based on incompatible displacements

An hybrid stress element formulation based on internal, incompatible displacements is used to develop efficient Mindlin plate elements. The 4-node quadrilateral Mindlin plate element is derived from a modified energy functional. Both displacements and stresses are defined in the natural co-ordinate interpolation system. The assumed stress field is obtained by tensor transformation and so chosen as to ensure that the element is co-ordinate invariant and stable. Shear locking is avoided through an appropriate identification of the internal, incompatible displacement field. The role played by incompatible displacements in the formulation of hybrid stress elements for thin and moderately thick plates is discussed. Numerical applications are presented to illustrate the accuracy and reliability of the suggested Mindlin plate element.     

A novel versatile multilayer hybrid stress solid-shell element

This paper presents a versatile multilayer locking free hybrid stress solid-shell element that can be readily employed for a wide range of geometrically linear elastic structural analyses, i.e. from shell-like isotropic structures to multilayer anisotropic composites. This solid-shell element has eight nodes with only displacement degrees of freedom and a few internal parameters that provide the locking free behavior and accurate interlaminar stress resolution through the element thickness. These elements can be stacked on top of each other to model multilayer structures, fulfilling the interlaminar stress continuity at the interlayer surfaces and zero traction conditions on the top and bottom surfaces of composite laminates. The element formulation is based on the modified form of the well-known Fraeijs de Veubeke–Hu–Washizu (FHW) multifield variational principle with enhanced assumed strains (EAS formulation) and assumed natural strains (ANS formulation) to alleviate the different types of locking phenomena in solid-shell elements. The distinct feature of the present formulation is its ability to accurately calculate the interlaminar stress field in multilayer structures, which is achieved by incorporating an assumed stress field in a standard EAS formulation based on the FHW principle. To assess the present formulation’s accuracy, a variety of popular numerical benchmark examples related to element patch tests, convergence, mesh distortion, shell and laminated composite analyses are investigated and the results are compared with those available in the literature. This assessment reveals that the proposed solid-shell formulation provides very accurate results for a wide range of structural analyses.