板料成形有限元模拟中自适应三角形网格的研究

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

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

板料成形CAE技术与其应用

板料成形CAE技术及应用 长期以来，困扰广大模具设计人员的主要问题就是较长的模具开发设计周期，特别是对于某些特殊复杂的板料成形零件，甚至制约了整个产品的开发进度，而板料成形CAE技术及分析软件的出现，有效地缩短模具设计周期，大大减少试模时间，帮助企业改进产品质量，降低生产成本，从根本上提高了企业的市场竞争力。 一、前言 计算机辅助设计技术以其强大的冲击力，影响和改变着工业的各个方面，甚至影响着社会的各个方面。它使传统的产品技术、工程技术发生了深刻的变革，极大地提高了产品质量，缩短了从设计到生产的周期，实现了设计的自动化。 板料成形是利用模具对金属板料的冲压加工，获得质量轻、表面光滑、造型美观的冲压件，具有节省材料、效率高和低成本等优点，在汽车、航空、模具等行业中占据着重要地位。由于板料成形是利用板材的变形得到所需的形状的，长期以来，困扰广大模具设计人员的主要问题就是较长的模具开发设计周期，特别是对于复杂的板料成形零件无法准确预测成形的结果，难以预防缺陷的产生，只能通过经验或类似零件的现有工艺资料，通过不断的试模、修模，才能成功。某些特殊复杂的板料成形零件甚至制约了整个产品的开发进度。 板料成形CAE技术及分析软件，可以在产品原型设计阶段进行工件坯料形状预示、产品可成形性分析以及工艺技术方案优化，从而有效地缩短模具设计周期，大大减少试模时间，帮助企业改进产品质量，降低生产成本，从根本上提高企业的市场竞争力。 板料成形CAE技术对传统开发模式的改进作用可以通过图1 和图2进行对比。

图1 传统板料成形模具开发模式 图2 CAE 技术模具开发方式 通过比较，就可发现板料成形CAE技术的主要优点。 (1)通过对工件的可成形工艺性分析，做出工件是否可制造的早期判断；通过对模具技术方案和冲压技术方案的模拟分析，及时调整修改模具结构，减少实际试模次数，缩短开发周期。 (2)通过缺陷预测来制定缺陷预防措施，改进产品设计和模具设计，增强模具结构设计以及冲压技术方案的可靠性，从而减少生产成本。 (3)通过CAE分析可以择优选择材料，可制造复杂的零件，并对各种成形参数进行优化,提高产品质量。 (4)通过CAE分析应用不仅可以弥补工艺人员在经验和应用工艺资料方面的不足，还可通过虚拟的冲压模拟，提高提高工艺人员的经验。 二、板料成形需要解决的问题 板料成形通过模具对板料施加压力，使板料产生永久性的塑性变形，以获得预期的产品形状。在这个过程中影响板材变形的因素非常多，要控制好变形的形状也非常困难。首先，金属受外力作用会发生变形，变形可分为弹性变形和塑性变形，弹性变形是可逆的，外力去除后变形体就会恢复成原来的形状；第二，材料的成分和组织对变形影响极大；第三，塑性变形有多种方式，再结晶温度下的塑性变形有晶内滑移和孪动、位错（位错分多种形式），再结晶温度上的塑性变形有晶间滑移、多晶体扩散和相变变形等；第四，变形温度、变形速度的影响；第五，变形体内部应力状态的影响；第六，摩擦与润滑的影响；第七，材料塑性变形后，当变形体内部各部分变形不一致时，

板料成形回弹问题研究新进展_朱东波

第7卷第1期2000年3月 塑性工程学报 JOU RN AL O F PLASTICITY EN GIN EERIN G V ol.7　No.1Ma r .　2000 板料成形回弹问题研究新进展 * (西安交通大学先进制造技术研究所　710049) 　 朱东波　孙　琨李涤尘　卢秉恒 摘　要:本文从回弹理论、回弹数值模拟分析、回弹控制三方面对弯曲成形、3-D 复杂浅拉深成形中回弹研究的历史和最新发展状况作了较全面的介绍。文章所引用的大量文献基本概括了前人在这些方面的主要研究方法和重要研究成果。 关键词:回弹;板料成形;模具 *国家“九·五”重点攻关资助项目(项目号: 85-951-19)。收稿日期: 1999-4-28 1　引　言 板料成形过程中普遍存在有回弹问题,特别在弯曲和浅拉深过程中回弹现象更为严重,对零件的尺寸精度和生产效率造成极大的影响,有必要对其进行深入的研究和有效的控制。零件的最后回弹形状是其整个成形历史的累积效应,而板料成形过程与模具几何形状、材料特性、摩擦接触等众多因素密切相关,所以板料成形的回弹问题非常复杂。半个多世纪来国内外许多学者对回弹问题进行了深入的研究和探讨,这些研究涵盖了从弯曲成形到复杂拉深成形、从理论分析到数值模拟、从回弹预测到回弹控制等诸多方面。本文从三个方面对前人的工作进行了概括性回顾,重点介绍了90年代回弹研究的一些新进展。 2　弯曲理论研究和回弹的解析分析方法 弯曲成形一般只涉及较为简单的几何形状和边界条件,所以有条件用解析方法对其进行深入的研究。50年代,R .H ill 、F .Proska 、F .J .Gardiner 等人的工作奠定了板料弯曲及回弹分析的理论基础 [1] ,后 来不断有学者对这些理论进行深化和发展。Huang ,etc [2] 在其文章中对50年代到80年代间诸多学者的 回弹研究工作做了较详细的回顾和评述。 回弹是弯曲卸载过程产生的反向弹性变形,板料回弹的经典计算公式为: Δk =1R -1R S =12M (1-ν2 ) Et 3 (1) 式中　Δk ——曲率变化量 R ——回弹前中面半径 R S ——回弹后中面半径E ——弹性模量ν ——泊松比t ——回弹前板料厚度 M ——回弹前板内弯矩 弯矩M 由截面纵向应力分布唯一确定。对同一弯曲过程,采用不同的弯曲模型(如是否考虑中性面内移,是否考虑材料强化、各向异性等)可得到不同的应力分布,从而由式(1)得到回弹量Δk 也就不同。所以在理论分析中,弯曲模型是否合理将直接影响回弹计算结果的准确程度。 弯曲的基本理论模型分为两大类。一类是以平截面假定和单向应力假定为基础的工程理论模型,该模型未考虑径向应力,认为弯曲过程中应力中性层、应变中性层始终和几何中面相重合;另一类是由H ill [3] 首先提出的精确理论模型,该模型考虑径向应力及中性层内移的影响,更接近板料弯曲的真实情况。从板料的外部受力状态和加载方式来看,弯曲过程可分为纯弯曲、拉伸弯曲、循环弯曲等几种典型情况。另外,材料模型对弯曲计算结果有很大的影响,常用的材料模型有刚塑性、理想弹塑性、刚性强化、弹性强化等多种形式。 以上基本模型、加载方式及材料模型的不同组合

