首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 93 毫秒
1.
基于哈密顿原理,得到水凝胶的化学-力学耦合控制方程的等效积分形式和有限元形式。在整体坐标系下推导出用形函数表示的化学-力学耦合应变矩阵和单元刚度矩阵,并且得到在局部坐标系下的离散化形式。结合ABAQUS软件,编制了用户单元子程序UEL,通过数值算例验证了所开发单元的正确性,为在ABAQUS软件中实现各种耦合问题的有限元UEL编程提供了参考依据。  相似文献   

2.
基于开放式结构有限元系统SiPESC.FEMS的单元计算模块的设计模式,研发设计一种通用的平板壳单元计算框架。考虑板壳单元的组合关系和程序编制过程中的重用性及灵活性等特点,采用了软件设计中的构造器(Builder)模式实现不同的组合单元。本框架具有很好的通用性和可扩展性,为有限元程序研发提供了一个新的方式;同时,系统能够处理复杂荷载和边界条件,并可灵活实现不同类型单元的组合分析。本文利用此方法构造五种平板壳单元,通过数值算例分析对比讨论其性能,为选取合适的平板壳单元类型进行结构数值分析提供参考。  相似文献   

3.
应用膜板比拟关系 ,可以避开 c1 连续性的困难 ,为板单元的构造提供了一种新的途径 ,并已成功地构造出一系列相应的板单元。本文构造了一个四节点二十四自由度的平板壳单元 ,该单元由平面四节点理性元 RQ4(膜部分 )和由膜板比拟理论构造的一个四节点十二自由度的板单元 (弯曲部分 )构成。该单元构造简单 ,数值结果表明具有很好的收敛性和精度。  相似文献   

4.
刘敏  李旭 《计算力学学报》2019,36(5):693-698
依据有限元理论,结合内聚力模型法则,推导出二维二次粘结界面单元在大位移情况下的数值格式,得到用形函数表示的单元位移模式、载荷向量和刚度矩阵,并进行了离散化。基于ABAQUS软件的自定义扩展模块,编制了相应的用户单元子程序UEL,通过数值算例验证了该程序的准确性和有效性。这一成果能为在ABAQUS软件中开展相关数值研究,以及开发其他类型的内聚力界面有限单元提供思路和参考。  相似文献   

5.
基于二维张量积区间B样条小波,构造了一种件能良好的小波平板壳单元.在小波单元的构造过程中,用二维区间B样条小波尺度函数取代传统多项式插值,在所构造的区间B样条平面弹性单元和平面Mindlin板单元的基础上组合而成.区间B样条小波单元同时具有B样条函数数值逼近精度高和多种用于结构分析的基函数的特点.数值算例表明:与传统有限元和解析解相比,构造的小波平板壳单元具有求解精度高,单元数量和自由度少等优点.  相似文献   

6.
本文从部分协调的三角形薄板弯曲单元出发,并假设横向剪切应变在单元内线性变化,提出了一种考虑横向剪切变形的具有15个自由度的三角平板弯曲单元。该单元应用于薄板和中等厚度板分析均有较高的精度,计算效率高,可用于工程中具有复杂形状的薄板和中等厚度板结构的分析。  相似文献   

7.
提出一类用于分析弹性板问题的算子自定义小波弹性板单元构造方法。该方法的优点在于根据工程问题的求解需要灵活构造具有解耦特性的算子自定义小波基,使得系统多尺度刚阵具有沿对角线的强稀疏性,从而实现了该算法在每个尺度上独立、快速求解,系统方程的求解效率得到较大提高。建立多分辨Lagrange有限元空间和多尺度计算理论,提出基于稳定完备法的算子自定义小波弹性板单元构造方法及解耦条件。依据两尺度相对误差估计,提出自适应算子自定义小波有限元算法。数值算例证明,算子自定义小波弹性板单元具有求解精度与计算效率高等特点。  相似文献   

8.
本文将板壳结构优化纳入了杆-膜-梁组合结构优化程序系统DDDU-2,从而发展成杆-膜-梁-壳组合结构优化的DDDU-3程序系统。该系统含有多种单元,可以处理组合结构的优化设计问题,具有良好的实用性。 在DDDU-3中,采用了基于离散Kirchhoff假定的三角形平板壳元。大量的数值试验和长期使用经验表明,这是一种有效的单元。在程序中,应力约束是以相当应力的形式给出  相似文献   

9.
板壳弹塑性屈曲的有限元分析   总被引:4,自引:0,他引:4  
1 拟协调双曲壳单元为解决板壳有限元分析的C~1连续性问题,本文给出三角形拟协调双曲扁壳单元的刚度矩阵.拟协调单元是基于域内假设应变场和边界拟协调位移场,应用最小势能原理所构  相似文献   

