第9章有限元建模方法

第9章有限元建模方法
第9章有限元建模方法

abaqus有限元分析过程

一、有限单元法的基本原理 有限单元法(The Finite Element Method)简称有限元(FEM),它是利用电子计算机进行的一种数值分析方法。它在工程技术领域中的应用十分广泛,几乎所有的弹塑性结构静力学和动力学问题都可用它求得满意的数值结果。 有限元方法的基本思路是:化整为零,积零为整。即应用有限元法求解任意连续体时,应把连续的求解区域分割成有限个单元,并在每个单元上指定有限个结点,假设一个简单的函数(称插值函数)近似地表示其位移分布规律,再利用弹塑性理论中的变分原理或其他方法,建立单元结点的力和位移之间的力学特性关系,得到一组以结点位移为未知量的代数方程组,从而求解结点的位移分量. 进而利用插值函数确定单元集合体上的场函数。由位移求出应变, 由应变求出应力 二、ABAQUS有限元分析过程 有限元分析过程可以分为以下几个阶段 1.建模阶段: 建模阶段是根据结构实际形状和实际工况条件建立有限元分析的计算模型――有限元模型,从而为有限元数值计算提供必要的输入数据。有限元建模的中心任务是结构离散,即划分网格。但是还是要处理许多与之相关的工作:如结构形式处理、集合模型建立、单元特性定义、单元质量检查、编号顺序以及模型边界条件的定义等。

2.计算阶段:计算阶段的任务是完成有限元方法有关的数值计算。 由于这一步运算量非常大,所以这部分工作由有限元分析软件控制并在计算机上自动完成 3.后处理阶段: 它的任务是对计算输出的结果惊醒必要的处理, 并按一定方式显示或打印出来,以便对结构性能的好坏或设计的合理性进行评估,并作为相应的改进或优化,这是惊醒结构有限元分析的目的所在。 下列的功能模块在ABAQUS/CAE操作整个过程中常常见到,这个表简明地描述了建立模型过程中要调用的每个功能模块。 “Part(部件) 用户在Part模块里生成单个部件,可以直接在ABAQUS/CAE环境下用图形工具生成部件的几何形状,也可以从其它的图形软件输入部件。 Property(特性) 截面(Section)的定义包括了部件特性或部件区域类信息,如区域的相关材料定义和横截面形状信息。在Property模块中,用户生成截面和材料定义,并把它们赋于(Assign)部件。 Assembly(装配件) 所生成的部件存在于自己的坐标系里,独立于模型中的其它部件。用户可使用Assembly模块生成部件的副本(instance),并且在整体坐标里把各部件的副本相互定位,从而生成一个装配件。 一个ABAQUS模型只包含一个装配件。

有限元网格划分的基本原则

有限元网格划分的基本原则 划分网格是建立有限元模型的一个重要环节,它要求考虑的问题较多,需要的工作量较大,所划分的网格形式对计算精度和计算规模将产生直接影响。为建立正确、合理的有限元模型,这里介绍划分网格时应考虑的一些基本原则。 1 网格数量 网格数量的多少将影响计算结果的精度和计算规模的大小。一般来讲,网格数量增加,计算精度会有所提高,但同时计算规模也会增加,所以在确定网格数量时应权衡两个因数综合考虑。图1中的曲线1表示结构中的位移随网格数量收敛的一般曲线,曲线2代表计算时间随网格数量的变化。可以看出,网格较少时增加网格数量可以使计算精度明显提高,而计算时间不会有大的增加。当网格数量增加到一定程度后,再继续增加网格时精度提高甚微,而计算时间却有大幅度增加。所以应注意增加网格的经济性。实际应用时可以比较两种网格划分的计算结果,如果两次计算结果相差较大,可以继续增加网格,相反则停止计算。 图1 位移精度和计算时间随网格数量的变化 在决定网格数量时应考虑分析数据的类型。在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。如果需要计算应力,则在精度要求相同的情况下应取相对较多的网格。同样在响应计算中,计算应力响应所取的网格数应比计算位移响应多。在计算结构固有动力特性时,若仅仅是计算少数低阶模态,可以选择较少的网格,如果计算的模态阶次较高,则应选择较多的网格。在热分析中,结构内部的温度梯度不大,不需要大量的内部单元,这时可划分较少的网格。 2 网格疏密 网格疏密是指在结构不同部位采用大小不同的网格,这是为了适应计算数据的分布特点。在计算数据变化梯度较大的部位(如应力集中处),为了较好地反映数据变化规律,需要采用比较密集的网格。而在计算数据变化梯度较小的部位,为减小模型规模,则应划分相对稀疏的网格。这样,整个结构便表现出疏密不同的网格划分形式。图2是中心带圆孔方板的四分之一模型,其网格反映了疏密不同的划分原则。小圆孔附近存在应力集中,采用了比较密的网格。板的四周应力梯度较小,网格分得较稀。其中图b中网格疏密相差更大,它比图a中的网格少48个,但计算出的孔缘最大应力相差1%,而计算时间却减小了36%。由此可见,采用疏密不同的网格划分,既可以保持相当的计算精度,又可使网格数量减

三维有限元建模方法的研究现状