有限元网格技巧

有限元网格技巧 1. 引言 有限元法是求解复杂工程问题的一种近似数值解法，现已广泛应用到力学、热学、电磁学等各个学科，主要分析工作环境下物体的线性和非线性静动态特性等性能。 有限元法求解问题的基本过程主要包括：分析对象的离散化?有限元求解?计算结果的处理三部分。 曾经有人做过统计：三个阶段所用的时间分别占总时间的40%~50%、5%及50%~55%。也就是说，当利用有限元分析对象时，主要时间是用于对象的离散及结果的处理。如果采用人工方法离散对象和处理计算结果，势必费力、费时且极易出错，尤其当分析模型复杂时，采用人工方法甚至很难进行，这将严重影响高级有限元分析程序的推广和使用。因此，开展自动离散对象及结果的计算机可视化显示的研究是一项重要而紧迫的任务。 可喜的是，随着计算机及计算技术的飞速发展，出现了开发对象的自动离散及有限元分析结果的计算机可视化显示的热潮，使有限元分析的“瓶颈”现象得以逐步解决，对象的离散从手工到半自动到全自动，从简单对象的单维单一网格到复杂对象的多维多种网格单元，从单材料到多种材料，从单纯的离散到自适应离散，从对象的性能校核到自动自适应动态设计/分析，这些重大发展使有限元分析摆脱了仅为性能校核工具的原始阶段，计算结果的计算机可视化显示从简单的应力、位移和温度等场的静动态显示、彩色调色显示一跃成为对受载对象可能出现缺陷(裂纹等)的位置、形状、大小及其可能波及区域的显示等，这种从抽象数据到计算机形象化显示的飞跃是现在甚至将来计算机集成设计/分析的重要组成部分。 2. 有限元分析对网格剖分的要求 有限元网格生成就是将工作环境下的物体离散成简单单元的过程，常用的简单单元包括：一维杆元及集中质量元、二维三角形、四边形元和三维四面体元、五面体元和六面体元。他们的边界形状主要有直线型、曲线型和曲面型。对于边界为曲线(面)型的单元，有限元分析要求各边或面上有若干点，这样，既可保证单元的形状，同时，又可提高求解精度、准确性及加快收敛速度。不同维数的同一物体可以剖分为由多种单元混合而成的网格。网格剖分应满足以下要求：合法性。一个单元的结点不能落入其他单元内部，在单元边界上的结点均应作为单元的结点，不可丢弃。相容性。单元必须落在待分区域内部，不能落入外部，且单元并集等于待分区域。逼近精确性。待分区域的顶点(包括特殊点)必须是单元的结点，待分区域的边界(包括特殊边及面)被单元边界所逼近。良好的单元形状。单元最佳形状是正多边形或正多面体。良好的剖分过渡性。单元之间过渡应相对平稳，否则，将影响计算结果的准确性甚至使有限元计算无法计算下去。网格剖分的自适应性。在几何尖角处、应力温度等变化大处网格应密，其他部位应较稀疏，这样可保证计算解精确可靠。 3. 现有有限元网格剖分方法 K. Ho-Le 对网格生成算法进行了系统分类，该分类方法可沿用至今，它们是拓扑分解法、结点连元法、网格模板法、映射法和几何分解法五种。目前，主要是上述方法的混合使用及现代技术的综合应用。

_基于ANSYS的有限元法网格划分浅析

