首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 187 毫秒
1.
考虑接触应力非线性分布的接触力元模式及其验证分析   总被引:1,自引:0,他引:1  
在作者提出的非连续变形计算力学模型中,采用接触力元模型描述多体接触界面上的接触特性.由于这种模型中假定接触应力沿接触界面为线性分布,从而得到的接触界面应力分布往往出现跳跃等非光滑性特征,该文对此进行了改进,采用具有高阶光滑性的非线性函数建立了能够考虑界面上接触应力非线性分布的接触力元模式,以期合理地揭示多体系统中界面的接触特性.对某一典型算例进行了数值计算,通过与大型通用非线性有限元结构分析软件ABASQUS的计算结果对比,验证了所建议计算模型的合理性与有效性.两种方法计算得到的界面接触对上的接触力基本相同;而由于采用的应力分布模式的假定不同,接触应力有所差别,由于在该文计算模型中接触对上的接触应力是按照未知量直接求得的,因此按照所建议的非线性接触力元模式所得到的接触应力更为合理.  相似文献   

2.
层合板界面层的弹簧界面元等效刚度计算模型   总被引:1,自引:1,他引:0  
赵鹏  石广玉 《计算力学学报》2011,28(Z1):131-135,146
弹簧界面元是一种最常用的层合结构界面脱层数值模拟的计算模型.但目前对应于给定的复合材料层合板界面层,弹簧界面元的长度和等效刚度的确定却没有严格的方法,所以时常因为弹簧元刚度定义的不当而导致脱层数值模拟结果的不稳定甚至计算困难.本文给出了一种基于层合板粘接层的实际厚度和材料特性确定弹簧界面元的长度和等效刚度的计算模型.通...  相似文献   

3.
变水深坝—库系统耦振分析的边界元—有限元混合法   总被引:2,自引:0,他引:2  
常用的混合元法解变水深坝-库系统的耦振,需要对变水深部分的流场进行域离散,计算工作量大,该文利用Friedman的算子函数理论,构造了势流问题在无限长带形域中的Green函数,从而使流场的边界元剖分只限于变水深区域的边界,关于坝体仍采用有限元离散,最后借助所导出的有限元-边界元格式对坝-库系统的实例作了数值计算,结果证明了它的有效性。  相似文献   

4.
在离散元-有限元耦合方法中,离散元和有限元交界面处的耦合方式对整体有砟道床的力学行为影响显著.采用基于球形单元的镶嵌单元或粘结单元模拟有砟道床时,由于球形单元和有限单元表面的自锁能力较差,使道砟层在列车载荷作用下容易产生侧向滑移,导致数值模型不稳定.此外,在实际铁路道床中,底部道砟均不同程度地嵌入路堤.为此,发展了一种嵌入式离散元-有限元耦合方法,通过设置一层嵌入地基有限元模型中的球形颗粒传递离散元域和有限元域间的力学参数,实现离散元和有限元方法的耦合.数值结果表明,嵌入式离散元-有限元耦合模型能够有效降低有砟道床的侧向位移,数值结果更加稳定,在处理与有砟道床类似的连续介质与散体介质的耦合问题时推荐采用嵌入式耦合算法.  相似文献   

5.
水利工程中混凝土重力坝在地震等作用下的三维破坏形式以及破坏过程对大坝安全影响重大。本文建立了固体材料的联合有限元-离散元(Combined Finite-Discrete Element)三维数学模型,并引入Single/Smeared断裂损伤模型,模拟混凝土材料的破坏。数学模型经验证后,用于地震作用下混凝土重力坝破坏的三维数值模拟,结果表明,联合有限元-离散元方法能够有效模拟坝体的破坏模式和破坏阈值条件,且具有将坝体连续变形与不连续断裂有机结合的优点。通过一系列数值模拟工况的对比,分析了模型网格分辨率对破坏条件、破坏模式以及程序计算效率的影响,研究结果能够为大坝的抗震设计提供有价值的参考。  相似文献   

6.
武文华  李锡夔 《力学学报》2007,39(4):473-478
提出了一个广义对流扩散方程的混合有限元方法,方程的基本变量及其空间梯度和流量在单 元内均作为独立变量分别插值. 基于胡海昌-Washizu三变量广义变分原理结合特征线法给 出了控制方程的单元弱形式. 混合元方法采用基于一点积分方案并结合可以滤掉虚假的 数值震荡的隐式特征线法. 数值结果证明了所提出的方法可以提供和四点积分同样的数 值计算结果,并能够提高计算效率.  相似文献   

7.
在文献[1]中建立的多孔介质中化学-热-渗流-力学(CTHM)本构模型基础上,针对文献[2]建立的非饱和多孔介质中热-渗流-力学耦合分析的混合有限元方法,发展了非饱和多孔介质中混合元的化学-热-渗流-力学(CTHM)耦合本构模拟算法。采用非关联流动多重屈服准则模拟非饱和多孔介质的材料非线性行为。推导了u-pw-pa-T形式的包含了耦合率本构方程积分的向后欧拉映射算法和一致性弹塑性切线模量矩阵(单元刚度矩阵)的混合元一致性算法。本文给出了临界状态线(CSL)和状态边界面(SBS)两个屈服准则的一致性算法。数值结果显示了本文所发展的混合元耦合本构模拟算法在模拟由热、化学、力学荷载共同引起的多孔介质中化学-热-渗流-力学(CTHM)耦合行为的能力和有效性。  相似文献   

8.
本文将摄动、边界元、有限元方法结合起来,提出一种求解线性蠕变问题的新方法。该方法不采用一般增量法中在一个时段内各物理量保持不变或作线性变化的假设,加大了计算步长提高了精度。文中构造了边界元摄动格式,构造了包含钢筋在内的边界元有限元耦合摄动格式,并给出了满意的数值结果。  相似文献   

