[整理]shell 单元厚度 复合材料.
nastran单元类型

nastran单元类型摘要:1.Nastran单元类型概述2.常见Nastran单元类型的特点及应用3.Nastran单元类型选择策略4.总结正文:astran是一种广泛应用于工程分析的有限元分析(FEA)软件。
Nastran 提供了多种单元类型,以适应各种工程领域的需求。
本文将介绍Nastran单元类型的概述、常见单元类型的特点及应用,以及如何选择合适的单元类型。
1.Nastran单元类型概述astran单元类型主要分为以下几类:(1)结构单元:用于模拟结构构件,如梁、柱、壳等。
(2)热传导单元:用于模拟热传导问题,如杆、板、壳等。
(3)热膨胀单元:用于模拟材料在温度变化下的膨胀或收缩。
(4)动力学单元:用于模拟振动、冲击等问题,如质量、弹簧、阻尼等。
(5)复合单元:用于模拟复合材料结构,如夹层板、网格结构等。
2.常见Nastran单元类型的特点及应用(1)BEAM(梁)单元:BEAM单元用于模拟简支梁、连续梁、悬臂梁等结构。
它具有弯曲、剪切、扭转等多种变形模式。
(2)PLATE(板)单元:PLATE单元用于模拟薄板结构,如平板、曲面板等。
它具有弯曲、剪切、拉伸等多种变形模式。
(3)SHELL(壳)单元:SHELL单元用于模拟壳状结构,如圆筒、球壳等。
它具有弯曲、剪切、拉伸等多种变形模式。
(4)COMPOSITE(复合)单元:复合单元用于模拟复合材料结构,如夹层板、网格结构等。
它可以模拟不同材料层的铺设角度、厚度等参数。
(5)THERMAL(热)单元:THERMAL单元用于模拟热传导问题。
它可以模拟一维、二维和三维热传导问题,以及材料的热膨胀特性。
3.Nastran单元类型选择策略(1)根据问题的实际情况选择合适的单元类型:例如,对于结构分析,可选择BEAM、PLATE或SHELL单元;对于热传导问题,选择THERMAL单元。
(2)考虑单元的阶次:高阶单元具有更高的精度和灵活性,但计算成本也较高。
单元类型

单元类型维数单元名称结构点单元1/2/3D MASS21结构线单元2D LINK13D LINK8、10、11、180J结构梁单元(beam)2D Beam3、23、543D BEAM4、24,44,188,189结构实体单元(solid)2D PLANE2,25,42,82,83,145,146,182,1833D SOLID45,64,65,92,95,147,148,185,186,187结构壳单元(shell) 2D SHELL51,613D SHELL28,41,42,63,93,143,150,181 结构管道单元(pipe) 3D PIPE16,17,18,20,59,60结构密封垫单元(gasket) 3D INTER192,193,194,195结构多点约束单元(rigidlink/beam)3D MPC184结构层复合材料单元(layered composite)3D SOLID46,91,99,191,显示动力学单元(explicit dynamic)Link160,beam161,plane162,shell163,s olid164,comni165,mass166,link167超弹实体单元(hyperelastic solid)3D Hyper56,58,74,84,86,158粘性实体单元(viscosolid)3D VISCO88,89,106,107,108热点单元MASS71热线单元LINK31,32,33,34热实体单元(solid)2D PLANE35,55,75,77,783D SOLID70,87,90热壳单元(shell)SHELL57,131,132热电单元(thermalelectric)PLANE67,6869,157流体单元(fluid)FLUID29,30,38,79,80,81,116,129,130,141,142电磁单元(magnetic electric)Plane53,solid96,solid97,Inter115,solid117,HF118,HF119,HF120,PLANE121,SO LID122,SOLID123,SOLID127,SOLI D128电路单元(electriccircuit)Sourc36,circu94,124,125机电转换单元(electric-mechanical)Trans109,126耦合场单元(coupled Solid5,plane13,solid62,solid98,rom14field) 4接触单元(contact)Contact12,26,48,49,52,TRAGE169,170CONTACT171,172,173,174,175,178矩阵单元(matrix)Matrix27,50无限边界单元(infinite)INFIN9,47,110,111表面效应单元(surface)SURF151,152,153,154Mesh200网格辅助单元(仅用于生成网格拓扑生成,不参与求解计算)初学ANSYS的人,通常会被ANSYS所提供的众多纷繁复杂的单元类型弄花了眼,如何选择正确的单元类型,也是新手学习时很头疼的问题。
Abaqus针对复合材料的优势