10.
虚拟裂纹闭合法在结构断裂分析中的应用   总被引:1,自引:0,他引:1  
基于三维虚拟裂纹技术(3DVCCT),利用ABAQUS用户单元子程序(UEL)编写裂纹界面单元,使3DVCCT集成于ABAQUS软件中,直接计算出裂纹的断裂参数.采用此方法对连杆杆身表面裂纹进行研究,得到了连杆裂纹的应力强度因子的分布规律.  相似文献   

11.
混凝土细观随机骨料结构与有限元网格剖分   总被引:23,自引:1,他引:22  
在细观层次上,混凝土被认为是一种由粗骨料、水泥砂浆及二者间的粘结带所组成三相非均质复合材料。本文首先基于蒙特卡罗随机抽样原理,用“取和放”方法在计算机上产生形状、尺寸和骨料颗粒分布与真实混凝土相似的随机骨料结构,再使用有限元分析软件ANSYS对骨料区域及砂浆区域分别划分网格,并编程在骨料和砂浆之间生成三角形三结点可控制厚度粘结单元,从而使三相网格缝合为一个整体,为混凝土非线性有限元分析提供可靠的细观计算模型。最后利用建立的模型进行混凝土轴心受拉和轴心受压的仿真模拟,在细观层次研究的基础上揭示出混凝土的宏观力学性能。  相似文献   

12.
This paper presents the user defined nine-node piezoelectric shell element and its implementation within a user element subroutine for modeling and analysis of thin-walled active composite structures. The element has been implemented within commercial finite element software ABAQUS. Application of the element covers modeling of arbitrary thin-walled structures also with complex geometries, whereby the automated mesh generation has been accomplished by developing a Python based interface for meshing procedure. In order to be able to perform the post-processing, a special adaptation of the user element had to be performed for visualization purposes. The implemented element regards the piezoelectric thin layers polarized in the thickness direction and it is based on the \(e_{31}\) piezoelectric effect. It has been also shown that this biquadratic nine-node element based on degenerated shell approach is less prone to locking effects and more suitable for implementation with curved structures. Through several examples the accuracy of the implemented user defined shell element as well as of the Python-based mesh extension has been demonstrated, along with the possibilities for post-processing. Meshing the structures with the nine-node user element is not third party software dependent.  相似文献   

13.
A simple and efficient triangular finite element for plate bending   总被引:1,自引:0,他引:1  
A simple and efficient triangular finite element is introduced for plate bending application. The element is a three-node triangular one with three basic degrees of freedom per node and two internal rotation degrees of freedom, using selective reduced integration. Numerical examples indicate that, despite its simplicity, the element is not only competitively accurate, but also useful as a thick/thin triangular plate bending element. It is also pointed out that this element using selective reduced integration is in fact a mixed element.  相似文献   

14.
基于重叠划分的自由网格四边形单元计算方法   总被引:1,自引:0,他引:1  
提出了一种基于重叠划分的自由网格四边形单元计算方法。这一方法将四边形单元引入到自由网格计算方法中,不仅提高了计算的精度,同时还保留了自由网格计算方法的特点。方法首先对分析域内自动生成的每一个节点建立一套临时三角形单元,利用这些临时三角形单元组合生成四边形单元,以节点为单位进行计算。由于各矩阵的计算与组集均以节点为中心进行处理,因而特别适合于并行计算环境。在详细介绍自由网格四边形单元计算方法的基础上,利用数值算例证实了这一方法改善计算精度方面的有效性。  相似文献   

15.
运用边光滑有限元法,研究分析了加筋板结构的静力和自由振动问题。在边光滑有限元法中,将基于边的应变光滑技术用于对原来的应变场进行光滑操作;由于应变光滑技术能够适当地软化原来过刚的有限元模型,从而能够得到更加接近于系统准确刚度的光滑有限元模型;鉴于三角形单元良好的适用性,选用三角形单元对模型进行网格划分;同时,为了解决低阶Reissner-Mindlin板单元弯曲过程中的横向剪切自锁问题,采用了一种新型的离散剪切间隙技术。算例的数值计算结果表明,与传统的有限元法相比,边光滑有限元法能够得到精度更高的计算结果,且收敛更快,计算效率更佳。  相似文献   

16.
The boundary-type finite element method has been investigated and applied to the Helmholz and mild-slope equations. Four types of interpolation function are examined based on trigonometric function series. Three-node triangular, four-node quadrilateral, six-node triangular and eight-node quadrilateral elements are tested; these are all non-conforming elements. Three types of numerical example show that the three-node triangular and four-node quadrilateral elements are useful for practical analysis.  相似文献   