9.
无限元方法及其应用   总被引:4,自引:0,他引:4  
限元是几何上趋于无穷的单元,它是一种特殊的有限元,也是对有限元在求解无界域 问题上的有效补充, 并可实现与有限元间的无缝连接.无限元分为映射无限元和非映射 无限元:映射无限元需要引入几何映射,在局部坐标系中构造插值形状函数,如Bettess 元和Astley元;非映射无限元则直接在整体坐标系中构造插值形状函数,如Burnett元. 本文评述求解无界域问题的无限元方法的研究现状和最新发展.首先介绍无限单元的概念 和无限元方法的特点;围绕求解以Helmholtz方程控制的波动问题,评述几种常规无限单 元的优劣,这些单元包括Bettess元、Astley元和Burnett元.然后介绍新近提出的广义 无限元方法,以及与常规无限元方法的区别与联系.最后对无限元方法在各种问题中的 应用做了总结.  相似文献   

10.
元胞单元法   总被引:4,自引:0,他引:4  
介绍一种力学计算方法——元胞单元法.它将结构的整体求解变成局域分析,通过力的局域间的不平衡传递达到最终的整体平衡.它对计算机容量要求低,在大规模计算方面具有良好的前景.同时,由于借用元胞自动机的思想,可形成一种高度并行的算法,适用于未来并行机的要求.是一种有发展潜力的数值方法,可望在大型结构和纳米力学的大规模摸拟方面得到应用.文中介绍了元胞单元法的基本原理,将其用于平面问题分析,进行了数值试验,给出了计算步骤和结果,并对其优点和存在的问题进行了讨论.  相似文献   

11.
A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.  相似文献   

12.
The element energy projection (EEP) method for computation of super- convergent resulting in a one-dimensional finite element method (FEM) is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton's method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation (ODE) of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.  相似文献   

13.
IntroductionManyinfiniteproblemscanbefoundincivilengineering ,suchastunnelconstruction ,structurefoundation ,etc ..Forotherengineeringproblems,whenthephysicaldimensionsofanobjectaresmallandthesurroundingmediaorstructuresaremuchbiggerthantheobject,thenumericalcalculatingmodelcanbetreatedasoneinaninfinitefield .Sofar,onlyafewanalyticalsolutionsforinfinitefieldproblemscanbefound[1- 3].ManyprojectsrelatedtoinfinitefieldproblemsaresolvedbytheFEM ,whereinfiniteelementmethodsareused[4 ,5 ].Sometime…  相似文献   

14.
采用格林公式和基本解推导出直接边界积分方程来求解渗流问题.边界积分方程数值离散基于格林元方法(Green element methond),改进了原方法中压力和压力导数的求解方法,命名为混合边界元方法(Mixed boundary element method).相较于格林元类方法,该方法显式考虑了求解节点的外法向流量值和压力值,并使求得的数值解在求解区域上能够连续,符合实际的物理过程,在不增加额外未知数的情况下提高了计算精度.分析了不同网格类型对模拟计算结果的影响,并对稳定渗流问题、非稳定(瞬态)渗流问题和非稳态问题进行了实例计算,结果显示改进方法提高了计算精度,并对各类渗流问题有较好的适应性.  相似文献   

15.
理性有限元     
提出了与常现有限元迥然不同的理性有限元列式。其方法论的差别在于理性有限元充分考虑了力学的微分方程,用方程的解来逼近单元内部场.即以力学的需求为主导,再用相应的数学方法推导。并用理住平面四边形元RQ4为典型予以表述。数值结果表明理性有限元的优越性质。  相似文献   

16.
Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.  相似文献   

17.
无网格Galerkin法(EFGM)处理不可压缩问题时不存在自锁现象,有限元方法(FEM)也常被用来与其耦合以方便地施加边界条件和提高计算效率。在有限元方法中使用等参元,EFGM与FEM的耦合方法在处理不可压缩问题时仍然存在自锁现象。本文在有限元方法中,采用非协调元,将无网格kGalerkin法与非协调元耦合,保留了耦合方法的优点,且避免了求解不可压缩问题时的自锁现象。算例显示本文方法在分析平面应变不可压缩问题时能得到合理的结果。  相似文献   

18.
A modified boundary element method (BEM) and the DEVSS‐G finite element method (FEM) are applied to model the deformation of a polymeric drop suspended in another fluid subjected to start‐up uniaxial extensional flow. The effects of viscoelasticity, via the Oldroyd‐B differential model, are considered for the drop phase using both FEM and BEM and for both the drop and matrix phases using FEM. Where possible, results are compared with the linear deformation theory. Consistent predictions are obtained among the BEM, FEM, and linear theory for purely Newtonian systems and between FEM and linear theory for fully viscoelastic systems. FEM and BEM predictions for viscoelastic drops in a Newtonian matrix agree very well at short times but differ at longer times, with worst agreement occurring as critical flow strength is approached. This suggests that the dominant computational advantages held by the BEM over the FEM for this and similar problems may diminish or even disappear when the issue of accuracy is appropriately considered. Fully viscoelastic problems, which are only feasible using the FEM formulation, shed new insight on the role of viscoelasticity of the matrix fluid in drop deformation. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

19.
This paper describes a numerical approximation scheme for the natural convection (NC) flow in a fluid-saturated porous medium. Our formulation of the problem is based on the mixed finite element method (FEM). Using the so-called consistent splitting scheme, a second-order accuracy in time and in space is verified by the numerical calculation. The resulting flow patterns and heat transfer for different Rayleigh numbers, Darcy numbers and porosities are numerically studied by the proposed scheme.  相似文献   

20.
A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.  相似文献   

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

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