四Abaqus在复合材料领域的优势4.1 复合材料介绍4.1.1 复合材料的应用复合材料有许多特性:1、制造工艺简单2、比强度高,比刚度大3、具有灵活的可设计性4、耐腐蚀,对疲劳不敏感5、热稳定性能、高温性能好由于复合材料的上述优点,在航空航天、汽车、船舶等领域,都有广泛的应用。
复合材料的大量应用对分析技术提出新的挑战。
4.1.2 复合材料的结构复合材料是一种至少由两种材料混合而成的宏观材料,其中的一种材料被称作基体,其它的材料称作纤维。
其中纤维可以包含很多不同的形式:离散的宏观粒子,任意方向的短纤维,规则排列的纤维和织物。
1)单向纤维层合板----冲击分析2)编织复合材料---- 挤压分析3)蜂窝夹心复合材料----不可见冲击损伤分析基体和纤维的存在形式以及材料属性对于复合材料的力学行为有着很大的影响。
改变纤维和基体的属性目的就是在于生成一种复合材料具有如下性质:1)低成本:原型,大规模生产,零件合并,维修,技术成熟。
2)期望的重量:轻重量,比重分配合理。
3)改进的强度和刚度:高强度/高刚度比。
4)改进的表面属性:良好的耐腐蚀性,表面抛光性好。
5)期望的热属性:较低的热传导性,热膨胀系数较低。
6)独特的电属性:具有较高的绝缘强度,无磁性。
7)空间适应性:大部件,特殊的几何构型。
4.1.4 复合材料的有限元模拟根据不同的分析目的,可以采用不同的复合材料模拟技术:1)微观模拟:将纤维和基体都分别模拟为可变形连续体。
2)宏观模拟:将复合材料模拟为一个正交各向异性体或是完全各向异性体。
3)混合模拟:将复合材料模拟为一系列离散、可见的纤维层合板。
4)离散纤维模拟:采用离散单元或是其它模拟工具进行模拟。
5)子模型模拟:对于研究加强纤维周围点的应力集中问题比较有效。
微观模拟:纤维-基体的单胞模拟混合模拟:层合板的混合模拟Abaqus中复合材料的单元技术Abaqus中复合材料的单元技术主要为三种:分层壳单元、分层实体单元以及实体壳单元。
ANSYS ACP复合材料前后处理软件