三维有限元建模方法的研究现状 作者:陈琼 作者单位:复旦大学附属华山医院口腔科,上海,200040 刊名: 口腔医学 英文刊名:STOMATOLOGY 年,卷(期):2006,26(2) 被引用次数:18次 参考文献(25条) 1.李青奕;董寅生;陈文静预加载"L"形曲力学行为的有限元分析[期刊论文]-口腔医学 2004(01) 2.Hirabayashi M;Motoyoshi M;Ishimarn T Stresses in mandibular cortical bone during mastication:biomechanical considerations using a three-dimensional finite element method 2002(01) 3.许文翠;陈文静;董寅生垂直曲的力学行为的研究[期刊论文]-口腔医学 2002(01) 4.周学军;赵志河;赵美英包括下颌骨的颞下颌关节三维有限元模型的建立[期刊论文]-实用口腔医学杂志 2000(01) 5.李玲;张睿;于力牛基于CT断层影像的下颌骨及下牙列三维几何学仿真[期刊论文]-上海口腔医学 2000(04) 6.于力牛;常伟;王成焘基于实体模型的牙颌组织三维有限元建模问题探讨[期刊论文]-机械设计与研究 2002(02) 7.张富强;魏斌;李玲牙颌组织及修复体三维几何学、有限元模型的设计[期刊论文]-上海口腔医学 2002(03) 8.陈剑虹一种基于断层测量的快速反求系统关键技术研究[学位论文] 2000 9.魏洪涛;张天夫;曾晨光牙颌三维有限元模型生成方法的探讨[期刊论文]-白求恩医科大学学报 2000(02) 10.朱静有限元分析方法在口腔临床中的应用进展[期刊论文]-上海生物医学工程 2003(03) 11.Huiskes R;Chao EY A survey of finite element analysis in orthopedic biomechanics:the first decade [外文期刊] 1983(06) 12.王宁;吴凤鸣;周小陆金属烤瓷冠瓷颈缘与金属颈缘的三维有限元应力分析[期刊论文]-口腔医学 2004(04) 13.龚璐璐口腔修复生物力学中三维有限元法应用的研究进展及展望[期刊论文]-医用生物力学 2002(02) 14.Aydin AK;Tekkaya AE Stresses induced by different loading around weak abutments[外文期刊] 1992(06) 15.Verdonschot N;Fennis WM;Kuijs R Generation of three-dimensional finite models of restored human teeth using micro-CT techniques 2001(04) 16.张富强;魏斌;于力牛个性化牙颌组织三维有限元模型库的建立[期刊论文]-上海口腔医学 2004(02) 17.于力牛;尚鹏;王成焘适用于口腔修复学的模块化牙列有限元建模[期刊论文]-上海交通大学学报 2002(08) 18.于力牛;张睿;李玲模块化牙列三维有限元模型的建立[期刊论文]-上海口腔医学 2000(04) 19.Nagasao T;Kobayashi M;Tsuchiya Y Finite element analysis of the stresses around endosseous implants in various reconstructed mandibular models 2002(03) 20.李玲上下颌三维重建及有限元建模[学位论文] 2001 21.李志华;陈天云;刘剑上颌第一磨牙的三维有限元模型的建立[期刊论文]-实用临床医学 2001(01) 22.张彤;刘洪臣;王延荣上颌骨复合体三维有限元模型的建立[期刊论文]-中华口腔医学杂志 2000(05) 23.高勃;王忠义;施长溪牙冠表面形状测量造型方法[期刊论文]-实用口腔医学杂志 1999(04) 24.牛晓明;李江;吴清文利用CAD/CAE技术进行骨骼的计算机模拟仿真[期刊论文]-光学精密工程 1999(06) 25.蒋孝煜有限元法基础 1992

专业ABAQUS有限元建模经验笔记

基于ABAQUS的有限元分析和应用 第一章绪论 1.有限元分析包括下列步骤: 2.为了将试验数据转换为输入文件,分析者必须清楚在程序中所应用的和由实验人员提供的材料数据的应力和应变的度量。 3.ABAQUS建模需注意以下内容: 4.对于许多包含过程仿真的大变形问题和破坏分析,选择合适的网格描述是非常重要的,需要认识网格畸变的影响,在选择网格时必须牢牢记住不同类型网格描述的优点。 第二章ABAQUS基础 1.一个分析模型至少要包含如下的信息:离散化的几何形体、单元截面属性、材料数据、载荷和边界条件、分析类型和输出要求。 ①离散化的几何形体:模型中所有的单元和节点的集合称为网格。 ②载荷和边界条件: 2.功能模块: (1)Assembly(装配):一个ABAQUS模型只能包含一个装配件。 (2)Interaction(相互作用):相互作用与分析步有关,这意味着用户必须规定相互作用是在哪些分析步中起作用。 (3)Load(载荷):载荷和边界条件与分析步有关,这意味着用户指定载荷和边界条件是在哪些分析步中起作用。 (4)Job(作业):多个模型和运算可以同时被提交并进行监控。 3.量纲系统 ABAQUS没有固定的量纲系统,所有的输入数据必须指定一致性的量纲系统,常用的一致性量纲系统如下:

