首页/文章/ 详情

UEL如何与接触非线性模型结合

1月前浏览419
   

概述

       

分别采用经典有限元法和比例边界有限法在abaqus中编制了C3D8单元的UEL程序,计算结果均与abaqus保持同样精度。

在此基础上将用户自定义单元与abaqus自带的接触非线性模型相结合,并完成了单元的位移、速度和加速度的可视化。

算例介绍

       

计算了滑块在斜面滑落算例。滑块分两个运动过程,分别为



受重力和摩擦力沿着斜面作直线匀加速运动离开斜面受重力作用作斜抛加速运动
 

采用三种方法计算,分别为




ABAQUS自带的C3D8单元采用UEL子程序,自编C3D8 UEL单元,且考虑B-BAR修正采用UEL子程序,依据比例边界有限元法自编C3D8 UEL单元
 

关于方法1,直接在ABAQUS中建立模型,给滑块施加重力荷载,在滑块和斜面之间建立接触即可。

关于方法2,首先编写C3D8 UEL程序,然后再给自定义单元世界重力荷载,并且在滑块和斜面之间建立接触。但是给用户自定义单元施加重力荷载不可以直接施加,需要借助虚假单元。在滑块的斜面建立接触过程中,因为用户自定义单元无法定义面,定义接触更无从谈起,所以也需要借助虚假单元才能定义SURFACE,最终完成SURFACE TO SURFACE接触定义

关于方法3,同方法2,首先依据SBFEM在UEL中完成程序编写,后面同样的操作,采用虚假单元给用户自定义单元施加重力荷载与定义接触。

帖子涉及到的内容:UEL二次开发、基于ABAUQS的SBFEM、用户自定义单元体力的施加、用户自定义单元面定义、用户自定义单元接触非线性分析

计算模型信息

       

滑块与斜面的几何尺寸

   

滑块与斜面的网格

   

该算例中的三种方法均采用动力隐式算法,计算总时长为10s,设置固定增量步长为0.01s,增量步数共计1000。重力加速度为9.81m/s²。具体参数见附件。

滑块两个运动过程受力分析为

   

上面左图为斜面内的受力分析,右图为斜面外的受力分析。

可知滑块沿斜面的分力与摩擦力大小比值有两种情况:

   

当滑块在斜面内运动的时,遵循如下控制方程:

   

当滑块做出斜面运动时,遵循如下控制方程:

   

通过上式,可求得位移和速度的解析解。

      滑块运动水平向位移时程对比如下图:

   

滑块运动竖直向位移时程对比如下图:

   

滑块运动水平向位速度时程对比如下图:

   

滑块运动竖直向速度时程对比如下图:

   

误差分析:观察四个时程曲线,发现滑块在斜面上运动的过程中,四种计算结果均吻合良好。而在滑块离开斜面的第二运动过程中,三个数值解吻合较好,但均与解析解有一定的差别。这是因为解析求解是建立在质点的假设上,但是数值解是一个小方块,因此方块第一个角点离开斜面的时候,方块并未完全离开斜面,所以误差会从方块的第一个角点离开斜面开始。

基于虚假单元的UEL可视化

       

在该算例中,完成荷载和接触定义的同时,附加的虚假单元兼具用户自定义单元的位移、速度和加速度可视化功能。

当虚假单元和用户自定义单元共存时,ABAQUS显示为

   

当只显示用户自定义单元的时候,ABAQUS显示如下图:

   

当只显示虚假单元的时候,ABAQUS显示为

   

可见虚假单元能够非常完美的完成单元位移、速度和加速度的可视化。在一定程度上解决了用户自定义单元在abaqus中无法可视化的问题

     

来源:有限元先生
Abaqus非线性二次开发控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-09-29
最近编辑:1月前
外太空土豆儿
硕士 我们穷极一生,究竟在追寻什么?
获赞 2粉丝 1文章 42课程 0
点赞
收藏
作者推荐

基于abaqus的sbfem及其在混凝土坝非线性地震响应分析中的应用

