基于ABAQUS的板材冲压仿真及回弹分析
华南理工大学模具研究室
1.5.4 Interaction
中定义接触： 在ABAQUS中定义接触： 中定义接触
定义接触面（从面和主面） 定义接触面（从面和主面） 接触相互作用 小滑动与有限滑动 摩擦
从面： 从面： 变形体 主面： 主面： 刚体
华南理工大学模具研究室
7
1.5.5 Load
• 施加边界条件 • 在ABAQUS中平移自由度和旋转自由度 中平移自由度和旋转自由度 的正方向规定如下： 的正方向规定如下：
华南理工大学模具研究室
1.5.6 Elements in ABAQUS
•单元类型（Family） 单元类型（ 单元类型 ） •A family of finite elements is the broadest category used to classify elements.
continuum (solid elements)
华南理工大学模具研究室
1.4ABAQUS/CAE的优势
• 将建模、分析、作业管理和结果评 将建模、分析、 估无缝集成。 估无缝集成。 • 可以为具体的应用系统定制界面。 可以为具体的应用系统定制界面。 • 与Ansys、Marc等通用有限元软件 、 等通用有限元软件 相比， 在单元种类、 相比，Abaqus在单元种类、材料 在单元种类 模型、 模型、丰富的用户子程序接口等方 面都有着显著的优势。 面都有着显著的优势。
基于ABAQUS的板材冲压仿真及回弹分析
华南理工大学模具研究室
主要内容
• 一、ABAQUS概述 概述
• What is ABAQUS? • ABAQUS有限元软件的功能 • ABAQUS/CAE的优势 • ABAQUS/CAE的基本流程
• 二、典型实例讲解
• 叶片板冲压回弹 • Springback of two-dimensional draw bending（二维弯曲回弹模拟）
华南理工大学模具研究室
3
1.3ABAQUS 有限元软件的功能
• 线性静力学, 动力学, 和热传导 线性静力学 动力学
应力, 振动, 声场, 地质力学, 压电效应, 例如 应力 振动 声场 地质力学 压电效应 等 汽车、飞机机身等的静力和动力学响应 结构刚度, 汽车、飞机机身等的静力和动力学响应, 结构刚度 等
•表2－1 常用量纲 表 －
华南理工大学模具研究室
1.5.2 定义部件(part)
• 可以创建模型的部件类型
• 可变形部件 • 任意形状的，可以包含不同维数的特征 （实体、表面、线）；在载荷作用下可以 变形。 任意形状的；在载荷作用下不可以变形。 只可以用直线、圆弧和抛物线创建的形状； 在载荷作用下不可以变形。 可变形部件: 轴机架
华南理工大学模具研究室
10
华南理工大学模具研究室
华南理工大学模具研究室
11
学习资料
• 帮助文件（Documentation) 帮助文件（ • 书籍： 书籍： 《ABAQUS在机械工程中的应用》 《ABAQUS有限元分析实例详解》 《ABAQUS非线性有限元分析与实例》 • 论坛： 论坛： 仿真科技论坛
shell elements
beam elements
•同类型单元有很多相 同的基本特。 Elements in the same family share many basic features. •同种类单元又有很多 变化：There are many variations within a family.
Visualization
华南理工大学模具研究室
5
1.5.1量纲 1.5.1量纲
• 在建立模型前，必须先决定量纲系统。 没有固定的量纲系统,但一 在建立模型前，必须先决定量纲系统。ABAQUS没有固定的量纲系统 但一 没有固定的量纲系统 个问题所有的输入数据只能用同一个量纲系统,下表为常用的量纲系统。 个问题所有的输入数据只能用同一个量纲系统 下表为常用的量纲系统。 下表为常用的量纲系统
• ABAQUS模型数据的全面视角 ABAQUS模型数据的全面视角
华南理工大学模具研究室
1.5ABAQUS/CAE的基本流程
Part
•创建部件的几 何模型（如果需 要，还需创建区 域）
Property
•定义材料 •定义部件区域的 截面属性 •为部件或区域定 义和分配截面属性
Assembly
•为初始构型定位部 件。如果模型中只包 含一个部件，装配件 将被自动创建。
ABAQUS/CAE main user interface
华南理工大学模具研究室
4
ABAQUS/CAE的基本特征优势
• 模型树 • 模型树使得对你的模型以及模型包含 的对象可以有一个图形上的直观概述
• 模型 查询/审核
• 例如, 显示截面, 材料, 当前定义的约束 的数目 • 显示某些特征/项目的状态 (无效的部件, 隐藏的特征)
• 三、常见问题讨论
华南理工大学模具研究室
1
ABAQUS的概述
华南理工大学模具研究室
1.1ABAQUS, Inc.?
• 创建于 1978年 年 • 2005年 和Dassault Systèmes 公司合并 年 • 总部位于 总部位于Providence, Rhode Island (USA) • ABAQUS 被广泛地认为是功能非常强的非线性有限元软件，可以 被广泛地认为是功能非常强的非线性有限元软件， 分析复杂的结构力学系统， 分析复杂的结构力学系统，尤其能够驾驭非常庞大复杂的问题和 模拟高度非线性问题。 模拟高度非线性问题。
•
华南理工大学模具研究室
6
1.5.3 Property
• 材料塑性的定义
• 在ABAQUS中必须用真实应力和真是应变定义 中必须用真实应力和真是应变定义 塑性。 塑性。大多数实验数据常常是名义应力和名义应变的 值给出的。这时，必须应用转换公式（见入门指南） 值给出的。这时，必须应用转换公式（见入门指南） 把塑性材料的数据从名义应力/应变的值转换为真实 把塑性材料的数据从名义应力 应变的值转换为真实 应力/应变值。 应力 应变值。 应变值
• 非线性和瞬态分析
接触, 塑性失效, 断裂和磨损, 复合材料, 接触 塑性失效 断裂和磨损 复合材料 超弹性 等 汽车碰撞, 电子器件跌落, 汽车碰撞 电子器件跌落 冲击和损毁等
• 多体动力学分析
同时结合刚体, 线性柔体, 同时结合刚体 线性柔体 和非线性柔体模拟各种连接件等 应用在：汽车运动 高速机械, 微机电系统MEMS, 应用在：汽车运动, 高速机械 微机电系统 航空航天机构, 医疗器械, 航空航天机构 医疗器械 等
• • • • • • 1. 方向1的平移：U1 2. 方向2的平移：U2 3. 方向3的平移：U3 4. 绕轴1的旋转：UR1 5. 绕轴2的旋转：UR2 6. 绕轴3的旋转：UR3
• 施加载荷 • 施加载荷
• 在加完约束后就应当在结构的底部加载荷。 在ABAQUS里，载荷通常是指从初始状态下使 结构响应发生变化的各种因素，如： • 集中力 • 压力 • 非零边界条件 • 温度（材料热膨胀特性被定义后） • ……
/forum/index.php 仿真在线 /bbs/ ……
华南理工大学模具研究室
Have fun with ABAQUS
Q&A
华南理工大学模具研究室
12
• •
离散刚体部件 • • 解析刚体部件
• 可以通过三种方法定义部件几何体： 可以通过三种方法定义部件几何体：
• • 直接由ABAQUS/CAE创建 创建 直接由 从CAD系统导入 系统导入 • 本地的CAD几何体 • 中性格式 作为ABAQUS (孤立 网格导入 孤立)网格导入 作为 孤立 离散刚体部件：任意形状的刚 体冲垫（带有可变形的毛坯） 解析刚体部件：旋转的刚性扎滚 （带有可变形的金属原料）
华南理工大学模具研究室
2
1.2What is ABAQUS?
• ABAQUS是一组有限元分析模块
华南理工大学模具研究室
What is ABAQUS?
• ABAQUS的核心是它的求解器模块，ABAQUS/Standard和 的核心是它的求解器模块， 的核心是它的求解器模块 和 ABAQUS/Explicit, 是互相补充的、集成的分析模块。 是互相补充的、 集成的分析模块。
Step
•定义分析步和输出需求In来自eraction•为区域、集定义接 触或其他形式的相 互作用，并在指定 的分析步应用它们
Load
•为区域或集指定 载荷、边界条件和 场，并在指定的分 析步应用它们
Mesh
•将装配件分割 成可以分网的 区域并划分网 格
Job
•提交、管理、并 监控分析作业 •检验结果
rigid elements
membrane elements
infinite elements
special-purpose elements like springs, dashpots, and masses
truss elements
华南理工大学模具研究室
8
二、典型实例讲解
华南理工大学模具研究室
ABAQUS软件冲压成形回弹模拟流程图 软件冲压成形回弹模拟流程图
ABAQUS软件冲压成形回弹模拟流程图
华南理工大学模具研究室
9
实例讲解
华南理工大学模具研究室
基于ABAQUS的冲裁仿真方法
Ｈｙ ｐ ｅ ｒ ｆ ｏ ｒ ｍ 模 块［ ８ ］ 。本文 使用 ＡＢ ＡＱＵＳHale Waihona Puke ／ Ｅ ｘ ｐ ｌ ｉ ｃ ｉ ｔ
计算结 果 相 近， 但计 算 代 价 较小， 所 以 使 用
ＭＳ Ｆ Ｌ Ｄ准 则 描述 材 料 的颈 缩 失效 ， 并 通 过指 定 成 形 极 限 曲线 的方 式定 义 ＭＳ Ｆ Ｌ Ｄ准则 。
中图 分类号 ： Ｔ Ｇ ３ ８ ６ ． ２
文献 标 志码 ： Ａ
文章 编号 ： ２ ０ ９ ５ — ４ ８ ５ ９ （ ２ ０ １ ６ ） ０ ３ — ０ ３ ３ ３ — ０ ４
冲裁是 使用安 装在 压力机 上 的模具 ， 使板 料 分 １ 本 构 方 程 离 以获 得零 件 的一 种 冲压工艺 ， 包 括落 料 、 冲孔 、 切 零件 的 材 料 使 用 铝 合 金 ， 材 料参 数 来 自 ＡＱＵＳ的帮助 文件 。在 冲 裁过 程 中 ， 零 件 材 料 边、 切断、 切 口、 剖切 、 冲槽 和修 边 等 ［ １ ］ 。冲裁 生 产 ＡＢ
对 冲 裁 力 和 断 口形 态 的 影 响 进行 了仿 真 。 进行 了 厚 板 冲 裁 的 仿 真 ， 研 究 了板 料 厚 度 对 冲 裁 力 的 影 响 。并 且 建 立 了精 密 冲 裁 模 型 ， 仿 真 了精 密 冲裁 中模 具 间 隙对 断 口形 态 的影 响 。
关 键 词 冲裁 ； 有 限元 ； 精密冲裁 ； 厚板 ； Ａ Ｂ Ａ Ｑ Ｕ Ｓ 仿真
是 一个剧 烈 的局部塑 性变形 过 程 ， 在狭 小 的剪切 带
异性的， 这里 将 材料 的弹 塑性行 为 近似作 为 各 向同 性处 理 ， 使用 线 弹性和 Ｍｉ ｓ ｅ ｓ 屈 服准 则 。
Abaqus钣金冲压成型例题讲解I
钣金成型例题讲解一、背景当前，制造行业加工工艺的趋势正朝着高新技术的方向发展。
由于新产品、新技术的开发成本太高、开发时间过长，加上开发成果没有保障，越来越多的公司在研发、制造过程中开始注重仿真技术的应用。
采用ABAQUS对加工工艺进行模拟有着诸多优点：1.数值模拟减少了耗时的原型实验，缩短了产品投放市场的时间；2.合理的参数设计可以降低对工件的损耗；3.合理的坯料设计，减少了飞边，也减少原材料的浪费；4.对模具的设计、加工提供合理建议；5.优化加工过程，提高产品成型质量；采用ABAQUS进行仿真模拟的目的：1.节约开发成本2.加快研发速度3.提高产品质量二、问题的描述本实例模拟油箱的冲压成型过程。
图1为实际的油箱形状，是由两个如图2所示的结构组成，考虑冲压成型过程中，它的结构的对称性，我们通过建立图3所示的结构，对其进行模拟分析，达到分析整个油箱成型分析的目的。
首先，我们将通过ABAQUS/CAE完成图4所示的装配图，其中平面铝板将被冲压成型为图3的结构。
三、建立模型3.1创建成型模具－阳模1、首先运行ABAQUS/CAE ，在出现的对话框内选择Create Model Database 。
2、在主菜单model 中命名新建模型为Forming example ，并保存文件为examle_forming.cae 。
3、从Module 列表中选择Part ，进入Part 模块。
4、选择Part→Create 来创建一个新的零件。
在提示区域会出现这样一图1图 2成型模具－阳模金属板压边条成型模具－阴模图4图3图5个信息。
5、CAE 弹出一个如图5的对话框。
将这个零件命名为punch ，确认Modeling Space 、Type 和BaseFeature 的选项如下图。
输入300作为Approximate size 的值。
点击Continue 。
ABAQUS/CAE 初始化草图，并显示格子。
6、在左侧工具条上点击 ，在提示栏中依次输入下表的坐标点，采用图标 连接1和2点、6和7点，采用连接图中2、3、4点和4、5、6点。
板料冲压成形及回弹有限元模拟分析
Ku=F
式中K为总刚度矩阵；u为位移向量；F为节点载荷向量。
在隐式算法中，对于第i个给定的加载增量，用Newton-Raphson迭代法，需要求解下面的方程：
板料冲压成形及回弹有限元模拟分析
摘要
回弹是板材冲压成形过程中不可避免的普遍现象，直接影响到冲压件的尺寸精度和零件最终形状。本文利用ANSYS/LS-DYNA有限元软件中的非线性动力的显式、隐式连续求解功能，模拟了板料冲压成形过程与卸载后板料回弹变形的全过程，得到了成形过程中任一时刻各处Von-Mises应力云图和应变值及卸载后板料的回弹结果，帮助我们更好的认识分析板料冲压成形以及回弹过程中物质内部的变化。
ANSYS／LS-DYNA是一个通用的显示非线性动力分析有限元程序，近年来开发的板料冲压成形数值模拟新功能，取得了很大成功。通过计算，可以观察板料冲压成形过程中的变形状态、应力状态和壁厚变化，了解可能出现的起皱和开裂现象。此外，ANSYS／LS-DYNA程序具有显式、隐式求解功能，用显式求解模拟动态成形过程，然后用隐式求解模拟线性回弹变形。这不仅能够模拟材料变形过程，而且也能较好地计算回弹现象，比较准确地得到材料最终成形状态。
关键词：板材冲压，回弹，非线性有限元分析，数值模拟
Sheet metal stamping and rebound finite element simulation analysis
Abstract
The rebound is inevitable common phenomenon in sheet metal forming process, a direct impact on the final shape to the dimensional accuracy of the stampings and parts. In this paper, the nonlinear dynamic finite element software ANSYS / LS-DYNA explicit, implicit sequential solution function to simulate thespringbackdeformation of the sheet after sheet metal stamping process and uninstall the whole process, forming process at any time throughout the Von-Misesstress cloud and strain and after unloading sheetspringbackresults, help us to a better understanding analysis sheet metal stamping and rebound process material internal changes.
基于ABAQUS软件的热冲压成形接触问题分析
基于ABAQUS软件的热冲压成形接触问题分析1. 引言1.1 背景介绍热冲压成形是一种重要的金属压制加工方法，广泛应用于汽车制造、航空航天等领域。
热冲压成形可以在金属板材加工过程中实现快速成形和高精度，同时也能改善材料的塑性变形性能。
热冲压成形接触问题是该工艺中的关键问题之一，直接影响产品质量和加工效率。
通过对热冲压成形接触问题的深入研究，可以优化工艺参数，提高成形质量和效率，降低生产成本。
目前，随着计算机仿真技术的发展，研究人员可以利用ABAQUS 等有限元软件对热冲压成形接触问题进行数值模拟分析。
通过建立合适的模型和边界条件，可以模拟实际加工过程中的力学行为和热传导过程，得到准确的结果。
基于ABAQUS软件的热冲压成形接触问题分析已经成为研究的重要方向，有助于揭示工艺中的关键问题，指导实际生产中的操作。
1.2 研究目的热冲压成形是一种先进的金属成形工艺，在汽车、航空航天和其他工业领域具有广泛应用。
研究旨在探讨基于ABAQUS软件的热冲压成形接触问题分析方法，以提高热冲压成形工艺的效率和成形质量。
具体目的包括：1. 分析热冲压成形中的接触问题，探讨接触压力分布、接触变形和接触面积等关键参数对成形质量的影响。
2. 建立热冲压成形接触问题的数值模拟方法，深入研究不同参数下的接触行为，为优化工艺提供理论依据。
3. 开展热冲压成形接触问题的实例分析，验证数值模拟方法的准确性和可靠性，为工程实践提供指导。
通过对热冲压成形接触问题的深入研究，旨在为提高热冲压成形工艺的稳定性、成形精度和生产效率提供理论支持，推动热冲压成形技术的进步和发展。
1.3 研究意义热冲压成形是一种在汽车工业、航空航天等领域广泛应用的先进制造工艺，其通过将金属材料在高温下进行塑性变形，以获得复杂形状和高强度的零部件。
研究热冲压成形接触问题对于优化成形工艺、提高产品质量和降低生产成本具有重要意义。
研究热冲压成形接触问题可以帮助了解成形工艺中材料之间的接触行为，验证成形过程中的接触压力分布和接触面积变化等参数。
CAE仿真经典文献系列之《板料冲压成形回弹的数值模拟》
（ i） （ i） 正向量； 为外力向量； 为内力向量 . P @ 位移增量由 （3） 式得到 . （ i +l） （ i） （ i） = !u + C !u 算法中的要点如下：
（3）
同的增量方法将所有接触节点力按比例卸载， 直 到接触节点力消失 . 具体的可总结如下： 首先给拉深结束后的板料一个假想的应力状 态! ， 该应力与拉深结束后板料的应力 !E 相平 衡， 可以使变形板料在进行回弹分析的开始处于 静态平衡状态 . 在分析过程中， 假想应力按照步长
图3
回弹后的零件截面线形状
回弹现象进行了模拟， 证明该算法还是比较稳定 可靠的 图 1 和图 2 给出了采用该算法对轨道件 拉深成形和回弹后形状的模 拟 结 果 其 中， E 2 6 8 MPa， 3， c 8 mm， f 1， ! " 4 216 （1 537 96 #） MPa 该 零 件 为 NUMISHEET’ 标准考题 图 3 为该标准考题给出的回弹后的零 件截面线形状， 该截面线过图 2 所示 E 、 F点 对 比图 1 和图 2 可以看出， 由于残余应力作用， 轨道 件拉深卸载后发生了较为明显的回弹现象 轨道 底部与侧壁间弯曲角变大， 侧壁向外倾斜； 侧壁与 凸缘间弯曲角变小， 凸缘向下倾斜 对比图 2 和 图 3 可知， 回弹模拟结果与标准考题参考结果基
［1］Joannic D，GeIin J C Accurate simuIation of springback in 3D sheet ［ A］ In：NUMIFORM’ ［ C］ 1995 729 metaI forming processes 95 ~ 734 ［2］ KarafiIIis A P，Boyce M C TooIing design in sheet metaI forming ［ J］ Int J Mech Sci，1992，34 （2 ） ： using springback caIcuIations 113 ~ 131 ［ 3］Wu Longwu Iterative FEM die surface design to compensate for spri［ A ］ In： ［ C］ NUMIFORM ’ 95 ngback in sheet metaI stampings 1995 637 ~ 641 ［4］宋 黎， 杨 坚， 黄天泽 板料弯曲成形的回弹分析与工程控 制综述 ［ J］ 锻压技术， （1） ：18 ~ 22 1996，21 ［5］ Shu Jawshi，Huang Chinghua Finite eIement anaIysis and optimiza“doubIe-bend”technigue ［ J］ Int J tion of springback reduction：the （4） ：423 ~ 434 Mech TooIs Manufacture，1996，36 ［6］ Mercer C D，NagtegaaI J D，RebeIo N，ec al Effective appIication
基于Abaqus的冲压模具仿真
板件冲压成型分析、背景描述当前，制造行业加工工艺的趋势正朝着高新技术的方向发展。
由于新产品、新技术的开发成本太高、开发时间过长，加上开发成果没有保障，越来越多的公司在研发、制造过程中开始注重仿真技术的应用。
采用ABAQUS 对加工工艺进行模拟有着诸多优点：1.数值模拟减少了耗时的原型实验，缩短了产品投放市场的时间；2.合理的参数设计可以降低对工件的损耗；3.合理的坯料设计，减少了飞边，也减少原材料的浪费；4.对模具的设计、加工提供合理建议；5.优化加工过程，提高产品成型质量；采用ABAQUS 进行仿真模拟的目的：1.节约开发成本2.加快研发速度3.提高产品质量二、问题描述1、本实例是关于板件的冲压成型仿真模拟，图 1 为拉延模具三维实体，图2为模具上模的刚性实体，图 3 为模具下模的刚性实体，图 4 为模具压边圈的刚性实体，图5 为毛坯的初始形状，其中定义上模、压边圈、下模为离散刚性的，即在模拟过程中假设这几个部件不发生变形，只有毛坯件发生变形，完成冲压，在分析过程中定义毛坯为可变性的。
图1拉延模三维模型图2模具上模刚性实体图3模具下模刚性实体图4压边圈刚性实体图5毛坯初始形状2、将上述模具四部分进行装配，其中上模在最上端，下模在最下端，压边圈在下模上，毛坯在压边圈下方下模上放着，在模具冲压的过程中，下模固定不动, 由上模向下运动来实现板料的冲压成型，压边圈在上模下压之前已将将板料压紧，防止板料在拉伸的过程毛坯串动，和提高拉伸后型面质量。
拉着板料的边缘, 使拉伸时板料不易起皱。
下图6模具为在abaqus中的装配体。
图6模具装配体3、定义毛坯的材料属性，进入属性模块，定义新材料名称为steel,该材料的密度为7.85E-9,杨氏模量为210000,泊松比为0.31,定义塑性属性为图7所示, True stress (MPa) Log plastic stra in910.01310.15910-21710.64910-22110.17710-12510.39510-12910.77610-13310.1393910.295图7材料的属性定义4、创建截面类型别壳，类型为均质的截面，同时将截面特性赋予给毛坯。
基于Abaqus的VRB板冲压成形数值模拟研究
图 1 板料柔性轧制
B
L1
L2
L3
L
T
t
（a） 模型装配图
（b） 坯料
注： 1— ——凸 模； 2—— —平板 压边 圈； 3—— —实 体 单 元 板 料 ； 4— ——壳单元板料 （具体用哪种依据仿真单元而定）；5 凹模。
图 2 有限元模型及坯料
材料
B210P B210P
表 1 材料性能参数
厚度 （mm）
2 1
杨氏模量 （GPa）
210 210
泊松比
0.3 0.3
屈服强度 抗拉强度 （MPa） （MPa）
240
421
265
412
n值
0.265 0.241
【作者简介】张波，男，广西人，本科，任职于上汽通用五菱汽车股份有限公司技术中心，研究方向：车身结构设计、车身轻量化设计。
76 qiyekejiyufazhan
2019 年第 6 期（总第 452 期） 过渡区材料的屈服变小，而靠近薄测的过渡区材料的屈服增 大，故在冲压过程中，由于靠近薄侧的过渡区材料抵抗塑性 变形的能力增加，所以薄测圆角区域的开裂趋势增加。
5.1.2 圆角区域减薄 减薄率云图显示最大减薄率在盒形件的 4 个圆角处 （如
图 10 所示）。 3 种离散程度情况下 4 个圆角处的减薄率数值如图 11 所
abaqus计算回弹的方法
Abaqus回弹计算过程回弹分析我倒是做过两个，说下简要步骤吧，同样是仅供参考啊1.首先用·explicit做成型过程的分析，加载方式选位移加载比较好，加载的幅值选smooth step（平滑变化）2.可适当的用质量放大来加快这一准静态分析的过程3.分析完成后可用standard观察工件的回弹，具体做法是：1.Model-Copy Model2.在新复制的模型中仅留下成型件，删除其他一切无关的边界条件以及上下模，包括在Explicit中定义的接触属性3.在step模块中创建predefine field request-others-initial state-last frame/last step（导入的job名称为之前做成型分析的那个job的名称）4.删除原来所有的后续分析步，并新建一个static，general的分析步5.创建一个新的作业提交分析，并观察回弹大致就是这样吧，希望对你有用！回弹分析，从explicit导入standard计算。
先copy explicit中模型进入standard模块，然后做一下改进，删除各个part、set和surface等，只留下需要回弹分析的变形体。
删除分析步，删除接触和属性。
然后在step中建立一个static分析步骤。
设置计算为非线性。
然后定义居于前面成形结果的回弹分析，在Model Tree中打开Predefined Fields，选择Initia 作为分析步，Other最为类别，选择Initial State，然后在视窗中选择需要分析的回弹体，然后点击done，然后Edit Predefined Field，选择你成形分析的job名字。
然后一致ok下去，对称的边界哦条件还要施加。
你可以在amplitude中设置，比如说你分析步设置时间为6s，然后在amplitude中设置0，0；4，1（也就是在4秒时冲头应景达到了要求的位移，也就是液晶冲完，那么剩下的2秒就是停留的时间了），然后在另外设置一个分析步把冲头往回移就可以了小弟这些天正好在做冲压回弹，刚做成功，从simwe论坛上学了很多东西。
abaqus计算回弹及方法
Abaqus回弹计算过程回弹分析我倒是做过两个，说下简要步骤吧，同样是仅供参考啊1.首先用·explicit做成型过程的分析，加载方式选位移加载比较好，加载的幅值选smooth step（平滑变化）2.可适当的用质量放大来加快这一准静态分析的过程3.分析完成后可用standard观察工件的回弹，具体做法是：1.Model-Copy Model2.在新复制的模型中仅留下成型件，删除其他一切无关的边界条件以及上下模，包括在Explicit中定义的接触属性3.在step模块中创建predefine field request-others-initial state-last frame/last step（导入的job名称为之前做成型分析的那个job的名称）4.删除原来所有的后续分析步，并新建一个static，general的分析步5.创建一个新的作业提交分析，并观察回弹大致就是这样吧，希望对你有用！回弹分析，从explicit导入standard计算。
先copy explicit中模型进入standard模块，然后做一下改进，删除各个part、set和surface等，只留下需要回弹分析的变形体。
删除分析步，删除接触和属性。
然后在step中建立一个static分析步骤。
设置计算为非线性。
然后定义居于前面成形结果的回弹分析，在Model Tree中打开Predefined Fields，选择Initia作为分析步，Other最为类别，选择Initial State，然后在视窗中选择需要分析的回弹体，然后点击done，然后Edit Predefined Field，选择你成形分析的job名字。
然后一致ok下去，对称的边界哦条件还要施加。
你可以在amplitude中设置，比如说你分析步设置时间为6s，然后在amplitude中设置0，0；4，1（也就是在4秒时冲头应景达到了要求的位移，也就是液晶冲完，那么剩下的2秒就是停留的时间了），然后在另外设置一个分析步把冲头往回移就可以了小弟这些天正好在做冲压回弹，刚做成功，从simwe论坛上学了很多东西。
基于ABAQUS的碳纤维复合材料板热冲压成形仿真
研究与开发合成纤维工业ꎬ２０１９ꎬ４２(２):１６ＣＨＩＮＡ㊀ＳＹＮＴＨＥＴＩＣ㊀ＦＩＢＥＲ㊀ＩＮＤＵＳＴＲＹ㊀㊀收稿日期:２０１８￣１０￣０２ꎻ修改稿收到日期:２０１９￣０１￣２５ꎮ作者简介:张华伟(１９８３ )ꎬ男ꎬ讲师ꎬ主要从事纤维增强复合材料的研究与教学工作ꎮＥ￣ｍａｉｌ:ｚｈａｎｇｈｗ＠ｎｅｕｑ.ｅｄｕ.ｃｎꎮ基金项目:国家自然科学基金项目(５１４７５０８６)ꎻ河北省自然科学基金项目(Ｅ２０１６５０１１１８)ꎻ中央高校基本科研业务费专项资金资助项目(Ｎ１７２３０４０３６)ꎻ河北省高等学校科学技术研究重点项目(ＺＤ２０１７３１５)ꎮ㊀∗通信联系人ꎮＥ￣ｍａｉｌ:２４７３２８２１０＠ｑｑ.ｃｏｍꎮ基于ＡＢＡＱＵＳ的碳纤维复合材料板热冲压成形仿真张华伟ꎬ李博宏∗(东北大学秦皇岛分校控制工程学院ꎬ河北秦皇岛０６６０００)摘㊀要:基于ＡＢＡＱＵＳ有限元分析软件对碳纤维增强聚醚醚酮(ＰＥＥＫ)复合材料板的热冲压成形工艺进行模拟ꎬ分析了复合材料在热冲压过程中的受力及变形特性ꎬ探讨了纤维铺层夹角㊁复合材料板层数和复合层类型对其成形性能的影响ꎮ结果表明:可以利用复合材料的工程常数建立材料的本构模型ꎻ在复合层类型的选择中ꎬ 三维实体￣连续壳 更符合实际情况ꎬ仿真效果更好ꎻ复合材料的纤维铺层夹角是其热冲压成形的主要影响因素ꎬ而且同其他的纤维铺层夹角(０ʎꎬ３０ʎꎬ４５ʎ)相比ꎬ复合材料在铺层夹角为９０ʎ时热冲压受力较好ꎬ但应变也较大ꎬ容易发生破坏ꎻ当总厚度一定时ꎬ复合材料板层数对材料的受力无明显影响ꎻ选择复合层类型为连续壳㊁９０ʎ夹角的复合材料板建立热冲压模型ꎬ仿真效果最好ꎮ关键词:碳纤维㊀复合材料㊀热冲压成形㊀有限元模型㊀模拟仿真中图分类号:ＴＱ３４２＋.７４㊀㊀文献标识码:Ａ㊀㊀文章编号:１００１￣００４２(２０１９)０２￣００１６￣０５㊀㊀复合材料是由两种或多种不同性质的材料用物理和化学方法在宏观尺度上组成的具有新性能的材料ꎮ通常来说ꎬ复合材料除了具有其组合材料的原始性能之外ꎬ还会产生一些原来组分材料所没有的新的性能ꎬ从而使得复合材料在强度㊁刚度㊁热力学等性能上有了一定幅度的改善[１]ꎮ碳纤维增强树脂基复合材料具有较高的强度ꎬ相对比金属材料来说ꎬ更加耐腐蚀㊁绝缘性好㊁制作方式简单ꎮ而且ꎬ碳纤维复合材料的密度小㊁比强度高ꎬ具有较好的轻量化效果和抗冲击性能ꎮ但这种材料也有其不足之处:材料各向异性严重㊁常温下成形性能差㊁生产成本较高等ꎮ张琦等[２]研究了碳纤维复合材料板的非等温模具热冲压ꎬ并以实验与有限元仿真结合的方式对纤维编织复合材料进行了分析ꎮ堵同亮等[３]建立了碳纤维编织复合材料的超弹性本构模型ꎬ研究了碳纤维编织复合材料在冲压成形带有双曲率曲面的结构件时的纤维重新排布和重新取向ꎮ但这些研究多数是以细观力学的分析方法来对复合材料冲压工艺仿真建模ꎬ对于从宏观的角度利用工程常数建模仿真的研究还相对较少[４－８]ꎮＡＢＡＱＵＳ软件是一款功能强大的有限元分析软件ꎮ在工程应用方面ꎬＡＢＡＱＵＳ基于丰富的单元库ꎬ可以用于模拟绝大多数的常见工程材料ꎬ如金属㊁聚合物㊁复合材料㊁橡胶等ꎮ作者利用ＡＢＡＱＵＳ软件对碳纤维复合材料板的宏观力学行为进行了研究ꎬ建立了碳纤维复合材料板在高温状态下的本构模型ꎬ在此基础上构建了碳纤维复合材料热冲压成形的有限元模型ꎬ完成了成形过程的有限元仿真ꎬ对比了不同的复合层类型㊁铺层角度㊁铺层层数对碳纤维复合材料板热冲压成形的影响ꎬ研究结果对于碳纤维复合材料零件在汽车轻量化上的应用[９]具有一定的指导意义ꎮ１㊀实验１.１㊀原材料连续碳纤维:牌号为Ｔ１１００Ｇꎬ日本东丽公司生产ꎻ聚醚醚酮(ＰＥＥＫ)树脂:牌号为ＫＴ￣８２０ＮＬꎬ熔点３４０ħꎬ索尔维集团产ꎮ１.２㊀碳纤维增强ＰＥＥＫ复合材料的制备对于碳纤维增强ＰＥＥＫ复合材料ꎬ工业上常采用热压法进行制备ꎮ其主要工艺流程为:首先将碳纤维和预处理过的ＰＥＥＫ薄膜交替铺层在清洗干净的模具中ꎬ然后将模具放入已经预热一段时间的平板硫化机内ꎬ将温度提升到一特定温度ꎬ同时加压至设定的压力ꎬ保压一段时间ꎮ待到温度自然冷却至ＰＥＥＫ树脂的玻璃化转变温度以下ꎬ完成脱模ꎬ则可以得到碳纤维增强ＰＥＥＫ复合材料板[１０－１１]ꎮ设置参数:热压温度３９５ħꎬ热压压力２.５ＭＰａꎮ１.３㊀复合材料板工程常数的拟定对于仿真建模来说ꎬ此处可以利用复合材料的工程常数来建立碳纤维复合材料板的本构模型ꎬ即从细观力学入手设置参数ꎬ最后结合宏观力学的方法来进行仿真模拟分析ꎮ复合材料工程弹性常数共有９个ꎬ即该材料的三个方向上的弹性模量(Ｅ１ꎬＥ２ꎬＥ３)㊁三个方向上的泊松比(ｖ１２ꎬｖ１３ꎬｖ２３)和三个方向上的剪切模量(Ｇ１２ꎬＧ１３ꎬＧ２３)ꎮ其中Ｅ１为主方向上的弹性模量ꎬ即本材料中的连续纤维方向ꎮ通过相关文献可以查得常温下碳纤维复合材料的Ｅ１为１５０~４５０ＧＰａ[１２]ꎮ横截面积为Ａꎬ连续纤维与基体的横截面积分别为ＡｆꎬＡｍꎬ于是有:σ１Ａ＝σｆＡｆ＋σｍＡｍ(１)式中:σ１为主方向上的应力ꎻσｆ为连续纤维截面应力ꎻσｍ为基体的应力ꎮ考虑到应力与模量的关系ꎬ则有:Ｅ１＝ＥｆＡｆＡ＋ＥｍＡｍＡ(２)式中:Ｅｆ为连续纤维截面的弹性模量ꎻＥｍ为基体的弹性模量ꎮ其他工程常数可选用ＰＥＥＫ材料的基本参数[１３－１５]ꎮ因此ꎬ拟定采用的复合材料工程常数如表１所示ꎮ表１㊀碳纤维增强复合板材工程常数Ｔａｂ.１㊀Ｅｎｇｉｎｅｅｒｉｎｇｃｏｎｓｔａｎｔｏｆｃａｒｂｏｎｆｉｂｅｒｒｅｉｎｆｏｒｃｅｄｃｏｍｐｏｓｉｔｅｐｌａｔｅ项目参数Ｅ１/ＭＰａ３２８００.０Ｅ２/ＭＰａ４８.３Ｅ３/ＭＰａ４８.３ｖ１２０.４ｖ１３０.４ｖ２３０.４Ｇ１２/ＭＰａ１１.４Ｇ１３/ＭＰａ１１.４Ｇ２３/ＭＰａ１１.４１.４㊀基于ＡＢＡＱＵＳ的冲压成形仿真建模方案在ＡＢＡＱＵＳ中进行几何建模ꎬ建立半球冲压仿真的模型ꎬ设定Ｚ方向为冲压方向ꎮ将冲头㊁压边框以及下模视为刚体ꎮ进行网格划分时ꎬ对于碳纤维复合材料板ꎬ选择了常规壳和连续壳两种复合层类型ꎬ其网格单元选择Ｓｈｅｌｌ中的Ｓ４Ｒ(常规壳)和ＣｏｎｔｉｎｕｕｍＳｈｅｌ１中的ＳＣ８Ｒ(连续壳)单元类型[１６－１９]ꎬ其余部件如冲头等采取仿真分析中常用的Ｃ３Ｄ１０Ｍ单元类型ꎮ为了保证精度ꎬ板材采用细化网格ꎬ同时为了提高运算效率ꎬ对冲头㊁压边框和下模采用粗网格划分ꎬ如图１ꎮ图１㊀热冲压模型示意Ｆｉｇ.１㊀Ｓｋｅｔｃｈｏｆｓｔａｍｐｉｎｇｍｏｄｅｌ㊀㊀设定碳纤维增强ＰＥＥＫ复合材料板处于２００ħ适宜冲压的状态ꎬ选用实验试件为厚度０.３ｍｍꎬ边长为５００ｍｍ的正方形材料板ꎬ用于成形热冲压深度为１００ｍｍ的半球形拉伸件ꎮ复合层纤维铺层夹角选取４种方式ꎬ即０ʎ㊁间隔３０ʎ㊁间隔４５ʎ( 米 字型铺层)与９０ʎ( 十 字型铺层)排布ꎮ复合层的层数与每层厚度由被冲压材料板的总厚度确定ꎬ即总厚度不变ꎬ单层材料厚度与层数成反比ꎮ冲压板总厚度为０.３ｍｍꎬ复合层层数设置为６ꎬ８ꎬ１２ꎬ２０ꎬ３０层５个等级ꎮ２㊀结果与讨论２.１㊀复合层的类型选择对仿真结果的影响针对复合层类型ꎬ进行了两种建模方式ꎬ即常规壳￣三维壳单元建模和连续壳￣三维实体单元建模ꎬ其应力云图分别见图２㊁图３ꎮ图２㊀４５ʎ复合层类型的应力云图对比Ｆｉｇ.２㊀Ｓｔｒｅｓｓｎｅｐｈｏｇｒａｍｓｃｏｍｐａｒｉｓｏｎｏｆｔｈｅｔｙｐｅｏｆ４５ʎｃｏｍｐｏｓｉｔｅｌａｙｅｒ７１第２期㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀张华伟等.基于ＡＢＡＱＵＳ的碳纤维复合材料板热冲压成形仿真图３㊀９０ʎ复合层类型的应力云图对比Ｆｉｇ.３㊀Ｓｔｒｅｓｓｎｅｐｈｏｇｒａｍｓｃｏｍｐａｒｉｓｏｎｏｆｔｈｅｔｙｐｅｏｆ９０ʎｃｏｍｐｏｓｉｔｅｌａｙｅｒ㊀㊀从两种建模方式的应力云图对比分析可以看出来ꎬ使用常规壳￣三维壳单元来进行建模的时候ꎬ与连续壳￣三维实体单元的应力数据几乎没有区别ꎬ但显示效果并没有连续壳￣三维实体单元的显示效果好ꎮ由此可知ꎬ在进行复合层建模的时候ꎬ常规壳对于碳纤维增强复合材料的连续性的因素影响没有连续壳好ꎬ应力分析是按照常规材料的应力计算方式进行的ꎬ所以在某些位置会出现没有任何过度的应力突变的情况ꎮ当改成连续壳￣三维实体单元建模的时候ꎬ由于考虑到了碳纤维复合材料的特性ꎬ在ＡＢＡＱＵＳ中有不同的计算方式ꎬ使得对于材料板的应力计算过度更加平滑ꎬ因此也不会出现类似于常规壳的那种应力突变点ꎮ因此ꎬ后文中复合层板间铺层夹角与板间层数的讨论均在建模类型为 连续壳￣三维实体单元 的基础上进行ꎮ２.２㊀复合材料的铺层夹角对冲压受力的影响碳纤维复合材料的主要受力元素是碳纤维ꎬ基体在受力过程中起到传递载荷的作用ꎮ碳纤维复合材料板冲压工艺主要受材料拉应力性能指标的影响ꎬ即可以通过分析比较最大拉应力来判断碳纤维受力情况ꎮ同时为了进一步探究复合材料层间夹角排布对其热冲压性能的影响ꎬ还设置了一组不同复合材料板层间夹角的真实应变对比ꎮ实验中以１２层复合材料板为例ꎬ其应力云图见图４ꎬ其应变云图见图５ꎮ结合图４ꎬ图５和表２ꎬ表３可分析出ꎬ当其他条件相同ꎬ当碳纤维复合材料板层间夹角不同的时候ꎬ其受到的冲压应力与应变是不同的ꎮ当层间夹角为９０ʎ叠层时ꎬ其碳纤维排列结构简单ꎬ层间滑移大ꎬ剪切柔度高ꎬ所以纤维夹角的可转角度也很大ꎻ当层间夹角为３０ʎ间隔与４５ʎ间隔分布时ꎬ由于纤维排布较为复杂ꎬ多个方向角都有限制截面剪切的纤维阻力ꎬ所以剪切柔度小ꎬ导致冲压深度相同的情况下ꎬ其纤维转角较小ꎬ即产生较大的应力ꎮ同时可观察得到ꎬ在层间夹角为３０ʎ或４５ʎ排布的时候ꎬ由于碳纤维材料在板材的几何中心排布最为密集ꎬ理论上在几何中心的受力应该是最大的ꎬ其次才是下模与板材接触的区域ꎬ即云图显示其主要受力区域在板材与冲头接触的中心位置ꎻ而间隔９０ʎ排布则不同ꎬ碳纤维在几何分布上是均匀的ꎬ不会出现某一区域的受力性能明显高于其他区域的情况ꎬ所以其主要受力区域为下模与板材的接触位置[２０－２２]ꎮ不难看出ꎬ无论对复合板中的纤维层夹角如何排布ꎬ都无法消除复合板冲压过程中的各向异性ꎮ图４㊀不同层间夹角的１２层复合材料板应力云图示意Ｆｉｇ.４㊀Ｓｔｒｅｓｓｎｅｐｈｏｇｒａｍｓｏｆ１２￣ｌａｙｅｒｃｏｍｐｏｓｉｔｅｓｈｅｅｔｗｉｔｈｄｉｆｆｅｒｅｎｔｉｎｔｅｒｌａｙｅｒａｎｇｌｅｓ图５㊀不同层间夹角的１２层复合材料板应变云图示意Ｆｉｇ.５㊀Ｓｔｒａｉｎｎｅｐｈｏｇｒａｍｓｏｆ１２￣ｌａｙｅｒｃｏｍｐｏｓｉｔｅｓｈｅｅｔｗｉｔｈｄｉｆｆｅｒｅｎｔｉｎｔｅｒｌａｙｅｒａｎｇｌｅｓ８１㊀合㊀成㊀纤㊀维㊀工㊀业㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０１９年第４２卷表２㊀不同层间夹角及层数的复合材料板的冲压应力Ｔａｂ.２㊀Ｓｔａｍｐｉｎｇｓｔｒｅｓｓｏｆｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｐｌａｔｅｓｗｉｔｈｄｉｆｆｅｒｅｎｔｉｎｔｅｒｌａｙｅｒａｎｇｌｅｓ层间夹角/(ʎ)冲击应力/ＭＰａ６层８层１２层２０层３０层０１０３.７１０６.１１１０.３８９.４８９.１３０８７２.４８７４.８９２０.２９９８.９９２３.６４５９０８.１９４９.０９３４.９９８２.９９６７.３９０１７０.３１０１.１１０８.５１１５.７１１４.６表３㊀不同层间夹角及层数的复合材料板的冲压应变Ｔａｂ.３㊀Ｓｔａｍｐｉｎｇｓｔｒａｉｎｏｆｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｐｌａｔｅｓｗｉｔｈｄｉｆｆｅｒｅｎｔｉｎｔｅｒｌａｙｅｒａｎｇｌｅｓａｎｄｌａｙｅｒｎｕｍｂｅｒｓ层间夹角/(ʎ)冲压应变６层８层１２层２０层３０层００.０９９２０.１００２０.１０１２０.０９９９０.０９９４３００.０６４２０.０３５９０.０３０００.０３１９０.０２９３４５０.０４２５０.０３０４０.０２９９０.０３０１０.０２９９９００.１００９０.１００２０.１０１５０.１０２５０.１０３４㊀㊀对比０ʎꎬ３０ʎꎬ４５ʎꎬ９０ʎ层间夹角的４种不同纤维排布方式的应变云图可以看出ꎬ 十 字型排布方式的应变与 米 字形排布方式的应变有着本质的区别ꎬ即其发生应变的区域完全不同ꎬ而其应变的大小也是差别很大ꎬ９０ʎ铺层的复合材料应变较大ꎬ达到了１０％的形变程度ꎻ而４５ʎ与３０ʎ铺层的复合材料应变较小ꎬ仅有３％的形变ꎮ这些现象都说明了复合材料板的不同板间夹角排布会给冲压工艺带来不同的受力情况ꎬ随着其层间夹角排布的不同ꎬ其纤维层之间剪切转动方式㊁复合材料板的宏观力学性能等均会产生较大的差异[２２]ꎮ一般说来ꎬ９０ʎ夹角的排布方式时其冲压受力与成形性能要好于另外两种ꎮ但这种排布方式伴随而来的是较大的应变幅度ꎬ这也意味着这种排布方式较另外两种来说更容易产生受力破坏[２３]ꎮ这个结论符合碳纤维复合材料的力学成形原理ꎮ除了这四种为工程上常用的铺层方式以外ꎬ本文还进行了单层０ʎ夹角的复合材料板的热冲压仿真ꎮ最后的现象也合理地解释了复合材料板冲压成形的原理ꎬ即当碳纤维受力后ꎬ其纤维方向上几乎没有拉伸形变ꎬ而是产生位移ꎮ所以在宏观上ꎬ板材显示出来的就是纤维方向上的板材的 收束 形变ꎮ２.３㊀复合材料的板间层数对冲压工艺的影响本研究设置了５组铺层层数ꎬ设定总厚度一定ꎬ即铺层层数越多ꎬ每一铺层就越薄ꎮ此处拟定最大铺层数为３０ꎬ当铺层数超过３０的时候ꎬ每一层的厚度均降至０.０１ｍｍ以下ꎬ已经不符合常规的生产要求ꎮ同时ꎬ由于本研究主要讨论的是碳纤维复合材料的纤维排布方向与复合材料板间层数对热冲压的影响ꎬ所以对各方向纤维占比受力的影响并未给予考虑ꎮ由１２层复合材料板冲压应力云图(图４)及表２㊁表３分析可知ꎬ当总厚度一定时ꎬ板间层数相对于铺层方式来说ꎬ对冲压性能的影响不大ꎮ其中细微差别可以通过纵向对比来发现:冲压过程中的应力一般随着板间铺层的层数上升而上升ꎮ其主要原理是由于铺层层数变多ꎬ板料的稳定性也会越来越好ꎬ而且由于层数增多ꎬ使得其总层间摩擦力增大ꎬ这也会对冲压受力的情况带来一定的影响ꎮ所以在相同的冲压位移的情况下ꎬ其受力变得更高ꎮ因此也可以认为板间层数对碳纤维增强复合材料板的冲压性能影响较小ꎬ与板间连续纤维排布角度相比ꎬ板间层数的影响可以忽略不记ꎮ同时ꎬ因为复合材料板板件的厚度对板材的力学性能有很大的影响ꎬ所以相比于复合材料板间的复合层层数ꎬ复合材料板的厚度对于热冲压的性能影响更大ꎮ３㊀结论ａ.以ＡＢＡＱＵＳ为建模工具ꎬ利用工程常数建模同宏观力学结合的方法构建了碳纤维复合材料高温本构模型ꎮｂ.在选择复合层类型的时候ꎬ 连续壳 类型的仿真效果要明显好于常规壳ꎬ更加符合实际的受力情况ꎮｃ.碳纤维增强复合材料的热冲压性能主要取决于复合板层间的夹角排布ꎬ当夹角排布的方式越松散ꎬ纤维之间剪切柔度越高ꎬ则其在高温状态下冲压受力情况就越好ꎻ反之ꎬ当夹角排布的方式越稳定ꎬ则其在高温状态下的冲压受力情况就越差ꎮ其中综合看来是９０ʎ夹角排布时冲压受力较好ꎮｄ.在同样的温度和夹角排布的情况下ꎬ碳纤维增强复合材料板的板间层数对于冲压受力性能的影响不大ꎮ参㊀考㊀文㊀献[１]㊀沈观林ꎬ胡更开ꎬ刘彬.复合材料力学[Ｍ].北京:清华大学出版社ꎬ２０１３:４－２０.９１第２期㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀张华伟等.基于ＡＢＡＱＵＳ的碳纤维复合材料板热冲压成形仿真[２]㊀张琦ꎬ高强ꎬ赵升吨.碳纤维复合材料板热冲压成形试验研究[Ｊ].机械工程学报ꎬ２０１２ꎬ４８(１８):７２－７４. [３]㊀堵同亮ꎬ彭雄奇ꎬ郭早阳ꎬ等.碳纤维编织复合材料冲压成形实验与仿真分析[Ｊ].功能材料ꎬ２０１２ꎬ１６(４４):２４０１－２４０４.[４]㊀代少俊.高性能纤维复合材料[Ｍ].上海:华东理工大学出版社ꎬ２０１３:４－９.[５]㊀罗益锋.碳纤维复合材料的研发方向和市场开发动向[Ｊ].高科技纤维与应用ꎬ２０１１ꎬ２２(１１):３８－４１. [６]㊀ＺｈａｎｇＱｉꎬＣａｉＪｉｎꎬＧａｏＱｉａｎｇ.Ｓｉｍｕｌａｔｉｏｎａｎｄｅｘｐｅｒｉｍｅｎｔａｌｓｔｕｄｙｏｎｔｈｅｒｍａｌｄｅｅｐｄｒａｗｉｎｇｏｆｃａｒｂｏｎｆｉｂｅｒｗｏｖｅｎｃｏｍｐｏｓｉｔｅｓ[Ｊ].ＪＭａｔｅｒＰｒｏｃＴｅｃｈꎬ２０１４ꎬ２１４(４):８０２－８１０. [７]㊀ＺｈａｎｇＱｉꎬＣａｉＪｉｎꎬＧａｏＱｉａｎｇ.Ｅｘｐｅｒｉｍｅｎｔａｌａｎｄｓｉｍｕｌａｔｉｏｎｒｅ￣ｓｅａｒｃｈｏｎｔｈｅｒｍａｌｓｔａｍｐｉｎｇｏｆｃａｒｂｏｎｆｉｂｅｒｃｏｍｐｏｓｉｔｅｓｈｅｅｔ[Ｊ].ＴｒａｎｓＮｏｎｆｅｒｒｏｕｓＭｅｔＳｏｃＣｈｉｎꎬ２０１４ꎬ２４(１):２１７－２２３.[８]㊀庄靖东.聚醚醚酮板材热成型性能研究[Ｄ].武汉:华中科技大学ꎬ２０１５.[９]㊀彭孟娜ꎬ马建伟.碳纤维及其在汽车轻量化中的应用[Ｊ].合成纤维工业ꎬ２０１８ꎬ４１(１):５３－５７.[１０]支建海ꎬ钱鑫ꎬ张永刚ꎬ等.国产碳纤维增强树脂基复合材料的界面结合性能研究[Ｊ].合成纤维工业ꎬ２０１８ꎬ４１(４):１４－１７.[１１]张照.碳纤维织物增强聚醚醚酮基(ＣＦＦ/ＰＥＥＫ)航空复合材料的制备及其界面改性[Ｄ].上海:东华大学ꎬ２０１７. [１２]唐见茂.高性能纤维及复合材料[Ｍ].北京:化学工业出版社ꎬ２０１２:１２－１９.[１３]张少实ꎬ庄茁.复合材料与粘弹性力学[Ｍ].北京:机械工业出版社ꎬ２０１１:５３－５５.[１４]庄靖东ꎬ黄志高ꎬ周华民.热成型条件下ＰＥＥＫ力学行为研究与建模[Ｊ].塑料工业ꎬ２０１５ꎬ４３(７):７３－７７.[１５]黄发荣ꎬ周燕.先进树脂基复合材料[Ｍ].北京:化学工业出版社ꎬ２００８:１６８－１６９.[１６]庄茁ꎬ由小川ꎬ廖剑晖ꎬ等.基于ＡＢＡＱＵＳ的有限元分析和应用[Ｍ].北京:清华大学出版社ꎬ２００９:１７－２４. [１７]刘展.ＡＢＡＱＵＳ有限元分析从入门到精通[Ｍ].北京:人民邮电出版社ꎬ２０１５:１９４－２０６.[１８]齐威.ＡＢＡＱＵＳ６.１４超级学习手册[Ｍ].北京:人民邮电出版社ꎬ２０１６:４９２－４９８.[１９]张建伟.ＡＢＡＱＵＳ有限元分析￣从入门到精通[Ｍ].北京:机械工业出版社ꎬ２０１５:２２－７９.[２０]ＣｈｅｎＱｉａｎｑｉａｎꎬＢｏｉｓｓｅＰꎬＰａｒｋＣＨꎬｅｔａｌ.Ｉｎｔｒａ/ｉｎｔｅｒ￣ｐｌｙｓｈｅａｒｂｅｈａｖｉｏｒｓｏｆｃｏｎｔｉｎｕｏｕｓｆｉｂｅｒｒｅｉｎｆｏｒｃｅｄｔｈｅｒｍｏｐｌａｓｔｉｃｃｏｍｐｏｓ￣ｉｔｅｓｉｎｔｈｅｒｍｏｆｏｒｍｉｎｇｐｒｏｃｅｓｓｅｓ[Ｊ].ＣｏｍｐＳｔｒｕｃｔꎬ２０１１ꎬ９３(７):１６９２－１７０３.[２１]㊀张衡ꎬ严飙ꎬ龚友坤ꎬ等.碳纤维机织物增强热塑性树脂复合材料热冲压叠层模型[Ｊ].复合材料学报ꎬ２０１７ꎬ３４(１２):２７４１－２７４６.[２２]丁纺纺ꎬ彭雄奇.复合材料用机织物非正交本构模型的半球形冲压成型验证[Ｊ].复合材料学报ꎬ２０１１ꎬ２８(１):１５６－１６０.[２３]韩宾ꎬ王宏ꎬ于杨惠文ꎬ等.碳纤维增强热塑性复合材料盒形件热冲压成型研究[Ｊ].航空制造技术ꎬ２０１７(１６):４０－４５.ＳｉｍｕｌａｔｉｏｎｏｆｔｈｅｒｍａｌｓｔａｍｐｉｎｇｆｏｒｍａｔｉｏｎｐｒｏｃｅｓｓｏｆｃａｒｂｏｎｆｉｂｅｒｃｏｍｐｏｓｉｔｅｓｈｅｅｔｂａｓｅｄｏｎＡＢＡＱＵＳＺｈａｎｇＨｕａｗｅｉꎬＬｉＢｏｈｏｎｇ(ＳｃｈｏｏｌｏｆＣｏｎｔｒｏｌＥｎｇｉｎｅｅｒｉｎｇꎬＮｏｒｔｈｅａｓｔｅｒｎＵｎｉｖｅｒｓｉｔｙａｔＱｉｎｈｕａｎｇｄａｏꎬＱｉｎｇｈｕａｎｇｄａｏ０６６０００)Ａｂｓｔｒａｃｔ:Ｔｈｅｔｈｅｒｍａｌｓｔａｍｐｉｎｇｐｒｏｃｅｓｓｏｆｃａｒｂｏｎｆｉｂｅｒｒｅｉｎｆｏｒｃｅｄｐｏｌｙｅｔｈｅｒｅｔｈｅｒｋｅｔｏｎｅ(ＰＥＥＫ)ｃｏｍｐｏｓｉｔｅｓｈｅｅｔｗａｓｓｉｍｕｌａ￣ｔｅｄｂａｓｅｄｏｎＡＢＡＱＵＳｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓｓｏｆｔｗａｒｅ.Ｔｈｅｍｅｃｈａｎｉｃａｌａｎｄｄｅｆｏｒｍａｔｉｏｎｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｓｗｅｒｅａｎａｌｙｚｅｄｄｕｒｉｎｇｔｈｅｒｍａｌｓｔａｍｐｉｎｇ.Ｔｈｅｅｆｆｅｃｔｓｏｆｔｈｅｆｉｂｅｒｉｎｔｅｒｌａｙｅｒａｎｇｌｅａｎｄｔｈｅｎｕｍｂｅｒａｎｄｔｈｅｔｙｐｅｏｆｃｏｍｐｏｓｉｔｅｌａｙｅｒｓｏｎｔｈｅｆｏｒｍａｂｉｌｉｔｙｗｅｒｅｄｉｓｃｕｓｓｅｄ.Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔｔｈｅｃｏｎｓｔｉｔｕｔｉｖｅｍｏｄｅｌｏｆｃｏｍｐｏｓｉｔｅｓｃｏｕｌｄｂｅｅｓｔａｂｌｉｓｈｅｄｂａｓｅｄｏｎｅｎ￣ｇｉｎｅｅｒｉｎｇｃｏｎｓｔａｎｔｓｏｆｃｏｍｐｏｓｉｔｅｓꎻｔｈｒｅｅ￣ｄｉｍｅｎｓｉｏｎａｌｓｏｌｉｄ￣ｃｏｎｔｉｎｕｏｕｓｓｈｅｌｌｗａｓａｍｏｒｅｒｅａｓｏｎａｂｌｅｃｏｍｐｏｓｉｔｅｓｈｅｅｔｔｙｐｅｉｎｌｉｎｅｗｉｔｈｔｈｅａｃｔｕａｌｓｉｔｕａｔｉｏｎꎬｐｒｏｖｉｄｉｎｇａｂｅｔｔｅｒｓｉｍｕｌａｔｉｏｎｅｆｆｅｃｔｔｈａｎｏｔｈｅｒｃｏｍｐｏｓｉｔｅｓｈｅｅｔｔｙｐｅｓꎻｔｈｅｆｉｂｅｒｉｎｔｅｒｌａｙｅｒａｎｇｌｅｏｆｃｏｍ￣ｐｏｓｉｔｅｓｗａｓｔｈｅｍａｉｎｆａｃｔｏｒａｆｆｅｃｔｉｎｇｔｈｅｔｈｅｒｍａｌｓｔａｍｐｉｎｇｐｒｏｃｅｓｓꎻｃｏｍｐａｒｅｄｗｉｔｈｏｔｈｅｒｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｓｗｉｔｈｔｈｅｆｉｂｅｒｉｎｔｅｒ￣ｌａｙｅｒａｎｇｌｅｓｏｆ０ʎꎬ３０ʎａｎｄ４５ʎꎬｔｈｅｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｗｉｔｈｔｈｅｉｎｔｅｒｌａｙｅｒａｎｇｌｅｏｆ９０ʎｓｈｏｗｅｄｂｅｔｔｅｒｓｔｒｅｓｓｐｅｒｆｏｒｍａｎｃｅｗｈｉｌｅｔｈｅｒｍａｌｓｔａｍｐｉｎｇａｎｄｗａｓｐｒｏｎｅｔｏｂｅｉｎｇｄａｍａｇｅｄｄｕｅｔｏｈｉｇｈｅｒｓｔｒａｉｎꎻｔｈｅｌａｙｅｒｎｕｍｂｅｒｏｆｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｈａｄｎｏｏｂｖｉｏｕｓｅｆｆｅｃｔｏｎｔｈｅｓｔｒｅｓｓａｔａｆｉｘｅｄｔｏｔａｌｔｈｉｃｋｎｅｓｓｏｆｔｈｅｃｏｍｐｏｓｉｔｅꎻａｎｄｔｈｅｓｉｍｕｌａｔｉｏｎｅｆｆｅｃｔｗａｓｏｐｔｉｍｉｚｅｄｗｈｅｎｔｈｅｔｈｅｒｍａｌｓｔａｍｐｉｎｇｍｏｄｅｌｏｆｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｗａｓｅｓｔａｂｌｉｓｈｅｄｗｉｔｈｃｏｎｔｉｎｕｏｕｓｓｈｅｌｌａｎｄｉｎｔｅｒｌａｙｅｒａｎｇｌｅｏｆ９０ʎ.Ｋｅｙｗｏｒｄｓ:ｃａｒｂｏｎｆｉｂｅｒꎻｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌꎻｔｈｅｒｍａｌｓｔａｍｐｉｎｇｆｏｒｍａｔｉｏｎꎻｆｉｎｉｔｅｅｌｅｍｅｎｔｍｏｄｅｌꎻｓｉｍｕｌａｔｉｏｎ０２㊀合㊀成㊀纤㊀维㊀工㊀业㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀２０１９年第４２卷。