4.建模要点 (1)创建部件:设定新部件的大致尺寸的原则必须是与最终模型的最大尺寸同一量级。(2)用户应当总是以一定的时间间隔保存模型数据(例如,在每次切换功能模块时)。(3)定义装配: 在模型视区左下角的三向坐标系标出了观察模型的方位。在视区中的第2个三向坐标系标出了坐标原点和整体坐标系的方向(X,Y和Z轴)。 (4)设置分析过程: (5)在模型上施加边界条件和荷载: 用户必须指定载荷和边界条件是在哪个或哪些分析步中起作用。 所有指定在初始步中的力学边界条件必须赋值为零,该条件是在ABAQUS/CAE中自动强加的。 在许多情况下,需要的约束方向并不一定与整体坐标方向对齐,此时用户可定义一个局部坐标系以施加边界条件。 在ABAQUS中,术语载荷通常代表从初始状态开始引起结构响应发生变化的各种因素,包括:集中力、压力、非零边界条件、体力、温度(与材料热膨胀同时定义)。

MSC Patran与LR ShipRight有限元建模技术的分析与比较

MSC.Patran与LR.ShipRight有限元建模技术的分析与比较 作者:江南造船集团朱彦 摘要:本文基于散货船CSR 探讨使用MSC.Patran 与LR.ShipRight 两款软件在进行有限元分析中的建模技术,并比较两款软件的特点以及相互联系。 关键字:Patran、ShipRight、散货船、CSR、有限元建模 1. 前言 在船舶详细设计阶段,对船体结构进行应力集中以及疲劳强度评估的一个有效的手段就是采用有限元分析。有限元分析的一般方法为选择有限元分析软件、确定单元形式、建立几何模型、网格划分、确定边界条件、判断载荷工况等,具体又可归纳为四个步骤: 1) 建立有限元模型; 2) 确定载荷及边界条件; 3) 进行详细应力应变评估(例如细化网格以评估高应力区域); 4) 对关键部位的结构进行疲劳强度评估。 在以上步骤中能否建立合理有效的有限元模型是前提条件,模型质量的好坏,特别是网格的类型与划分方法,直接影响后续的分析结果。目前常用的有限元分析软件主要有 MSC.Patran\Nastran、LR.ShipRight、基于Patran 的CCS.TOOLS、DNV.Sesam 等,本文以散货船CSR 有限元建模为例,探讨Patran 与ShipRight 两种软件的建模技术和异同点。 2. Patran 与ShipRight 的简介 MSC.Patran 作为一个优秀的前后置处理器,具有高度的集成能力和良好的适用性,模型处理智能化、自动有限元建模、分析的集成、用户自主开发新功能、分析结果的可视化处理等等是其典型的特征,它提供了功能全面、方便灵活的可满足各种精度要求的复杂有限元的建模功能,其综合全面先进的网格划分技术,为用户根据不同的几何模型提供了多种不同的生成和定义的有限元模型工具。 ShipRight 是LR 自主开发的一款基于CSR 的有限元分析应用软件,具有很强的针对性,

有限元分析过程

有限元分析过程可以分为以下三个阶段: 1.建模阶段: 建模阶段是根据结构实际形状和实际工况条件建立有限元分析的计算模型——有限元模型,从而为有限元数值计算提供必要的输入数据。有限元建模的中心任务是结构离散,即划分网格。但是还是要处理许多与之相关的工作:如结构形式处理、集合模型建立、单元特性定义、单元质量检查、编号顺序以及模型边界条件的定义等。 2.计算阶段: 计算阶段的任务是完成有限元方法有关的数值计算。由于这一步运算量非常大,所以这部分工作由有限元分析软件控制并在计算机上自动完成。 3.后处理阶段: 它的任务是对计算输出的结果惊醒必要的处理,并按一定方式显示或打印出来,以便对结构性能的好坏或设计的合理性进行评估,并作为相应的改进或优化,这是惊醒结构有限元分析的目的所在。 注意:在上述三个阶段中,建立有限元模型是整个有限分析过程的关键。首先,有限元模型为计算提供所以原始数据,这些输入数据的误差将直接决定计算结果的精度;其次,有限元模型的形式将对计算过程产生很大的影响,合理的模型既能保证计算结构的精度,又不致使计算量太大和对计算机存储容量的要求太高;再次,由于结构形状和工况条件的复杂性,要建立一个符合实际的有限元模型并非易事,它要考虑的综合因素很多,对分析人员提出了较高的要求;最后,建模所花费的时间在整个分析过程中占有相当大的比重,约占整个分析时间的70%,因此,把主要精力放在模型的建立上以及提高建模速度是缩短整个分析周期的关键。 原始数据的计算模型,模型中一般包括以下三类数据: 1.节点数据: 包括每个节点的编号、坐标值等; 2.单元数据: a.单元编号和组成单元的节点编号;b.单元材料特性,如弹性模量、泊松比、密度等;c.单元物理特征值,如弹簧单元的刚度系数、单元厚度、曲率半径等;d.一维单元的截面特征值,如截面面积、惯性矩等;e.相关几何数据 3.边界条件数据:a.位移约束数据;b.载荷条件数据;c.热边界条件数据;d.其他边界数据. 建立有限元模型的一般过程: 1.分析问题定义 在进行有限元分析之前,首先应对结果的形状、尺寸、工况条件等进行仔细分析,只有正确掌握了分析结构的具体特征才能建立合理的几何模型。总的来说,要定义一个有限元分析问题时,应明确以下几点: a.结构类型; b.分析类型; c.分析内容; d.计算精度要求; e.模型规模; f.计算数据的大致规律 2.几何模型建立 几何模型是从结构实际形状中抽象出来的,并不是完全照搬结构的实际形状,而是需要根据结构的具体特征对结构进行必要的简化、变化和处理,以适应有限元分析的特点。 3.单元类型选择 划分网格前首先要确定采用哪种类型的单元,包括单元的形状和阶次。单元类型选择应根据结构的类型、形状特征、应力和变形特点、精度要求和硬件条件等因素综合进行考虑。