ANSYS Composite PrepPost Modeling Composites the Simple Way安世亚太科技股份有限公司演讲人:刘程伟复合材料简介ANSYS复合材料分析技术 ACP介绍ACP特色功能ACP新功能复合材料简介ANSYS复合材料分析技术 ACP介绍ACP特色功能ACP新功能•复合材料(Composite materials),是由两种或两种以上不同性质的材料,通过物理或化学的方法,在宏观上组成具有新性能的材料。
•各种材料在性能上取长补短复合材料的分类A.颗粒增强复合材料B.纤维增强复合材料•短纤维增强•长纤维增强复合材料薄片- 由基体和增强体组成的单层结构。
材料属性采用等效的方式。
基体- 一般采用各向同性材料,起到包裹增强体的作用。
纤维增强体- 被基体材料包裹,是复合材料各向异性力学行为的主要原因。
LT叠层组合- 多组复合材料薄片组成,纤维的方向一般不同应用领域•航空航天•风能•运动和娱乐•建筑行业•汽车行业•船舶行业•国防领域•……目录复合材料简介ANSYS复合材料分析技术 ACP介绍ACP特色功能ACP新功能ANSYS 复合材料分析技术•复合材料单元:•Mechanical APDL包含多种复合单元•不同截面属性的3D 梁单元–BEAM188, BEAM189, ELBOW290 elements•2D 轴对称壳单元–SHELL208, SHELL209 elements•3D 铺层壳单元–SHELL181, SHELL281 elements (SHELL131 and SHELL132 are thermal shells)•3D 铺层实体单元–SOLID185, SOLID186, SOLSH190,SOLID278,SOLID279 elementsANSYS 复合材料分析技术•厚度方向仅采用一个单元,模拟复合材料的铺层力学性能–不需要单独的对每层划分网格•可以采用多种失效准则,对复合材料的强度进行评估–最大应变准则–最大应力准则–Tsai-Wu 准则–用户自定义准则•更多的内容参看ANSYS–“Mechanical APDL (formerly ANSYS) > Structural Analysis Guide > Ch. 13 Composites”–“Mechanical APDL (formerly ANSYS) > Structural Analysis Guide > Ch. 16 Beam Analysis and Cross Sections”–“Mechanical APDL (formerly ANSYS) > Structural Analysis Guide > Ch. 17 Shell Analysis and Cross Sections”ET,1,SHELL181! LAYERS PROPERTIESSECTYPE,1,SHELL,,MON_STRATIFIE SECDATA,0.025,1, 0SECDATA,0.500,1, 45SECDATA,0.500,1, -45 SECDATA,0.025,1, 0! ORTHOTROPIC MATERIAL PROPERTIES MP,EX,1,25E6MP,EY,1,1E6MP,EZ,1,1E6MP,GXY,1,5E5MP,GYZ,1,2E5MP,GXZ,1,5E5MP,PRXY,1,0.25MP,PRYZ,1,0.01 SECPLOT,1 (or LAYPLOT command)/PSYMB,ESYS,1THETA = Angle (in degrees) withrespect to element coordinatesystem (ESYS)/ESHAPE,1/EFACET,2EPLOTANSYS 复合材料分析技术目录复合材料简介ANSYS复合材料分析技术 ACP介绍ACP特色功能ACP新功能ANSYS Composite PrepPost介绍ANSYS ACP的主要功能•创建各种形式的复合材料模型•定义复合材料的铺层设置•定义复合材料的纤维方向•对复合材料结构进行评估•评估每层的应力情况•计算复合材料的失效ANSYS ACP分析流程建立几何模型施加边界条件ACP复合材料前处理:铺层信息定义(每层材料属性、厚度、铺层方向角等)ANSYS求解ACP后处理为复合材料加工制作提供必要数据ACP与ANSYS数据传递ANSYS WorkbenchSOLVER定义几何以及边界条件定义铺层信息Write an acp.cdbFileRead the acp.cdbFileImport the results (rstfiles) for postprocessingWrite a *.cdb file withdefined lay-upsLaunch solver其他方式后处理ANSYS Mechanical APDL前处理后处理求解定义材料属性•材料属性的定义同样在ANSYS材料库中进行•定义材料属性的过程中需要考虑纤维的影响•必须的材料属性设置•x, y, z三个方向的杨氏模型•xy, yz, xz三个方向的剪切模量•xy, yz, xz三个方向的泊松比1, x 2, y3, z•失效准则需要定义应力极限值和应变极限值•应力极限值和应变极限值在拉伸和压缩方向一般是不同的应力极限应变极限拉伸X, Y 和 Z 拉伸X, Y 和 Z压缩X, Y 和 Z 压缩X, Y 和 Z剪切XY, YZ 和 XZ 剪切XY, YZ 和 XZ定义材料属性•网格划分过程与普通静力学通用•Mesh中所有的尺寸和选项控制都可以使用•复合材料的建模分析,以面体的网格为起点•ACP同样能够进行实体复合材料的分析•单元集通过Named Selections进行定义•基于Named Selections定义复合材料的铺层•更改几何模型时候需要检查Named Selections的定义单元集Named SelectionsANSYS Composite PrepPost 分析过程•开始复合材料的分析ANSYS Composite PrepPost 界面标准视图模式和视图设置模型树后处理显示设置铺层设置•复合材料铺层设置•采用的纤维层(fabric)•确定需要铺层的位置•铺层的参考方向•增强纤维的方向Static Structural 分析求解ACP (Post) 后处理•失效准则•Max. Strain & Max. Stress •Tsai-Wu•Tsai-Hill•Hashin•Puck•LaRC•Cuntze•Face Sheet Wrinkling•Core Failure•可以在后处理中查看•失效准则•失效模式•关键层•临界载荷步s2t(5) s2t(5)s2t(5)s2t(5)ACP基本功能总结•完全集成于WorkBench平台•直观的创建复合材料铺层•依据制造的过程创建模型•简单高效的修改铺层设置•包含目前通用的失效准则•高效的后处理过程目录复合材料简介ANSYS复合材料分析技术 ACP介绍ACP特色功能ACP新功能Oriented Element Sets利用“方向化单元集”来定义复杂的铺层方向(OES :Oriented Element Sets) the OES normal the OES reference direction the ply angle直观的铺层定义(Build the laminate lay-up)铺层材料铺层方向角铺层数量查看铺层截面,确认铺层正确无误提供复合材料“Draping”功能Draping coefficients提供多种失效模式定义复合材料失效通常有以下两种:●层间失效(由于剪切或拉伸力作用引起)●层内的纤维或基体破坏失效failure patternout-of-phase micro buckling in-phasemicro bucklingshear failurefiber tensionfiber compression fiber ruptureεm > εfmatrix failureεf > εmtransverse crack interfacial failureshear failure with fiber rupture shear w/o fiber rupturematrix tension matrix compressionACP 强大的后处理功能 Agarwal/Broutman (1990) Agarwal/Broutman (1990) delaminationmatrix crack fiber rupturedebondingmicro-buckling / shear failure实体模型•复合材料模型通常是薄壳类结构•壳体厚度方向应力为零实体模型•壳体厚度方向应力不为零•厚度方向应力产生的效果与垂直于纤维方向应力效果相同,很小的应力会造成复合材料的失效•以下情况建议采用实体复合材料模型•多层复合材料•大变形•厚度方向存在载荷•实体模型的分析流程•实体复合材料模型创建过程类似于制造过程。
【ANSYS培训】单元库及常用单元、材料库、高级有限元模型技术