零概述壹基于abaqus的sbfem采用abaqus提供的uel和uexternaldb子程序二次开发平台将sbfem嵌入到abaqus中,其中子程序uel用于实现sbfem中关键矩阵的计算,包括amartix、rhs等等,下面是基于sbfem的uel子程序流程图子程序uexternaldb用于将任意多面体单元以外部文件的形式导入进abaqus中,下面展示的是2个任意多面体的单元-节点信息上述的单元-节点信息无法仅仅从节点顺序重现单元的几何拓扑关系,因此需要准备外部txt文件将上述单元-节点信息导入到abaqus中,下面的数据是上面2个单元的信息采用Fortran程序编写完成uel之后,还需要在inp文件中定义用户自定义单元,目前的abaqus只支持手动修改,对于上述单元,在inp文件中对应的关键字为贰耦合sbfem-femabaqus有许许多多的非线性模型,比如材料模型和界面模型,因此发展一种耦合方法,将用户自定义单元与abaqus自带的非线性模型相结合,以便同时发挥abaqus和用户自定义单元的优势,下面简要介绍两种耦合算法。首先是将自定义单元与abaqus的接触非线性模型结合,需要在用于自定义单元的界面进行处理,如其对应的公式为如果要对用户自定义单元进行可视化,需要在原有的自定义单元基础上再嵌套一层abaqus自带的单元,并将自带单元的材料属性设置为无限小,如肆悬臂梁受动荷载设计悬臂梁受动荷载算例,验证子程序求解刚度矩阵、阻尼矩阵和质量矩阵的正确性,模型的边界条件和受荷载为荷载表达式为采用三种单元对该几何模型进行离散,三种单元分别为六面体单元四面体单元棱柱体单元统计悬臂端关键点的位移等数据,最终计算的结果基本一致。伍koyna混凝土坝地震响应分析以koyna混凝土重力坝为研究对象,将文章提出的方法应用到实际工程的地震响应分析,koyna混凝土坝的截面为采用两种单元对该模型进行离散,分别为六面体单元棱柱体单元对上述网格设置三个计算工况分别为case1:六面体单元+fem计算case2:六面体单元+sbfem计算case3:棱柱体单元+sbfem计算首先对上述计算模型进行自振频率分析,统计前20阶自振频率数值对比为可以发现自振频率分析计算结果保持一致。然后对koyna混凝土坝进行地震动力响应分析,提取了坝顶-坝锺的相对数据,分别为竖直相对位移水平相对位移竖直相对速度水平相对速度计算结果均保持一致。计算时长比较为可以发现sbfem的计算时长均比fem有缩短,其中case3的计算时长缩短最明显,说明棱柱体单元有较好的计算性能。陆NG5拱坝地震响应分析选取NG5拱坝为研究对象,考虑坝体横缝接触非线性与土-结构动力相互作用模型,探究文章提出的方法在拱坝地震响应分析中的作用。NG5拱坝的平面布置图及其横缝为坝体-地基-库水的几何模型为采用六面体单元离散几何模型,最终的计算模型为对上述计算模型设置三种计算工况土-结构动力相互作用采用黏弹性边界模拟,黏弹性边界是计算模型边界节点上一系列的弹簧-阻尼器元件坝体横缝接触非线性采用abaqus自带的指数接触模型输入地震动的位移和速度为坝体横缝附近的关键点为首先统计了坝体横的开度时程曲线统计了关键点的动水压力时程曲线,c点数据为d点数据为e点数据为统计了15s时刻坝面的动水压力云图结合以上数据,可以发现,考虑了土-结构动力相互作用以后,混凝土坝整体的响应,包括相对位移、横缝开度和坝面动水压力等数据都有不同程度的降幅,说明了土-结构动力相互作用在混凝土坝抗震计算中应予以重点考虑。柒结论Thescaledboundaryfiniteelementmethod(SBFEM)isimplementedintheABAQUSthroughtheusersubroutineinterface.WiththecoupledSBFEMandstandardfiniteelementmethod(FEM),thedynamicinteractionmodelofconcretedam–reservoir–foundationconsideringthegeometricalnonlinearityoftransversejointsisestablishedusingtheviscous–springabsorbingboundaryconditions.Theoverlayingelementtechniqueisappliedtodefinethecontactinterfaceoftransversejoints.ThedynamiccharacteristicsandlinearseismicresponseoftheKoynagravitydamhavebeenanalyzedusingpolyhedralelements.Finally,nonlinearseismicanalysisoftheNG5archdamhasbeenconducted.Thefollowingconclusionscanbedrawn:(1)Accordingtoseveralbenchmarkexamples,theproposedmethodhasgoodaccuracyandefficiency.(2)Theoverlayingelementtechnique,whichdoesnotrequiretheintroductionofextraDOFs,completesthecombinationoftheSBFEMwiththebuilt-inelements,models,andotherfunctionsofABAQUS.(3)BasedonthecoupledSBFEM-FEM,thedynamicinteractionmodelofdam,reservoirandfoundation,thenonlinearityofjointcontact,andtheseismicinputcanbefullyconsidered.ThecalculatedresultsareconsistentwiththoseofABAQUSbuilt-inelementswithahighcomputationalefficiency.Itwillprovideanovelmethodforthedynamicanalysisofcomplexengineeringstructures.来源:有限元先生

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