随机地质建模技术方法简介

随机地质建模技术方法简介 李 燕 (胜利油田物探研究院,山东东营 257000) 摘 要:随机建模是指以已知的信息为基础,以随机函数为理论,应用随机模拟方法,产生可选的、等概率的储层模型方法。该方法承认控制点以外的储层参数具有一定的不确定性,即具有一定的随机性。Deautch等根据模拟单元的特征,将随机模型分为基于目标的随机模型和基于象元的随机模型。 关键词:随机建模;克立金方程;地质统计学;储层结构 地下储层本身是确定的,在每一个位置点都具有确定的性质和特征。但是,地下储层又是复杂的,它是许多复杂地质过程(沉积作用、成岩作用和构造作用)综合作用的结果,具有复杂的储层结构(储层相)空间配置及储层参数的空间变化。在现有资料不完善的条件下,人们对它的认识总会存在一些不确定的因素,难于掌握任意尺度下储层的真实特征或性质。特别是对于连续性较差且非均质性强的陆相储层来说,更难于精确表征储层的特征。从而认为储层描述便具有不确定性即随机性。 1 随机建模技术的产生和发展 在地质统计学技术的形成和发展中,法国枫丹白露地质统计学与数学形态学中心起了重要的作用,其核心人物M atheron是地质统计学的创始人。他的许多学生(如Journel,David等)后来都成了该领域的继承者和发展者。在随机建模的发展中, Jo urnel领导的斯坦福大学油藏预测中心则是令人起敬的先锋。他们研制的GSLIB是公认的较完整、先进的地质统计学软件包。近年来研制了许多随机建模的算法,并做了应用研究。另外加拿大的David、原英国BP公司的H aldorsen、加拿大FSS International公司的Srivastava、美国斯坦福大学的Deutsch以及科罗拉多矿业学院、得克萨斯大学澳斯万分校、挪威计算中心、澳大利亚新南威尔士大学等处的一些学者都在这一领域有很高的造诣。 地质统计学创建于本世纪60年代初期,当时人们基本上把克里金作为地质统计学的同义词。70年代末,Jo urnel(1978)在所著的《Minging Geostatistics》一书中,介绍了随机建模的基本思想。80年代中后期,尤其在90年代,随着克里金方法不但被用作插值方法,越来越多的被用来建立数据的条件累积分布函数(CCDF),随机建模得到了飞速发展。出于对解决不同问题的需要以及对时间、经费、人力和软硬件的考虑,发展了种类繁多、功能不同的随机建模方法和算法。 地质统计学引入我国较晚,早期都把克里金认为是地质统计学。随机建模仅在近几年才得到重视,并引入油藏勘探开发研究中。西安石油学院张团峰、王家华等人(1995a,b)在引进国外资料的基础上,研制了一套储层地质统计分析系统(GASOR2.0),可用于建立储层模型。北京石油勘探开发科学研究院刘明新等人在“八五”期间利用分形理论进行了储层建模研究。胜利油田“八五”期间在其研制的油藏描述软件中也加进了随机建模内容。一些青年学者在利用随机建模解决油田问题方面做了有益的工作;石油大学纪发华(1994)在其博士论文中利用随机建模技术对油藏特征做了研究,利用序贯指示模拟、模拟退火研究了渗透率的空间分布。文键(1995)在其博士论文讨论了随机建模技术应用中的几个问题: 统计特征量与储层空间分布非均质性特征的关系; 储层空间分布不确定性对开发可行性研究的影响;统计特征量与样本间距、容量的关系;得出了很有价值的经验(诸如岩性指示变差函数与砂岩面密度结合和表征砂体连续性特征),同时还利用序 收稿日期:2009-07-28 作者简介:李燕(1973—),女,现从事岩石物理反演工作。

薄壁组合结构的有限元快速建模技术