文章编号:1003-0794(2005)01-0038-02 基于ANSYS的有限元法网格划分浅析 杨小兰,刘极峰,陈 旋 (南京工程学院,南京210013) 摘要:为提高有限元数值的计算精度和对复杂结构力学分析的准确性,针对不同分析类型采用了不同的网格划分方法,结合实例阐述了ANSYS有限元网格划分的方法和技巧,指出了采用ANSYS有限元软件在网格划分时应注意的技术问题。 关键词:ANSYS;有限元;网格;计算精度 中图号:O241 82;TP391 7文献标识码:A 1 引言 ANSYS有限元分析程序是著名的C AE供应商美国ANSYS公司的产品,主要用于结构、热、流体和电磁四大物理场独立或耦合分析的CAE应用,功能强大,应用广泛,是一个便于学习和使用的优秀有限元分析程序。在ANSYS得到广泛应用的同时,许多技术人员对ANSYS程序的了解和认识还不够系统全面,在工作和研究中存在许多隐患和障碍,尤为突出的是有限元网格划分技术。本文结合工程实例,就如何合理地进行网格划分作一浅析。 2 网格划分对有限元法求解的影响 有限元法的基本思想是把复杂的形体拆分为若干个形状简单的单元,利用单元节点变量对单元内部变量进行插值来实现对总体结构的分析,将连续体进行离散化即称网格划分,离散而成的有限元集合将替代原来的弹性连续体,所有的计算分析都将在这个模型上进行。因此,网格划分将关系到有限元分析的规模、速度和精度以及计算的成败。实验表明:随着网格数量的增加,计算精确度逐渐提高,计算时间增加不多;但当网格数量增加到一定程度后,再继续增加网格数量,计算精确度提高甚微,而计算时间却大大增加。在进行网格划分时,应注意网格划分的有效性和合理性。 3 网格划分的有效性和合理性 (1)根据分析数据的类型选择合理的网格划分数量 在决定网格数量时应考虑分析数据的类型。在静力分析时,如果仅仅是计算结构的变形,网格数量可以少一些。如果需要计算应力,则在精度要求相同的情况下取相对较多的网格。同样在响应计算中,计算应力响应所取的网格数应比计算位移响应多。在计算结构固有动力特性时,若仅仅是计算少数低阶模态,可以选择较少的网格。如果计算的模态阶次较高,则应选择较多的网格。在热分析中,结构内部的温度梯度不大,不需要大量的内部单元,可划分较少的网格。 (2)根据分析数据的分布特点选择合理的网格疏密度 在决定网格疏密度时应考虑计算数据的分布特点,在计算固有特性时,因为固有频率和振型主要取决于结构质量分布和刚度分布,采用均匀网格可使结构刚度矩阵和质量矩阵的元素不致相差很大,可减小数值计算误差。同样,在结构温度场计算中也趋于采用均匀的网格形式。在计算数据变化梯度较大的部位时,为了更好地反映数据变化规律,需要采用比较密集的网格,而在计算数据变化梯度较小的部位,为了减小模型规模,则应划分相对稀疏的网格,这样整个结构就表现出疏密不同的网格划分形式。 以齿轮轮齿的有限元分析模型为例,由于分析的目的是求出齿轮啮合传动过程中齿根部分的弯曲应力,因此,分析计算时并不需要对整个齿轮进行计算,可根据圣文男原理将整个区域缩小到直接参与啮合的轮齿。虽然实际上参与啮合的齿数总大于1,但考虑到真正起作用的是单齿,通常只取一个轮齿作为分析对象,这样作可以大大节省计算机内存。考虑到轮齿应力在齿根过渡圆角和靠近齿面处变化较大,网格可划分得密一些。在进行疏密不同网格划分操作时可采用ANSYS提供的网格细化工具调整网格的疏密,也可采用分块建模法设置网格疏密度。 图1所示即为采用分块建模法进行网格划分。图1(a)为内燃机中重要运动零件连杆的有限元应力分析图,由于连杆结构对称于其摆动的中间平面,其厚度方向的尺寸远小于长度方向的尺寸,且载荷沿厚度方向近似均匀分布,故可按平面应力分析处 38 煤 矿 机 械 2005年第1期

板料成形CAE分析

板料成形CAE分析 实验报告 班级： 学号： 姓名：

板料成形CAE分析 一、实验目的和要求： 通过本实验的教学，使学生基本掌握有限元技术在板料塑性成形领域的应用情况，拓宽学生的知识面，开阔视野，使学生对塑性成形过程的数值模拟技术有深刻的理解，预测板料弯曲成形的性能。 二、教学基本要求： 学会使用Dynaform数值模拟软件进行板料弯曲成形过程的仿真模拟，对模拟结果具有一定的分析和处理能力。 三、实验内容提要： 掌握前处理的关键参数设置，如零件定义、网格划分、模型检查、工具定义、坯料定义、工具定位和移动、工具动画、运行分析。了解后处理模块对模拟结果的分析，如读入d3plot 文件、动画显示变形和生成动画文件、成形极限图分析、坯料厚度变化分析等。 四、实验步骤 1、导入零件模型，保存文件 打开下拉菜单File->Import，如图2所示，在F:\dynaform\BLANK_CAE目录下分别导入文件punch.igs,binder.igs,die.igs和blank.igs。 图1 导入文件窗口

3、更改零件层名 打开下拉菜单Parts－>Edit,对应不同的零件更改层名，改好层名后保存文件。 图2 修改层名窗口 4、进行网格划分 以blinder为例进行说明。 （1）、点击，只选择binder1（红色），点击OK退出。 图3

（2）、选择Preprocess—>Element进入如图3界面。选择，在surf mesh中将max size 改为5. 图4 图5 （3）、依次选select surfaces—>displayed surf-->0k-->apply，然后依次退出各个页面。网格化后的零件如图6所示。 图6网格化后的零件 4、检查工具。 仍然以binder为例。

CATIA有限元高级划分网格教程