梁网格划分
• 在单元绘图中显示横截面形状, 激活显示单元形状:
– Utility Menu > PlotCtrls > Style > Size and Shape… – 或使用命令 /ESHAPE,1
ANSYS TRAINING
梁网格划分
• 网格划分完成后,接着施加荷载并求解
ANSYS TRAINING
ANSYS TRAINING
• 建筑框架
练习
ANSYS TRAINING
壳的高级使用
ANSYS TRAINING
壳用于薄面板或曲面的模型,壳分析 应用的基本原则是每块面板的主尺寸不低 于其厚度的10倍。
ANSYS TRAINING
壳的使用
薄壳 (Love/Kirchhoff)
– 假设最初垂直于壳的中面的横截面 在加载过程中保持平直并垂直于中 性轴,该假设不包括剪切变形。
• 对非线性分析,不同的壳单元在面内和面外有不同的积分 点数。例如, 具有 URI 的SHELL181 在面内有1个积分点 ,在厚度方向上有5 个。特定单元的详细说明可参阅理论 手册第14 章。注意SHELL181 在厚度上可以有用户定义 的积分点。
ANSYS TRAINING
SHELL181
• SHELL181 属于 18x 系列单元,由于其丰富的特 点,对非线性应用建议选择该单元。
梁属性
• “自定义”横截面形状步骤: 建立面 指定线分割数 保存截面 编辑截面(若需要,如定义复合材料截面等) 读入截面(使用时)
举例:有圆孔的方形梁 截面上的不同材料
ANSYS TRAINING
梁属性
材料属性 • 线性和非线性材料属性均可 • 所有梁的属性定义好以后,下一步是对几何模型进
abaqus单元属性大总结