薄壁组合结构的有限元快速建模技术 作者:北京航空制造工程研究所岳中第高光波 摘要:在并行设计过程中,面向数字化、无纸化的技术要求,如何实现复杂薄壁组合结构的有限元快速建模,是工程设计人员最为关心的重要技术问题之一。本文基于MSC的开放平台,综述了薄壁组合结构三维实体设计中FEA建模的特点与复杂性,探索了实现CAD/CAE集成与二次开发、FEA结构模型简化及快速建模的技术途径与方法。应用测试表明,它能有效地缩短薄壁组合结构的有限元建模周期,使有限元模型在简化上更能精确于CAD 数字化模型。 1 薄壁组合结构FEA建模的特点与复杂性 飞机与舰船都是复杂的薄壁组合结构,是MSC/NASTRAN系统的传统工程应用领域。如何快速有效地为MSC/NASTRAN系统建立有限元分析(FEA)模型,往往成为这些行业工程技术人员的日常工作。今天通过电子商务,我们能够容易地采购到现代化的数字化技术的基础平台,例如CAD/CAE/CAM/PDM等商用软件,包括IBM/CATIA、MSC/NASTRAN 等产品。在这方面,我们与国外同类公司相比,相差无几;但是,在应用这些基础平台解决复杂薄壁组合结构的分析方面,我们与国外同行确有不小的差距。例如,某些军民机要建立整机的有限元分析模型,需要集中近20多人,准备一年以上,而国外同行解决同样的问题所用的人力与时间却要少得多;在方案设计阶段要实现整机有限元分析模型,国外需要几个星期,而我们加班加点却还需要几个月。同样的软硬件技术环境,却有如此不同的结果,其原因在于我们缺乏在通用数字化平台基础上的实用开发技术与个性化的技术支持。 CAD软件商向我们推荐的整体解决方案,即三维CAD实体模型自动转换为三维有限元分析模型,在零件实体一级也许可以使用(例如,CATIA的FEA功能),但在复杂薄壁组合结构的部件级或整机级,这种方案往往导致几十万至几百万节点的FEA模型,而最终所得到的结果却无法让专家们相信它的正确性。在这种解决方案里,薄壁组合结构视为三维连续实体结构,CAD/3D模型直接转换为FEA/3D模型,“组合”变成“熔合”,“薄壁”变成“体元”,使结构的传力与力学特性受到扭曲。 正因为如此,复杂薄壁组合结构的FEA建模应该有它固有的力学方法和途径。我国工程技术人员利用薄壁组合结构(板杆梁)理论简化飞行器结构模型的基本实践,已经在航空航天领域极大地普及。自“九五”以来,我们按照传统的薄壁组合结构的模型简化方法,探索

GOCAD 软件三维地质建模方法

GOCAD 软件三维地质建模方法 1建模方法 GOCAD 三维地质建模主要包括两类:一类是构造模型(structural modeling)建模,一类是三维储层栅格结构(3D Reservoir Grid Construction)建模。 (1)构造模型(structural modeling)建模建立地质体构造模型具有非常重要的意义。通过建立构造模型能够模拟地层面、断层面的形态、位置和相互关系;结合反映地质体的各种属性模型的可视化图形,还能够用于辅助设计钻井轨迹。此外,构造模型还是地震勘探过程中地震反演的重要手段。 (2)三维储层栅格结构(3D Reservoir Grid Construction)建模根据建立的构造模型,在3D Reservoir Grid Construction 中可以建立其体模型;同时地质体含有多种反映岩层岩性、资源分布等特性的参数,如岩层的孔隙度、渗透率等,可对这些物性参数进行计算和综合分析,得到地质体的物性参数模型。 当采样值在地质体内密集、规则分布时,可以直接建立采样值到应用模型的映射关系,把对采样值的处理转化为对物性参数的处理,这样可以充分利用计算机的存储量大、计算速度快的特点。 当采样值呈散乱分布,并且数据量有限时,需要采用数学插值方法,拟合出连续的数据分布,充分利用由采样值所隐含的数据场的内部联系,精确的模拟模型中属性场的分布。 图1-1孔隙度参数模型分布图 2 建模流程 2.1数据分析 (1)钻孔、测井分布及数据分析 支持三维建模的数据主要为钻孔和测井。由于对区域范围和建立三维地质建模的精度要求不同,得对所得到的钻孔、测井的分布和根据其取得的数据进行分析和处理是的必要。根据钻孔、测井的分布范围和稠密程度可以大致确定地层的分布界限,对钻孔较少区域采取补充钻探或者采用其它方法进行处理。 (2)地质剖面

ANSYS有限元网格划分的基本原则

ANSYS有限元网格划分的基本原则 引言 ANSYS中有两种建立有限元模型的方法:实体建模和直接生成。使用实体建模,首先生成能描述模型的几何形状的几何模型,然后由ANSYS程序按照指定的单元大小和形状对几何体进行网格划分产生节点和单元。对于直接生成法,需要手工定义每个节点的位置和单元的连接关系。 一般来说对于规模较小的问题才适于采用直接生成法,常见的问题都需要先通过实体建模生成几何模型,然后再对其划分网格生成有限元模型。随着计算机性能的提高,分析模型的复杂性和规模都越来越大,而直接生成法也因其自身的局限性逐渐的被淘汰,所以正确的理解划分网格的目的和掌握划分网格的方法不论是对ANSYS的学习还是对二次开发都有重要的作用,尤其是当模型复杂度大,对模型的某些部分网格需要特殊处理时,这种对划分网格深度的理解作用更加明显。 2 常用高级网格划分方法 随着ANSYS功能的越来越强大和计算机性能的飞速提高,有限元分析向着大型化、复杂化的方向发展,而划分网格的观念也需要逐渐从二维模型向三维模型上上转变。这里主要描述三种常见的高级划分网格的方法,正确的理解和掌握这些划分网格的思想对于二次开发者来说非常的重要。 1)延伸网格划分 延伸网格划分是指将一个二维网格延伸生成一个三维网格;三维网格生成后去掉二维网格,延伸网格划分的步骤大体包括:先生成横截面、指定网格密度并对面进行网格划分、拖拉面网格生成体网格、指定单元属性、拖拉、完成体网格划分、释放已选的平面单元。 这里通过一个延伸网格划分的简单例子来加深对这种网格划分的理解。 图1 延伸网格划分举例 建立如图1所示的三维模型并划分网格,我们可以先建立z方向的端面,然后划分网格,通过拖拉的方法在z方向按照图中所示尺寸要求的三维模型,只需

