1、打开已有的铺层顺序设计
把相关区域定义成
set,然后在Region
中选择set。
右击,选择write to file
输入文件名
打开相应的文件看一下,可
以改角度等。Picked需要改
成相应的set名
新建一个复合材料铺层设计。
右击,导入文件,
注意必须选第一个
单元格,“对号”
选择要导入的文件
双击。
闲话不说,直接上链接 https://www.doczj.com/doc/1116158426.html,/s/122MWa abaqus 6.13对操作系统的新要求: 自abaqus 6.13版本开始,将不再支持windows 的32位操作平台; 同时,也不再支持windows xp和 windows vista操作系统; 安装: 1. Run "Install Abagus Product & Licensing" 2. In SIMULIA FLEXnet License Server window select "Just install the license utilities" NOTE: If you already have SIMULIA FLEXnet License Server for ABAQUS 6.12-3 installed and running you can use it for 6.13-1 too 3. After finishing License Utilities setup copy files "ABAQUS.lic" and "ABAQUS.log" to
四Abaqus在复合材料领域的优势 4.1 复合材料介绍 4.1.1 复合材料的应用 复合材料有许多特性: 1、制造工艺简单 2、比强度高,比刚度大 3、具有灵活的可设计性 4、耐腐蚀,对疲劳不敏感 5、热稳定性能、高温性能好 由于复合材料的上述优点,在航空航天、汽车、船舶等领域,都有广泛的应用。复合材料的大量应用对分析技术提出新的挑战。
4.1.2 复合材料的结构 复合材料是一种至少由两种材料混合而成的宏观材料,其中的一种材料被称作基体,其它的材料称作纤维。其中纤维可以包含很多不同的 形式:离散的宏观粒子,任意方向的短纤维,规则排列的纤维和织物。 4.1.3 典型的复合材料 1)单向纤维层合板----冲击分析
2)编织复合材料---- 挤压分析 3)蜂窝夹心复合材料----不可见冲击损伤分析
基体和纤维的存在形式以及材料属性对于复合材料的力学行为有 着很大的影响。改变纤维和基体的属性目的就是在于生成一种复合材料具有如下性质: 1)低成本:原型,大规模生产,零件合并,维修,技术成熟。 2)期望的重量:轻重量,比重分配合理。 3)改进的强度和刚度:高强度/高刚度比。 4)改进的表面属性:良好的耐腐蚀性,表面抛光性好。 5)期望的热属性:较低的热传导性,热膨胀系数较低。 6)独特的电属性:具有较高的绝缘强度,无磁性。 7)空间适应性:大部件,特殊的几何构型。 4.1.4 复合材料的有限元模拟 根据不同的分析目的,可以采用不同的复合材料模拟技术: 1)微观模拟:将纤维和基体都分别模拟为可变形连续体。 2)宏观模拟:将复合材料模拟为一个正交各向异性体或是完全各向
第二章 ABAQUS 基本使用方法 [2](pp15)快捷键:Ctrl+Alt+左键来缩放模型;Ctrl+Alt+中键来平移模型;Ctrl+Alt+右键来旋转模型。 ②(pp16)ABAQUS/CAE 不会自动保存模型数据,用户应当每隔一段时间自己保存模型以避免意外丢失。 [3](pp17)平面应力问题的截面属性类型是Solid(实心体)而不是Shell(壳)。 ABAQUS/CAE 推荐的建模方法是把整个数值模型(如材料、边界条件、载荷等)都直接定义在几何模型上。载荷类型Pressure 的含义是单位面积上的力,正值表示压力,负值表示拉力。 [4](pp22)对于应力集中问题,使用二次单元可以提高应力结果的精度。 [5](pp23)Dismiss 和Cancel 按钮的作用都是关闭当前对话框,其区别在于:前者出现在包含只读数 据的对话框中;后者出现在允许作出修改的对话框中,点击Cancel 按钮可关闭对话框,而不保存 所修改的内容。 [6](pp26)每个模型中只能有一个装配件,它是由一个或多个实体组成的,所谓的“实体”(instance) 是部件(part)在装配件中的一种映射,一个部件可以对应多个实体。材料和截面属性定义在部件上,相互作用(interaction)、边界条件、载荷等定义在实体上,网格可以定义在部件上或实体上,对求解过程和输出结果的控制参数定义在整个模型上。 [7](pp26) ABAQUS/CAE 中的部件有两种:几何部件(native part)和网格部件(orphan mesh part)。 创建几何部件有两种方法:(1)使用Part 功能模块中的拉伸、旋转、扫掠、倒角和放样等特征来直 接创建几何部件。(2)导入已有的CAD 模型文件,方法是:点击主菜单File→Import→Part。网格部件不包含特征,只包含节点、单元、面、集合的信息。创建网格部件有三种方法:(1)导入ODB 文件中的网格。(2)导入INP 文件中的网格。(3)把几何部件转化为网格部件,方法是:进入Mesh 功能模块,点击主菜单Mesh→Create Mesh Part。 [8](pp31)初始分析步只有一个,名称是initial,它不能被编辑、重命名、替换、复制或删除。在初始分析步之后,需要创建一个或多个后续分析步,主要有两大类:(1)通用分析步(general analysis step)可以用于线性或非线性分析。常用的通用分析步包含以下类型: —Static, General: ABAQUS/Standard 静力分析 —Dynamics, Implicit: ABAQUS/Standard 隐式动力分析 —Dynamics, Explicit: ABAQUS/ Explicit 显式动态分析 (2)线性摄动分析步(linear perturbation step)只能用来分析线性问题。在ABAQUS/Explicit 中 不能使用线性摄动分析步。在ABAQUS/Standard 中以下分析类型总是采用线性摄动分析步。 —Buckle: 线性特征值屈曲。 —Frequency: 频率提取分析。 —Modal dynamics: 瞬时模态动态分析。 —Random response: 随机响应分析。 —Response spectrum: 反应谱分析。 —Steady-state dynamics: 稳态动态分析。 [9](pp33)在静态分析中,如果模型中不含阻尼或与速率相关的材料性质,“时间”就没有实际的物 理意义。为方便起见,一般都把分析步时间设为默认的 1。每创建一个分析步,ABAQUS/CAE 就会自动生成一个该分析步的输出要求。 [10] (pp34)自适应网格主要用于ABAQUS/Explicit 以及ABAQUS/Standard 中的表面磨损过程 模拟。在一般的ABAQUS/Standard 分析中,尽管也可设定自适应网格,但不会起到明显的作用。 Step 功能模块中,主菜单Other→Adaptive Mesh Domain 和Other→Adaptive Mesh Controls 分别 设置划分区域和参数。 [11](pp37)使用主菜单Field 可以定义场变量(包括初始速度场和温度场变量)。有些场变量与分析步有关,也有些仅仅作用于分析的开始阶段。使用主菜单Load Case 可以定义载荷状况。载荷状况由一系列的载荷和边界条件组成,用于静力摄动分析和稳态动力分析。
超详细A b a q u s安装教 程 Pleasure Group Office【T985AB-B866SYT-B182C-BS682T-STT18】
Abaqus2017 安装教程 将安装镜像加载至虚拟光驱 以管理员身份运行J:\1下的,保持默认点击下一步。 首先安装的是拓展产品文档,点击“下一步”,选择安装目标,并点击“下一步” 选择文档需要包含的部件,并点击“下一步” 提示程序安装信息足够,点击“安装” 提示SIMULIA 2017文档安装成功,点击“关闭”。 接着会自动弹出Abaqus Simulation Services,修改安装目标地址,并点击“下一步” 选择您需要的部件,并点击“下一步” 检测安装信息足够,点击“安装”, 提示安装成功够点击“关闭” 接来下会自动安装Abaqus Simulation Services CAA API,点击“下一步”, 选择您需要的安装的部件,点击下一步 安装完成后“关闭” 接下来是自动安装 Abaqus CAE 找到安装包里的_SolidSQUAD_文件夹,将里面的License文件复制到Abaqus安装目录里的SIMULIA文件夹里。然后打开License文件夹,改.lic许可证文件的计算机名,同时再新建一个.log日志文件。 粘贴到D:\SIMULIA 下 在License下新建一个文本文档,重命名为 用记事本打开 使用计算机名替换this_host 保存后退出。 以管理员方式打开,点击Config Services,按如下配置 点击Save Service。切换到Start/Stop/Reread选项卡,点击Stop Server,再点击Start Server。下面提示Server Start Successful则配置成功。关闭窗口。 回到Abaqus CAE的安装界面。 在License Server 1 中输入27011@DESKTOP-Q8CNNLR 注:DESKTOP- Q8CNNLR 是计算机的用户名。点击下一步。 点击下一步 修改目录,点击下一步 设置工作空间路径,建议选择较大的硬盘分区。点击“下一步” 信息确认无误开始安装。 点击Continue 提示Abaqus CAE安装完成,点击关闭。 下面是Tosca 2017的安装。根据需求选择安装。 根据需求选择接口 若安装按ANSA可选择路径,没有则直接下一步。 若没有安装Fluent ,取消勾选FLUENT 剩下步骤类似上面。最后点击安装。 最后是 Isight 2017的安装。步骤后之前相同,一直点下一步记忆可以了
复合材料不只是几种材料的混合物。它具有普通材料所没有的一些特性。它在潮湿和高温环境,冲击,电化学腐蚀,雷电和电磁屏蔽环境中具有与普通材料不同的特性。 复合材料的结构形式包括层压板,三明治结构,微模型,编织预成型件等。 复合材料的结构和材料具有同一性,并且可以在结构形成时同时确定材料分布。它的性能与制造过程密切相关,但是制造过程很复杂。由于复合结构不同层的材料特性不同,复合结构在复杂载荷作用下的破坏模式和破坏准则是多种多样的。 在ABAQUS中,复合材料的分析方法如下 1,造型 它的结构形式决定了它的建模方法,并且可以使用基于连续体的壳单元和常规壳单元。复合材料被广泛使用,但是复合材料的建模是一个困难。铺设复杂的结构光需要一个月 2,材料
使用薄片类型(层材料)建立材料参数。材料参数可以工程参数的形式给出,或者材料强度数据可以通过子选项给出。这种材料仅使用平面应力问题。 ABAQUS可以通过两种方式定义层压板:复合截面定义和复合层压板定义 复合截面定义对每个区域使用相同的图层属性。这样,我们只需要建立壳体组合即可将截面属性分配给二维(在网格中定义的常规壳体元素)或三维(三维的大小应与壳体中给定的厚度一致)。基于网格中定义的连续体的壳单元) ABAQUS复合材料分析方法介绍 复合叠加定义是由复合布局管理器定义的,它主要用于在模型的不同区域中构造不同的层。因此,应在定义之前对区域进行划分,并且应将不同的层分配给不同的区域。可以根据常规外壳的元素和属性进行定义。 传统的壳单元定义了每个层的厚度,并将其分配给二维模型。应该给基于连续体的壳单元或实体单元提供3D模型(厚度是相对于单元长度的系数,因此厚度方向可以分为一层单元)。
复合材料模型建模与分析 1. Cohesive单元建模方法 1.1 几何模型 使用内聚力模型(cohesive zone)模拟裂纹的产生和扩展,需要在预计产生裂纹的区域加入cohesive层。建立cohesive层的方法主要有: 方法一、建立完整的结构(如图1(a)所示),然后在上面切割出一个薄层来模拟cohesive 单元,用这种方法建立的cohesive单元与其他单元公用节点,并以此传递力和位移。 方法二、分别建立cohesive层和其他结构部件的实体模型,通过“tie”绑定约束,使得cohesive单元两侧的单元位移和应力协调,如图1(b)所示。 (a)cohesive单元与其他单元公用节点(b)独立的网格通过“tie”绑定 图1.建模方法 上述两种方法都可以用来模拟复合材料的分层失效,第一种方法划分网格比较复杂;第二种方法赋材料属性简单,划分网格也方便,但是装配及“tie”很繁琐;因此在实际建模中我们应根据实际结构选取较简单的方法。 1.2 材料属性 应用cohesive单元模拟复合材料失效,包括两种模型:一种是基于traction-separation 描述;另一种是基于连续体描述。其中基于traction-separation描述的方法应用更加广泛。 而在基于traction-separation描述的方法中,最常用的本构模型为图2所示的双线性本构模型。它给出了材料达到强度极限前的线弹性段和材料达到强度极限后的刚度线性降低软化阶段。注意图中纵坐标为应力,而横坐标为位移,因此线弹性段的斜率代表的实际是cohesive单元的刚度。曲线下的面积即为材料断裂时的能量释放率。因此在定义cohesive的力学性能时,实际就是要确定上述本构模型的具体形状:包括刚度、极限强度、以及临界断裂能量释放率,或者最终失效时单元的位移。常用的定义方法是给定上述参数中的前三项,也就确定了cohesive的本构模型。Cohesive单元可理解为一种准二维单元,可以将它看作被一个厚度隔开的两个面,这两个面分别和其他实体单元连接。Cohesive单元只考虑面外的力,包括法向的正应力以及XZ,YZ两个方向的剪应力。 下文对cohesive单元的参数进行阐述,并介绍参数的选择方法。
1,建立模型Part Module :类型三维,solid,旋转;按尺寸绘图,done,设置旋转角此处为360度。 2,建立参考面,将圆筒分成两半 3,Assembly Module :类型Independent 分区partition截面 4,Mesh module : 点击remove空二,选择cells消隐分区 X Select entities to remove: Cells Undo 撒种子时,需要分几层就在边缘上撒多少个种子,在每条边上尽量都撒相同数量的种子, 生成结构网格,生成的网格才比较规整。 (注意,此处的mesh,对象为assembly,而不是part) 生成网格后,Mesh: Create Mesh Part Module I- Mesh * Model:j Model-1 abject: * Awembly Part「 4,Job Module : Create Job,例如job-007-01,运行生成job-007-01.inp 文件,保存成007-01.cae 文件。 5,File: New打开新窗口
6,File: Import : Model 选择job-007-01.inp 打开 7,Mesh Module: Tools: Surface manager: create: by angle 定义surface 集合 Tools: Set manager: create: Element: by angle 定义Element 集合 用以下三个命令操作,选择恰当的面。 丄i Select the Entity Closest to the Screen, ---- Select From Exterior En tities '包i 一 J Select From Interior Entities (左键点击第二个图标不放拖出即可) 注:定义Element集合时,可以从外到内,定以一层后,在display中--- -:把定义的那层remove掉再定义下面一层。 8,Mesh: Edit :Mesh : Mesh Offset (create solid layers): Surfaces (选择相应的面):Total thickness定义厚度,生成cohesive单元,把其之前定义的几层surface,都生成cohesive单丿元。 9,Mesh: Element type :对cohesive 单元,Family 选择Cohesive,对其他单元,Family 选择3D Stress;对于静态运算,Element Library选择Standard,对于动态(显式)运算,Element Library 选择Explicit。 10,Property: Create Material: jiti (材料名字):Mechanical : Elastic: Type: Isotropic =tdrt Matetial 邑 M<)terial-jiti Description; NLrnnb?r of field v-arid4)l?:0 ' Moduli tme scale [forvi&ctwlKlicrty^ Long-term No compr-eision 3 Nc Datia Voungi'i P鈕1刖n1* 1 4D0C Create Material: xianwei (材料名字):Mechanical : Elastic : Type : Isotropic
《现代机械设计方法》课程结业论文 ( 2011 级) 题目:ABAQUS 实例分析 学生姓名XXXX 学号XXXXX 专业机械工程 学院名称机电工程与自动化学院 指导老师XX 2013 年 5 月8 日
目录 第一章Abaqus 简介 (1) 一、Abaqus 总体介绍 (1) 二、Abaqus 基本使用方法 (2) 1.2.1Abaqus 分析步骤 (2) 1.2.2Abaqus/CAE 界面 (3) 1.2.3Abaqus/CAE 的功能模块 (3) 第二章基于Abaqus 的通孔端盖分析实例 (4) 一、工作任务的明确 (4) 二、具体步骤 (4) 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) 第三章课程学习心得与作业体会 (23)
第一章: Abaqus 简介 一、 Abaqus 总体介绍 Abaqus 是功能强大的有限元分析软件,可以分析复杂的固体力学和结构力学系统,模拟非常庞大的模型,处理高度非线性问题。Abaqus 不但可以做单一零件的力学和多物理场的分析,同时还可以完成系统级的分析和研究。 Abaqus 使用起来十分简便,可以很容易的为复杂问题建立模型。Abaqus 具备十分丰富的单元库,可以模拟任意几何形状,其丰富的材料模型库可以模拟大多数典型工程材料的性能,包括金属、橡胶、聚合物、复合材料、钢筋混泥土、可压缩的弹性泡沫以及地质材料(例如土壤、岩石)等。 Abaqus 主要具有以下分析功能: 1.静态应力/位移分析 2.动态分析 3.非线性动态应力/位移分析 4.粘弹性/粘塑性响应分析 5.热传导分析 6.退火成形过程分析 7.质量扩散分析 8.准静态分析 9.耦合分析 10.海洋工程结构分析 11.瞬态温度/位移耦合分析 12.疲劳分析 13.水下冲击分析 14.设计灵敏度分析 二、 Abaqus 基本使用方法 1.2.1Abaqus 分析步骤 有限元分析包括以下三个步骤: 1.前处理(Abaqus/CAE):在前期处理阶段需要定义物理问题的模型,并生 成一个 Abaqus 输入文件。提交给 Abaqus/Standard 或 Abaqus/Explicit。 2.分析计算(Abaqus/Standard 或 Abaqus/Explicit):在分析计算阶段, 使用 Abaqus/Standard 或Abaqus/Explicit 求解输入文件中所定义的
一、定义材料的刚度矩阵 从弹性力学理论可以知道,各向异性材料的刚度矩阵由于有对称性,刚度系数有最初的36个减少到21个,如下图: 在实际应用中,大多数工程材料都有对称的内部结构,因此材料具有弹性对称性,这种对称性可以进一步简化上述的刚度矩阵。 1、有一个弹性对称面的材料(如结晶学中的单斜体) 例如取x-y平面为对称面,则D1112= D1113= D2212= D2213= D3312= D3313= D1223= D1323=0,刚度系数又减少8个,剩下13个。 2、有两个正交(相互垂直)弹性对称面的材料 例如进一步取x-z平面为对称面,则D1123= D2223= D3323= D1213=0,刚度系数又减少4个,剩下9个,如下图: 在Abaqus编辑材料中进行个刚度系数的设定。
3、有三个正交弹性对称面的材料 如果材料有三个相互垂直的弹性对称面,没有新的刚度系数为零,也只有9个。 4、横观各项同性材料 若经过弹性体材料一轴线,在垂直该轴线的平面内,各点的弹性性能在各方向上都相同,我们称此材料横观各向同性材料,如单向复合材料。对于这种材料最终的刚度系数只剩下D1111,D1122,D1133,D3333,D1212五项,其余各项均为零。在复合材料中,经常遇到正交各项异性和横观各项同性两种材料。 二、定义材料工程弹性常数 通过指定工程弹性常数定义线弹性正交各向异性材料是最便捷的一种方法,根据复合材料力学理论,用工程弹性常数表示的柔度矩阵表示如下:
其中,γij/Ei=γji/Ej,所以用9个独立弹性常数可以表征材料属性,即三个材料主 方向上的弹性模量E1,E2,E3,三个泊松比γ12,γ13,γ23,三个平面内的剪切弹性模量G12,G13,G23。 例如测得复合材料一组材料数据为:E1=39GPa,E2=8.4GPa,E3=5.2GPa,γ12=0.26,γ13=0.3,γ23=0.28,G12=4.2GPa,G13=3.6GPa,G23=2.4GPa (随便给出的)。在Abaqus编辑材料对话框中输入对应数据,完成正交各向异性材料的定义。 对于横观各向同性材料,E1=E2,γ13=γ23,γ31=γ32,G13=G23,弹性常数
材料本构模型及编程-ABAQUS-UMAT 材料本构模型及编程实现:简介 1、什么时候用用户定义材料(User-defined material, UMAT)? 很简单,当ABAQUS没有提供我们需要的材料模型时。所以,在决定自己定义一种新的材料模型之前,最好对ABAQUS已经提供的模型心中有数,并且尽量使用现有的模型,因为这些模型已经经过详细的验证,并被广泛接受。 2、好学吗?需要哪些基础知识? 先看一下ABAQUS手册(ABAQUS Analysis User's Manual)里的一段话: Warning: The use of this option generally requires considerable expertise. The user is cautioned that the imple mentation of any realistic constitutive model requires extensive development and testing. Initial testing on a s ingle element model with prescribed traction loading is strongly recommended. 但这并不意味着非力学专业,或者力学基础知识不很丰富者就只能望洋兴叹,因为我们的任务不是开发一套完整的有限元软件,而只是提供一个描述材料力学性能的本构方程(Constitutive equation)而已。当然,最基本的一些概念和知识还是要具备的,比如 应力(stress),应变(strain)及其分量;volumetric part和deviatoric part;模量(modulus)、泊松比(Poisson’s ratio)、拉美常数(Lame constant);矩阵的加减乘除甚至求逆;还有一些高等数学知识如积分、微分等。 3、UMAT的基本任务? 我们知道,有限元计算(增量方法)的基本问题是: 已知第n步的结果(应力,应变等),;然后给出一个应变增量, 计算新的应力。UMAT要完成这一计算,并要计算Jacobian矩阵DDSDDE(I,J) =。是应力增量矩阵(张量或许更合适),是应变增量矩阵。DDSDDE(I,J) 定义了第J个应变分量的微小变化对第I 个应力分量带来的变化。该矩阵只影响收敛速度,不影响计算结果的准确性(当然,不收敛自然得不到结果)。 4、怎样建立自己的材料模型? 本构方程就是描述材料应力应变(增量)关系的数学公式,不是凭空想象出来的,而是根据实验结果作出的合理归纳。比如对弹性材料,实验发现应力和应变同步线性增长,所以用一个简单的数学公式描述。为了解释弹塑性材料的实验现象,又提出了一些弹塑性模型,并用数学公式表示出来。 对各向同性材料(Isotropic material),经常采用的办法是先研究材料单向应力-应变规律(如单向拉伸、压缩试验),并用一数学公式加以描述,然后把讲该规律推广到各应力分量。这叫做“泛化“(generalization)。 5、一个完整的例子及解释 下面这个UMAT取自ABAQUS手册,是一个用于大变形下的弹塑性材料模型。希望我的注释能帮助初学者理解。需要了解J2理论。SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT, 1 DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED, 2 CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT, 3 PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC) STRESS--应力矩阵,在增量步的开始,保存并作为已知量传入UMAT ;在增量步的结束应该保存更新的应力; STRAN--当前应变,已知。 DSTRAN—应变增量,已知。 STATEV--状态变量矩阵,用来保存用户自己定义的一些变量,如累计塑性应变,粘弹性应变等等。增量步开始时作为已知量传入,增量步结束应该更新; DDSDDE=。需要更新 DTIME—时间增量dt。已知。 NDI—正应力、应变个数,对三维问题、轴对称问题自然是3(11,22,33),平面问题是2(11,22);已知。 NSHR —剪应力、应变个数,三维问题时3(12,13,23),轴对称问题是1(12);已知。
ABAQUS6.13百度网盘!附中英文安装教程,win8.1测试通 过,需要的人拿去! 作者: 909933622(站内联系TA)发布: 2013-10-11 闲话不说,直接上链接 https://www.doczj.com/doc/1116158426.html,/s/122MWa abaqus 6.13对操作系统的新要求: 自abaqus 6.13版本开始,将不再支持windows 的32位操作平台; 同时,也不再支持windows xp和windows vista操作系统; 安装: 1. Run "Install Abagus Product & Licensing" 2. In SIMULIA FLEXnet License Server window select "Just install the license utilities " NOTE: If you already have SIMULIA FLEXnet License Server for ABAQUS 6.12-3 inst alled and running you can use it for 6.13-1 too 3. After finishing License Utilities setup copy files "ABAQUS.lic" and "ABAQUS.log" to
基于ABAQUS复合材料薄壁圆筒的屈曲分析 由于玻璃钢复合材料的薄壁圆筒结构具有强度高、重量轻、刚度大、耐腐蚀,电绝缘及透微波等优点,目前已广泛应用于航空航天和民用领域中。工程中广泛使用的这些薄壁圆筒,当它们受压缩、剪切、弯曲和扭转等荷载作用时,最常见的失效模式为屈曲。因此,为了保证结构的安全,需要进行屈曲分析。 对结构进行屈曲分析,涉及到较复杂的弹(塑)性理论和数学计算,要通过求解高阶偏微分方程组,才能求解失稳临界荷载,而且只有少数简单结构才能求得精确的解析解。因此,只能采用能量法、数值方法和有限元方法等近似的分析方法进行分析。近20年来,随着计算机和有限元方法的迅猛发展,形成了许多的实用分析程序,提高了对复杂结构进行屈曲分析的能力和设计水平。ABAQUS 就是其中的杰出代表。 1.屈曲有限元理论 有限元方法中,对结构的屈曲失稳问题的分析方法主要有两类:一类是通过特征值分析计算屈曲载荷,另一类是利用结合Newton—Raphson迭代的弧长法来确定加载方向,追踪失稳路径的几何非线性分析方法,能有效分析高度非线性屈曲和后屈曲问题。 1.1线性屈曲 假设结构受到的外载荷模式为。,幅值大小为,结构内力为Q,则静力平衡方程应为 进一步考察结构在载荷作用下的平衡方程,得到 由于结构达到保持稳定的临界载荷时有,代入上式得 该方程对应的特征值问题为 如果忽略几何刚度增量的影响,屈曲分析的方程又可进一步简化为 该方程即为求解线性屈曲的特征值方程。为屈曲失稳载荷因子,为结构失稳形态的特征向量。
1.2非线性屈曲 非线性屈曲分析方法多采用弧长法进行分步迭代计算,在增量非线性有限元分析中,沿着平衡路径迭代位移增量的大小(也叫弧长)和方向,确定载荷增量的自动加载方案,可用于高度非线性的屈曲失稳问题。与提取特征值的线性屈曲分析相比,弧长法不仅考虑刚度奇异的失稳点附近的平衡,而且通过追踪整个失稳过程中实际的载荷、位移关系,获得结构失稳前后的全部信息,适合于高度非线性的屈曲失稳问题。 2.ABAQUS的线性屈曲分析 ABAQUS中提供两种分析方法来确定结构的临界荷载和结构发生屈曲响应的特征形状:线性屈曲分析(特征值屈曲分析)、非线性屈曲分析。 线性屈曲分析用于预测一个理想的弹性结构的理论屈曲强度。它是预期的线性屈曲荷载的上限,可以作为非线性屈曲分析的给定荷载,在渐进加载达到此荷载前,非线性求解必然发散;它还可以作为施加初始缺陷或扰动荷载的依据。所以预先进行特征值屈曲分析有助于非线性屈曲分析,进行特征值屈曲分析是必要的。 3.算例 3.1问题概述 图3-1 实例模型 如图所示两端开口的复合材料薄壁圆筒,底端固支,顶端作用有均匀分布的轴压边载。半径R=152mm,高度300mm,厚度t=0.804mm,对称铺层[±45,0]s,
复合材料模型建模与分析 1.Cohesive单元建模方法 1.1 几何模型 使用内聚力模型(cohesive zone)模拟裂纹的产生和扩展,需要在预计产生裂纹的区域加入cohesive层。建立cohesive层的方法主要有: 方法一、建立完整的结构(如图1(a)所示),然后在上面切割出一个薄层来模拟cohesive 单元,用这种方法建立的cohesive单元与其他单元公用节点,并以此传递力和位移。 方法二、分别建立cohesive层和其他结构部件的实体模型,通过“tie”绑定约束,使得cohesive单元两侧的单元位移和应力协调,如图1(b)所示。 (a)cohesive单元与其他单元公用节点(b)独立的网格通过“tie”绑定 图1.建模方法 上述两种方法都可以用来模拟复合材料的分层失效,第一种方法划分网格比较复杂;第二种方法赋材料属性简单,划分网格也方便,但是装配及“tie”很繁琐;因此在实际建模中我们应根据实际结构选取较简单的方法。 1.2 材料属性 应用cohesive单元模拟复合材料失效,包括两种模型:一种是基于traction-separation描述;另一种是基于连续体描述。其中基于traction-separation描述的方法应用更加广泛。 而在基于traction-separation描述的方法中,最常用的本构模型为图2所示的双线性本构模型。它给出了材料达到强度极限前的线弹性段和材料达到强度极限后的刚度线性降低软化阶段。注意图中纵坐标为应力,而横坐标为位移,因此线弹性段的斜率代表的实际是cohesive单元的刚度。曲线下的面积即为材料断裂时的能量释放率。因此在定义cohesive 的力学性能时,实际就是要确定上述本构模型的具体形状:包括刚度、极限强度、以及临界断裂能量释放率,或者最终失效时单元的位移。常用的定义方法是给定上述参数中的前三项,也就确定了cohesive的本构模型。Cohesive单元可理解为一种准二维单元,可以将它看作被一个厚度隔开的两个面,这两个面分别和其他实体单元连接。Cohesive单元只考虑面外的力,包括法向的正应力以及XZ,YZ两个方向的剪应力。 下文对cohesive单元的参数进行阐述,并介绍参数的选择方法。
纤维增强材料的累积损伤与失效:Abaqus拥有纤维增强材料的各向异性损伤的建模功能(纤维增强材料的损伤与失效概论,19.3.1节)。假设未损伤材料为线弹性材料。因为该材料在损伤的初始阶段没有大量的塑性变形,所以用来预测纤维增强材料的损伤行为。Hashin标准最开始用来预测损伤的产生,而损伤演化规律基于损伤过程和线性材料软化过程中的能量耗散理论。 另外,Abaqus也提供混凝土损伤模型,动态失效模型和在粘着单元以及连接单元中进行损伤与失效建模的专业功能。 本章节给出了累积损伤与失效的概论和损伤产生与演变规律的概念简介,并且仅限于塑性金属材料和纤维增强材料的损伤模型。 损伤与失效模型的通用框架 Abaqus提供材料失效模型的通用建模框架,其中允许同一种的材料应用多种失效机制。材料失效就是由材料刚度的逐渐减弱而引起的材料承担载荷的能力完全丧失。刚度逐渐减弱的过程采用损伤力学建模。 为了更好的了解Abaqus中失效建模的功能,考虑简单拉伸测试中的典型金属样品的变形。如图19.1.1-1中所示,应力应变图显示出明确的划分阶段。材料变形的初始阶段是线弹性变形(a-b段),之后随着应变的加强,材料进入塑性屈服阶段(b-c段)。超过c点后,材料的承载能力显著下降直到断裂(c-d段)。最后阶段的变形仅发生在样品变窄的区域。C点表明材料损伤的开始,也被称为损伤开始的标准。超过这一点之后,应力-应变曲线(c-d)由局部变形区域刚度减弱进展决定。根据损伤力学可知,曲线c-d可以看成曲线c-d‘的衰减,曲线c-d‘是在没有损伤的情况下,材料应该遵循的应力-应变规律曲线。
图19.1.1-1 金属样品典型的轴向应力-应变曲线 因此,在Abaqus中失效机制的详细说明里包括四个明显的部分: ●材料无损伤阶段的定义(如图19.1.1-1中曲线a-b-c-d‘) ●损伤开始的标准(如图19.1.1-1中c点) ●损伤发展演变的规律(如图19.1.1-1中曲线c-d) ●单元的选择性删除,因为一旦材料的刚度完全减退就会有单元从计算中移除(如图19.1.1-1 中的d点)。 关于这几部分的内容,我们会对金属塑性材料(金属塑性材料的损伤与失效概论,19.2.1节)和纤维增强材料(纤维增强符合材料的损伤与失效概论,19.3.1节)进行分开讨论。 网格依赖性 在连续介质力学中,通常是根据应力-应变关系建立材料本构模型。当材料表现出导致应变局部化的应变软化行为时,有限元分析的结果带有强烈的网格依赖性,能量的耗散程度取决于网格的精简程度。在Abaqus中所有可使用损伤演化模型都使用减轻网格依赖性的公式。这是通过在公式中引入特征长度来实现的,特征长度作为一个应力-位移关系可以表达本构关系中软化部分,它与单元尺寸有关系。在此情况下,损伤过程中耗散的能量不是由每个单位体积衡量,而是由每个单位面积衡量。这个能量值作为另外一个材料参数,用来计算材料发生完全损伤时的位移。这是与材料断裂力学中临界能量释放率的概念一致的。此公式确保了合适能量的耗散以及最大程度减轻网格的依赖。
复合材料Abaqus仿真分析——精讲版 本文以一个非常简单的复合材料层合板为例,应用Abaqus/CAE对其进行线性静态分析。一块边长为254mm的方形两层层合板,两层厚度均为2.54mm,第一层铺层角45°,第二层铺层角-45°;板的四边完全固支,板的上表面受到689.4kpa的压强。各单层的材料相同,材料属性如下: E1=276GPa,E2=6.9GPa,E3=5.2GPa,γ12=0.25,G12=3.4GPa,G13=3.4GPa,G23=3.4G。 定义模型的几何形状 创建一个具有平面壳体单元基本特征的三维变形体,在草图环境绘制板的几何形状如下图:
定义材料属性和局部材料方向 Create coordinate system
定义局部坐标系,对于像本例这样的简单几何体,本可以不用另外建立局部坐标系,但笔者还是在本例中用了局部坐标系,主要是考虑到以后再复杂问题中会经常用到这一方法。 创建铺层 或者使用菜单栏
此处使用全局坐标系
使用用户自定义坐标系 Rotation angle depends on the coordinate system defined by user. Par example, if x-axe in the user defined system is parallel to the direction of fiber; we should replace the angles by 0 and 90. 使用全局坐标系和局部坐标系的区别在下面这一步可以查看 如果使用全局坐标系,会有方向指示,如果使用用户自定义坐标系,在层中没有方向指示可以通过’工具——查询’来检查铺层(Tool ---- Q uery----ply stack plot) Case 1 全局坐标系