S3/S3R 单元可以作为通用壳单元使用。
由于单元中的常应变近似,需要划分较细的网格来模拟弯曲变形或高应变梯度。
S4R 单元性能稳定,适用围很广对于复合材料,为模拟剪切变形的影响,应使用适于厚壳的单元(例如S4、S4R、S3、S3R、S8R),并要注意检查截面是否保持平面。
对于几何非线性分析,在ABAQUS/Standard中的小应变壳单元(S4R5, S8R, S8R5, S8RT, S9R5, STRI3, 和STRI65)使用总体拉格朗日应变算法,应力应变可以相对于参考构型的材料方向改定。
垫片单元是小应变小位移单元,默认情况下其应力应变值也是以初始参考构型定义的行为方向输出。
对于有限膜应变单元(所有的膜单元以及S3/S3R, S4, S4R, SAX,和 SAXA单元)和在ABAQUS/Explicit中的小应变单元,其材料方向是随着曲面的平均刚性旋转运动而变以形成当前构型的材料方向。
此时这些单元的应力应变则是根据当前的参考构型中的材料方向给出的。
(更详细地说明可以参考ABAQUS相关手册)。
用户可以决定与*section print和*section file相关的局部坐标系统是固定不动还是随着曲面的平均刚性运动而旋转。
C3D20RP 20-node brick, triquadratic displacement, trilinear pore pressure, reducedintegration20节点实体,三重二次位移,三线孔隙压力,缩减积分C3D20R PH 20-node brick, triquadratic displacement, trilinear pore pressure, hybrid, linearpressure, reduced integration20节点实体,三重二次位移,三线孔隙压力,混合动力,线压力,缩减积分。
圆锥罩聚能射流LSDYNA仿真计算(图解)
*MAT_HIGH_EXPLOSIVE_BURN //注:该关键字用于炸药参数设定,须通过编辑、修改 DYNA 的数据输入文件 Jobname.K 来进行加入。
MID, RO, D, PCJ, BETA, K, G, SIGY
其中:
MID 为材料模型的 ID 号;
RO 为材料的密度;
D 为材料的爆速;
选择 Main Menu > Preprocessor > Material Props > Material Models 命令,弹 出 Define Material Model Behavior 对话框。左侧框图中,显示为材料模型代号,有多 种材料时,点击菜单 Material > New Model… 命令,弹出 Define Material ID 对话框, 输入材料 ID,即可新建材料。右侧框图的 Material Models Available 选择栏对应的是
QQ:759195361
第 4 页 共 44 页 Email:qingzhouwang@
ANSYS/LS-DYNA
在对话框的 Library of Element Types 下拉列表框中选择需要的单元类型,单击 OK 退 出即可。图示如下:
相应关键字: *SECTION_SOLID(_ALE) SECID, ELFORM, AET SECID 为截面单元 ID 号; ELFORM 为实体单元的算法公式,有以下选项: 0:用于模拟蜂窝材料; 11:单点 Euler/ALE 多物质单元(一个单元最多含 10 种材料)。
QQ:759195361
第 3 页 共 44 页 Email:qingzhouwang@
ANSYS/LS-DYNA
ansys关于薄板、厚板、壳单元的特性区别
一、板壳弯曲理论简介1. 板壳分类按板面内特征尺寸与厚度之比划分:当L/h < (5~8) 时为厚板,应采用实体单元。
当(5~8) < L/h < (80~100) 时为薄板,可选2D 实体或壳单元当L/h > (80~100) 时为薄膜,可采用薄膜单元。
壳类结构按曲率半径与壳厚度之比划分:当R/h >= 20 时为薄壳结构,可选择薄壳单元。
当6 < R/h < 20 时为中厚壳结构,选择中厚壳单元。
当R/h <= 6 时为厚壳结构。
上述各式中h 为板壳厚度,L 为平板面内特征尺度,R 为壳体中面的曲率半径。
2. 薄板理论的基本假定薄板所受外力有如下三种情况:①外力为作用于中面内的面内荷载。
弹性力学平面应力问题。
②外力为垂直于中面的侧向荷载。
薄板弯曲问题。
③面内荷载与侧向荷载共同作用。
所谓薄板理论即板的厚度远小于中面的最小尺寸,而挠度又远小于板厚的情况,也称为古典薄板理论。
薄板通常采用Kirchhoff-Love 基本假定:①平行于板中面的各层互不挤压,即σz = 0。
②直法线假定:该假定忽略了剪应力和所引起的剪切变形,且认为板弯曲时沿板厚方向各点的挠度相等。
③中面内各点都无平行于中面的位移。
薄板小挠度理论在板的边界附近、开孔板、复合材料板等情况中,其结果不够精确。
3. 中厚板理论的基本假定考虑横向剪切变形的板理论,一般称为中厚板理论或Reissner(瑞斯纳)理论。
该理论不再采用直法线假定,而是采用直线假定,同时板内各点的挠度不等于中面挠度。
自Reissner 提出考虑横向剪切变形的平板弯曲理论后,又出现了许多精化理论。
但大致分为两类,如Mindlin(明特林)等人的理论和Власов(符拉索夫)等人的理论。
厚板理论是平板弯曲的精确理论,即从3D 弹性力学出发研究弹性曲面的精确表达式。
4. 薄壳理论的基本假定也称为Kirchhoff-Love(克希霍夫-勒夫)假定:①薄壳变形前与中曲面垂直的直线,变形后仍然位于已变形中曲面的垂直线上,且其长度保持不变。
复合材料整体成型大后掠机翼设计与验证研究
复合材料整体成型大后掠机翼设计与验证研究蒋建军1何利军2何建3赵琛41.陆军装备部驻北京地区航空军事代表室 北京 100012;2.海鹰航空通用装备有限公司 北京 100018;3.陆军装备部驻株洲地区航空军事代表室 湖南株洲 412000;4.陆军装备部驻哈尔滨地区航空军事代表室 黑龙江哈尔滨 150060摘要:根据机翼设计及材料工程力学性能要求,基于给定的机翼外形设计并制造了一种满足工程应用要求的全复合材料整体模压成型机翼结构。
通过有限元仿真分析和工程静力学试验方法对该机翼结构进行了全面验证,有限元仿真分析结果与工程静力学试验结果吻合良好,满足复合材料工程力学许用值要求及机翼最大变形不大于半翼展长8%的刚度设计需求,该机翼的力学承载性能得到了充分验证。
关键词:复合材料 大后掠机翼 整体成型 有限元分析 力学试验中图分类号:V279文献标识码:A 文章编号:1672-3791(2024)01-0098-04 Research of the Design and Verification of Highly Swept-BackWings Based on the Integral Molding of Composite MaterialsJIANG Jianjun1HE Lijun2HE Jian3ZHAO Chen41.Aviation Military Representative Office of Army Equipment Department in Beijing Area, Beijing, 100012 China;2. Seahawk General Aviation Equipment Co., Ltd., Beijing, 100018 China;3.Aviation Military RepresentativeOffice of Army Equipment Department in Zhuzhou Area, Zhuzhou, Hunan Province, 412000 China;4.Aviation Military Representative Office of Army Equipment Department in Harbin Area, Harbin,Heilongjiang Province, 150060 ChinaAbstract:According to the requirements of wing design and the engineering mechanical properties of materials, this paper designs and manufactures a kind of wing structure made of the integral molding of all-composite materials which meets the requirements of engineering application based on the shape of the given wing. This paper com‐prehensively verifies the structure of the wing by the methods of finite element simulation analysis and engineering statics tests. The results of finite element simulation analysis are well consistent with the results of engineering statics tests, which satisfies the requirement of the engineering mechanical allowable value of composite materials and the stiffness design requirement that the maximum deformation of the wing is not more than the 8% of the span length of the half wing. The mechanical bearing capacity of the wing has been fully verified.Key Words: Composite materials; Highly swept-back wing; Integral molding; Finite element method; Mechanical test复合材料相较于传统的金属材料具有比强度高、比刚度高、耐腐蚀、可设计性等诸多优点,在航空航天领域中得到了广泛的应用[1-3]。
复合材料薄壁圆管压溃吸能机理分析及层叠壳建模方法研究
复合材料薄壁圆管压溃吸能机理分析及层叠壳建模方法研究冯振宇;周建;张雪晗;马骢瑶;解江;牟浩蕾【摘要】In order to study the energy-absorbing mechanism of T700/3234 composite thin-walled circulartubes,the material properties were measured and the quasi-static axial crush tests of the thin-walled circular tubes were performed.The energy-absorbing indexes,such as,peak load (Fmax),and specific energy absorption (Es),et al of composite thin-walled circular tubes were observed.The effects of induced mechanism and ply orientation of composite material on failure modes and energy-absorbing characteristics of the tubes were analyzed.The test results showed that the initial crushing peak load can be significantly reduced by setting 45-degree outside chamfer at the crush top end of the circular tubes;different fiber ply orientations can lead to different failure modes,and then affect the energy-absorbing characteristics of the tubes accordingly.The finite element modeling method for thin walled-tube laminated shell was developed based on Puck2000 and Yamada Sun failure criteria,the FE model was used to study the energy-absorbing characteristics of T700/3234 composite thin-walled circular tubes.Its accuracy was verified by comparing the simulated results with test ones.The comparison showed that the initial crushing peak load and the specific energy absorption simulated agree well with test ones.%为了研究T700/3234复合材料薄壁圆管的吸能机理,先通过试验测试得到复合材料层板的力学性能,然后对复合材料薄壁圆管进行轴向准静态压缩试验研究.通过测得的峰值载荷(Fmax)、比吸能(Es)等指标参数,分析了引发形式及铺层角度对复合材料薄壁圆管破坏模式和吸能特性的影响.结果表明,在圆管顶端设置45°外倒角引发形式,能够很大程度降低压溃初始峰值.此外,不同纤维铺层角度会导致不同破坏模式,进而影响其吸能特性.基于Puck2000和Yamada Sun失效准则,发展了一种薄壁圆管层叠壳有限元建模方法,用于研究T700/3234复合材料薄壁圆管的吸能特性,对比试验数据,仿真结果的初始峰值载荷、比吸能等指标均与试验结果吻合性较好.【期刊名称】《振动与冲击》【年(卷),期】2017(036)023【总页数】8页(P268-275)【关键词】复合材料薄壁圆管;引发形式;铺层角度;层叠壳;破坏模式;吸能特性【作者】冯振宇;周建;张雪晗;马骢瑶;解江;牟浩蕾【作者单位】中国民航大学天津市民用航空器适航与维修重点实验室,天津300300;中国民航大学天津市民用航空器适航与维修重点实验室,天津300300;中国民航大学天津市民用航空器适航与维修重点实验室,天津300300;中国民航大学天津市民用航空器适航与维修重点实验室,天津300300;中国民航大学天津市民用航空器适航与维修重点实验室,天津300300;中国民航大学天津市民用航空器适航与维修重点实验室,天津300300【正文语种】中文【中图分类】TH212;TH213.3对于普遍意义的运载工具和交通工具来说,无论是车辆还是航空器,都需要对乘客具有一定冲击碰撞安全保护能力,这就要求其结构具有一定的耐撞性[1]。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
------------- ------------- 6.1 复合材料结构分析基本过程 6.1.1 概述 复合材料是由两种或两种以上物理或化学性质不同的材料复合在一起而形成的一种多相固体材料,其主要优点是具有很高的比刚度(刚度与密度之比)。复合材料作为结构材料应用已有很长的历史。目前,复合材料的应用已非常普遍,其应用范围涉及航空、航天、军事、民用等诸多领域。 ANSYS程序提供了一种特殊的单元、层单元来模拟复合材料,利用这些单元就可以进行任意的复合材料结构分析。 复合材料结构分析也包括建模、加载求解及后处理3个基本步骤,其中加载求解及后处理基本同于一般的结构分析过程,建模部分具有其特殊性,下面主要对其建模部分进行详细讨论。 6.1.2 建立复合材料模型 与一般的各向同性材料相比,复合材料的建模过程要相对复杂。由于各层材料性能为任意正交各向异性,材料性能与材料主轴取向有关,所以在定义各层材料性能和方向时要特别注意。在本节中,主要探讨以下4个问题。 1.选择适当的单元类型 用于建立复合材料模型的单元有SHELL99、SHELL91、SHELL181、SOLSH190、SOLID46、SOLID186和SOLID191 7种单元。单元类型的选择主要根据具体的应用和所需计算的结果类型来确定。 (1)SHELL99单元 SHELL99是一种8节点3D壳单元,每个节点有6个自由度。该类单元主要适用于薄到中等厚度的板和壳结构,一般要求结构宽厚比大于10。对于宽厚比小于10的结构,则应考虑选择SOLID46单元建模。SHELL99允许有多达250层的等厚度材料层,或者是125层的厚度在单元面内成双线性变化的不等厚度材料层。如果材料层大于250,用户可通过输入自定义的材料矩阵来建立模型。SHELL99单元可进行失效分析。另外,该类型单元可以将单元节点偏置到结构的表层或底层。 (2)SHELL91单元 SHELL91和SHELL99相类似,只是它允许的复合材料最多有100层,而且用户不能输入自定义的材料性能矩阵。但是,SHELL91支持塑性、大应变行为以及具有一个特殊的“三明治”选项,而SHELL99则无此功能。另外,SHELL91更适用于大变形的情况。 (3)SHELL181单元 SHELL181单元是一种4节点3D壳单元,每个节点有6个自由度。该单元具有包含大应变的完全的非线性性能,最多允许有255层的复合材料结构。各层的信息可通过横截面相关命令输入。 (4)SOLSH190单元 SOLSH190单元是一种4节点3D壳单元,每个节点有3个自由度。该单元具有包含大应变的完全的非线性性能,最多允许有255层的复合材料结构,允许沿厚度方向的变形斜率可以不连续,各层的信息可通过横截面相关命令输入,可以使用FC命令输入失效准则。 (5)SOLID46单元 SOILD46单元是8节点3D实体单元SOLID45的一种叠层形式,它的每个节点有3个自由度。可以用来建立叠层壳或实体的有限元模型,每个单元最多允许有250层的等厚度材料层,同样允许125层的厚度在单元面内成双线性变化的不等厚度材料层。该单元的另一个优点是可以用几个单元叠加的方式来对多于250层的复合材料建立模型,并允许沿厚度方向的变形斜率可以不连续,而且用户也可以输入自定义的本构矩阵。SOLID46单元有一个等效的------------- ------------- 横向刚度,允许在横向上存在非零应力、应变和位移。它可以指定失效准则。与8节点壳单元相比较,SOLID46单元的阶次要低。因此,在壳结构中要得到与SHELL91或SHELL99单元相同的求解结果,需要更大的网格密度。 (6)SOLID186单元 SOILD186单元是20节点3D实体单元,它的每个节点有3个自由度。可以用来建立叠层壳或实体的有限元模型,每个单元最多允许有250层的等厚度材料层,允许沿厚度方向的变形斜率可以不连续,支持材料的非线性行为和大变形。 (7)SOLID191单元 SOILD191单元是20节点3D实体单元SOLID95的一种叠层形式,它的每个节点有3个自由度。可以用来建立叠层壳或实体的有限元模型,每个单元最多允许有100层的等厚度材料层,允许沿厚度方向的变形斜率可以不连续。该单元有一选项允许材料横向上常应力的存在。SOLID191单元不支持材料的非线性行为和大变形。 2.定义材料的叠层结构 复合材料最重要的特征就是叠层结构,每层材料都可能由不同的正交各向异性材料构成,并且其主轴方向也有可能不同。对于叠层复合材料,纤维的方向决定了层的主方向。 下面两种方法可以定义材料层的配置: 1)通过定义各层的材料性质。 2)通过定义表征宏观力、力矩与宏观应变、曲率之间相互关系的本构矩阵(仅适用于SOLID46和SHELL99单元) 这两种方法都通过单元实常数来保存层的信息,每个单元都通过其REAL属性来获得材料层的结构信息。 (1)定义各层材料性质 这种方法由下到上一层一层定义材料层的配置,底层为第一层,后续层沿单元坐标系的正Z轴方向自底向上叠加。如果叠层结构是对称的,SOLID46和SHELL99允许只定义一半的材料层。 有时,某些层可能只延续到模型的一部分,为了建立连续的层,可以把这些中断的层的厚度设置为0。图6.1所示为一个5层模型,其中第2层在某处中断。
通过如下命令定义材料层性质: Command:R、RMORE和RMODIF GUI:Main Menu︱Preprocessor︱Real Constants︱Add/Edit/Delete 所定义的材料层性质包括材料性质、层的定向角和层厚度。 (2)定义本构矩阵 这是定义各层材料性质的另一种方式,适用于SOLID46和SHELL99单元。该矩阵表征了单元力、力矩和应变、曲率的关系,必须在ANSYS外进行定义,可参阅ANSYS理论手册,这种方法的主要优点是: 1)允许用户对聚合物基复合材料的性质进行合并。 ------------- ------------- 2)支持热载荷向量。 3)可表征层数无限制的材料。 矩阵项通过实常数来定义,通过定义单元平均密度可以将质量影响考虑进去。如果使用这种方法,就无法得到每层材料的详细信息。 (3)夹层(“三明治”)结构和多层结构 夹层结构是由两个薄的面板和一个厚的但相对较软的夹心层(至少为总厚度的一半)构成,并假定夹心层承受了所有的横向剪切载荷,而面板则承受了所有的(或几乎所有的)弯曲载荷。图6.2所示为一夹层结构。
夹层结构可以用SHELL63、SHELL91或SHELL181单元来进行建模。SHELL63单元只能有1层,但可以通过实常数选项来模拟夹层结构,即通过修改有效弯曲惯性矩和中面到外层纤维的距离来考虑夹心层的影响。 SHELL91单元可用于夹层结构并允许面板和夹心具有不同的性质,将该单元的KEYOPT(9)设置为1即可激活“夹层”选项。 SHELL181单元采用能量平衡方法描述结构的横向剪切变形,夹层选项在此便失去了意义。 (4)节点偏置 对于SHELL181单元而言,在定义横截面的过程中通过使用SECOFFSET命令设置节点;对于SHELL91和SHELL99单元而言,使用其节点设置选项KEYOPT(11)可将单元的节点设置在壳的底面、中面或顶面上。如果壳是由不同厚度的几部分构成的(台阶状),则节点设置应遵循以下原则: 1)若各部分壳的底面在同一平面上,则将单元节点设置到壳的底面上(参见图6.3a所示)。 2)若各部分壳的中面在同一平面上,则将单元节点设置到壳的中面上(参见图6.3b所示)。 3)若各部分壳的顶面在同一平面上,则将单元节点设置到壳的顶面上(参见图6.3c所示)。 -------------
------------- 3.定义失效准则 失效准则用于获知在所加载荷条件下,各层是否会失效。用户可以从3种已定义的失效准则中进行选择,也可以自定义多达6种的失效准则。以下是3种已定义的失效准则。 1)最大应变失效准则;允许有9个失效应变。 2)最大应力失效准则;允许有9个失效应力。 3)Tsai-Wu失效准则;允许有9个失效应力和3个附加的耦合系数。 定义失效准则的命令如下: (1)Command:FC(FCLIST、FCDELE) GUI:Main Menu︱Preprocessor︱Material Props︱Material Models︱Structural︱Nonlinear︱Inelastic︱Non-Metal Plasticity︱Failure Criteria GUI:Main Menu︱General Postproc︱Failure Criteria 定义失效准则的典型的命令流如下: FC,1,TEMP,,100,200 FC,1,S,XTEN,1500,1200 FC,1,S,YTEN,400,500 FC,1,S,ZTEN,10000,8000 FC,1,S,XY ,200,200 FC,1,S,YZ ,10000,8000 FC,1,S,XZ ,10000,8000 FCLIST,,100 FCLIST,,150 FCLIST,,200 PRNSOL,S,FAIL (2)Command:TB、TBTEMP、TBDATA GUI:Main Menu︱Preprocessor︱Material Props︱Failure Criteria 定义失效准则的典型的命令流如下: TB,FAIL,1,2 TBTEMP,,CRIT TBDATA,2,1 TBTEMP,100 TBDATA,10,1500,,400,,10000