有限元建模基本原则

?确保精度 ?控制规模 ?确保精 度: 表格1:误差分析及处理 即使采用较少的单元和较低的差值函数阶次,也能获得较满意的离散精度。例如,假设场函数在整个结构内的分布是二次函数,则用一个二次单元离散就能得到场函数的精确解。如果场函数是线性或接近于线性分布,则用线性单元离散也能得到很好的离散精度。但实际问题的场函数往往很复杂(如存在应力集中),在整个结构内很难遵循某一种函数规律,某些部位可能按高阶函数规律分布,某些部位又可能接近低阶函数的性质。故,在划网格时,结构内的不同部位可能采用不同密度和阶次的网格形式。 综上所述:提高精度的措施: 1?提高单元阶次(单元插值函数完全多项式的最高次数) 阶次越高,插值函数越能逼近复杂的真实场函数,物理离散精度越高。 其次,高阶单元的边界可以是曲线或曲面,因此在离散具有曲线或曲面边界 的结构时,几何离散误差也较线性单元小。所以当结构的场函数和形状较复杂时,可以采用这种方法来提高精度。 单元的阶次越高,收敛速度越快。 2?增加单元数量 等同于减小单元尺寸,尺寸减小时,单元的插值函数和边界能够逼近结构的 实际的场函数和实际边界,物理和几何离散误差都将减小。当模型规模不太大时, 可以采用这种方法提高精度。 但是值得注意的是:精度随着单元数量增加是有限的,当数量增加到一定程

度后,继续增加单元数量,精度却提高甚微,再采用这种方法就不经济了。实际操作时可以比较两种单元数量的计算结果,如果两次计算的差别较大,可以继续增加单元数量,否则停止增加。 3.划分规则的单元形状 单元形状的好坏将影响模型的局部精度,如果模型中存在较多的形状较差的单元,则会影响整个模型的精度。 直观上看,单元各条棱边或各个内角相差不大的形状是较好的形状。 4.建立与实际相符的边界条件 如果模型边界条件与实际工况相差较大,计算结果就会出现较大的误差,这 种误差有时甚至会超过有限元法本身带来的原理性误差。 可采用组合结构模型法,这种方法可以较好地考虑影响较大的结构间的相互作用,避免人为设置边界条件带来的误差。或采用一些测试结果,将计算值与测试值进行比较,以逐步将边界条件调整合理。 5.减少模型规模 计算误差与运算次数有关,运算次数越多,误差累计就可能越大,所以采取适当的措施降低模型规模,减少运算次数,也可能提高计算精度。 模型规模直观上可以用节点数和单元数来衡量,一般讲,节点数和单元数越多,模型规模越大,反之则越小。 在估计模型规模时,除考虑节点的多少外,还应考虑节点的自由度数,总刚度矩阵的阶次等于节点数与其自由度数的乘积,即结构的总自由度数。 减小模型规模的方法: (1)对模型进行处理:建立几何模型时,并不总是照搬结构的原有形状和尺寸,有时要做适当的简化和变换处理。合理的近似和变换可以降低模型规模,而仍然保持一定的工程精度要求。几何模型的处理方法有:降维处理、细节简化、等效变化、对称性利用和划分局部结构等。 此处很重要,参考《有限元法-原理、建模及应用》第二版.杜平安编著154 页.左下角 (2)采用子结构法:将一个复杂的结构从几何上分割为一定数量的相对简单的子结构,首先对每个子结构进行分析,然后将每个子结构的计算结果组集成整体结构的有限元模型。这种模型比直接离散结构所得到的模型要相对简单的多,从而使模型规模得到控制。这种方法适用于静力分析和动力分析。还有三种方法,不适合初级学者,待续… 看abaqus视频时了解到,对于三角形单元,一般要用二阶单元来提高精度,二阶单元会增加自由度数;但对于四边形或六面体单元,一般一阶单元已有很好的精度,不必使用二阶单元。

地质体三维建模方法与技术指南

内容简介 本书系统分析了目前国内外地质体三维模拟技术和应用软件开发的现状,由此提出了不同领域地质体 三维建模的数据需求、技术流程和主要建模软件的数据接口;详细阐述了Micmmine、surpac、Mapgis、3D-Grid等三维地质体模拟软件在矿山、地下水、城市地质等领域的应用实践和示范工作,以及提交的相 应三维模型成果;并对今后如何展开相关工作提出了建议。 本书可作为开展三维地质建模工作的指导用书,同时亦可作为地质及相关专业学生的专业参考书。 【节选】 (一)地下水三维地质建模所需数据类型 在地下水三维地质建模中,会涉及的地质现象主要有:地貌(或地形)、地层、褶 皱、断裂、透镜体及侵人体等,为刻画这些地质现象,就需要用到地表数字高程模型数据 (DEM)、遥感影像数据、地理信息数据、钻孔数据及剖面数据等。具体来说,为刻画三 维模型中的各种地质现象,需要的相关数据包括以下几种: 1.地表数字高程模型(DEM)数据 地表数学高程模型数据用于生成三维地质结构模型顶面(地表面),此部分数据可以 从测绘主管部门获取或向国家测绘局基础地理信息中心购买,从基础地理信息中心购买的 数据属于标准数据,数据以ARCINFO数据格式存放。DEM数据比例尺有多种,其中,全 国的1:25万数据库在空间上包含816幅地形图数据,覆盖整个国土范围,国外部分沿国 界外延25公里采集数据。地貌统一在TERLK层中存放,包括等高线、等深线、冲沟等, DEM等高线的等高距,在全国范围内共分40 m、50 m、100 m三种,使用时可参照等分 布图确定。对于标准数据,可以根据需要进行数据格式转换、比例变换、投影变换等多种 处理。 另外,如果不能获取现成的DEM数据,也可以自己使用专门的地理信息系统软件用 地形图生产。即把纸质地形图数字化及几何纠正校准,然后进行高程信息的提取——对等 高线进行屏幕矢量跟踪并对等高线标赋高程值,同时编辑、检查、拼接以生成各种拓扑关 系,最后用软件进行内插值、裁剪生成DEM数据。 2.遥感影像数据