CATIA有限元高级网格划分教程 盛选禹李明志 1.1进入高级网格划分工作台 （1）打开例题中的文件Sample01.CATPart。 （2）点击主菜单中的【开始】→【分析与模拟】→【Advanced Meshing Tools】(高级网格划分工具），就进入【Advanced Meshing Tools】（高级网格划分工具）工作台，如图1－1所示。进入工作台后，生成一个新的分析文件，并且显示一个【New Analysis Case】（新分析算题）对话框，如图1－2所示。 图1－1【开始】→【分析与模拟】→【Advanced Meshing Tools】(高级网格划分工具）（3）在【New Analysis Case】（新分析算题）对话框内选择【Static Analysis】（静力分析）选项。如果以后打开该对话框的时候均希望是计算静力分析，可以把对话框内的【Keep as default starting analysis case】（在开始时保持为默认选项）勾选。这样，下次进入本工作台时，将自动选择静力分析。 （4）点击【新分析算题】对话框内的【确定】按钮，关闭对话框。 1.2定义曲面网格划分参数 本节说明如何定义一个曲面零件的网格类型和全局参数。 （1）点击【Meshing Method】(网格划分方法)工具栏内的【高级曲面划分】按钮

，如图1－3所示。需要在【Meshing Method】(网格划分方法)工具栏内点击中间按钮的下拉箭头才能够显示出【高级曲 面划分】按钮。 图1－2【New Analysis Case】（新分析算题）对话框图1－3【高级曲面划分】按钮

板料冲压件螺纹底孔冲压成形技术

板料冲压件螺纹底孔冲压成形技术 摘要：在板料冲压件上，按其料厚不同分别采用精冲小孔、变薄翻边、冷冲挤等工艺方法，成形螺纹底孔。本文论述了上述螺纹冲压成形工艺、冲模结构及其设计与制造技术。 主题词：冲件螺纹底孔冲小孔变薄翻边冷冲挤成形技术 螺纹联接结构，尤其紧螺纹联接结构，是各种机电与家电产品中零部件最主要的联接结构型式。薄板冲压件进行紧螺纹联接，需要有大于料厚的联接螺纹长度，以确保其联接可靠性，增强其负载能力，才能达到使薄板冲件联接牢靠、重量小的目的，从而使其成为结实、轻巧、紧凑的理想结构零件。 在仪器仪表、电子电器、各类家电、家用器具、玩具等产品的板料冲压件上，经常采用M2-M10的小螺纹紧联接结构。为提高效率并满足大量生产的需求，采用精冲小孔、变薄翻边、冷冲挤等工艺方法，冲压成形这些小螺纹底孔，不仅能以冲压制孔取代钻孔而大幅度提高生产效率，同时能获得尺寸精确、一致性好的底孔，并可使螺纹联接有足够的长度，从而确保其联接可靠性及设计要求的承载能力。所以，用冲压成形技术加工小螺纹底孔，具有优质高产的效果，也是一种成熟而值得推广的工艺技术。 1 螺纹底孔的计算 合适螺纹底孔的大小，不仅取决于螺纹直径，而且与其螺距有着密切的关系，通常可按下式计算： 当t L≤1时，取：d Z=d-t L

当t L>1时，取：d Z=d-～t L (2) 式中 t L-螺距，mm d z-螺纹底孔直径，mm d-螺纹直径，mm 表1 螺纹底孔直径的合理值(mm) 螺纹直径d 螺 距 t L 底 孔 直 径d z M1 M2 M3 M4 M5 M6 M8 M10 M12 1 5

M14 M16 M18 M20 M22 M24 M27 M302 2 3 3 2 冲制螺纹底孔的基本工艺方法 用冷冲压冲制板料冲压件上螺纹底孔的主要工艺方法有如下几种： (1)厚料冲小孔与精冲孔 当冲件厚t可以满足螺纹联接所需长度时，可用冲压制孔工艺解决。通常在这种情况下，多为厚料冲小孔，即冲制螺纹底孔的直径dz≤t或稍大于t，见表2。螺纹联接的最小有效长度取决于螺纹直径、螺距并与联接件的材料种类密切相关。

有限元网格划分

有限元网格划分 摘要:总结近十年有限元网格划分技术发展状况。首先,研究和分析有限元网格划分的基本原则；其次,对当前典型网格划分方法进行科学地分类,结合实例,系统地分析各种网格划分方法的机理、特点及其适用范围,如映射法、基于栅格法、节点连元法、拓扑分解法、几何分解法和扫描法等；再次,阐述当前网格划分的研究热点,综述六面体网格和曲面网格划分技术；最后,展望有限元网格划分的发展趋势。 关键词:有限元网格划分；映射法；节点连元法;拓扑分解法；几何分解法；扫描法；六面体网格 1 引言 有限元网格划分是进行有限元数值模拟分析至关重要的一步，它直接影响着后续数值计算分析结果的精确性。网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。在有限元数值求解中，单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成，连续体单元以及壳、板、梁单元的面内均采用高斯（Gauss）积分，而壳、板、梁单元的厚度方向采用辛普生（Simpson）积分。 2 有限元网格划分的基本原则 有限元方法的基本思想是将结构离散化,即对连续体进行离散化,利用简化几何单元来近似逼近连续体,然后根据变形协调条件综合求解。所以有限元网格的划分一方面要考虑对各物体几何形状的准确描述,另一方面也要考虑变形梯度的准确描述。为正确、合理地建立有限元模型,这里介绍划分网格时应考虑的一些基本原则。 2.1 网格数量 网格数量直接影响计算精度和计算时耗,网格数量增加会提高计

算精度,但同时计算时耗也会增加。当网格数量较少时增加网格，计算精度可明显提高,但计算时耗不会有明显增加;当网格数量增加到一定程度后,再继续增加网格时精度提高就很小,而计算时耗却大幅度增加。所以在确定网格数量时应权衡这两个因素综合考虑。 2.2 网格密度 为了适应应力等计算数据的分布特点,在结构不同部位需要采用大小不同的网格。在孔的附近有集中应力,因此网格需要加密;周边应力梯度相对较小,网格划分较稀。由此反映了疏密不同的网格划分原则:在计算数据变化梯度较大的部位,为了较好地反映数据变化规律,需要采用比较密集的网格;而在计算数据变化梯度较小的部位,为减小模型规模,网格则应相对稀疏。 2.3 单元阶次 单元阶次与有限元的计算精度有着密切的关联,单元一般具有线性、二次和三次等形式,其中二次和三次形式的单元称为高阶单元。高阶单元的曲线或曲面边界能够更好地逼近结构的曲线和曲面边界,且高次插值函数可更高精度地逼近复杂场函数,所以增加单元阶次可提高计算精度。但增加单元阶次的同时网格的节点数也会随之增加,在网格数量相同的情况下由高阶单元组成的模型规模相对较大,因此在使用时应权衡考虑计算精度和时耗。 2.4 单元形状 网格单元形状的好坏对计算精度有着很大的影响,单元形状太差的网格甚至会中止计算。单元形状评价一般有以下几个指标: (1)单元的边长比、面积比或体积比以正三角形、正四面体、正六面体为参考基准。 (2)扭曲度:单元面内的扭转和面外的翘曲程度。 (3)节点编号:节点编号对于求解过程中总刚矩阵的带宽和波前因数有较大的影响,从而影响计算时耗和存储容量的大小 2.5 单元协调性 单元协调是指单元上的力和力矩能够通过节点传递给相邻单元。为保证单元协调,必须满足的条件是: (1)一个单元的节点必须同时也是相邻点,而不应是内点或边界

带孔板的建模及有限元分析Word版

基于SolidWorks带孔板的建模及有限元分析 李军 摘要：利用SolidWorks对带孔矩形板进行虚拟建模，通过赋予板材材质、载荷后进行网格划分，进而进行有限元分析，得出其应力、应变和位移的分布图，并对结果进行分析研究对板材安全性的影响。 关键词：SolidWorks；带孔板；建模；有限元分析 0 SolidWorks简介 Solidworks是一款优秀的三维设计软件，具有十分强大的零件设计功能及装配模块，同时也拥有丰富的后置处理模块。由于其功能强大，新手上手快，应用领域广，所以成为了主流的三维造型软件。经过17年的发展，在全球已经拥有30多万的客户，最新版本为SolidWorks 2011版。在中国SolidWorks在计算机辅助设计、计算机辅助工程、计算机辅助制造、计算机辅助工艺、数据管理等方面为企业提供了强大的动力，使企业在管理、设计和制造方面有了很大的提升。 1 带孔板的模型建立 矩形板材的尺寸为300*180*10mm，孔位于中心，直径为50mm，模型如图1。 图1 带孔矩形板模型 2前置处理 2.1在Command Manager中点击SIMULATION选项，建立新算例，名称默认，确认。 2.2赋予板材材料属性 材料为AISI304，材料属性如表1

表1 材料的属性 模型参考属性零部件 名称:AISI 304 模型类型:线性弹性同向性 默认失败准则:最大von Mises 应力屈服强度: 2.06807e+008 N/m^2 张力强度: 5.17017e+008 N/m^2 弹性模量: 1.9e+011 N/m^2 泊松比:0.29 质量密度:8000 kg/m^3 抗剪模量:7.5e+010 N/m^2 热扩张系数: 1.8e-005 /Kelvin SolidBody 1(凸台-拉伸1)(aisi304带孔矩形钢板静力分析) 曲线数据:N/A 2.3网格生成 在SIMULATION选项中选择“运行”中的“生成网格”，使用默认网格划分。网格 信息如表2，网格信息细节如表3，网格划分后的模型如图2。 表2 网格信息 网格类型实体网格 所用网格器: 基于曲率的网格 雅可比点 4 点 最大单元大小7.44196 mm 最小单元大小7.44196 mm 网格品质高 表3 网格信息细节 节点总数23523 单元总数13612 最大高宽比例 3.9347 单元(%)，其高宽比例< 3 99.7 单元(%)，其高宽比例> 10 0 扭曲单元(雅可比)的% 0 完成网格的时间(时;分;秒): 00:00:03 计算机名: PC-201009062016

板料成形CAE技术

板料成形CAE技术 贵州风华机器厂童春桥 一、前言 计算机辅助设计技术以其强大的冲击力，影响和改变着工业的各个方面，甚至影响着社会的各个方面。它使传统的产品技术、工程技术发生了深刻的变革，极大地提高了产品质量，缩短了从设计到生产的周期，实现了设计的自动化。 板料成形是利用模具对金属板料的冲压加工，获得质量轻、表面光滑、造型美观的冲压件，具有节省材料、效率高和低成本等优点，在汽车、航空、模具等行业中占据着重要地位。由于板料成形是利用板材的变形得到所需的形状的，长期以来，困扰广大模具设计人员的主要问题就是较长的模具开发设计周期，特别是对于复杂的板料成形零件无法准确预测成形的结果，难以预防缺陷的产生，只能通过经验或类似零件的现有工艺资料，通过不断的试模、修模，才能成功。某些特殊复杂的板料成形零件甚至制约了整个产品的开发进度。 板料成形CAE技术及分析软件，可以在产品原型设计阶段进行工件坯料形状预示、产品可成形性分析以及工艺方案优化，从而有效地缩短模具设计周期，大大减少试模时间，帮助企业改进产品质量，降低生产成本，从根本上提高企业的市场竞争力。 板料成形CAE技术对传统开发模式的改进作用可以通过图1 和图2进行对比。

图1 传统板料成形模具开发模式 图2 CAE 技术模具开发方式 通过比较，就可发现板料成形CAE技术的主要优点。 (1)通过对工件的可成形工艺性分析，做出工件是否可制造的早期判断；通过对模具方案和冲压方案的模拟分析，及时调整修改模具结构，减少实际试模次数，缩短开发周期。 (2)通过缺陷预测来制定缺陷预防措施，改进产品设计和模具设计，增强模具结构设计以及冲压方案的可靠性，从而减少生产成本。 (3)通过CAE分析可以择优选择材料，可制造复杂的零件，并对各种成形参数进行优化,提高产品质量。 (4)通过CAE分析应用不仅可以弥补工艺人员在经验和应用工艺资料方面的不足，还可通过虚拟的冲压模拟，提高提高工艺人员的经验。 二、板料成形需要解决的问题 板料成形通过模具对板料施加压力，使板料产生永久性的塑性变形，以获得预期的产品形状。在这个过程中影响板材变形的因素非常多，要控制好变形的形状也非常困难。首先，金属受外力作用会发生变形，变形可分为弹性变形和塑性变形，弹性变形是可逆的，外力去除后变形体就会恢复成原来的形状；第二，材料的成分和组织对变形影响极大；第三，塑性变形有多种方式，再结晶温度下的塑性变形有晶内滑移和孪动、位错(位错分多种形式)，再结晶温度上的塑性变形有晶间滑移、多晶体扩散和相变变形等；第四，变形温度、变形速度的影响；第五，变形体内部应力状态的影响；第六，摩擦与润滑的影响；第七，材料塑

板料成形回弹特征及其控制技术

板料成形回弹特征及其控制技术 1 前言 回弹是板材冲压成形过程的主要缺陷之一．严重影响着威形件的威形质量和尺寸精度，是实际工艺中很难有效克服的成形缺陷之一，它不仅降低了产品质量和生产效率．还制约了自动化装配生产线的实施，是我 国汽车制造工业中亟待解决的关键性问题。 从理论上说，板材冲压成形过程可以被看作是板材经过塑性变形变为想要获得的形状的过程。然而实际上．板料尺寸．材料特性和环境条件使冲压成形过程的预测性和可重复性变得困难。以韧性金属板材为主的冲压成形件从模具上取出后，必然产生一定量的回弹。回弹是板材冲压成形的3种主要缺陷(起皱．破裂和回弹)中最难控制的一种，因为它涉及到对回弹量的准确预示．不同的材料和尺寸的零件其回弹规律大不相同，单凭经验和工艺过程类比是很难进行准确的回弹补偿的．这就使得一个模具设计的周期变长．因此在板材冲压成形中回弹变形是使模具设计明显变复杂的一个基本参数。在大多数板材冲压成形中．强烈的非线性变形过程致使板料产生很大的弹性应变能．在模具与板料动态接触过程中存在于板料中的这种弹性应变能会随着接触压力的消除而自动释放掉，回弹的驱动力一般是朝着板料原始形状变形。因此，冲压成形中的最终产品形状不但依赖于凹模形状．而且依赖于成形后存储在板料中的弹性应变能。弹性应变能与许多诸如材料特性．接触载荷等参数有关，因此在成形过程中预测回弹变得很复杂．这也就给那些必须精 确评估回弹量的设计者提出了很重要的问题。 近40年来，有许多研究人员一直在对回弹行为进行着研究．并提出了很多解决方法和计算机仿真算法．发表了大量相关论文。就有限元仿真方法而言．在众多仿真算法模拟应用中，采用显式算法模拟成形过程．用隐式算法模拟回弹过程的方法最多；其次是冲压成形和卸载回弹过程都采用隐式算法。而G．Y-L．等学者提出一种新算法，冲压成形和回弹过程全部采用显式算法。U．Abdelsalam等学者还提出了采用一步成形算法模拟冲压成形过程，再用隐式算法计算卸载回弹过程．并应用该算法模拟了3个复杂冲压件的卸载回弹过程．这种算法的模拟精度虽然不高．但计算速度很快．可以为模具在设计阶段提供一个定性的参考方案。T-C．Hsu等学者采用隐式TL(Total Lagrangian)算法，引入Hill--次方屈服函数模拟了轴对称问题的冲压成形和回弹过程。M．Kawka等学者采用静态显式有限元(实际上也是隐式算法)算法软件ITAS3D模拟了轿车顶盖和轮毂的多阶段成形过程，以及卸载回弹和切边回弹过程．并与试验结果进行了比较。 以上这些对于回弹的研究只限于理论方面．其与实际试验的对比验证还鲜有涉及。对于如何补偿所产生的

有限元网格剖分方法概述

有限元网格剖分方法概述 在采用有限元法进行结构分析时，首先必须对结构进行离散，形成有限元网格，并给出与此网格相应的各种信息，如单元信息、节点坐标、材料信息、约束信息和荷载信息等等，是一项十分复杂、艰巨的工作。如果采用人工方法离散对象和处理计算结果，势必费力、费时且极易出错，尤其当分析模型复杂时，采用人工方法甚至很难进行，这将严重影响高级有限元分析程序的推广和使用。因此，开展自动离散对象及结果的计算机可视化显示的研究是一项重要而紧迫的任务。 有限元网格生成技术发展到现在, 已经出现了大量的不同实现方法，列举如下： 映射法 映射法是一种半自动网格生成方法，根据映射函数的不同，主要可分为超限映射和等参映射。因前一种映射在几何逼近精度上比后一种高，故被广泛采用。映射法的基本思想是：在简单区域内采用某种映射函数构造简单区域的边界点和内点，并按某种规则连接结点构成网格单元。也就是根据形体边界的参数方程，利用映射函数，把参数空间内单元正方形或单元三角形(对于三维问题是单元立方体或单元四面体)的网格映射到欧氏空间，从而生成实际的网格。这种方法的主要步骤是，首先人为地把分析域分成一个个简单可映射的子域，每个子域为三角形或四边形，然后根据网格密度的需要，定义每个子域边界上的节点数，再根据这些信息，利用映射函数划分网格。 这种网格控制机理有以下几个缺点： （1）它不是完全面向几何特征的，很难完成自动化，尤其是对于3D区域。 （2）它是通过低维点来生成高维单元。例如，在2D问题中，先定义映射边界上的点数，然后形成平面单元。这对于单元的定位，尤其是对于远离映射边界的单元的定位，是十分困难的，使得对局部的控制能力下降。 （3）各映射块之间的网格密度相互影响程度很大。也就是说，改变某一映射块的网格密度，其它各映射块的网格都要做相应的调整。 其优点是：由于概念明确，方法简单，单元性能较好，对规则均一的区域，适用性很强，因此得到了较大的发展，并在一些商用软件如ANSYS等得到应用。 2 。拓扑分解法 拓扑分解法较其它方法发展较晚, 它首先是由Wordenwaber提出来的。该方法假设最后网格顶点全部由目标边界顶点组成, 那么可以用一种三角化算法将目标用尽量少的三角形完全分割覆盖。这些三角形主要是由目标的拓扑结构决定, 这样目标的复杂拓扑结构被分解成简单的三角形拓扑结构。该方法生成的网格一般相当粗糙, 必须与其它方法相结合, 通过网格加密等过程, 才能生成合适的网格。该方法后来被发展为普遍使用的目标初始三角化算法, 用来实现从实体表述到初始三角化表述的自动化转换。 单一的拓扑分解法因只依赖于几何体的拓扑结构使网格剖分不理想，有时甚至很差。 3．连接节点法 这类方法一般包括二步：区域内布点及其三角化。早期的方法通常是先在区域内布点, 然后再将它们联成三角形或四面体, 在三角化过程中, 对所生成的单元形状难于控制。随着Delaunay三角化(简称为DT ) 方法的出现, 该类方法已成为目前三大最流行的全自动网格生成方法之一。 DT法的基本原理：任意给定N个平面点Pi(i=1,2,…,N)构成的点集为S，称满足下列条件的点集Vi为Voronoi多边形。其中，Vi满足下列条件： Vi ={ X：|X- Pi|(|X- Pj|,X(R2,i(j,j=1,2,…,N }Vi为凸多边形，称{ Vi}mi=1为Dirichlet Tesselation

ANSYS 网格划分方法总结

(1) 网格划分定义：实体模型是无法直接用来进行有限元计算得，故需对它进行网格划分以生成有限元模型。有限元模型是实际结构和物质的数学表示方法。 在ANSYS中，可以用单元来对实体模型进行划分，以产生有限元模型，这个过程称作实体模型的网格化。本质上对实体模型进行网格划分也就是用一个个单元将实体模型划分成众多子区域。这些子区域（单元），是有属性的，也就是前面设置的单元属性。 另外也可以直接利用单元和节点生成有限元模型。 实体模型进行网格划分就是用一个个单元将实体模型划分成众多子区域（单元）。 （2）为什么我选用plane55这个四边形单元后，仍可以把实体模型划分成三角 形区域集合？？？ 答案：ansys为面模型的划分只提供三角形单元和四边形单元，为体单元只提供四面体单元和六面体单元。不管你选择的单元是多少个节点，只要是2D单元，肯定构成一个四边形或者是三角形，绝对没有五、六边形等特殊形状。网格划分也就是用所选单元将实体模型划分成众多三角形单元和四边形子区域。 见下面的plane77/78/55都是节点数目大于4的，但都是通过各种插值或者是合并的方式形成一个四边形或者三角形。 所以不管你选择什么单元，只要是对面的划分，meshtool上的划分类型设置就只有tri和quad两种选择。 如果这个单元只构成三角形，例如plane35，则无论你在meshtool上划分设置时tri还是quad，划分出的结果都是三角形。

所以在选用plane55单元，而划分的是采用tri划分时，就会把两个点合并为一个点。如上图的plane55,下面是plane单元的节点组成，可见每一个单元上都有两个节点标号相同，表明两个节点是重合的。 。 同样在采用plane77 单元，进行tri划分时，会有三个节点重合。这里不再一一列出。（3）如何使用在线帮助： 点击对话框中的help，例如你想了解plane35的相关属性，你可以

先进板料成形技术与性能

板料成形有限元分析的发展综述 摘要：在参阅和分析大量有关文献的基础上，对有限元法的产生和弹塑性有限元的发展进行了总结，特别是对当前应用广泛的板料成形有限元数值模拟在国内外的发展概况和发展趋势进行了详尽的剖析，为深入了解板料成形有限元的发展提供了有益的参考。 关键词：板料成形；数值模拟；有限元法；有限元分析；弹塑性 引言 有限单元法是工程计算领域的一种主要的数值计算方法，其基本思想就是将连续区域上的物理力学关系近似地转化为离散规则区域上的物理力学方程。它是一种将连续介质力学理论、计算数学和计算机技术相结合的一种数值分析方法。此方法由于其灵活、快捷和有效，已迅速发展成为板料冲压成形中求解数理方程的一种通用的数值计算方法。 有限元法源于40年代提出的结构力学的矩阵算法。“有限元法”这一术语是R.W.Clough于1960年在论文“The finite element method in plane stress analysis”中首次提出来的，他用这种方法首次求解了弹性力学的二维平面应力问题。1963年，Besseling证明了有限元法是基于变分原理的Ritz法的另一种形式，从而使Ritz分析的所有理论基础都适用于有限元法，确认了有限元法是处理连续介质问题的一种普遍方法。 板料成形数值模拟涉及到连续介质力学中材料非线性、几何非线性、边界条件非线性等三非线性问题的计算，难度很大。随着非线性连续介质力学理论、有限元法和计算机技术的发展，通过高精度的数值计算来模拟板料成形过程已成为可能。从70年代后期开始，经过近二十年的发展，板料成形数值模拟逐渐走向成熟，并开始在汽车、飞机等工业领域得到实际的应用。 1 弹塑性有限元分析研究发展概况 有限元法建立之初，只能处理弹性力学问题，无法应用于金属塑性成形分析。1965年Marcal提出了弹塑性小变形的有限元列式求解弹塑性变形问题，揭开了有限元在塑性加工领域应用的序幕。1968年日本东京大学的Yamada推导了弹塑性小变形本构的显式表达式，为小变形弹塑性有限元法奠定了基础。但小变形理论不适于板料冲压成形这样的大变形弹塑性成形问题，因此人们开始致力于研究大变形弹塑性有限元法。1970年美国学者Hibbitt等首次利用有限变形理论建立了基于Lagrange格式（T.L格式）的弹塑性大变形有限元列式。1973年Lee 和Kabayashi提出了刚塑性有限元法。1973年Oden等建立了热-弹粘塑性大变形有限元列式。1975年Mcmeeking建立了更新Lagrange格式（U.L格式）的弹塑性大变形有限元列式。1978年Zienkiewicz等提出了热耦合的刚塑性有限元法。1980年Owen出版了第一本塑性力学有限元的专著，全面系统地论述了材料非线性和几何非线性的问题。至此，大变形弹塑性有限元理论系统地建立起来了。 2 板料成形有限元数值模拟国内外研究发展概况

ANSYS结构有限元分析中的网格划分技术及其应用实例

一、前言 有限元网格划分是进行有限元数值模拟分析至关重要的一步，它直接影响着后续数值计算分析结果的精确性。网格划分涉及单元的形状及其拓扑类型、单元类型、网格生成器的选择、网格的密度、单元的编号以及几何体素。从几何表达上讲，梁和杆是相同的，从物理和数值求解上讲则是有区别的。同理，平面应力和平面应变情况设计的单元求解方程也不相同。在有限元数值求解中，单元的等效节点力、刚度矩阵、质量矩阵等均用数值积分生成，连续体单元以及壳、板、梁单元的面内均采用高斯（Gauss）积分，而壳、板、梁单元的厚度方向采用辛普生（Simpson）积分。辛普生积分点的间隔是一定的，沿厚度分成奇数积分点。由于不同单元的刚度矩阵不同，采用数值积分的求解方式不同，因此实际应用中，一定要采用合理的单元来模拟求解。 CAD软件中流行的实体建模包括基于特征的参数化建模和空间自由曲面混合造型两种 方法。Pro/E和SoildWorks是特征参数化造型的代表，而CATIA与Unigraphics等则将特征参数化和空间自由曲面混合造型有机的结合起来。现有CAD软件对表面形态的表示法已经大大超过了CAE软件，因此，在将CAD实体模型导入CAE软件的过程中，必须将CAD 模型中其他表示法的表面形态转换到CAE软件的表示法上，转换精度的高低取决于接口程序的好坏。在转换过程中，程序需要解决好几何图形（曲线与曲面的空间位置）和拓扑关系（各图形数据的逻辑关系）两个关键问题。其中几何图形的传递相对容易实现，而图形间的拓扑关系容易出现传递失败的情况。数据传递面临的一个重大挑战是，将导入CAE程序的CAD模型改造成适合有限元分析的网格模型。在很多情况下，导入CAE程序的模型可能包含许多设计细节，如细小的孔、狭窄的槽，甚至是建模过程中形成的小曲面等。这些细节往往不是基于结构的考虑，保留这些细节，单元数量势必增加，甚至会掩盖问题的主要矛盾，对分析结果造成负面影响。 CAD模型的“完整性”问题是困扰网格剖分的障碍之一。对于同一接口程序，数据传递的品质取决于CAD模型的精度。部分CAD模型对制造检测来说具备足够的精度，但对有限元网格剖分来说却不能满足要求。值得庆幸的是，这种问题通常可通过CAD软件的“完整性检查”来修正。改造模型可取的办法是回到CAD系统中按照分析的要求修改模型。一方面检查模型的完整性，另一方面剔除对分析无用的细节特征。但在很多情况下，这种“回归”很难实现，模型的改造只有依靠CAE软件自身。CAE中最直接的办法是依靠软件具有的“重构”功能，即剔除细部特征、缝补面和将小面“融入”大曲面等。有些专用接口在模型传递过程中甚至允许自动完成这种工作，并且通过网格剖分器检验模型的“完整性”，如发现“完整性”不能满足要求，接口程序可自动进行“完整性”修复。当几何模型距CAE分析的要求相差太大时，还可利用CAE程序的造型功能修正几何模型。“布尔运算”是切除细节和修理非完整特征的有效工具之一。 目前数据传递一般可通过专用数据接口，CAE程序可与CAD程序“交流”后生成与CAE 程序兼容的数据格式。另一种方式是通过标准图形格式如IGES、SAT和ParaSolid传递。现有的CAD平台与通用有限元平台一般通过IGES、STL、Step、Parasolid等格式来数据