17.
杜成斌  黄文仓  江守燕 《力学学报》2022,54(4):1026-1039
混凝土是一种被广泛应用于土木和水利工程中的准脆性材料, 在各种内外部因素的作用下, 开裂是混凝土结构最为普遍的破坏形式, 准确模拟结构的开裂过程, 对于结构的安全评估至关重要. 将比例边界有限元与非局部宏微观损伤模型相结合提出一种准脆性材料开裂模拟新方法. 以比例边界有限元子域的比例中心作为物质点, 通过两比例中心(物质点对)之间的物质键的正伸长率来定义微细观损伤, 将某点影响域内物质键的微细观损伤加权平均得到该点的宏观拓扑损伤. 再引入能量退化函数, 将宏观拓扑损伤嵌入到比例边界有限元的基本框架中. 充分利用比例边界有限元网格允许存在悬挂节点的优势, 采用四叉树网格离散技术进行快速、高质量的多级网格划分与过渡. 通过一个I型开裂与一个混合型开裂的两个典型算例, 验证了该方法可捕获结构裂纹扩展路径与荷载变形曲线. 与现有的方法相比, 本文的损伤模型可得到更准确的局部开裂损伤带, 结果更为合理, 且具有更高的计算精度和计算效率. 当损伤过程区网格尺寸小于影响域半径的1/5时, 计算结果不存在网格敏感性问题.   相似文献   

18.
In this work, a new plane stress element is proposed for the nonlinear static and dynamic analysis of plane stress/plane strain problems. The four node quadrilateral element formulation for the elastic case is extended by introducing a novel hysteretic constitutive relation, based on the Bouc–Wen model of hysteresis. The hysteretic model introduced is directly derived from the governing equations of classical plasticity based on the flow rule and specific hardening law. The stiffness matrix of the element is formulated using the principle of virtual displacements, where the elastic stress–strain relation is substituted by the hysteretic relation proposed. The derived stiffness matrix is expressed as a smooth function of the internal stress field both in the elastic and inelastic regime. The efficiency of the proposed element in the simulation of the cyclic behavior in plane structures is presented through illustrative examples.  相似文献   

19.
Finite elements with different orders can be used in the analysis of constrained deformable bodies that undergo large rigid body displacements. The constrained mode shapes resulting from the use of finite elements with different orders differ in the way the stiffness of the body bending and extension are defined. The constrained modes also depend on the selection of the boundary conditions. Using the same type of finite element, different sets of boundary conditions lead to different sets of constrained modes. In this investigation, the effect of the order of the element as well as the selection of the constrained mode shapes is examined numerically. To this end, the constant strain three node triangular element and the quadratic six node triangular element are used. The results obtained using the three node triangular element are compared with the higher order six node triangular element. The equations of motion for the three and six node triangular elements are formulated from assumed linear and quadratic displacement fields, respectively. Both assumed displacement fields can describe large rigid body translational and rotational displacements. Consequently, the dynamic formulation presented in this investigation can also be used in the large deformation analysis. Using the finite element displacement field, the mass, stiffness, and inertia invariants of the three and six-node triangular elements are formulated. Standard finite element assembly techniques are used to formulate the differential equations of motion for mechanical systems consisting of interconnected deformable bodies. Using a multibody four bar mechanism, numerical results of the different elements and their respective performance are presented. These results indicate that the three node triangular element does not perform well in bending modes of deformation.  相似文献   

20.
In this paper, new planar isoparametric triangular finite elements (FE) based on the absolute nodal coordinate formulation (ANCF) are developed. The proposed ANCF elements have six coordinates per node: two position coordinates that define the absolute position vector of the node and four gradient coordinates that define vectors tangent to coordinate lines (parameters) at the same node. To shed light on the importance of the element geometry and to facilitate the development of some of the new elements presented in this paper, two different parametric definitions of the gradient vectors are used. The first parametrization, called area parameterization, is based on coordinate lines along the sides of the element in the reference configuration, while the second parameterization, called Cartesian parameterization, employs coordinate lines defined along the axes of the structure (body) coordinate system. The fundamental differences between the ANCF parameterizations used in this investigation and the parametrizations used for conventional finite elements are highlighted. The Cartesian parameterization serves as a unique standard for the triangular FE assembly. To this end, a transformation matrix that defines the relationship between the area and the Cartesian parameterizations is introduced for each element in order to allow for the use of standard FE assembly procedure and define the structure (body) inertia and elastic forces. Using Bezier geometry and a linear mapping, cubic displacement fields of the new ANCF triangular elements are systematically developed. Specifically, two new ANCF triangular finite elements are developed in this investigation, namely four-node mixed-coordinate and three-node ANCF triangles. The performance of the proposed new ANCF elements is evaluated by comparison with the conventional linear and quadratic triangular elements as well as previously developed ANCF rectangular and triangular elements. The results obtained in this investigation show that in the case of small and large deformations as well as finite rotations, all the elements considered can produce correct results, which are in a good agreement if appropriate mesh sizes are used.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号