ansys有限元建模与分析实例-详细步骤

《有限元法及其应用》课程作业ANSYS应用分析 学号: 姓名: 专业:建筑与土木工程

角托架的有限元建模与分析 一 、模型介绍 本模型是关于一个角托架的简单加载,线性静态结构分析问题,托架的具体形状和尺寸如图所示。托架左上方的销孔被焊接完全固定,其右下角的销孔受到锥形压力载荷,角托架材料为Q235A 优质钢。角托架材料参数为:弹性模量366E e psi =;泊松比0.27ν= 托架图(厚度:0.5) 二、问题分析 因为角托架在Z 方向尺寸相对于其在X,Y 方向的尺寸来说很小,并且压力荷载仅作用在X,Y 平面上,因此可以认为这个分析为平面应力状态。 三、模型建立 3.1 指定工作文件名和分析标题 (1)选择菜单栏Utility Menu → 命令.系统将弹出Jobname(修改文件名)对话框,输入bracket (2)定义分析标题 GUI :Utility Menu>Preprocess>Element Type>Add/Edit/Delete 执行命令后,弹出对话框,输入stress in a bracket 作为ANSYS 图形显示时的标题。 3.2设置计算类型 Main Menu: Preferences … →select Structural → OK 3.3定义单元类型 PLANE82 GUI :Main Menu →Preprocessor →Element Type →Add/Edit/Delete 命令,系统将弹出Element Types 对话框。单击Add 按钮,在对话框左边的下拉列表中单击Structural Solid →Quad 8node 82,选择8节点平面单元PLANE82。单击ok ,Element Types 对话框,单击Option ,在Element behavior 后面窗口中选取Plane strs w/thk 后单击ok 完成定义单元类型。 3.4定义单元实常数 GUI :Main Menu: Preprocessor →Real Constants →Add/Edit/Delete ,弹出定义实常数对话框,单击Add ,弹出要定义实常数单元对话框,选中PLANE82单元后,单击OK →定义单元厚度对话框,在THK 中输入0.5.

__地质建模型操作步骤62.

