ABAQUS实例分析论文
悬臂梁—有限元ABAQUS线性静力学分析实例
线性静力学分析实例——以悬臂梁为例线性静力学问题是简单且常见的有限元分析类型，不涉及任何非线性（材料非线性、几何非线性、接触等），也不考虑惯性及时间相关的材料属性。
在ABAQUS中，该类问题通常采用静态通用（Static，General）分析步或静态线性摄动（Static，Linear perturbation）分析步进行分析。
线性静力学问题很容易求解，往往用户更关系的是计算效率和求解效率，希望在获得较高精度的前提下尽量缩短计算时间，特别是大型模型。
这主要取决于网格的划分，包括种子的设置、网格控制和单元类型的选取。
在一般的分析中，应尽量选用精度和效率都较高的二次四边形/六面体单元，在主要的分析部位设置较密的种子；若主要分析部位的网格没有大的扭曲，使用非协调单元（如CPS4I、C3D8I）的性价比很高。
对于复杂模型，可以采用分割模型的方法划分二次四边形/六面体单元；有时分割过程过于繁琐，用户可以采用精度较高的二次三角形/四面体单元进行网格划分。
悬臂梁的线性静力学分析1.1 问题的描述一悬臂梁左端受固定约束，右端自由，结构尺寸如图1-1所示，求梁受载后的Mises应力、位移分布。
ν材料性质：弹性模量3=E=，泊松比3.02e均布载荷：F=103N图1-1 悬臂梁受均布载荷图1.2 启动ABAQUS启动ABAQUS有两种方法，用户可以任选一种。
（1）在Windows操作系统中单击“开始”--“程序”--ABAQUS 6.10 --ABAQUS/CAE。
（2）在操作系统的DOS窗口中输入命令：abaqus cae。
启动ABAQUS/CAE后，在出现的Start Section（开始任务）对话框中选择Create Model Database。
1.3 创建部件在ABAQUS/CAE顶部的环境栏中，可以看到模块列表：Module：Part，这表示当前处在Part（部件）模块，在这个模块中可以定义模型各部分的几何形体。
基于abaqus的ujoint有限元分析有限元分析论文大学论文
有限元分析课程论文课程名称：有限元分析论文题目：ujoint有限元分析学生班级；学生姓名：任课教师：学位类别：评分标准及分值选题与参阅资料（分值）论文内容（分值）论文表述（分值）创新性（分值）评分论文评语：总评分评阅教师: 评阅时间年月日注：此表为每个学生的论文封面，请任课教师填写分项分值基于abaqus的ujoint有限元分析摘要：万向传动装置在汽车中起到了传递扭矩的关键作用,在abaqus中导入ujoint实体模型，之后对其进行坐标系建立，wire 建立，以及各部件之间的连接关系的建立，最后对该模型施加边界条件，令其运动。
关键词：abaqus、有限元、ujoint一问题的描述对导入的ujoint在所有步骤完成后，施加力：在stepinitial：均设为0；step SPIN：doundary1：限制除UR2的所有，且把UR2值设为：pi。
在boundary2 中，限制UR1和UR3自由度。
二在abaqus中导入ujoint实体模型启动abaqus CAE，在文件下拉菜单中选择：import ，选择最终文件位置or 输入ws_connector_ujoint.py.inp打开文件ujoint。
（如下图所示）2.1 创建坐标系单机操作界面中的tool，从下拉菜单中选择datum，再出来的窗口中选择coordinate，3points。
首先选择origin，在选择x正方向，Y正方向、z正方向。
创建完成。
2.2创建VERT和CROSS之间的2坐标系。
根据 2.1所述操作步骤创建坐标系V-C 和V-G （VERT和GROUND）。
Notice：1、创建过程中为了清晰分辨，可将IN的suppress，创建完成后再将其resume。
其他同样2、在V-C和I-C中，x轴与cross转动所绕轴平行。
根据2.1所属步骤创建I-C 和I-G. 结果如图；2.3 定义connector geometry1. 2.3.1 创建disjoint型wire在选项中选择interaction，在所出现窗口中点击Create Wire Feature tool.，在所出现的窗口中选择Disjointwires，单机添加要成wire的点。
abaqus-铝合金A357切削加工有限元模拟
.2.
论文主要内容 1.铝合金A357切削加工有限元模型
刀具的儿何参数； 影 响 因 素 装夹条件； 切削参； 假 设 条 件
刀具是刚体且锋利，只考虑刀具的温度传导; 忽略加工过程中，由于温度变化引起的金相组 织及其它的化学变化; 被加工对象的材料是各向同性的； 不考虑刀具、工件的振动； 由于刀具和工件的切削厚度方向上，切削工程 中层厚不变，所以按平面应变来模拟；
1.1.3 A357与刀具材料参数
A357铝合金，密度ρ=2680Kg/m3,弹性模量E=79GP，泊松比 μ=0.33其他参数如下表：
.4.
论文主要内容 1.铝合金A357切削加工有限元模型
1.1.3 A357与刀具材料参数
A357铝合金，密度ρ=2680Kg/m3,弹性模量E=79GP，泊松比μ=0.33其他 参数如下表：
Gf pl L yd ydu
0
o
fpl
pl
uf
pl
pl
表达式中的 u pl 为等效塑性位移，它描述了当损伤开始之后裂纹 pl 变化的屈服应力，在损伤开始之前 u pl=0.在损伤开始之后 u =L pl ，L 为与积分点相关的单元特征长度，单元特征长度的定义基于单元的集 合形状，平面单元长度为积分点区域面积的平方根，而立方体单元长 度为积分点区域体积的立方根。
.4.
论文主要内容 1.铝合金A357切削加工有限元模型
1.1.2 材料失效准则
实现切屑从工件分离，本文采用的是剪切失效模型。剪切 失效模型是基于等效塑性应变在积分点的值，当损伤参数达到 1时，单元即失效，失效参数定义如下：
0 f
pl
pl
pl
式中： 为失效参数， 为等效塑性应变初始值， 为等效塑 性应变增量， 为失效应变。失效应变 设定以来于以下几个方 面：依据塑性应变率，无量纲压应力与偏应力之比p/q(p为压应 力，q为Mises应力),温度，预定义域变量。这里采用 Johnson357切削加工有限元模型
abaqus1典型实例分析
1.应用背景概述随着科学技术的发展，汽车已经成为人们生活中必不可少的交通工具。
但当今由于交通事故造成的损失日益剧增，研究汽车的碰撞安全性能，提高其耐撞性成为各国汽车行业研究的重要课题。
目前国内外许多著名大学、研究机构以及汽车生产厂商都在大力研究节省成本的汽车安全检测方法，而汽车碰撞理论以及模拟技术随之迅速发展，其中运用有限元方法来研究车辆碰撞模拟得到了相当的重视。
而本案例就是取材于汽车碰撞模拟分析中的一个小案例------------------------------------------------------- 保险杠撞击刚性墙。
2.问题描述该案例选取的几何模型是通过导入已有的*」GS文件来生成的（已经通过Solidworks 软件建好模型的），共包括刚性墙（PART-wall ）、保险杠（PART-bumpe）、平板（PART-plane）以及横梁（PART-rail ）四个部件，该分析案例的关注要点就是主要吸能部件（保险杠）的变形模拟，即发生车体碰撞时其是否能够对车体有足够的保护能力这里根据具体车体模型建立了保险杠撞击刚性墙的有限元分析模型，为了节省计算资源和时间成本这里也对保险杠的对称模型进行了简化，详细的撞击模型请参照图1所示，撞击时保险杠分析模型以2000mm/s的速度撞击刚性墙，其中分析模型中的保险杠与平板之间、平板与横梁之间不定义接触，采用焊接进行连接，对于保险杠和刚性墙之间的接触采用接触对算法来定义。
1.横梁（rail）2.平板（plane）3保险杠（bumper）4.刚性墙（wall）图碰撞模型的SolidWorks图为了使模拟结果尽可能真实，通过查阅相关资料，定义了在碰撞过程中相关的数据以及各部件的材料属性。
其中，刚性墙的材料密度为X 10-9，弹性模量为X 105,泊松比为；保险杠、平板以及横梁的材料密度为X 10-9，弹性模量为X 105,泊松比为，塑形应力-应变数据如表所示。
abaqus回弹分析实例：在分析步之间传递数据
Abaqus Example Problems Guide1.5.1 Springback of two-dimensional draw bendingProducts: Abaqus/Standard Abaqus/Explicit是用explicit进行成形，然后使用standard分析回弹Problem description这个例子描述了在93年成形模拟数值会议中报道的基准测试。
这个基准包括使用三种材料描述六种问题以及两种不同的夹持力。
六个问题中的一个描述如下。
原文见Taylor等人在93年的论文。
坯料初始尺寸350*35，厚度0.78.问题本质上是个平面应变问题（垂直于模型平面的尺寸是35mm）。
夹持力是2.45kN，夹具的质量是5kg。
摩擦系数0.144坯料为低碳钢，材料为弹塑性材料，弹性为各向同性，对于塑性使用Hill 异性屈服准则。
材料的性质如下：Young's modulus = 206.0 GPaPoisson's ratio = 0.3Density = 7800.Yield stress = 167.0 MPaAnisotropic yield criterion: =1.0, =1.0402, =1.24897, =1.07895, =1.0, =1.0此例是对称问题，只取一半建模。
坯料使用一行175个一阶壳单元建模。
对称的边界条件施加到对称面上。
边界条件施加到了坯料所有的节点上以模拟平面应变的条面；建模平面外的尺寸是5mm；因此，坯料的加持力经过了粗略缩放。
成形的过程由explicit中的两个分析步完成。
坯料的加持力在第一个分析步施加。
加载使用平滑分析步，以将惯性效应最小化。
第二个分析步，通过设定冲头刚体参考点的速度冲头下行70mm。
速度使用triangular smooth step amplitude 功能施加，初始速度和最终速度都是0，峰值速度在这个过程之间。
【完整版】Abaqus课程设计论文
学校代码：10128学号：************课程论文题目：球体接触问题的有限元分析****：***学院：理学院班级：工程力学12-1班指导教师：韦广梅、周承恩2015年9月26日内蒙古工业大学课程设计（论文）任务书课程名称：有限单元法课程设计学院：理学院班级：力学12-1班 _ 学生姓名：李劲波学号：201220907039 指导教师：韦广梅、周承恩摘要许多工程问题都涉及两个或多个部件之间的接触，在这些问题中，当两个物体彼此接触时，垂直于接触面上的力作用在两个物体上。
如果在接触面之间存在摩擦，可能产生剪力以阻止物体的切向运动。
本文主要是运用ABAQUS软件对椭球体与刚性半无限大体接触问题进行探讨，因为椭球体是前后左右都对称的三维实体，所以对它取四分之一分析，这样不仅计算方便而且结果直观。
先将材料定义为弹性材料通过前处理建模、计算、后处理来分析椭球体与大体之间的接触，运用多种方法讨论接触问题。
结果得到的数据远超弹性屈服极限，后将材料重新赋值为弹塑性，仿照弹性的分析方法对弹塑性进行分析，得到的结果完全正确。
关键词：接触；ABAQUS；弹塑性；屈服极限AbstractMany engineering problems involve contact between two or more components, in these problems, when two objects come into contact with each other, on the surface of the perpendicular to the contact force between two objects. If there is friction between contact surface, shear may be produced to prevent the tangential motion of the object. This article mainly using ABAQUS software to ellipsoid and rigid semi-infinite discusses the general contact problem, because the ellipsoid is before and after the left and right sides is symmetrical three-dimensional entity, so take a quarter analysis to it, that not only the calculation results of convenient and intuitive. Before the adoption of the first material is defined as elastic material modeling, calculation, after processing to analyze between ellipsoid and general contact, using a variety of methods to discuss contact problem. Than elastic yield limit, the resulting data will reopen assignment for elastic-plastic materials, imitates the elastic analysis method of elastic-plastic analysis, get the right result.Key words: Contact; ABAQUS; Elastic-plastic; Yield limit目录第一章绪论1.1 有限单元法课程设计目的与任务 (1)1.2 接触问题概述 (1)1.3 有限元法概述 (1)第二章有限元工具ABAQUS介绍及应用方法 (3)2.1 ABAQUS简介 (3)2.2 ABAQUS功能简介 (3)第三章椭球—平面接触问题有限元分析 (7)3.1 平面－球接触问题描述 (7)3.2 用ABAQUS建立有限元模型 (7)3.2.1 建立实体 (7)3.2.2 赋值材料与装配部件 (8)3.2.3 创建分析步 (8)3.2.4 定义接触和边界条件 (9)3.2.5 划分网格 (9)3.2.6 计算及后处理 (10)3.2.7 理想弹塑性分析 (11)第四章结论与收获 (14)参考文献 (15)第一章绪论1.1有限单元法课程设计目的与任务目的：《有限单元法课程设计》是“有限单元法”课程的拓展部分。
ABAQUS有限元分析毕业论文
1.1.1振动与冲击对电子设备的危害
在电子设备所处的机械环境中，各种机械力和干扰形式都有可能对设 备的可 靠性造成危害，其中危害最大的是振动和冲击。它们造成的危害主
要有两种：
(1)设备在某一激振频率下导致振幅很大的共振现象，最终因振动加速
度超 过设备所能承受的极限加速度而破坏；或者由于冲击所产生的冲击力
Studying the general analysis method and solving process of the phone'flip con tact and impact duing dropp ing with fin ite eleme nt an alysis simulati on tech no logy based on the con tact dyn amics will be the set up for non li near simulati on an alysis」t is be provied that simulati on method with complex structure,duri ng dropp ing impact load and comb ined with con tact dyn amics and fin ite eleme nt an alysis method.
Drop and impact are ofen see n in engin eeri ng practice and recog ni zed as a con tact collision problem,and the problem shows highly characteristics of nonlinearity. Normally it'effective method to solve dropping issue complex structure with finite element method and actual dropping test data,and it will improve the precision of finite element simulation.So the analysis and research for drop phenomena with finite element method is important in both theoretical and engineering fields.
(完整word版)Abaqus弹塑性分析简单实例
Abaqus弹塑性分析简单实例
ABAQUS默认的塑性材料特性应用金属材料的经典塑性理论，采用MISES屈服面来定义各向屈服。
金属材料的弹塑性行为可以简述如下：在小应变时，材料性质基本为线弹性，弹性模量E为常数；应力超过屈服应力后，刚度会显著下降，此时材料的应变包括塑性应变和弹性应变两部分；在卸载后，弹性应变消失，而塑性应变是不可恢复的；如果再次加载，材料的屈服应力会提高，即所谓的加工硬化。
在abaqus中，等效塑性应变PEEQ大于0表明材料发生了屈服。
在工程结构中，等效塑性应变一般不应超过材料的破坏应变。
对于金属成形等大变形问题，应根据生产工艺要求来确定许可的等效塑性应变量。
需要注意的是在比例加载时，大多数材料的PEMAG和PEEQ相等。
这两个量的区别在于，PEMAG描述的是变形过程中某一时刻的塑性应变，与加载历史无关，而PEEQ是整个变形过程中塑性应变的累积结果。
下面我们以单向压缩过程的模拟来演示ABAQUS弹塑性仿真设置。
模型如图所示，压头用解析刚体来模拟，试样用SHELL来模拟。
采用轴对称模型。
试样的截面属性设置如下图所示，注意塑性应变必须从0开始。
在压头与试样之间定义无摩擦的接触。
固定对称轴
上的径向位移U1和底边的轴向位移U2。
压头是轴对称刚体，U2边界条件需要施加在压头的参考点上。
设定两个分析步，第一个分析步让压头与试样建立平稳的接触，设置压头下移-5.001mm。
第二个分析步，设定压头下移20mm。
具体如下图所示：
提交分析，结果如下图所示：有限元在线因为专注所以卓越。
ABAQUS实例分析论文
目录第一章Abaqus简介1一、Abaqus总体介绍1二、Abaqus基本使用方法21.2.1 Abaqus分析步骤21.2.2 Abaqus/CAE界面31.2.3 Abaqus/CAE的功能模块3第二章基于Abaqus的通孔端盖分析实例4一、工作任务的明确5二、具体步骤52.2.1 启动Abaqus/CAE42.2.2 导入零件52.2.3 创建材料和截面属性62.2.4 定义装配件72.2.5 定义接触和绑定约束（tie）102.2.6 定义分析步142.2.7 划分网格152.2.8 施加载荷192.2.9 定义边界条件202.2.10提交分析作业212.2.11 后处理22第三章课程学习心得与作业体会22第一章：Abaqus简介一、Abaqus总体介绍Abaqus是功能强大的有限元分析软件，可以分析复杂的固体力学和结构力学系统，模拟非常庞大的模型，处理高度非线性问题。
Abaqus不但可以做单一零件的力学和多物理场的分析，同时还可以完成系统级的分析和研究。
Abaqus使用起来十分简便，可以很容易的为复杂问题建立模型。
Abaqus具备十分丰富的单元库，可以模拟任意几何形状，其丰富的材料模型库可以模拟大多数典型工程材料的性能，包括金属、橡胶、聚合物、复合材料、钢筋混泥土、可压缩的弹性泡沫以和地质材料（例如土壤、岩石）等。
Abaqus主要具有以下分析功能：1.静态应力/位移分析2.动态分析3.非线性动态应力/位移分析4.粘弹性/粘塑性响应分析5.热传导分析6.退火成形过程分析7.质量扩散分析8.准静态分析9.耦合分析10.海洋工程结构分析11.瞬态温度/位移耦合分析12.疲劳分析13.水下冲击分析14.设计灵敏度分析二、Abaqus基本使用方法1.2.1 Abaqus分析步骤有限元分析包括以下三个步骤：1.前处理（Abaqus/CAE）:在前期处理阶段需要定义物理问题的模型，并生成一个Abaqus输入文件。
基于abaqus仿真的电动机振动特性分析及振动故障诊断
ABSTRACT
Vibration problem of the motor is the common problem of the motor industry, in motor design phase of the vibration characteristic analysis of forecast, and the vibration problem of events can put forward reasonable solution, is the motor enterprises urgently need to solve problems in production. How to accurately calculate vibration frequency conversion motor, so as to find out the reasonable structure design, and thus improve the quality of motor design and level, strengthen the competitiveness of the enterprise, has become one of the key technology of motor manufacturing enterprises.
the Vibration Fault Diagnosis
领 域： 研 究 生： 指导教师： 企业导师：
软件工程 任鸣 王秀锦 吴德忠
天津大学软件学院 二零一三年十一月
独创性声明
本人声明所呈交的学位论文是本人在导师指导下进行的研究工作和取得的 研究成果，除了文中特别加以标注和致谢之处外，论文中不包含其他人已经发表
abaqus案例
abaqus案例Abaqus是一款由法国达索系统公司开发的有限元分析软件。
它提供了强大的分析工具和高度可定制的建模环境，使工程师能够在各种领域进行准确的仿真分析。
在本文中，我们将介绍一些使用Abaqus的真实案例，以便更好地了解该软件的功能和应用。
1. 案例一：汽车碰撞分析汽车碰撞是交通事故中最常见的一种，也是最危险的一种。
利用Abaqus进行汽车碰撞仿真分析可以帮助工程师更好地理解碰撞过程和车辆的结构变形情况。
在这个案例中，我们将对一辆小型轿车进行碰撞测试。
首先，我们需要建立一个准确的车辆模型。
通过CAD软件，我们可以设计出车辆的外形，并将其导入Abaqus中。
接下来，我们需要添加材料属性和约束条件，以便对车辆进行仿真分析。
在这个案例中，我们使用了钢材作为车辆的材料，并设置了车辆的边界条件和碰撞速度。
通过Abaqus进行汽车碰撞仿真分析后，我们可以得到车辆在碰撞过程中的应力、应变和变形情况。
此外，我们还可以观察到车辆的安全性能和结构强度是否符合标准要求。
这些结果可以帮助汽车制造商更好地设计和改进汽车结构，以提高车辆的安全性能和耐用性。
2. 案例二：桥梁结构分析桥梁是重要的基础设施之一，其结构的稳定性和安全性对人们的出行和生活具有重要影响。
利用Abaqus进行桥梁结构分析可以帮助工程师更好地了解桥梁的结构特性和受力状态，以便更好地设计和改进桥梁结构。
在这个案例中，我们将对一座悬索桥进行分析。
首先，我们需要建立一个准确的桥梁模型，并添加材料属性和约束条件。
接下来，我们需要对桥梁进行荷载分析，以便了解桥梁在不同荷载下的受力状态和变形情况。
通过Abaqus进行桥梁结构分析后，我们可以得到桥梁在不同荷载下的应力、应变和变形情况。
此外，我们还可以观察到桥梁的结构稳定性和安全性能是否符合标准要求。
这些结果可以帮助桥梁工程师更好地设计和改进桥梁结构，以提高其安全性能和耐用性。
3. 案例三：电子设备热分析电子设备的热管理是一个重要的问题，因为过热会导致设备故障或损坏。
ABAQUS裂纹分析例子
Analysis of a Deep Single Edge Notch Crack in a PlateThis example demonstrates the use of the "deep crack" capability in ZENCRACK with both sides of the crack modelled. In this type of analysis, in addition to inserting crack-blocks into the mesh, ZENCRACK separates other standard elements to form parts of the crack face.The example consists of a simple mesh of a square plate under prescribed end loading (displacement boundary condition). The full plate is modelled even though the geometry and loading is symmetric. A schematic of the geometry and loading is shown in Figure QA16-1.Figure QA16-1 - Square plate with single edge notchF.E. Mesh for the Uncracked ComponentThe finite element model for the uncracked component is shown in Figure QA16-2. The mesh consists of a 4x4x1 grid of elements. The crack is to be inserted by replacing elements 7 and 11. The nodes used to define crack orientation for each of these elements are shown in Figure QA16-3.To extend the crack faces outside the crack-blocks, it is required that elements 8 and 12 are split.ABAQUS Input Deck for the Uncracked ComponentThe ABAQUS input deck for the uncracked component is given in listing 1. The following points should be noted :1.PropertiesThe material property name ZCRMAT has been used toallow ZENCRACK to extract Young’s modulus andPoisson ratio from the uncracked mesh.ZENCRACK Input and OutputThe ZENCRACK input file is shown in listing 2. The analysis output file is shown in listing 3.The following points are specific to the "deep crack" definition :Input record 3 : 2 2 2 0.25 1 0 1Item 7, ISPLIT, is set to 1 to indicate that one split set is to be defined.Input record 3c : 1This indicates that there is 1 element pair in the split set.Input record 3d : 8 12This indicates that elements 8 and 12 are to be split.Figure QA16-2 - Uncracked mesh - element numbersFigure QA16-3 - Uncracked mesh - element 7 and 11 Nodes for crack front orientationFigure QA16-4 - Cracked meshFigure QA16-5 - Cracked mesh - displaced plotFigure QA16-6 - Cracked mesh - "lower left" region - displaced plot LISTING 1File: a16uc01.inpDescription: ABAQUS input file for uncracked component*HEADINGQA16 - 4x4x1 mesh*PREPRINT,ECHO=NO,HISTORY=NO,MODEL=NO*RESTART,WRITE*ELEMENT,TYPE=C3D20,ELSET=EALL1 12 7 6 26 27 32 31 126 127 128 129 130 131 132133 134 135 136 1372 23 8 7 27 28 33 32 384 376 385 127 386 380 387131 135 381 383 1363 34 9 8 28 29 34 33 373 374 375 376 377 378 379380 381 382 371 3834 45 10 9 29 30 35 34 388 389 365 374 390 391 368378 382 392 372 3715 26 27 32 31 51 52 57 56 130 131 132 133 302 213 224303 402 403 404 4056 27 28 33 32 52 53 58 57 386 380 387 131 211 203 212213 403 419 421 4047 28 29 34 33 53 54 59 58 377 378 379 380 200 201 202203 419 420 417 4218 29 30 35 34 54 55 60 59 390 391 368 378 219 220 192201 420 422 418 4179 51 52 57 56 76 77 82 81 302 213 224 303 238 216 227239 304 217 218 23010 52 53 58 57 77 78 83 82 211 203 212 213 214 207 215216 217 208 210 21811 53 54 59 58 78 79 84 83 200 201 202 203 204 205 206207 208 209 198 21012 54 55 60 59 79 80 85 84 219 220 192 201 221 222 195205 209 223 199 19813 76 77 82 81 101 102 107 106 238 216 227 239 240 241 242243 244 245 246 24714 77 78 83 82 102 103 108 107 214 207 215 216 294 290 295241 245 291 293 24615 78 79 84 83 103 104 109 108 204 205 206 207 287 288 289290 291 292 285 29316 79 80 85 84 104 105 110 109 221 222 195 205 296 297 282288 292 298 286 285*NODE1 0.0000 0.0000 0.00002 2.5000 0.0000 0.00003 5.0000 0.0000 0.00004 7.5000 0.0000 0.00005 10.0000 0.0000 0.00006 0.0000 2.5000 0.00007 2.5000 2.5000 0.00008 5.0000 2.5000 0.00009 7.5000 2.5000 0.000010 10.0000 2.5000 0.000026 0.0000 0.0000 2.500027 2.5000 0.0000 2.500028 5.0000 0.0000 2.500029 7.5000 0.0000 2.500030 10.0000 0.0000 2.500031 0.0000 2.5000 2.500032 2.5000 2.5000 2.500033 5.0000 2.5000 2.500034 7.5000 2.5000 2.500035 10.0000 2.5000 2.500051 0.0000 0.0000 5.000052 2.5000 0.0000 5.000053 5.0000 0.0000 5.000054 7.5000 0.0000 5.000055 10.0000 0.0000 5.000056 0.0000 2.5000 5.000057 2.5000 2.5000 5.000058 5.0000 2.5000 5.000059 7.5000 2.5000 5.000060 10.0000 2.5000 5.000076 0.0000 0.0000 7.500077 2.5000 0.0000 7.500078 5.0000 0.0000 7.500079 7.5000 0.0000 7.500080 10.0000 0.0000 7.500081 0.0000 2.5000 7.500082 2.5000 2.5000 7.500083 5.0000 2.5000 7.500084 7.5000 2.5000 7.500085 10.0000 2.5000 7.5000 101 0.0000 0.0000 10.0000 102 2.5000 0.0000 10.0000 103 5.0000 0.0000 10.0000 104 7.5000 0.0000 10.0000 105 10.0000 0.0000 10.0000 106 0.0000 2.5000 10.0000 107 2.5000 2.5000 10.0000 108 5.0000 2.5000 10.0000 109 7.5000 2.5000 10.0000 110 10.0000 2.5000 10.0000 126 1.2500 0.0000 0.0000 127 2.5000 1.2500 0.0000 128 1.2500 2.5000 0.0000129 0.0000 1.2500 0.0000 130 1.2500 0.0000 2.5000 131 2.5000 1.2500 2.5000 132 1.2500 2.5000 2.5000 133 0.0000 1.2500 2.5000 134 0.0000 0.0000 1.2500 135 2.5000 0.0000 1.2500 136 2.5000 2.5000 1.2500 137 0.0000 2.5000 1.2500 192 8.7500 2.5000 5.0000 195 8.7500 2.5000 7.5000 198 7.5000 2.5000 6.2500 199 10.0000 2.5000 6.2500 200 6.2500 0.0000 5.0000 201 7.5000 1.2500 5.0000 202 6.2500 2.5000 5.0000 203 5.0000 1.2500 5.0000 204 6.2500 0.0000 7.5000 205 7.5000 1.2500 7.5000 206 6.2500 2.5000 7.5000 207 5.0000 1.2500 7.5000 208 5.0000 0.0000 6.2500 209 7.5000 0.0000 6.2500 210 5.0000 2.5000 6.2500 211 3.7500 0.0000 5.0000 212 3.7500 2.5000 5.0000 213 2.5000 1.2500 5.0000 214 3.7500 0.0000 7.5000 215 3.7500 2.5000 7.5000 216 2.5000 1.2500 7.5000 217 2.5000 0.0000 6.2500 218 2.5000 2.5000 6.2500 219 8.7500 0.0000 5.0000 220 10.0000 1.2500 5.0000 221 8.7500 0.0000 7.5000 222 10.0000 1.2500 7.5000 223 10.0000 0.0000 6.2500 224 1.2500 2.5000 5.0000 227 1.2500 2.5000 7.5000 230 0.0000 2.5000 6.2500 238 1.2500 0.0000 7.5000 239 0.0000 1.2500 7.5000 240 1.2500 0.0000 10.0000 241 2.5000 1.2500 10.0000 242 1.2500 2.5000 10.0000243 0.0000 1.2500 10.0000 244 0.0000 0.0000 8.7500 245 2.5000 0.0000 8.7500 246 2.5000 2.5000 8.7500 247 0.0000 2.5000 8.7500 282 8.7500 2.5000 10.0000 285 7.5000 2.5000 8.7500 286 10.0000 2.5000 8.7500 287 6.2500 0.0000 10.0000 288 7.5000 1.2500 10.0000 289 6.2500 2.5000 10.0000 290 5.0000 1.2500 10.0000 291 5.0000 0.0000 8.7500 292 7.5000 0.0000 8.7500 293 5.0000 2.5000 8.7500 294 3.7500 0.0000 10.0000 295 3.7500 2.5000 10.0000 296 8.7500 0.0000 10.0000 297 10.0000 1.2500 10.0000 298 10.0000 0.0000 8.7500 302 1.2500 0.0000 5.0000 303 0.0000 1.2500 5.0000 304 0.0000 0.0000 6.2500 365 8.7500 2.5000 0.0000 368 8.7500 2.5000 2.5000 371 7.5000 2.5000 1.2500 372 10.0000 2.5000 1.2500 373 6.2500 0.0000 0.0000 374 7.5000 1.2500 0.0000 375 6.2500 2.5000 0.0000 376 5.0000 1.2500 0.0000 377 6.2500 0.0000 2.5000 378 7.5000 1.2500 2.5000 379 6.2500 2.5000 2.5000 380 5.0000 1.2500 2.5000 381 5.0000 0.0000 1.2500 382 7.5000 0.0000 1.2500 383 5.0000 2.5000 1.2500 384 3.7500 0.0000 0.0000 385 3.7500 2.5000 0.0000 386 3.7500 0.0000 2.5000 387 3.7500 2.5000 2.5000 388 8.7500 0.0000 0.0000 389 10.0000 1.2500 0.0000 390 8.7500 0.0000 2.5000391 10.0000 1.2500 2.5000392 10.0000 0.0000 1.2500402 0.0000 0.0000 3.7500403 2.5000 0.0000 3.7500404 2.5000 2.5000 3.7500405 0.0000 2.5000 3.7500417 7.5000 2.5000 3.7500418 10.0000 2.5000 3.7500419 5.0000 0.0000 3.7500420 7.5000 0.0000 3.7500421 5.0000 2.5000 3.7500422 10.0000 0.0000 3.7500*NSET,NSET=BASE1 2 3 4 5 6 7 8 9 10 126 127 128 129 365 373374 375 376 384 385 388 389*NSET,NSET=TOP101 102 103 104 105 106 107 108 109 110 240 241 242 243 282 287288 289 290 294 295 296 297*NSET,NSET=FIX151 56*NSET,NSET=FIX251*NSET,NSET=NOUTBASETOPFIX1FIX2*ELSET,ELSET=EOUT1*MATERIAL,NAME=ZCRMAT*ELASTIC200000,0.3*SOLID SECTION,ELSET=EALL,MATERIAL=ZCRMAT*BOUNDARYFIX1,1FIX2,2*STEP*STATIC*BOUNDARYTOP,3,3,0.01BASE,3,3,-0.01*EL PRINT,ELSET=EOUT*NODE PRINT,NSET=NOUT,TOTAL=YESRF*END STEPLISTING 2File: a16c05.zcrDescription: Input file for analysis of the initial configuration## ZENCRACK Q.A. Example #16## Uncracked mesh:a16uc01.inp# Job id:a16c05r0 1 0 3 2 0 02 2 2 0.25 1 0 1st111x5.sup7 54 53st111x5.sup11 59 5818 120.5 0.50.5 0.510 11LISTING 3File: a16c05r.repDescription: Output file for analysis of the initial configuration***************************************************************************** ** ZENCRACK is proprietary software of: ** ** ZENTECH INTERNATIONAL LTD., ** 103 Mytchett Road, ** Mytchett, ** Camberley, ** Surrey, ** GU16 6ES, ** GREAT BRITAIN ** Tel. (+44) 1252 376388 ** Fax. (+44) 1252 376389 ***********************.uk ** Web ** ****************************************************************************** ** The following system has been licensed: ** ** Program : ZENCRACK ** Version : 6.0 ** Release date : 09-NOV-1998 ** License type : Commercial ** Interface to : MARC & ABAQUS **** This system is licensed to: ** ** Company : Zentech International Limited ** Address : 103 Mytchett Road ** Camberley ** Surrey GU16 6ES ** U.K. ** Phone : +44 1252 376388 ** Fax : +44 1252 376389 ** Computer : PC ** ******************************************************************************** .zcr input file for analysis : ** E:\zencrack\qa16\a16c05 ** ** Date of analysis (D-M-Y) : 21-11-1998 ** Start time (H-M-S) : 15:18:40 ** *****************************************************************************INPUT RECORD 1***************UNCRACKED MESH : a16uc01.inp JOBNAME IS : E:\zencrack\qa16\a16c05rINPUT RECORD 2***************NUMBER OF DEBUG MODULES = 0WARNING FLAG = 1DISTORTION FLAG = 0F.E. CODE FLAG = 3RUN TYPE = 2READ RESTART = 0WRITE RESTART = 0PROGRAM WILL NOT STOP IF A WARNING IS ISSUEDELEMENT DISTORTION PARAMETERS WILL NOT BE CALCULATEDF.E. ANALYSIS USING ABAQUS 5.4 OR LATERNO DATA WILL BE READ FROM A RESTART FILENO DATA WILL BE WRITTEN TO RESTART FILESRUN TYPE 2 INDICATES :- ANALYSIS OF THE INITIAL CRACK ONLY OR FULL CRACK GROWTH ANALYSIS- EACH FE ANALYSIS SHOULD BE LINEAR ELASTIC- LOADING IS PROPORTIONAL(i.e. ENERGY RELEASE RATES CAN BE PROCESSED FROM ONLY THE FIRSTANALYSIS INCREMENT)INPUT RECORD 3***************NUMBER OF CRACK-BLOCKS TO BE INCLUDED = 2FLAG FOR INTERMEDIATE STEP DURING MAPPING OF CRACK-BLOCKS INTO THEMESHICONFO = 2 FLAGS TWO STAGE MAPPING OF CRACK-BLOCKS INTO THE MESHFLAG FOR COLLAPSED ELEMENT MODIFICATION AFTER MAPPINGIFIXER = 2 FLAGS LOCAL MODIFICATION OF CRACK FRONT ELEMENTS AFTERMAPPING"QUARTER POINT NODE" PARAMETER. VALID RANGE IS 0.25 TO 0.5,WITH 0.25 FOR QUARTER-POINT NODES AND 0.5 FOR MIDSIDE NODES.VALUE READ FOR QPR IS 0.2500VALUE USED FOR QUARTER POINT PARAMETER IS 0.2500FLAG FOR GENERATING EXTRA CRACK FRONT ELEMENT SETSICFSET HAS BEEN SET TO 1 - EXTRA ELEMENT SETS WILL BE GENERATED. SETSARE :CFACEn FOR ELEMENTS ON THE FACE(S) OF CRACK FRONT nCFRONTn FOR CRACK FRONT ELEMENTS OF CRACK FRONT nFLAG FOR GENERATING EXTRA CRACK RELATED NODE SETSINDSET HAS NOT BEEN SET TO 1,2 OR 12 - NO EXTRA NODE SETS WILL BEGENERATEDNUMBER OF SETS OF ELEMENTS TO BE SPLIT FOR CRACK FACES OF DEEP CRACKS= 1BIAS FACTOR NOT RELEVANT - ICONFO NOT SET TO 3USING DEFAULT VALUE FOR CONMIDF : 0.90000***WARNINGVALUE OF QUARTER POINT NODE PARAMETER HAS BEEN RESET TO 0.26TO AVOID POSSIBLE NUMERICAL PROBLEMS (INSIDE OUT ELEMENTS) WHENEXACTLY 0.25 IS USED.MINIMUM VALUE ALLOWED WITHOUT RESETTING IS 0.252.USE -0.25 TO FORCE A VALUE OF EXACTLY 0.25.CRACK-BLOCK NUMBER 1READ FROM FILE p:\009-zencrack\temptest\install\crack\st111x5.sup THIS CRACK-BLOCK REPLACES ORIGINAL ELEMENT 7CRACK PLANE IN MESH DEFINED BY CORNER NODES 54 53CRACK-BLOCK NUMBER 2READ FROM FILE p:\009-zencrack\temptest\install\crack\st111x5.sup THIS CRACK-BLOCK REPLACES ORIGINAL ELEMENT 11CRACK PLANE IN MESH DEFINED BY CORNER NODES 59 58SETS FOR SPLITTING ELEMENTS ON DEEP CRACK FACES1 PAIR OF ELEMENTS DEFINED FOR SET 1ELEMENT PAIR FOR THIS SET IS :8 & 12INPUT RECORD 4***************CRACK-BLOCK NUMBER 1RATIO OF CRACK LENGTH TO CRACK-BLOCK EDGE 1 LENGTH IS .50000E+00RATIO OF CRACK LENGTH TO CRACK-BLOCK EDGE 2 LENGTH IS .50000E+00RING SIZE NOT RELEVANT - ICONFO NOT SET TO 3CRACK-BLOCK NUMBER 2RATIO OF CRACK LENGTH TO CRACK-BLOCK EDGE 1 LENGTH IS .50000E+00RATIO OF CRACK LENGTH TO CRACK-BLOCK EDGE 2 LENGTH IS .50000E+00RING SIZE NOT RELEVANT - ICONFO NOT SET TO 3INPUT RECORD 5***************VALUE OF IPROG = 1NO USER SUBROUTINE REQUIREDINPUT RECORD 6***************NO FATIGUE CRACK GROWTH PREDICTIONS ARE REQUESTEDINPUT RECORD 10***************FACTORS ON APPLIED LOAD TO GIVE CYCLIC LOAD RANGEFACTOR TO GIVE MINIMUM LOAD = 0.00000E+00FACTOR TO GIVE MAXIMUM LOAD = 0.10000E+01INPUT RECORD 12***************FLAG FOR ABAQUS CONTOUR SELECTION FOR PROCESSINGNUCONT = 1 : USE CONTOUR 1 TO CALCULATE GmaxABAQUS CONTOUR INTEGRAL OPTIONJITYPE = 0 : *CONTOUR INTEGRAL OPTION SPECIFIES J-INTEGRALCALCULATIONSJVCE = 0 : *CONTOUR INTEGRAL WILL SPECIFY CRACK EXTENSIONDIRECTIONSNUMBER OF CONTOURS TO BE EVALUATED BY ABAQUSNCONT IS SET TO THE DEFAULT VALUE - 3 CONTOURS*NORMAL AT ENDS OF CRACK FRONTSJSNORM = 0 : *NORMAL AT ENDS OF CRACKS WILL BE USED***WARNINGUSE OF CONTOUR 1 VALUES IS NOT RECOMMENDEDYOU SHOULD COMPARE CONTOUR 1 VALUES WITH RESULTS FOR HIGHER CONTOURSFINAL INPUT RECORD******************ELEMENT DISTORTION PARAMETERS WILL NOT BE CALCULATED (IDIST=0 ON CARD2)--------------------------------------------------------------------------------END OF READ FROM .ZCR FILE--------------------------------------------------------------------------------***INFORMATION : WARNINGS WHILE READING .ZCR INPUT DATA2 warning messages generated while reading .ZCR data***INFORMATION : MATERIAL PROPERTY DATAZCRMAT material property defined on *MATERIAL input :*MATERIAL,NAME=ZCRMAT*ELASTICYoungs modulus= 0.200000E+06 Poisson ratio= 0.30000Temperature= 0.00***INFORMATION : CRACK FRONT DEFINITIONCrack front 1 Side 1 Side 2C.B. Element C.B. Element1 72 11***INFORMATION : NODES FOR SPLIT SET 1ORIGINAL NODE NODE ON OPPOSITE FACE INITIAL COORDINATES (X,Y,Z)54 1633 0.75000E+01 0.00000E+000.50000E+0155 1634 0.10000E+02 0.00000E+000.50000E+0159 1635 0.75000E+01 0.25000E+010.50000E+0160 1636 0.10000E+02 0.25000E+010.50000E+01192 1637 0.87500E+01 0.25000E+010.50000E+01201 1638 0.75000E+01 0.12500E+010.50000E+01219 1639 0.87500E+01 0.00000E+000.50000E+01220 1640 0.10000E+02 0.12500E+010.50000E+01***INFORMATION : CRACK FRONT NODE LISTSCrack front number 1472 541 588 636 673 724 777 816 861908 973***INFORMATION : CONTOURS AND Gmax EVALUATION FOR EACH CRACK FRONTCrack front 1 has 3 contours : Gmax evaluated using contour 1***INFORMATION : ANALYSIS DESCRIPTIONSingle finite element analysis (IRUN=2) with no fatigue crack growth--------------------------------------------------------------------------------ZENCRACK 6.0 : Zentech International Ltd. Date of analysis :21-11-1998Licensed to : Zentech International LimitedLicense type : CommercialInput ref. : E:\zencrack\qa16\a16c05 Output ref. : E:\zencrack\qa16\a16c05r --------------------------------------------------------------------------------RESULTS FOR FINITE ELEMENT ANALYSIS NUMBER 1--------------------------------------------------------------------------------NEW MESH WRITTEN.SYSTEM TIME IS (H,M,S) 15:18:47ELAPSED TIME IS (H,M,S) 0: 0: 7J-INTEGRAL VALUES FROM ABAQUS OUTPUT FILE, SET 1**************************************************ANALYSIS USED ABAQUS VERSION 5.8-1CRACK SET NODE CONTOUR 1 CONTOUR 2 CONTOUR 31 1 1 0.16840E+01 0.16330E+01 0.16450E+011 12 0.19440E+01 0.19600E+01 0.19550E+011 1 3 0.20310E+01 0.19690E+01 0.19740E+011 1 4 0.19990E+01 0.20210E+01 0.20200E+011 1 5 0.21040E+01 0.20360E+01 0.20370E+011 1 6 0.20210E+01 0.20440E+01 0.20430E+011 1 7 0.21030E+01 0.20360E+01 0.20370E+011 1 8 0.19990E+01 0.20210E+01 0.20200E+011 1 9 0.20300E+01 0.19690E+01 0.19740E+011 1 10 0.19450E+01 0.19600E+01 0.19550E+011 1 11 0.16830E+01 0.16340E+01 0.16450E+011 2 1 0.25790E+01 0.25020E+01 0.25200E+011 2 2 0.29790E+01 0.30030E+01 0.29950E+011 2 3 0.31110E+01 0.30170E+01 0.30240E+011 2 4 0.30620E+01 0.30960E+01 0.30950E+011 2 5 0.32230E+01 0.31200E+01 0.31210E+011 2 6 0.30970E+01 0.31310E+01 0.31300E+011 2 7 0.32220E+01 0.31200E+01 0.31210E+011 2 8 0.30630E+01 0.30960E+01 0.30950E+011 2 9 0.31100E+01 0.30170E+01 0.30240E+011 2 10 0.29800E+01 0.30030E+01 0.29950E+011 2 11 0.25790E+01 0.25030E+01 0.25210E+011 3 1 0.31630E+01 0.30700E+01 0.30910E+011 32 0.36550E+01 0.36840E+01 0.36740E+011 3 3 0.38150E+01 0.37000E+01 0.37100E+011 3 4 0.37570E+01 0.37980E+01 0.37970E+011 3 5 0.39530E+01 0.38270E+01 0.38280E+01 1 3 6 0.37990E+01 0.38410E+01 0.38390E+01 1 3 7 0.39530E+01 0.38270E+01 0.38280E+01 1 3 8 0.37570E+01 0.37980E+01 0.37970E+01 1 3 9 0.38150E+01 0.37010E+01 0.37100E+01 1 3 10 0.36550E+01 0.36840E+01 0.36740E+01 1 3 11 0.31630E+01 0.30710E+01 0.30920E+01 1 4 1 0.33650E+01 0.32670E+01 0.32900E+01 1 42 0.38900E+01 0.39200E+01 0.39100E+01 1 43 0.40600E+01 0.39380E+01 0.39480E+01 14 4 0.39980E+01 0.40420E+01 0.40410E+01 1 45 0.42060E+01 0.40720E+01 0.40740E+01 1 46 0.40430E+01 0.40870E+01 0.40860E+01 1 47 0.42060E+01 0.40720E+01 0.40740E+01 1 48 0.39980E+01 0.40420E+01 0.40410E+01 1 49 0.40600E+01 0.39380E+01 0.39480E+01 1 4 10 0.38900E+01 0.39200E+01 0.39100E+01 1 4 11 0.33660E+01 0.32680E+01 0.32910E+01 1 5 1 0.31620E+01 0.30700E+01 0.30910E+01 1 5 2 0.36550E+01 0.36840E+01 0.36740E+01 1 5 3 0.38150E+01 0.37000E+01 0.37100E+01 1 5 4 0.37570E+01 0.37980E+01 0.37970E+01 1 5 5 0.39520E+01 0.38270E+01 0.38280E+01 1 5 6 0.38000E+01 0.38410E+01 0.38390E+01 1 5 7 0.39520E+01 0.38270E+01 0.38280E+01 1 5 8 0.37570E+01 0.37980E+01 0.37970E+01 1 5 9 0.38150E+01 0.37010E+01 0.37100E+01 1 5 10 0.36550E+01 0.36840E+01 0.36740E+01 1 5 11 0.31630E+01 0.30710E+01 0.30920E+01 1 6 1 0.25770E+01 0.25020E+01 0.25200E+01 1 6 2 0.29800E+01 0.30030E+01 0.29950E+01 1 6 3 0.31100E+01 0.30170E+01 0.30240E+01 1 6 4 0.30630E+01 0.30960E+01 0.30950E+01 1 6 5 0.32220E+01 0.31200E+01 0.31210E+01 1 6 6 0.30980E+01 0.31310E+01 0.31300E+01 1 6 7 0.32220E+01 0.31200E+01 0.31210E+01 1 6 8 0.30620E+01 0.30960E+01 0.30950E+01 1 6 9 0.31100E+01 0.30170E+01 0.30240E+01 1 6 10 0.29800E+01 0.30030E+01 0.29950E+01 1 6 11 0.25790E+01 0.25030E+01 0.25210E+01 1 7 1 0.16810E+01 0.16330E+01 0.16450E+01 1 7 2 0.19450E+01 0.19600E+01 0.19550E+01 1 7 3 0.20290E+01 0.19690E+01 0.19740E+01 1 7 4 0.19990E+01 0.20210E+01 0.20200E+01 1 7 5 0.21020E+01 0.20360E+01 0.20370E+011 7 6 0.20220E+01 0.20430E+01 0.20430E+011 7 7 0.21030E+01 0.20360E+01 0.20370E+011 7 8 0.19990E+01 0.20210E+01 0.20200E+011 7 9 0.20300E+01 0.19690E+01 0.19740E+011 7 10 0.19450E+01 0.19600E+01 0.19550E+011 7 11 0.16830E+01 0.16340E+01 0.16450E+01DATA FOR CRACK FRONT 1***********************ENERGY RELEASE RATE SET 1 ALONG CRACK FRONT 1 - from contour 1 N G-LOCAL NX NY NZ1 0.16840E+01 -0.50000E+00 0.00000E+00 0.86603E+002 0.20310E+01 -0.50000E+00 0.00000E+00 0.86603E+003 0.21040E+01 -0.50000E+00 0.00000E+00 0.86603E+004 0.21030E+01 -0.50000E+00 0.95367E-06 0.86603E+005 0.20300E+01 -0.50000E+00 -0.95370E-06 0.86603E+006 0.16830E+01 -0.50000E+00 0.00000E+00 0.86603E+00ENERGY RELEASE RATE SET 2 ALONG CRACK FRONT 1 - from contour 1 N G-LOCAL NX NY NZ1 0.25790E+01 -0.76604E+00 0.00000E+00 0.64279E+002 0.31110E+01 -0.76604E+00 0.00000E+00 0.64279E+003 0.32230E+01 -0.76604E+00 0.00000E+00 0.64279E+004 0.32220E+01 -0.76604E+00 0.14611E-05 0.64279E+005 0.31100E+01 -0.76604E+00 -0.14612E-05 0.64279E+006 0.25790E+01 -0.76604E+00 0.00000E+00 0.64279E+00ENERGY RELEASE RATE SET 3 ALONG CRACK FRONT 1 - from contour 1 N G-LOCAL NX NY NZ1 0.31630E+01 -0.93969E+00 0.00000E+00 0.34202E+002 0.38150E+01 -0.93969E+00 0.00000E+00 0.34202E+003 0.39530E+01 -0.93969E+00 0.00000E+00 0.34202E+004 0.39530E+01 -0.93969E+00 0.17923E-05 0.34202E+005 0.38150E+01 -0.93969E+00 -0.17924E-05 0.34202E+006 0.31630E+01 -0.93969E+00 0.00000E+00 0.34202E+00ENERGY RELEASE RATE SET 4 ALONG CRACK FRONT 1 - from contour 1 N G-LOCAL NX NY NZ1 0.33650E+01 -0.10000E+01 0.00000E+00 0.38147E-062 0.40600E+01 -0.10000E+01 0.00000E+00 0.00000E+003 0.42060E+01 -0.10000E+01 0.00000E+00 0.00000E+004 0.42060E+01 -0.10000E+01 0.19073E-05 0.38147E-065 0.40600E+01 -0.10000E+01 -0.19074E-05 0.00000E+006 0.33660E+01 -0.10000E+01 0.00000E+00 0.00000E+00ENERGY RELEASE RATE SET 5 ALONG CRACK FRONT 1 - from contour 1 N G-LOCAL NX NY NZ1 0.31620E+01 -0.93969E+00 0.00000E+00 -0.34202E+002 0.38150E+01 -0.93969E+00 0.00000E+00 -0.34202E+003 0.39520E+01 -0.93969E+00 0.00000E+00 -0.34202E+004 0.39520E+01 -0.93969E+00 0.17923E-05 -0.34202E+005 0.38150E+01 -0.93969E+00 -0.17924E-05 -0.34202E+006 0.31630E+01 -0.93969E+00 0.00000E+00 -0.34202E+00ENERGY RELEASE RATE SET 6 ALONG CRACK FRONT 1 - from contour 1 N G-LOCAL NX NY NZ1 0.25770E+01 -0.76604E+00 0.00000E+00 -0.64279E+002 0.31100E+01 -0.76604E+00 0.00000E+00 -0.64279E+003 0.32220E+01 -0.76604E+00 0.00000E+00 -0.64279E+004 0.32220E+01 -0.76604E+00 0.14611E-05 -0.64279E+005 0.31100E+01 -0.76604E+00 -0.14612E-05 -0.64279E+006 0.25790E+01 -0.76604E+00 0.00000E+00 -0.64279E+00ENERGY RELEASE RATE SET 7 ALONG CRACK FRONT 1 - from contour 1 N G-LOCAL NX NY NZ1 0.16810E+01 -0.50000E+00 0.00000E+00 -0.86603E+002 0.20290E+01 -0.50000E+00 0.00000E+00 -0.86603E+003 0.21020E+01 -0.50000E+00 0.00000E+00 -0.86603E+004 0.21030E+01 -0.50000E+00 0.95367E-06 -0.86603E+005 0.20300E+01 -0.50000E+00 -0.95370E-06 -0.86603E+006 0.16830E+01 -0.50000E+00 0.00000E+00 -0.86603E+00GMAX DISTRIBUTION FOR CRACK FRONT 1N GMAX VALUES DIRECTION (X,Y,Z)1 0.33653E+01 -0.10000E+01 0.00000E+00 0.49112E-032 0.40600E+01 -0.10000E+01 0.00000E+00 0.22922E-033 0.42062E+01 -0.10000E+01 0.00000E+00 0.25253E-034 0.42061E+01 -0.10000E+01 0.19073E-05 0.31615E-045 0.40599E+01 -0.10000E+01 -0.19074E-05 -0.13273E-076 0.33662E+01 -0.10000E+01 0.00000E+00 -0.20204E-07G SUMMARY : MODE I CONVERSION FROM G TO K USING: K =SQRT(ExG/(1-NuxNu))THESE ARE PLANE STRESS AND PLANE STRAIN K VALUESN GMAX VALUES K, Nu=0.000 K, Nu1 0.33653E+01 0.82040E+03 0.86002E+032 0.40600E+01 0.90111E+03 0.94463E+033 0.42062E+01 0.91720E+03 0.96148E+034 0.42061E+01 0.91718E+03 0.96146E+035 0.40599E+01 0.90110E+03 0.94461E+036 0.33662E+01 0.82051E+03 0.86013E+03------------ ------------ ------------MEAN: 0.38773E+01 0.87958E+03 0.92205E+03------------ ------------ ------------MIN.: 0.33653E+01 0.82040E+03 0.86002E+03 AT N = 1MAX.: 0.42062E+01 0.91720E+03 0.96148E+03 AT N = 3N NODE NO. CURRENT COORDS (AT WHICH Gmax IS CALCULATED)1 472 0.62500E+01 0.12482E-04 0.50000E+012 588 0.62500E+01 0.50000E+00 0.50000E+013 673 0.62500E+01 0.10000E+01 0.50000E+014 777 0.62500E+01 0.15000E+01 0.50000E+015 861 0.62500E+01 0.20000E+01 0.50000E+016 973 0.62500E+01 0.25000E+01 0.50000E+01SUMMARY OF LOADING******************REPORTS OF MAXIMUM AND MINIMUM ENERGY RELEASE RATE VALUES,DELTA sqrt(G), DELTA K AND CRACK GROWTH DIRECTION VECTORS(DELTA K0 = PLANE STRESS VALUE, DELTA Knu = PLANE STRA IN VALUE)CRACK FRONT 1MAXIMUM Gmax MINIMUM Gmax DELTA sqrt(G) DELTA K0 DELTA Knu1 0.33653E+01 0.00000E+00 0.18345E+01 0.82040E+03 0.86002E+032 0.40600E+01 0.00000E+00 0.20150E+01 0.90111E+03 0.94463E+033 0.42062E+01 0.00000E+00 0.20509E+01 0.91720E+03 0.96148E+034 0.42061E+01 0.00000E+00 0.20509E+01 0.91718E+03 0.96146E+035 0.40599E+01 0.00000E+00 0.20149E+01 0.90110E+03 0.94461E+036 0.33662E+01 0.00000E+00 0.18347E+01 0.82051E+03 0.86013E+03DIRECTION (X,Y,Z)1 -0.10000E+01 0.00000E+00 0.49112E-032 -0.10000E+01 0.00000E+00 0.22922E-033 -0.10000E+01 0.00000E+00 0.25253E-034 -0.10000E+01 0.19073E-05 0.31615E-045 -0.10000E+01 -0.19074E-05 -0.13273E-076 -0.10000E+01 0.00000E+00 -0.20204E-07***ANALYSIS COMPLETE***NO CRACK GROWTH REQUESTED。
abaqus论文
《从ABAQUS在土木工程中的一个应用实例简要说下我对ABAQUS软件的认识》专业：班级：学好：姓名：一、ABAQUS的发展作为一种功能强大的有限元分析软件，ABAQUS在商业有限元软件中占有极其重要的位置。
从简单的线弹性问题到复杂的几何非线性和材料非线性问题均获得了广泛应用。
其有效性不论在工程应用还是在科学研究方面均得到了验证。
ABAQUS包含了丰富的单元库和材料库，能够模拟各种材料受力和变形行为，特别是其提供了UMA T和UEL，为研究者提供了开发平台。
二、ABAQUS功能简介以ABAUQS6. 5为例，ABAQUS安装完成后，在程序菜单出现的AEAQUS 工具条下有以下几个选项: ABAQUSCAE ，ABAQUS Command，ABAQUS Documentation，ABAQUS Licensing，ABAQUS V erification，ABAQUS Viewer，My ABAQUS，Uninstall ABAQUS。
经常用到的主要是前三项，为了使ABAQUS 计算所涉及的文件均存储到同一目录下，运行ABAQUS Command比较方便，在ABAQUS Command的Dos窗口下，运行批处理文件abq651.bat，该文件存放在C：\ABAQUS\Command目录下(假定ABAQUS安装在C盘根目录下)，为了方便可将该文件另外保存为aba.bat。
该文件所包含的内容为@echo off C：\ABAQUS\6.5.1\exec\abq651. exe%，接下来就可以运行ABAQUS不同模块了。
如Aba cae进人CAE界面；ABAQUS在土木工程中的应用Aba viewer进人后处理；Aba job= wang interactive交互式运行wang. inp文件；Aba job= wang datacheck interactive交互运行wang. inp文件，且仅对其进行数据检查；Aba fetch job= terzaphi_cpe8p. inp将terzaphi_cpe8p. inp解压释放到当前目录下，因为文件名字比较长，可以仿前面定义批处理文件的方式进行类似处理，如定义文件名为 f. bat，其内容为Aba fetch job= terzaphi_cpe8p. inp，在当前目录下键人f并回车，即可完成文件的释放存储。
用ABAQUS分析工程中的多场耦合问题
[K ]{Cδ } − [L]{Cu } = {P} −{I} [B]T {v} + [Hˆ ]{u} = {Q}
(14)
-3-
Back
ABAQUS 软件 2003 年度用户论文集
[B]T 与[L]为耦合矩阵。 引入差分算子
{δ t+∆t } -{δ t } + ∆t[(1-ξ ){vt } + ξ{vt+∆t }]
防渗墙底部渗流速度矢量图 碎石土心墙上游坡脚速度矢量图
-7-
Back
ABAQUS 软件 2003 年度用户论文集
首先进行方案比较分析，分析如下三种方案： ①碎石土心墙＋防渗墙：在碎石土心墙上游面至碎石土心墙底部与刺墙交界处铺设土工 膜，防渗墙位于碎石土心墙的底部，采用倒悬挂式，封闭至覆盖层 4 的底部，心墙下游的坝底 设一层反滤层，下游坝坡脚处设排水沟；②混凝土面板＋防渗墙：堆石坝上有铺设混凝土面板， 防渗墙位于上游坝坡脚处于趾板相联接，采用倒悬挂式，封闭至覆盖层 4 的底部，坝底设一层 反滤层，下游坝坡脚处设排水沟；③坝面复合土工膜＋防渗墙：在坝上游面至趾板铺设土工膜， 防渗墙位于上游坝坡脚处于趾板相联接，采用倒悬挂式，封闭至覆盖层 4 的底部，坝底设一层 反滤层，下游坝坡脚处设排水沟。这三种防渗方案的防渗墙深度约为 37.1m。 经分析比较后，最终采取坝面复合土工膜＋碎石土心墙＋防渗墙防渗方案：在上游坝面至 碎石土心墙底部与刺墙交界处铺设土工膜，将碎石土心墙前移至上游坝坡脚，防渗墙位于心墙 下方，采用倒悬挂式、封闭至覆盖层 2 的底部，如图 1 和图 2 所示，该方案防渗墙的最大深度 约为：82.33m。 分析结果如下： 1．对于深厚覆盖层的基础，在现有施工技术条件下，防渗墙不能完全封闭覆盖层时，仍 具有筑坝的条件。 2．在深厚覆盖层的基础上修筑堆石坝，若覆盖层没有完全被防渗墙封闭时，以下几个问 题值得注意：①未封闭的覆盖层成为渗流的主要路径，工程的渗流量是否满足要求，在于未被 封闭的覆盖层上的渗流量；②防渗体下游的水位较高，使得防渗体下游的坝底存在一定的水头。 因此，在这种情况下，下游坝底的反滤层的铺设和排水沟的设置尤为重要；③在防渗墙的底部 存在着一个很小的区域，这个区域的水力坡降较大，可能造成该处的覆盖层不能满足渗透稳定 的要求。减小这个区域的水力坡降有效办法是增大渗径，加深防渗墙的深度是减小该处的水力 坡降的有效措施之一。 3．土工膜对坝体的防渗效果非常明显。土工膜的失效使得碎石土心墙的局部水力坡降超 出允许值。 4．防渗墙的局部开裂开叉对工程的防渗效果不会造成重大影响。 5．对于深厚覆盖层的坝基，在覆盖层没有被防渗墙完全封闭的情况下，由于防渗体下游 的水位较高，在作渗流场分析时，应将防渗体下游的堆石区考虑在分析的范围内，以确定反滤 层及排水沟的效果。计算也表明，位于防渗体下游的堆石区上的渗流满足 Darcy 定律适用范 围，可以参与渗流场计算。 6．由二维与三维模型比较可见，在最高坝段，二维位移－渗流场计算的水力坡降较三维 的大。三维分析表明，渗流场有明显的绕渗现象，二维分析不能表现这一现象，进行三维渗流 场分析是必要的。 （二）掺 MgO 混凝土拱坝的施工／运行仿真分析（TCM 分析）
清华牛人的ABAQUS总结[5篇模版]
![清华牛人的ABAQUS总结[5篇模版]](https://img.taocdn.com/s3/m/cb530f2c974bcf84b9d528ea81c758f5f61f2902.png)
清华牛人的ABAQUS总结[5篇模版]第一篇：清华牛人的ABAQUS总结准确的说，应该是谈谈我对发在这里的帖子质量的看法。
因为专业的原因，我对有限元及其软件是有很深感情的。
又因为一直对清华深有好感，可惜因为一些原因最终没有来清华深造，所以比较关注清华的BBS。
写这篇文章要耗费我差不多一个完整的下午，但是我愿意。
我知道学有限元其实不是一件容易的事情，我把我的想法说出来，希望对初学者有所裨益。
坦率的说，我认为这里有限元板块的质量是不高的。
之所以如此，是因为在这里很多人问的问题是太简单而且对自己不负责任的。
这不是版主的错，是因为我们许多人还没有养成良好的专业素养和严谨的精神。
请不要轻易的否认我的这个评价，好吗？至少，现在请不要。
就是对我的话不屑一顾，也应该是在看完我的话之后吧。
我用我的思维方式来说话，并不是每一个人都会习惯，请见谅！我还要声明的是，我本人的水平一般般，自己也对自己有很多的不满，所以在这里说的很可能很幼稚或者有错误。
请大家指教！我们应该有一个良好的讨论气氛。
有限元对许多工科的人而言，其必要性和重要性不言而喻。
问题在于，应该怎样的学习它呢？学习它，至少不用它到处害人也害己的话，我觉得至少要在下面四个方面有些基本知识：1、有限元基本理论及其求解基本步骤（数学基础）；2、有限元专业英语（英语基础）；3、你自己所属专业的东东（专业基础）；4、几何造型及拓扑学知识（建模基础）。
这个排序是由重到轻的。
接下来，我首先说一说上面四个方面的意义和作用；之后谈一下为什么我认为在这里问的相当一部分问题是太简单而且对自己不负责任的。
1、做专业就要有做专业的样子。
咱们理工科的学生，没有辛苦的付出是不可能有真正收获的。
收获和付出在这里成正比。
常常有人觉得有限元的软件很难，不好学，不好用，很多东西搞不懂，一提就头痛。
其实这里面相当的一部分是有限元基本理论可以解决的问题，而不是软件的设计思想不好。
现在的商用有限元软件，比如我用过的abaqus，ansys，adina以及algor，应该说它们的界面已经很友好了，包括帮助文档等等都不错。
ABAQUS静力结构分析实例
非线性静力分析总论simwe会员：LXCAD一，原创：大开大变形的悬臂压杆大变形例子*STEP, NLGEOM杆长2540mm，正方形截面h=12.7mm，弹性模量E=207.E3，临界FCR=171.5N 端点给-323.098N压力，求UX、UY。
*Heading** Job name: ansys3_3 Model name: Model-1*Preprint, echo=NO, model=NO, history=NO, contact=NO**** PAR TS***Part, name=Part-1*Node1, 0., 2540., 0.2, 0., 2286., 0.3, 0., 2032., 0.4, 0., 1778., 0.5, 0., 1524., 0.6, 0., 1270., 0.7, 0., 1016., 0.8, 0., 762., 0.9, 0., 508., 0.10, 0., 254., 0.11, 0., 0., 0.12, 0., 2413., 0.13, 0., 2159., 0.14, 0., 1905., 0.15, 0., 1651., 0.16, 0., 1397., 0.17, 0., 1143., 0.18, 0., 889., 0.19, 0., 635., 0.20, 0., 381., 0.21, 0., 127., 0.*Element, type=B321, 2, 12, 12, 3, 13, 23, 4, 14, 34, 5, 15, 45, 6, 16, 56, 7, 17, 67, 8, 18, 79, 10, 20, 910, 11, 21, 10*Nset, nset=_PickedSet2, internal, generate1, 21, 1*Elset, elset=_PickedSet2, inter nal, generate1, 10, 1*Nset, nset=_PickedSet3, internal, generate1, 21, 1*Elset, elset=_PickedSet3, inter nal, generate1, 10, 1*Nset, nset=_PickedSet4, internal, generate1, 21, 1*Elset, elset=_PickedSet4, inter nal, generate1, 10, 1** Region: (Section-1:Picked), (Beam Orientation:Picked)*Elset, elset=_PickedSet2, inter nal, generate1, 10, 1** Section: Section-1 Profile: Profile-1*Beam Section, elset=_PickedSet2, mater ial=steel, temperature=GRADIEN TS, section=REC T 12.7, 12.70.,0.,-1.*End Part****** ASSEMBL Y***Assembly, name=Assembly***Instance, name=Part-1-1, part=Part-1*End Instance***Nset, nset=_PickedSet4, internal, instance=Part-1-111,*Nset, nset=_PickedSet5, internal, instance=Part-1-11,*Nset, nset=_PickedSet8, internal, instance=Part-1-11,*End Assembly**** MA TERIALS***Material, name=steel*Elastic** ----------------------------------------------------------------**** STEP: Step-1***Step, name=Step-1, nlgeom=YES*Static1., 1., 1e-05, 1.**** BOUNDARY CONDITION S**** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre *Boundary_PickedSet4, ENCASTRE**** LOAD S**** Name: Load-1 Type: Concentrated force*Cload_PickedSet5, 1, 10._PickedSet5, 2, 0.**** OUTPUT REQUESTS***Restart, write, frequency=0**** FIELD OUTPUT: F-Output-1***Output, field, variable=PRESELECT**** HISTOR Y OU TPU T: H-Output-1***Output, history, variable=PRESELEC T*End Step** ----------------------------------------------------------------**** STEP: Step-2***Step, name=Step-2, nlgeom=YES*Static1., 1., 1e-05, 1.**** LOAD S**** Name: Load-1 Type: Concentrated force*Cload, op=N EW** Name: Load-2 Type: Concentrated force*Cload, op=N EW_PickedSet8, 2, -323.098**** OUTPUT REQUESTS***Restart, write, frequency=0**** FIELD OUTPUT: F-Output-1***Output, field, variable=PRESELECT**** HISTOR Y OU TPU T: H-Output-1***Output, history, variable=PRESELEC T*End Step要点：打开大变形，加一水平干扰力。
ABAQUS在岩土工程中的应用
ABAQUS在岩土工程中的应用ABAQUS是一款功能强大的工程仿真软件，广泛应用于各种领域，包括岩土工程。
在岩土工程中，ABAQUS可以帮助工程师模拟复杂的物理现象，如沉降、位移、应力等，为工程设计和优化提供强有力的支持。
岩土工程涉及到许多关键问题，如沉降、位移和应力等。
这些问题是岩土工程设计的核心，也是工程师们的重点。
借助ABAQUS软件，我们可以对这些复杂的问题进行有效的模拟和分析。
在ABAQUS中，对岩土工程问题进行模拟分析的流程大致可以分为以下步骤：模型建立：首先需要建立反映实际工程问题的三维模型。
ABAQUS提供了强大的建模功能，允许用户根据需要创建复杂的几何形状。
材料设置：在模型建立完成后，需要定义材料的属性，如弹性模量、泊松比、密度等。
ABAQUS提供了多种材料模型，可以根据实际情况选择合适的模型。
边界条件：定义模型的边界条件，如固定边界、自由边界等。
这些边界条件对于模拟结果的准确性至关重要。
网格划分：对模型进行离散化，将其划分为一系列小的单元，以便进行数值计算。
ABAQUS提供了多种网格划分方式，可以根据实际问题进行选择。
求解设置：设置求解器和求解参数，如迭代次数、收敛准则等。
求解过程：运行求解器，得出模型的位移、应力等结果。
结果后处理：对求解结果进行可视化处理，如生成云图、动画等，以便进行深入分析和优化。
以一个实际的应用实例来说明ABAQUS在岩土工程中的应用。
某桥梁基础设计面临着复杂的土壤环境，为了确保桥梁的安全性和稳定性，我们需要对土壤进行详细的模拟和分析。
我们使用ABAQUS建立桥梁和土壤的三维模型，并定义材料的属性。
然后，我们根据实际工程情况设置边界条件，如固定桥梁的底部和侧向边界，施加桥梁的重量等。
接着，我们对模型进行网格划分，并设置求解器和求解参数。
在求解过程中，ABAQUS通过迭代计算得出位移和应力的分布情况。
我们进行结果后处理，将位移和应力云图进行可视化，以便于工程师们进行深入的分析和优化设计。
abaqus案例
abaqus案例
Abaqus是一种广泛使用的有限元分析软件，用于模拟各种工程问题。
以下是一个关于抗风塔设计的Abaqus案例。
针对抗风塔的设计，需要进行结构分析来确定其承受最大风荷载
的能力。
这个案例使用了Abaqus这一强大的有限元分析工具，其中包
括了以下步骤：
首先，建立三维模型。
这个模型代表了一个矩形筒体形的抗风塔，其尺寸和几何形状与实际设计一致。
其次，定义材料性能。
结构材料
是铝合金，使用Abaqus的材料数据库来定义材料的应力-应变曲线。
接下来，设置边界条件。
在这个案例中，地面为固定边界，风荷
载通过侧面分布在塔身上。
对于风荷载，使用预定义的速度和压力分布，这相当于对风场进行了划分并进行平均计算。
完成了这些步骤后，可以进行结构稳态解的计算。
在这个阶段，Abaqus将使用所定义的材料和边界条件，对模型进行处理，来计算塔
结构的稳态响应和位移。
这个过程可以帮助我们确定抗风塔当前的承载能力及其弹性变形情况。
最后，我们可以将模型加入最大风荷载，以模拟超载情况下的反应。
这可以帮助确定抗风塔的耐久性，以及可能的破坏形式。
整个过程都可以通过Abaqus快速灵活地完成，它可以模拟出各种不同的条件下的结构响应，并且提供了极为详尽的计算结果。
这可以帮助设计人员更好地优化结构设计，并确保其能够承受最大的负载。
Abaqus课程设计报告论文
学校代码：10128学号：201220907039课程论文题目：球体接触问题的有限元分析学生：劲波学院：理学院班级：工程力学12-1班指导教师：韦广梅、周承恩2015年9月26日工业大学课程设计（论文）任务书课程名称：有限单元法课程设计学院：理学院班级：力学12-1班_ 学生：劲波学号：9 指导教师：韦广梅、周承恩摘要许多工程问题都涉及两个或多个部件之间的接触，在这些问题中，当两个物体彼此接触时，垂直于接触面上的力作用在两个物体上。
如果在接触面之间存在摩擦，可能产生剪力以阻止物体的切向运动。
本文主要是运用ABAQUS软件对椭球体与刚性半无限大体接触问题进行探讨，因为椭球体是前后左右都对称的三维实体，所以对它取四分之一分析，这样不仅计算方便而且结果直观。
先将材料定义为弹性材料通过前处理建模、计算、后处理来分析椭球体与大体之间的接触，运用多种方法讨论接触问题。
结果得到的数据远超弹性屈服极限，后将材料重新赋值为弹塑性，仿照弹性的分析方法对弹塑性进行分析，得到的结果完全正确。
关键词：接触；ABAQUS；弹塑性；屈服极限AbstractMany engineering problems involve contact between two or more components, in these problems, when two objects come into contact with each other, on the surface of the perpendicular to the contact force between two objects. If there is friction between contact surface, shear may be produced to prevent the tangential motion of the object. This article mainly using ABAQUS software to ellipsoid and rigid semi-infinite discusses the general contact problem, because the ellipsoid is before and after the left and right sides is symmetrical three-dimensional entity, so take a quarter analysis to it, that not only the calculation results of convenient and intuitive. Before the adoption of the first material is defined as elastic material modeling, calculation, after processing to analyze between ellipsoid and general contact, using a variety of methods to discuss contact problem. Than elastic yield limit, the resulting data will reopen assignment for elastic-plastic materials, imitates the elastic analysis method of elastic-plastic analysis, get the right result.Key words: Contact; ABAQUS; Elastic-plastic; Yield limit目录第一章绪论1.1 有限单元法课程设计目的与任务 (1)1.2 接触问题概述 (1)1.3 有限元法概述 (1)第二章有限元工具ABAQUS介绍及应用方法 (3)2.1 ABAQUS简介 (3)2.2 ABAQUS功能简介 (3)第三章椭球—平面接触问题有限元分析 (7)3.1 平面－球接触问题描述 (7)3.2 用ABAQUS建立有限元模型 (7)3.2.1 建立实体 (7)3.2.2 赋值材料与装配部件 (8)3.2.3 创建分析步 (8)3.2.4 定义接触和边界条件 (9)3.2.5 划分网格 (9)3.2.6 计算及后处理 (10)3.2.7 理想弹塑性分析 (11)第四章结论与收获 (14)参考文献 (15)第一章绪论1.1有限单元法课程设计目的与任务目的：《有限单元法课程设计》是“有限单元法”课程的拓展部分。
ABAQUS热分析毕业设计论文
摘要复合材料因具有高的比强度、比刚度、耐腐蚀性以及材料的可设计性等优良性能，已成为新一代飞机中必不可少的材料。
目前，金属工装因其热膨胀系数较大，成型的复合材料制件以不能满足要求，采用复合材料作为复材制件成型工装的材料是今后发展的趋势。
由于传统工装设计以经验为主，造成工装结构不合理，降低了制件的质量，增加了生产成本。
因此，优化工装结构，稳定产品质量、降低成本已成为复材制件生产中继续解决的问题。
本文提出利用Python语言在ABAQUS平台上开发复材工装分析系统。
在此基础上，对一套实际复材工装进行热应力分析，得到工装的整体应变云图结果。
这些结果对工装的进一步结构优化提供了定量参考。
关键词：复材工装、ABAQUS、Python、热应力分析Thermally analyzed of composited material processequipment of aircraft based on ABAQUS simulationAbstractWith high specific strength，high specific stiffness ，erode-resistant and material designable, polymer composite material become essential in the new generation of aircraft. nowadays ,the composite parts formed using metal frock with a large coefficient of thermal expansion can not meet the requirements. therefore, the uses of composite material frock forming the composite part are the development trend in the future. the traditional design of frock depends on experiences mainly. resulting in the structure of frock is unreasonable. Affecting the quality of the products and increasing the production costs. Consequently it is problem need to be solved urgently in forming composite parts such as optimization of frock structure .in this paper, composite material frock design system is programmed in ABAQUS using python language. Base on these, an actual composite material frock is thermally analyzed. These results provide quantificational reference for more optimization of frock. Keywords: composite material process equipment, ABAQUS, python, thermal analysis毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
目录第一章Abaqus简介 (1)一、Abaqus总体介绍 (1)二、Abaqus基本使用方法 (2)1.2.1 Abaqus分析步骤 (2)1.2.2 Abaqus/CAE界面 (3)1.2.3 Abaqus/CAE的功能模块 (3)第二章基于Abaqus的通孔端盖分析实例 (4)一、工作任务的明确 (5)二、具体步骤 (5)2.2.1 启动Abaqus/CAE (4)2.2.2 导入零件 (5)2.2.3 创建材料和截面属性 (6)2.2.4 定义装配件 (7)2.2.5 定义接触和绑定约束（tie） (10)2.2.6 定义分析步 (14)2.2.7 划分网格 (15)2.2.8 施加载荷 (19)2.2.9 定义边界条件 (20)2.2.10 提交分析作业 (21)2.2.11 后处理 (22)第三章课程学习心得与作业体会 (22)第一章： Abaqus简介一、Abaqus总体介绍Abaqus是功能强大的有限元分析软件，可以分析复杂的固体力学和结构力学系统，模拟非常庞大的模型，处理高度非线性问题。
Abaqus不但可以做单一零件的力学和多物理场的分析，同时还可以完成系统级的分析和研究。
Abaqus使用起来十分简便，可以很容易的为复杂问题建立模型。
Abaqus具备十分丰富的单元库，可以模拟任意几何形状，其丰富的材料模型库可以模拟大多数典型工程材料的性能，包括金属、橡胶、聚合物、复合材料、钢筋混泥土、可压缩的弹性泡沫以及地质材料（例如土壤、岩石）等。
Abaqus主要具有以下分析功能：1.静态应力/位移分析2.动态分析3.非线性动态应力/位移分析4.粘弹性/粘塑性响应分析5.热传导分析6.退火成形过程分析7.质量扩散分析8.准静态分析9.耦合分析10.海洋工程结构分析11.瞬态温度/位移耦合分析12.疲劳分析13.水下冲击分析14.设计灵敏度分析二、Abaqus基本使用方法1.2.1 Abaqus分析步骤有限元分析包括以下三个步骤：1.前处理（Abaqus/CAE）:在前期处理阶段需要定义物理问题的模型，并生成一个Abaqus输入文件。
提交给Abaqus/Standard或Abaqus/Explicit。
2.分析计算（Abaqus/Standard或Abaqus/Explicit）:在分析计算阶段，使用Abaqus/Standard或Abaqus/Explicit求解输入文件中所定义的数值模型，通常以后台方式运行，分析结果保存在二进制文件中，以便于后处理。
3.后处理（Abaqus/CAE或Abaqus/Viewer）: Abaqus/CAE的后处理部分又称为Abaqus/Viewer，可用来读入分析结果数据，以多种方法显示分析结果，包括彩色云纹图、动画、变形图和XY曲线图等。
1.2.2 Abaqus/CAE界面图1-1 Abaqus/CAE的主窗口1.2.3 Abaqus/CAE的功能模块一般使用Abaqus/CAE进行应力分析有以下几个经典的步骤：1、问题的描述→2、启动Abaqus/CAE→3、创建部件→4、创建材料和截面属性→5、定义装配件→6、设置分析步→7、定义边界条件和载荷→8、划分网格→9、提交分析作业→10、后处理→11、退出Abaqus/CAEAbaqus/CAE包括一系列的功能模块（module）,每个模块均包含其特定的工具，在Module（模块）列表中可以选择各个功能模块（如图1-2），图1-2 选择功能模块这些模块的次序同时也是Abaqus/CAE所推荐的模型创建顺序（如图1-3），当然也可以首先划分网格（如图1-4），这样做的好处是，往往在划分网格的过程中，会发现部件的几何模型需要进一步修改，而经过这些修改后，已经定义好的边界条件、载荷和接触等可能变为无效的，需要再重新定义。
图1-3 Abaqus/CAE所推荐的模型创建顺序图1-4 根据需要来选择适当的建模顺序第二章基于Abaqus的通孔端盖分析实例一、工作任务的明确在安装文件夹中有一个simple文件夹，里面有三个.x_t文件，如图2-1所示。
它们是三个零件，该文章目的是演示将这三个文件通过Abaqus装配好并达到划分网格的目的，最终生成一个.ace文件。
图2-1二、具体步骤2.2.1 启动Abaqus/CAE在Windows操作系统中：[开始]→[程序]→[Abaqus 6.8-1]→[ABAQUS CAE]。
启动ABAQUS CAE后，在出现的Start Session(开始任务)对话框中选择Create Model Database(创建新模型数据库)。
如图2-2所示。
图2-22.2.2 导入零件如图2-3左图，依次单击主菜单中的File→Import→Part，弹出图2-3右图对话框，选择simple文件夹后可看到三个.x_t文件，它们是通孔端盖的三个组成部件，分别是螺钉、安装板和固定圈。
图2-3 导入零件导入后可看到零件如图2-4所示。
此时零件尚未定义材料，所以表面是白色，待定义了材料后零件表面显示蓝色。
零件zdc 零件LD 零件azb图2-42.2.3 创建材料和截面属性在窗口左上角的Module(模块)列表中选择Property(特性)功能模块，按照以下步骤来定义材料。
1)创建材料点击左侧工具区域中的（Create Material）,也可以直接双击左侧模型树中的Material来完成此项操作。
在Name(材料名称)后面输入Steel，点击此对话框中的Mechanical(力学特性)→Elasticity(弹性)→Elastic。
在数据表中设置Young’s Modulus(弹性模量)为210000，Poisson’s Ratio(泊松比)为0.3，其余参数不需改变（如图2-5所示），点击OK。
按上述同样的方法创建ZM-6合金材料。
图2-5 创建材料2)创建截面属性点击左侧工具区中的（Create Section），在弹出的对话框中的Name中输入Section-Steel，点击continue，在接着弹出的对话框中将material下拉菜单中的Steel选中，单击OK，此时已定义好Section-Steel。
按照同样的方法定义Section-ZM，如图2-6所示。
图2-6创建截面属性3)给部件赋予截面属性以安装板为例，点击左侧工具区的（Assign Section），点击视图区中的安装板模型，Abaqus/CAE以红色高亮度显示被选中的实体边界，如图2-7所示。
在视图区的空白点击鼠标中键，弹出Edit Section Assignment对话框，将Section设置为Section-Steel，点击OK。
另外两个零件的赋予截面属性方法与此类似。
赋予截面属性后，各个部件表面变成绿色。
图2-72.2.4 定义装配件1)添加实体进入Assembly功能模块，点击，选中全部部件，然后点击OK，如图2-8所示。
图2-8 添加装配实体2)使固定圈下表面与安装板上表面相接触长按按钮，弹出接触面约束类型选项，将鼠标移动至（Face to Face）然后放开，如图2-9所示。
图2-9首先点击固定圈下表面，再点击安装板的上表面，在所点击的位置会显示出面的方向。
在窗口底部的提示区中点击Flip来更正方向，然后点击OK。
如图2-10左图所示。
提示区中显示的两个面的默认距离为0.0，按回车键确认。
得到的模型如图2-10右图所示。
图2-103)使安装板和固定圈的螺钉孔对齐重复图2-9所示的步骤，将鼠标移动至（Coaxial）放开，先点击固定圈上的某个螺纹孔的内表面，再点击安装板的某个螺纹孔的内表面，点击OK。
图2-11如图2-11所示，此时只定义了一对螺纹孔，还需要定义一对螺纹孔才能把固定圈装配到正确的位置，再次点击，先后选择需对齐的两个螺纹孔的内表面，点击OK，此时安装板与固定圈就已经装配到正确的位置，如图2-12表示。
图2-124)安装螺钉安装螺钉分为三个部分，先要使螺钉螺纹柱与螺纹孔对齐，然后是安装板的底面与螺钉头部的底面相接触，最后对这个螺钉进行阵列。
首先单击，先后选择螺钉的螺纹柱和安装板的螺纹孔，点击Flip选好方向后点击OK确认，螺钉的位置如图2-13所示。
图2-13单击，分别选择螺钉头部的底面和安装板的底面，使螺钉装配到正确的位置，这里要注意箭头的方向，如果方向相反则装配会出现错误，点击Flip改正方向后点击OK确认。
正确的螺钉位置如图2-14所示。
图2-14点击工具区中的（环形阵列），在视图区点击螺钉，在提示区中点击Done,弹出阵列对话框（如图2-15左图），将阵列数量改成8，角度改为360度，点击Axis然后选择Y轴，待看到阵列出的螺钉处于正确位置后点击OK退出,阵列后如图2-15右图所示。
图2-15 阵列螺钉2.2.5 定义接触和绑定约束（tie）下面将在安装板和固定环之间定义接触，在螺钉和底座的螺纹处建立绑定约束。
1)定义接触和绑定约束所要用到的各个面在定义这些面之前，由于已经完成assembly过程；造成平面选择困难，此过程成可以使用view-assembly display options；如图2-16所示，分别勾选不同的零部件进行显示，从而正确快速地确定出对应的面。
图2-16进入Interaction功能模块，双击模型树中的Surface,新建一个接触面，先给这个面命名，确认后开始在模型上选择要定义的面，选好后点击提示区的Done。
注意：在选择阵列过的螺钉上的接触面时要先把安装板设为Suppress状态。
依次定义如图2-17所示的接触面（螺钉只选一个作解释，其余的都一样）：图2-17 定义接触面2)在螺纹处定义绑定约束在模型树中双击Constraint，出现如图2-18左图所示对话框，按图中设置好名称和类型后点击Continue，点击窗口底部提示区右侧的Surface按钮，在弹出的Region Selection对话框中，选中Surf-zdc-ld1来作为绑定约束的主面，点击Continue，再次点击提示区中的Surface按钮,选中Surf-ld1-zdc来作为绑定约束的从面。
在弹出的Edit Constraint对话框中（如图2-18中图），将Position Tolerance设为Specify distance:0.02，点击OK退出。
按同样的方法设置如图2-18右图中的各个绑定约束。
图2-183)定义带库伦摩擦的接触属性和不带摩擦的接触属性点击（Create Interaction Property），在Name后面输入IntProp-friction-015，点击Continue。
点击Mechanical→Tangential Behavior，把摩擦类型改为Friction formulation:Penalty，在Friction Coeff下面输入0.15，然后点击OK。