地质建模培训操作步骤 Gemcom国际矿业软件公司中国办事处 目录 第1 章地质数据库结构 (1 1.1 分析实例地质数据 (1

第2 章创建一个Surpac地质数据库 (5 2.1 创建数据库结构 (5 2.2 编辑转换表 (8 2.3 导入数据 (9 第3 章三维空间显示钻孔 (14 3.1 设置钻孔显示风格 (14 3.2 显示钻孔 (18 第4 章地质解译 (22 4.1 创建剖面 (22 4.2 地质解译(圈定矿体 (25 第5 章形成实体模型 (38 5.1 按矿体保存线文件 (38 5.2 形成实体模型 (39 第6 章数据提取 (50 6.1 提取钻孔平面图/剖面图 (50 6.2 提取样品表中数据,获得取样长度 (52 6.3 统计样长,确定组合样长 (54 第7 章组合 (57 7.1 样品组合原理 (57

7.2 根据勘探工程组合 (58 第8 章基本统计与分析 (66 8.1 基本统计 (66 8.2 处理特高品位 (71 第9 章创建块体模型 (73 9.1 创建块模型 (73 9.2 增加属性及背景值 (77 9.3 约束 (79 9.4 块体模型的显示 (80 第10 章块体模型赋值 (84 10.1 直接赋值法 (84 10.2 QV1矿体的距离幂次反比法 (87 10.2.1 QV1第一次估值 (87 10.2.2 根据属性为模型着色 (91 10.2.3 赋QV1矿体第一次估值次数 (93 10.2.4 QV1矿体第二次估值 (94 10.2.5 赋QV1矿体第二次估值次数 (97 10.2.6 判断QV1矿体内的块是否已经全部估值完成 (99 10.3 QV2矿体的距离幂次反比法 (99

有限元分析步骤

有限元建模与分析 有限元分析(FEA)是一种预测结构的偏移与其它应力影响的过程,有限元建模(FEM)将这个结构分割成单元网格以形成实际结构的模型,每个单元具有简单形态(如正方形或三角形)。这样有限元程序就有了可写出在刚度矩阵结构中控制方程方面的信息。每个单元上的未知量就是在节点上的位移,这个点就是单元元的连接点。有限元程序将这些单个单元的刚度矩阵组合起来以形成整个模型的总刚度矩阵,并给予已知力和边界条件来求解该刚度矩阵以得出未知位移,从节点上位移的变化就可以计算出每个单元中的应力。 有限单元由假定的应变方程式导出,有些单元可假设其应变是常量,而另外一些可采用更高阶的函数。利用给定单元的这些方程和实际几何体,则可以写出外力和节点位移之间的平衡方程。对于单元的每个节点来说,每个自由度就有一个方程,这些方程被十分便利地写成矩阵的形式以用于计算机的演算中,这个系数的矩阵就变成了一个显示出力对位移的关系的刚度矩阵:{F}=[K]、{d} 尽管求知量处于离散的自由度,内部方程仍被写成表述为连续集的应变函数。这就意味着如果选择了正确单元的话,纵然这个有限元模型有一组离散的方程,只要用有限的节点和单元也可以收敛出正确的答案。 有限元模型是解决全部结构问题的完全理想的模型。这些问题包括节点的定位,单元,物理的和材料的特性,载荷和边界条件,根据分析类型的不同,如静态结构载荷,动态的或热力分析,这个模型就确定得不同。 一个有限元模型常常由不止一种单元类型来建立,有限元模型是以结构的偏移来建立成数学模型,而不只是在外观上象原结构。也许某个零件用梁单元最好,而另外的零件则可能用薄壳单元最理想。 对于给定的问题来讲,求解结果的准确性将取决于结构建模的好坏,负载和边界条件的确定,以及所用单元的精度。 一般来讲,如模型细分更小的单元,则求解将更准确。了解你在最终的求解结果上有充分收敛的唯一确信的方法是用更细网格的单元来建立更多的模型,以检查求解结果的收敛性。 新的有限元用户经常产生想象上的错误,即建立一个有限元模型的目的是建立一个看起来象这种结构的模型。有限元建模的目的是建立一个从数学意义是“相似”的模型,而不是一个外观相似的模型。一个有经验的使用者学会了怎样选择单元的正确类型,和在模型的不同区域中怎样来细分网格。 一个经常忽略的错误根源是在一个模型中的负载和边界条件上进行了错误的假设。同时也很轻易地相信一个有限元模型的每个十进位的结果。以及忘掉了在负载和边界条件上粗糙的假设。如果有一个关于怎样建立边界条件模型的问题的话,宁可用你的模型以不同的方法去测试其灵敏度,而不是仅遵循一种方法,得出一种答案,

abaqus有限元建模小例子

问题一: 工字梁弯曲 1.1 问题描述: 在<<材料力学实验>>中,弯曲实验測定了工字梁弯曲应变大小及其分布,以验证弯曲正应力公式。在这里,採用ABAQUS/CAE建立试验件的有限元模型,ABAQUS/Standard模块进行分析求解,得到应力、应变分布,对比其与理论公式计算值及实验測量值的差別。 弯曲实验的相关数据: 材料:铝合金E=70GPa 泊松比0.3 实验装置结构简图如图所示: 结构尺寸测量值:H=50(+/-0.5mm) h=46(+/-0.5mm) B=40(+/-0.5mm) b=2(+/-0.02mm) a=300(+/-1mm) F1=30N Fmax=300N N ? F100 = 1.2 ABAQUS有限元建模及分析 一对象: 工字型截面铝合金梁 梁的结构简图如图1所示,結构尺寸、载荷、約束根据1.1设定,L取1600mm,两端各伸出100mm。 二用ABAQUS/CAE建立实验件的有限元模型,效果图如下: 边界条件简化: 左侧固定铰支座简化为下表面左参考点处的约束U1=U2=U3=0

右侧活动铰支座简化为下表面右参考点处的约束U1=U2=UR3=0 几何模型

有限元模型 三ABAQUS有限元分析結果 ①应力云图(Z方向正应力分量):施加载荷前 F=300N

②应变(Z方向分量): 中间竖直平面的厚度方向应变分布图: F=100N F=200N

F=300N 由上图可以看出应变沿着厚度方向呈线性比例趋势变化,与实验测得的应变值变化趋势相同。中性轴处应变均接近零值,应变与距离中性轴位移基本为正比关系。 1.3分析结果: 中间竖直截面上下边缘轴向应力数值对比:*10^-6 MPa 距中性轴距ABAQUS模拟实验测量值平均理论值 1/2H -96.182*70000 -97*70000 -6.9165=-70000*98.807 -1/2H 95.789*70000 92*70000 6.9165

petrel地质建模软件的饱和度计算方法

Sensitivity on imported Saturation Functions Background This small workflow creates a new simulation run for each saturation function set that are in the project. It can use both sets created in Petrel and sets that has been generated outside Petrel and imported. Result Produces several simulation runs, one for each set of saturation functions. Prerequisites This workflow requires an ECLIPSE and/or Frontsim model that runs and has been built in Petrel, and a set of different saturation functions. The workflow creates a new folder on the cases tree and puts the simulation cases in there once they are run. It then loops over all the saturation functions, and creates and run an ELCIPSE and/or Frontsim case for each saturation function set. Having different saturation functions will change the dynamic simulation results. From the workflow editor the Variable A is dropped into the Saturation functions table, and the case is exported using all the different sets in the saturation functions folder. Workflow Specifics Workflow Specifics Petrel Version 2009.1 Domain Reservoir Engineering Modules Reservoir Engineering Core or Combined Core (optional Frontsim Locked) Expert Level Advanced

相关文档
最新文档