ANSYS-复合材料实例2

ANSYS-复合材料实例2
ANSYS-复合材料实例2

MSC/NASTRAN 113 Exercise Workbook

1a-1

WORKSHOP PROBLEM 1a

Uniaxial Loading of a Laminar

Composite Plate (Part I)

X

Y Z

1”

1”

F

t = 0.

0108”

Objectives:

s Create composite material definition.s Create model.s Specify loads.

s Create a MSC/NASTRAN input file directly or by using

MSC/PATRAN.s Run the analysis using MSC/NASTRAN.s Review deformed shape.

1a-2MSC/NASTRAN 113 Exercise Workbook

WORKSHOP 1a

Uniaxial Loading - Part I

MSC/NASTRAN 113 Exercise Workbook 1a-3

Model Description:

The figure below shows a 2-ply composite plate with uniaxial loading.

Figure 1a-1

The plies are a typical graphite/epoxy tape with the following properties:

Table 1a-1: Material Properties

Elastic Modulus, 1-120 x 106 psi Elastic Modulus, 1-2 2 x 106 psi Poisson Ratio 0.35Shear Modulus 1 x 106 psi Layer Thickness (in)

.0054 in

X

Y Z

1”

1”

F

t = 0.

0108”

Suggested Exercise Steps:

s Create the elements.

s Define the ply material (MAT8) and composite layer

element properties (PCOMP) and apply them to the model.

s Apply uniform loads and constraints.

s Prepare the model for a linear static analysis (SOL 101).

s Submit it for a linear static analysis.

s Review results.

1a-4MSC/NASTRAN 113 Exercise Workbook

WORKSHOP 1a

Uniaxial Loading - Part I

MSC/NASTRAN 113 Exercise Workbook

1a-5

Exercise Procedure:

https://www.360docs.net/doc/6018642581.html,ers who are not utilitizing MSC/PATRAN for generating an input file should go to Step 12, otherwise,proceed to step

2.

2.Create a new database called prob1a.db .

In the New Model Preferences form, set the following: 3. Create geometry for the finite element model.

For the next step, turn on entity labels using the icon:

4.Create distributed load to the surface at Surface (Edge) 1.3.

File/New...

New Database Name:prob1a

OK

Tolerance:x Default

Analysis Code:MSC/NASTRAN Analysis Type:Structural

OK

x Geometry

Action: Create Object: Surface Method:

XYZ Vector Coordinate List:< 1, 1, 0 >Origin Coordinate List:[ 0, 0, 0 ]

Apply

x Loads/BCs

Action:

Create

1a-6

MSC/NASTRAN 113 Exercise Workbook

5.Create boundary conditions for Surface (Edge) 1.1 and Point 1.

Constraining Surface 1.1 translations in x and z directions and rotation in y-direction:

Object: Distributed Load Type: Element Uniform New Set Name: uniform_load Target Element Type: 2D

Input Data...

Edge Distr Load : <0 -15 0>

OK

Select Application Region...Geometry Filter:q Geometry Select Surface Edges: Surface 1.3

Add OK Apply

x Loads/BCs

Action: Create Object: Displacement Type: Nodal

New Set Name: simple_constraint

Input Data...

Translations : <0, ,0>Rotations : < ,0, >

OK

Select Application Region...Select Geometric Entities:

Surface 1.1

WORKSHOP 1a

Uniaxial Loading - Part I

MSC/NASTRAN 113 Exercise Workbook 1a-7

(Remember...“

Curve or Edge ” icon to select the edge.)

Constraining Point 1 in y-direction:

The model should appear as follows:

Add OK Apply

x Loads/BCs

Action: Create Object: Displacement Type: Nodal New Set Name: y_constraint

Input Data...

Translations : < ,0, >Rotations : < >OK

Select Application Region...Select Geometry Entities:Point 1Add OK Apply

1a-8MSC/NASTRAN 113 Exercise Workbook

Figure 1a-2: Geometry with loads and boundary conditions.

6.Define nodes and element connectivities by generating a mesh over

the geometry.

View geometry:

The completed model should appear as follows:

x Finite Elements

Action: Create Object: Mesh Type:

Surface Global Edge Length 0.5Surface List Surface 1

Apply

23

4

1

215.00

15.00

135

135

1X

Y

Z

WORKSHOP 1a

Uniaxial Loading - Part I

MSC/NASTRAN 113 Exercise Workbook 1a-9

Figure 1a-3: The finite element model.

7.Next, define the ply material describing the constituent layers of the

composite using the specified modulus of elasticity, shear modulus,and poisson ratio.

x Materials

Action: Create Object: 2d Orthotropic Method: Manual Input Material Name:ply_mat Input Properties ...

Constitutive Model: Linear Elastic Elastic Modulus 11 =20e6Elastic Modulus 22 =2e6Poisson Ratio 12 =.35Shear Modulus 12 =1e6

Apply Cancel

23412

3

4

56

78

9

1

2

3

4

215.00

15.00

135

135

1

X

Y Z

1a-10MSC/NASTRAN 113 Exercise Workbook

7a.

Ply direction or angle is the direction of the 1 axis of the ply coordinate system. Since E 11 equals 10 times E 22, it is inferred that the fibers are going in the direction of the 1 axis of the ply coordinate system.

7b.

NASTRAN defines Ply 1 as being the most negative ply in the z direction of the element coordinate system (Z e ). The element z axis can be displayed using the normal display as shown below:

Hide Labels:

Figure 1a-4: Normal vectors displayed.

Looking back, Figure 1a-1 shows that Ply 1 has fibers lying in the direction of the global y axis (Y g ) and layer 2 has fibers in the direction of the global x axis (X g ).7c.

NASTRAN defines the zero ply orientation angle as being in the direction of the x axis of the material coordinate system (X m ).

x Finite Elements

Action: Verify Object: Element Test:

Normals

Display Control:q Draw Normal Vectors Test Control: Display Only

Apply

WORKSHOP 1a

Uniaxial Loading - Part I

MSC/NASTRAN 113 Exercise Workbook 1a-11

X m

will be defined by the projection of the x axis of a user defined coordinate system (X u ) onto the surface of the elements.We will define a material coordinate system with an x axis in the direction of Y g as shown below:

(Remember...“Node ” i con t o s elect n ode.)

Figure 1a-5: Material coordinate system created.

By projecting X u onto the elements, a ply with fibers going in the direction of Y g has a ply orientation angle of 0 degrees, such as Ply 1(most negative Z e ), and plies in the direction of X g will have a ply orientation angle of 90 degrees.

https://www.360docs.net/doc/6018642581.html,ing the material ply_mat , create a composite.

x Geometry

Action: Create Object: Coord Method:

3Point Point on Plane 1-3:

Node 4

Apply

x Materials

Action:

Create

1a-12

MSC/NASTRAN 113 Exercise Workbook

To enter material layers for the composite: Click on ply_mat in the

“Existing Materials” databox twice.

To enter layer thickness: (1) Click “Thickness For All Layers of...”databox in the Laminated Composite form, (2) type .0054, and (3) hit Enter .

To enter layer orientations: (1) Under Orientation , click cell on the first row. (2) Click Overwrite Orientations databox, (3) type 0, Enter ,and 90, Enter . Then click on the following: The form should look like the following:

9.Next, create an element property for referencing the composite created in the previous step.

Object: Composite Method: Laminate Material Name:

composite1

Load Text Into Spreadsheet

Material Name Thickness Orientation

1ply_mat .005402ply_mat

.0054

90

Apply

x Properties

Action: Create Dimension: 2 D Type:

Shell Property Set Name:plate Option(s): (pull-down menu) Laminate

Input Properties ...Material Name m:composite1Material Orientation:

Coord 1

WORKSHOP 1a

Uniaxial Loading - Part I

MSC/NASTRAN 113 Exercise Workbook 1a-13

9a.

In order to check the X m direction (orientation angle), do the following:

Figure 1a-6: Material orientation angle shown.

OK

Select Members:Surface 1

Add Apply x Properties

Action: Show

Select Property: Orientation Angle Display Method: Vector Plot Select Groups:

q Current Viewport default_group

Apply

1.000

1.000

1.000

1.000

X Y

Z

X

Y Z

1a-14MSC/NASTRAN 113 Exercise Workbook

10.Submit model for linear static analysis.

Click on the Analysis radio button on the Top Menu Bar and complete the entries as shown here:

An MSC/NASTRAN input file called prob1a.bdf will be generated.This process of translating your model into an input file is called the Forward Translation. The Forward Translation is complete when the Heartbeat turns green.

x Analysis

Action: Analyze Object: Entire Model Method: Analysis Deck Job Name:prob1a

Apply

WORKSHOP 1a Uniaxial Loading - Part I

Generating an input file for MSC/NASTRAN Users: 11.The generated input file (prob1a.bdf) should be similar to the following:

(Type “more prob1a.bdf” at UNIX shell window to compare.)

$ NASTRAN input file created by the MSC MSC/NASTRAN input file

$ translator ( MSC/PATRAN Version 7.5 ) on January 15, 1998 at

$ 19:10:09.

ASSIGN OUTPUT2 = ’prob1a.op2’, UNIT = 12

$ Direct Text Input for File Management Section

$ Linear Static Analysis, Database

SOL 101

TIME 600

$ Direct Text Input for Executive Control

CEND

SEALL = ALL

SUPER = ALL

TITLE = MSC/NASTRAN job created on 15-Jan-98 at 19:09:52

ECHO = NONE

MAXLINES = 999999999

$ Direct Text Input for Global Case Control Data

SUBCASE 1

$ Subcase name : Default

SUBTITLE=Default

SPC = 2

LOAD = 2

DISPLACEMENT(SORT1,REAL)=ALL

SPCFORCES(SORT1,REAL)=ALL

STRESS(SORT1,REAL,VONMISES,BILIN)=ALL

BEGIN BULK

PARAM POST -1

PARAM PATVER 3.

PARAM AUTOSPC YES

PARAM INREL 0

PARAM ALTRED NO

PARAM COUPMASS -1

PARAM K6ROT 0.

PARAM WTMASS 1.

PARAM,NOCOMPS,-1

PARAM PRTMAXIM YES

$ Direct Text Input for Bulk Data

$ Elements and Element Properties for region : plate

$ Composite Property Record created from P3/PATRAN composite material $ record : composite1

$ Composite Material Description :

PCOMP 1 0. 0. + A

+ A 1 .0054 0. YES 1 .0054 90. YES CQUAD4 1 1 1 2 5 4 1

CQUAD4 2 1 2 3 6 5 1

CQUAD4 3 1 4 5 8 7 1

CQUAD4 4 1 5 6 9 8 1

$ Referenced Material Records

MSC/NASTRAN 113 Exercise Workbook1a-15

$ Material Record : ply_mat

$ Description of Material : Date: 15-Jan-98 Time: 17:40:54 MAT8 1 2.+7 2.+6 .35 1.+6

$ Nodes of the Entire Model

GRID 1 0. 0. 0.

GRID 2 .5 0. 0.

GRID 3 1. 0. 0.

GRID 4 0. .5 0.

GRID 5 .5 .5 0.

GRID 6 1. .5 0.

GRID 7 0. 1. 0.

GRID 8 .5 1. 0.

GRID 9 1. 1. 0.

$ Loads for Load Case : Default

SPCADD 2 1 3

LOAD 2 1. 1. 1

$ Displacement Constraints of Load Set : simple_constraint

SPC1 1 135 1 4 7

$ Displacement Constraints of Load Set : y_constraint

SPC1 3 2 1

$ Distributed Loads of Load Set : uniform_load

FORCE 1 3 3.75 1. 0. 0.

FORCE 1 6 3.75 1. 0. 0.

FORCE 1 6 3.75 1. 0. 0.

FORCE 1 9 3.75 1. 0. 0.

$ Referenced Coordinate Frames

CORD2R 1 0. 0. 0. 0. 0. 1. + B

+ B 0. 1. 0.

ENDDATA f64cb716

1a-16MSC/NASTRAN 113 Exercise Workbook

WORKSHOP 1a

Uniaxial Loading - Part I

MSC/NASTRAN 113 Exercise Workbook

1a-17

SUBMITTING THE INPUT FILE FOR MSC/NASTRAN and MSC/PATRAN USERS:

12.Submit the input file to MSC/NASTRAN for analysis by finding an available UNIX shell window. At the command prompt enter nastran prob1a.bdf scr=yes . Monitor the run using the UNIX ps command. 13.Proceed with the Reverse Translation process, that is, importing the prob1a.op2 results file into MSC/PATRAN. To do this, return to the Analysis form and proceed as follows:

Before postprocessing the results, clear the LBC markers from the screen by selecting the following main menu icon:

When the translation is complete and the Heartbeat turns green, bring up the Results form. 14.Create Fringe plot.

Click Select Results icon:

x Analysis

Action: Read Output 2Object: Result Entities Method:

Translate

Select Results File...Selected Results File prob1a.op2OK Apply

x Results

Action: Create Object:

Fringe

Select Result Cases:

Default, Static Subcase

Reset Graphics

1a-18

MSC/NASTRAN 113 Exercise Workbook

Click Target Entities icon: (Note: Target Entity allows you to view fringe plots of entities of your choice.)

Click Display Attributes icon: (Note: Display Attributes form allows you to change the displayed graphics of fringe plots.)

Now click Plot Options icon.

You can press Apply whenever you choose; after working on each icon menu or, as we have done here, after working on all of them.15.Create Deformation plot.

Click Select Results icon:

Select Fringe Result:

Displacements,Translational

Position...((NON-LAYERED))Quantity:

Z Component

Target Entity:

Current Viewport

Style: Discrete/Smooth Display:

Free Edges

Coordinate Transfomation: None Filter Values: None

Apply

x Results

Action: Create Object:

Deformation

Select Result Cases:Default, Static Subcase Select Deformation Result:Displacements,Translational

Position...((NON-LAYERED))Show As:

Resultant

Target Entities Plot Options

WORKSHOP 1a

Uniaxial Loading - Part I

MSC/NASTRAN 113 Exercise Workbook 1a-19

Go through all the icon menus as we have done with the Fringe plot (if you wish), and then:

Figure 1a-7: Fringe and Deformation plots of our model.

Notice that the z displacement of the model dominates the deformation. This is due to the effects of the ply orientation. The first layer’s fibers (the bottom layer) are oriented at 90 degrees to the load, while the second layer’s fibers (top) are oriented at 0degrees. The ply modulus in the primary direction (0 degrees orientation) is ten times that of the secondary direction. Therefore,under tension, the bottom layer will translate more than the top layer and “bow” the material upwards.

Quit MSC/PATRAN when you have completed this exercise.

Apply

1a-20MSC/NASTRAN 113 Exercise Workbook

Ansys在复合材料结构优化设计中的应用

A一13玻璃钢学会第十六届玻璃钢/复合材料学术年会论文集2006年 Amys在复合材料结构优化设计中的应用 覃海艺,邓京兰 (武汉理工大学材料科学与工程学院,武汉430070) 摘要:优化设计方法在复合材料结构设计中起着十分重要的作用。本文详细介绍了Ansys两种优化设计方法.目标函数最优设计和拓扑优化设计的过程,并运用目标函数最优设计方法对复合材料夹层结构进行了最优结构层合设计和运用拓扑优化设计方'法对玻璃钢圆凳进行了最佳形状设计。结果证明Ansys优化设计方法在复合材料结构设计中的有效性。 关键词:Ansys;优化设计方法;目标函数最优设计;拓扑优化设计;复合材料 l前言 复合材料是由两种或多种性质不同的材料组成,具有比强度、比刚度高、耐疲劳性能好及材料与性能可设计强等特点,广泛应用于汽车、建筑、航空、卫生等领域。复合材料通过各相组分性能的互补和关联获得优异的性能,因此复合材料各组分之间及材料整体结构的合理布置,充分发挥复合材料的性能已成为设计的关键所在…。Ansys软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。优化设计是一种寻找确定最优设计方案的技术,Ansys强大的优化设计功能已广泛地应用于复合材料制品的结构设计心J。 2Ansys中的优化设计方法【3娟j 2.1目标函数最优设计 “最优设计”是指满足所有的设计要求,而且所需(如重量、面积、体积、应力、费用等)的方案最小,即目标函数值最小。也就是说,最优设计方案是一个最有效率的方案。在Ansys中设计方案的任何方面都是可以优化的,如尺寸(如厚度)、形状(如过渡圆角的大小)、支撑位置、制造费用、自然频率、材料特性等。实际上,所有可以参数化的Ansys选项都可以作优化设计。目标函数最优设计是通过改变设计变量(自变量)的数值,使状态变量(设计变量的函数,因变量)在满足一定条件时,目标函数(因设计变量的改变而有所改变)的值最小。 目标函数最优设计的一般步骤为①生成循环所用的分析文件,该文件须包括整个分析的过程,并满足以下条件:参数化建立模型(PREIy7),对模型进行初次求解(SOLUTION),对初次求解的结果提取并指定状态变量和目标函数(POSTl/POST26);②在Ansys数据库里建立与分析文件中变量相对应的参数,这一步是标准的做法,但不是必须的(BEGIN或OPT);③进入OPT优化处理器,指定要进行优化设计循环的分析文件(oPT);④声明优化变量:指定哪些参数是设计变量,哪些参数是状态变量,哪个参数是目标函数;⑤选择优化工具或优化算法:优化算法是使单个函数(目标函数)在控制条件下达到最小值的传统算法,包括零阶算法和一阶算法;⑥指定优化循环控制方式,每种优化方法和工具都有相应的循环控制参数,比如最大迭代次数等;⑦进行优化分析;⑧查看设计序列结果(OPT)和后处理(POSTl/POST26)。 2.2拓扑优化设计 拓扑优化是指形状优化,有时也称为外型优化。拓扑优化的目标是寻找承受单载荷或多作者简介:覃海艺(1980?),男,在读硕士。 49

ANSYS复合材料仿真分析及其在航空领域的应用

ANSYS复合材料仿真分析及其在航空领域的应用 复合材料,是由两种或两种以上性质不同的材料组成。主要组分是增强材料和基体材料。复合材料不仅保持了增强材料和基体材料本身的优点,而且通过各相组分性能的互补和关联,获得优异的性能。复合材料具有比强度大、比刚度高、抗疲劳性能好、各向异性、以及材料性能可设计的特点,应用于航空领域中,可以获得显著的减重效益,并改善结构性能。目前,复合材料技术已成为影响飞机发展的关键技术之一,逐渐应用于飞机等结构的主承力构件中,西方先进战斗机上复合材料使用量已达结构总重量的25%以上。飞机结构中,复合材料最常见的结构形式有板壳、实体、夹层、杆梁等结构。板壳结构如机翼蒙皮,实体结构如结构连接件,夹层结构如某些薄翼型和楔型结构,杆梁结构如梁、肋、壁板。此外,采用缠绕工艺制造的筒身结构也可视为层合结构的一种形式。一.复合材料设计分析与有限元方法复合材料层合结构的设计,就是对铺层层数、铺层厚度及铺层角的设计。采用传统的等代设计(等刚度、等强度)、准网络设计等设计方法,复合材料的优异性能难以充分发挥。在复合材料结构分析中,已经广泛采用有限元数值仿真分析,其基本原理在本质上与各向同性材料相同,只是离散方法和本构矩阵不同。复合材 料有限元法中的离散化是双重的,包括了对结构的离散和每一铺层的离散。这样的离散可以使铺层的力学性能、铺层方向、铺层形式直接体现在刚度矩阵中。有限元分析软件,均把增强材料和基体复合在一起,讨论结构的宏观力学行为,因此可以忽略复合材料的多相性导致的微观力学行为,以每一铺层为分析单元。二.ANSYS复合材料仿真技术及其在航空领域应用复合材料具有各向异性、耦合效应、层间剪切等特殊性质,因此复合材料结构的精确仿真,已成为现代航空结构的迫切需求。许多CAE程序都可以进行复合材料的分析,但是大多程序并没有提供完备的功能,使复合材料的精确仿真难以完成。如有些程序不提供非线性分析能力,有些不提供层间剪切应力的求解能力,有些不提供考虑材料失效破坏继续计算能力等等。ANSYS作为一款著名的商业化大型通用有限元软件,广泛应用于航空航天领域,为飞机结构中的复合材料层合结构分析提供了完整精确的解决方案。1.复合材料的有限元模型建立针对飞机结构中的复合材料层合板、梁、实体以及加筋板等结构类型,ANSYS提供一种特殊的复合材料单元———层单元,以模拟各种复合材料,铺层数可达250层以上,并提供一系列技术模拟各种复杂层合结构。复合材料层单元支持非线性、振动特性、热应力、疲劳断裂等各种结构和热的分析功能和算法。2.复合材料的层合结构定义:■铺层结构:ANSYS对于每一铺层可先定义材料性质、铺层角、铺层厚度,然后通过由下到上的顺序逐层叠加组合为复合材料层合结构;也可以通过直接输入材料本构矩阵来定义复合材料性质。■板壳和梁单元截面形状:ANSYS利用截面形状工具可定义矩形、I型、槽型等各种形式;还可以定义各种函数曲线以模拟变厚度截面。3.特殊层合结构的模拟:?变厚度板壳铺层切断:将切断的某铺层厚度定义为零,即可模拟铺层切断前后的板壳实际形状。(图1上)?不同铺层板壳的节点协调:ANSYS板壳层单元的节点均可偏置到任意位置,使不同铺层数板壳的节点在中面或顶面、底面对齐。(图1下)?蜂窝/泡沫夹层结构:ANSYS通过板壳层单元来模拟夹层结构的特性,夹层面板和芯子可以是不同材料。(图2)?板-梁-实体组合结构:ANSYS将实体、板壳与梁等不同类型单元通过MPC技术相联系,各类单元的节点不需要重合并协调,便于飞机等复杂结构模型的处理。4.复合材料有限元模型的检查:复合材料结构模型建立后,可以将板壳和梁单元显示为实际形状,还可以通过图形显示和列表直观地观察铺层厚度、铺层角度和铺层组合形式,方便模型的检查及校对。(图3)5.复合材料层合结构分析ANSYS层单元支持各种静强度刚度、非线性、稳定性、疲劳断裂和振动特性等结构分析。完成分析后,可以图形显示或输出每个铺层及层间的应力和应变等结果(虽然一个单元包含许多铺层),根据这些结果可以判断结构是否失效破坏和满足设计要求。6.复合材料失效准则ANSYS已经预定义了三种复合材料破坏准则来评价复合材料结构安全性,包括最大应变/应力失效准则,蔡-吴(Tsai-Wu)准则。每种强度准则均可定义与温度相关,考虑不同温度下的材料性能。另外,用户也可自定义最多达六种的

ANSYS结构分析指 复合材料

ANSYS结构分析指南第五章复合材料 5.1 复合材料的相关概念 复合材料作为结构应用已有相当长的历史。在现代,复合材料构件已被大量应用于飞行器结构、汽车、体育器材及许多消费产品中。 复合材料由一种以上具有不同结构性质的材料构成,它的主要优点是具有很高的比刚度(刚度与重量之比)。在工程应用中,典型复合材料有纤维和叠层型材料,如玻璃纤维、玻璃环氧树脂、石墨环氧树脂、硼环氧树脂等。 ANSYS程序中提供一种特殊单元--层单元来模拟复合材料。利用这些单元就可以作任意的结构分析了(包括非线性如大挠度和应力刚化等问题)。对于热、磁、电场分析,目前尚未提供层单元。 5.2 建立复合材料模型 与铁或钢等各向同性材料相比,建立复合材料的模型要复杂一些。由于各层材料性能为任意正交各向异性,材料性能与材料主轴取向有关,在定义各层材料的材料性能和方向时要特别注意。本节主要探讨如下问题: 选择合适的单元类型; 定义材料层; 确定失效准则; 应遵循的建模和后处理规则。 5.2.1 选择合适的单元类型 用于建立复合材料模型的单元类型有SHELL99、SHELL91、SHELL181、SOLID46和SOLID191 五种单元。但 ANSYS/Professional 只能使用 SHELL99 和 SHELL46 单元。具体应选择哪一类单元要根据具体应用和所需计算结果类型等来确定。所有的层单元允许失效准则计算。 1、SHELL99--线性层状结构壳单元 SHELL99 是一种八节点三维壳单元,每个节点有六个自由度。该单元主要适用于薄到中等厚度的板和壳结构,一般要求宽厚比应大于10。对于宽厚比小于10的结构,则应考虑选用 SOLID46 来建立模型。SHELL99 允许有多达 250 层的等厚材料层,或者 125 层厚度在单元面内呈现双线性变化的不等材料层。如果材料层大于 250,用户可通过输入自己的材料矩阵形式来建立模型。还可以通过一个选项将单元节点偏置到结构的表层或底层。 2、SHELL91--非线性层状结构壳单元 SHELL91 与 SHELL99 有些类似,只是它允许复合材料最多只有 100 层,而且用户不能输入自己的材料性能矩阵。但是,SHELL91 支持塑性、大应变行为

ANSYS命令流学习笔记圆柱形shell单元的复合材料分析

A N S Y S命令流学习笔记圆柱形s h e l l单元的 复合材料分析 Document serial number【LGGKGB-LGG98YT-LGGT8CB-LGUT-

! ANSYS命令流学习笔记15-圆柱形shell单元的复合材料分析 !学习重点: !1、熟悉单元坐标系下的铺层 当零件形状为规则圆筒时,如何进行铺层建立局部的柱坐标系,将需要铺层单元坐标设置为局部坐标系,进行铺层即可。 譬如圆筒铺层的单元坐标系要建立局部圆柱坐标系。如果还使用笛卡尔坐标系,铺层也能进行,但是铺层方向有很大不同,求解结果也会异常。所以划分网格时,确认单元坐标系选择,划分网格之后,检查单元坐标系情况。确认铺层方向符合预期要求。 本例中要特别注意横向(即Y向)是否符合要求。 !2、熟悉圆面的建模和局部坐标系建立 不解释 !3、熟悉利用MPC施加扭矩 APDL如何对一个圆周施加扭矩在圆心处建立一个节点,然后用MPC单元连接圆心节点和圆周节点,然后在圆心节点上施加一个扭矩即可。 注意将MPC单元的属性改为刚性梁。 注意这里MPC单元的利用,也是自己的一些理解。很多细节也不知道如何在APDL 实现。 !问题描述 ! 传动轴长度为1m,壁厚,直径,铺层共十层,角度为-45/45/-45/45/-45/45/- 45/45/-45/45。一端固定,一端圆周施加扭矩M=2000N·m。 复合材料为横向正交各向异性Ex,Ey,Ez,Vxy,Vyz,Vxz,Gxy,Gyz,Gxz分别为195e9Pa, 35e9Pa, 35e9Pa,, , , 15e9Pa, , 15e9Pa。 应力失效参数:+X:767E6Pa; -X:392E6Pa; +Y:20E6Pa; -Y:70E6Pa; +Z:30E6Pa; -Z:55E6Pa; Sxy: 41E6Pa; Syz: 30E6Pa; Sxz: 41E6Pa。

ANSYS命令流学习笔记15-圆柱形shell单元的复合材料分析

! ANSYS命令流学习笔记15-圆柱形shell单元的复合材料分析 !学习重点: !1、熟悉单元坐标系下的铺层 当零件形状为规则圆筒时,如何进行铺层?建立局部的柱坐标系,将需要铺层单元坐标设置为局部坐标系,进行铺层即可。 譬如圆筒铺层的单元坐标系要建立局部圆柱坐标系。如果还使用笛卡尔坐标系,铺层也能进行,但是铺层方向有很大不同,求解结果也会异常。所以划分网格时,确认单元坐标系选择,划分网格之后,检查单元坐标系情况。确认铺层方向符合预期要求。 本例中要特别注意横向(即Y向)是否符合要求。 !2、熟悉圆面的建模和局部坐标系建立 不解释 !3、熟悉利用MPC施加扭矩 APDL如何对一个圆周施加扭矩?在圆心处建立一个节点,然后用MPC单元连接圆心节点和圆周节点,然后在圆心节点上施加一个扭矩即可。 注意将MPC单元的属性改为刚性梁。 注意这里MPC单元的利用,也是自己的一些理解。很多细节也不知道如何在APDL实现。 !问题描述 ! 传动轴长度为1m,壁厚0.003m,直径0.08m,铺层共十层,角度为-45/45/-45/45/-45/45/-45/45/-45/45。一端固定,一端圆周施加扭矩M=2000N·m。 复合材料为横向正交各向异性Ex,Ey,Ez,Vxy,Vyz,Vxz,Gxy,Gyz,Gxz分别为195e9Pa, 35e9Pa, 35e9Pa,0.28, 0.3, 0.3, 15e9Pa, 3.78e9Pa, 15e9Pa。 应力失效参数:+X:767E6Pa; -X:392E6Pa; +Y:20E6Pa; -Y:70E6Pa; +Z:30E6Pa; -Z:55E6Pa; Sxy: 41E6Pa; Syz: 30E6Pa; Sxz: 41E6Pa。 应变失效参数:+X:0.05; -X:0.045; +Y:0.08; -Y:0.06; +Z:0.04; -Z:0.045; Sxy: 0.035; Syz: 0.042; Sxz:0.025。

ANSYS命令流学习笔记14-shell单元的铺层复合材料分析

! ANSYS命令流学习笔记14-shell单元的铺层复合材料分析!学习重点: !1、熟悉复合材料的材料特点 工程应用中典型的复合材料为纤维增强复合材料。玻璃纤维增强塑料(玻璃钢)、碳纤维、石墨纤维、硼纤维等高强度和高模量纤维。复合材料各层为正交各向异性材料(Orthotropic)或者横向各向异性材料(Transversal Isotropic),材料的性能与材料主轴的取向有关。 各向异性Anisotropic,一般的各项同性材料需要两个材料参数弹性模量E和泊松比v。而各向异性在XYZ有着不同的材料属性,而且拉伸行为和剪切行为互相关联。定义其几何方程需要21个参数。 正交各向异性orthotropic,在XYZ有着不同的材料属性,而且拉伸行为和剪切行为无关,定义材料需要9个参数:Ex,Ey,Ez,Vxy,Vyz,Vxz,Gxy,Gyz,Gxz。 横向各向异性Transversal Isotropic,属于各向异性材料,但是在某个平面上表现出二维上的各向同性。

!2、熟悉复合材料分析所用的ANSYS单元 复合材料单元关键在于能够实现铺层。不同截面属性的梁单元(beam188, beam189, elbow290),2D对称壳单元(shell208, shell209),3D铺层壳单元(shell181, shell281, shell131, shell132),3D铺层实体单元(solid185, solid186, solsh190, solid278, solid279),均能实现复合材料的搭建。其中Beam单元和2D对称壳单元很少使用。SHELL91、SHELL99、SOLID46、SOLID191用于一些以前的分析教程中,但是现在这些单元已经被淘汰,最好选择下列单元区替代他们。用越来越少的单元做越来越多的事情也是趋势。 Shell208和shell209,2D对称壳单元 前者为2节点3自由度单元,后者为3节点3自由度单元,均能用于薄板和中厚板结构(L/h > 5-8)。能够用于复合材料铺层,三明治结构建模。 shell181和shell281, 3D铺层壳单元 前者为4节点6自由度单元,后者为8节点6自由度单元,均能用于薄板和中厚板结构(L/h > 5-8)。能够用于复合材料铺层,三明治结构建模。复合材料计算精度由一阶剪切变形理论决定。shell131, shell132为热分析单元,单元类型分别类似于shell181,shell281。 [注:经典变形理论假设变形后的中位线仍然垂直于中面,且长度不变。一阶变形理论假设变形后的法线仍然为直线且长度不变。三阶阶变形理论假设变形后的法线为三阶曲线。] solid185和solid186, 3D铺层实体单元 前者为8节点3自由度单元,后者为20节点3自由度单元,用于厚板和实体的复合材料分析,均为六面体单元,均可退化为六棱柱单元。Solid278, solid279为热分析单元,单元类型分别类似于solid185,solid186。 Solsh190,3D铺层实体壳单元 8节点3自由度单元,类似实体单元,但是用于薄板和中厚度板的壳结构分析,其结构行为遵循一阶剪切变形理论。 !3、熟悉复合材料的失效准则

ANSYS-复合材料实例2

MSC/NASTRAN 113 Exercise Workbook 1a-1 WORKSHOP PROBLEM 1a Uniaxial Loading of a Laminar Composite Plate (Part I) X Y Z 1” 1” F t = 0. 0108” Objectives: s Create composite material definition.s Create model.s Specify loads. s Create a MSC/NASTRAN input file directly or by using MSC/PATRAN.s Run the analysis using MSC/NASTRAN.s Review deformed shape.

1a-2MSC/NASTRAN 113 Exercise Workbook

WORKSHOP 1a Uniaxial Loading - Part I MSC/NASTRAN 113 Exercise Workbook 1a-3 Model Description: The figure below shows a 2-ply composite plate with uniaxial loading. Figure 1a-1 The plies are a typical graphite/epoxy tape with the following properties: Table 1a-1: Material Properties Elastic Modulus, 1-120 x 106 psi Elastic Modulus, 1-2 2 x 106 psi Poisson Ratio 0.35Shear Modulus 1 x 106 psi Layer Thickness (in) .0054 in X Y Z 1” 1” F t = 0. 0108”

基于ANSYS的大型复合材料风力机叶片结构分析

国 防 科 技 大 学 学 报 第32卷第2期 JOURNA L OF NA TIONA L UNIVERSITY OF DEFE NSE TECHNO LOGY V ol.32N o.22010文章编号:1001-2486(2010)02-0046-05 基于ANSYS的大型复合材料风力机叶片结构分析Ξ 周鹏展1,2,3,肖加余1,曾竟成1,王 进2,杨 军2 (1.国防科技大学航天与材料工程学院,湖南长沙 410073; 2.株洲时代新材料科技股份有限公司,湖南株洲 412007; 3.长沙理工大学能源与动力工程学院,湖南长沙 410076) 摘 要:基于ANSY S软件,对某款应用于G L3A风场的1500kW大型复合材料风力机叶片进行了结构分析。分析结果表明:该叶片的振型以一阶挥舞和一阶摆振为主,其频率分别为0186H z和1159H z;在极限挥舞 载荷作用下,该叶片有限元模型计算得到的叶尖挠度为81445m,而该叶片全尺寸静力试验得到的极限挥舞载 荷作用下的叶尖挠度为8112m,计算值与试验值的误差只有318%;另外,该叶片的最大计算拉应力和压应力 分别为228MPa和201MPa,而该叶片玻纤Π环氧复合材料实测拉伸强度和实测压缩失稳强度分别为720MPa和 380MPa,其计算最大应力只有对应实测极限强度的3117%和5219%。 关键词:复合材料;风力机叶片;结构分析;极限挥舞载荷 中图分类号:TK8 文献标识码:A Structural Analysis of Large2scale Composite Wind Turbine B lade B ased on ANSYS ZH OU Peng2zhan1,2,3,XI AO Jia2yu1,ZE NGJing2cheng1,W ANGJin2,Y ANGJun2 (1.C ollege of Aerospace and M aterial Engineering,National Univ.of Defense T echnology,Changsha410073,China; 2.Zhuzhou T imes New M aterial T echnology C o.Ltd.,Zhuzhou412007,China; 3.C ollege of Energy and P ower Engineering,Changsha Univ.of Science&T echnology,Changsha410076,China) Abstract:Based on the ANSY S s oftware,the structural analysis of a kind of1500kW large2scale com posite wind turbine blade which applied in G L3A wind farm was carried out.The analysis results show that the vibration m odes of this blade are mainly presented as first flapwise m ode and first edgewise m ode,the frequencies of the vibration are respectively0.86H z and1.59H z.At the action of ultimate flapwise loads,the FE M analysis results show that the blade tip deformation is8.445m,while the blade tip deformation of the full scale blade under static test is8.12m,s o the deviation between the calculated and tested value of the blade tip deformation is only 3.8%.M oreover,the calculated maximum tensile stress and the com pressive stress are228MPa and201MPa,while the tested tensile strength and com pressive buckling strength of the glass2fiberΠepoxy com posite are720MPa and380MPa,respectively.C onsequently,the percentages of the calculated maximum stress and the tested ultimate strength are respectively31.7%and52.9%. K ey w ords:com posite;wind turbine blade;structural analysis;ultimate flapwise load 风力机叶片是风力发电机组的关键部件之一,随着世界风力发电机组向大功率方向发展,风力机叶片的长度越来越长,目前世界最长的复合材料风力机叶片是丹麦LM公司生产的,其长度已达6115m,单片重约18t,从而对叶片结构的强度、刚度、重量等的设计提出了更高的要求[1-3]。复合材料具有比强度高、比刚度高、重量轻、可设计性强、承力性能好等特点[4-5],因而在大型风力机叶片中获得了广泛应用。风力机叶片的结构分析作为风力机叶片结构设计的技术基础之一,开始在大功率风力机叶片结构的校核与优化设计中发挥着日益重要的作用。 由于大型复合材料风力机叶片的外形结构和铺层结构都非常复杂,其外形由不同翼型构建而成,属Ξ收稿日期:2009-09-22 基金项目:国家863计划资助项目(2007AA03Z563);中国博士后科学基金资助项目(20070420832);湖南省科技资助项目(2008RS4033) 作者简介:周鹏展(1973—),男,博士后。

Ansys复合材料结构分析操作指导书

Ansys10.0 复合材料结构分析操作指导书

第一章概述 复合材料是两种或两种以上物理或化学性质不同的材料复合在一起而形成的一种多相固体材料,具有很高的比刚度和比强度(刚度和强度与密度的比值),因而应用相当广泛,其应用即涉及航空、航天等高科技领域,也包括游艇、风电叶片等诸多民用领域。由于复合材料结构复杂,材料性质特殊,对其结构进行分析需要借助数值模拟的方法,众多数值模拟软件中Ansys是个不错的选择。 Ansys软件由美国ANSYS公司开发,是目前世界上唯一一款通过ISO9001质量体系认证的分析设计软件,有着近40年的发展历史,经过多次升级和收购其它CAE(Computer Aided Engineering )软件,目前已经发展成集结构力学、流体力学、电磁学、声学和热学分析于一体的大型通用有限元分析软件,是一款不可多得的工程分析软件。Ansys在做复合材料结构分析方面也有不俗的表现,此书将介绍如何使用该款软件进行复合材料结构分析。在开始之前有以下几点需要说明,希望大家能对有限元法有大体的认识,以及Ansys软件有哪些改进,最后给出一些学习Ansys软件的建议。 1、有限元分析方法应用简介 有限元法(Finite Element Method,简称FEM)是建立在严格数学分析理论上的一种数值分析方法。该方法的基本思想是离散化模型,将求解目标离散成有限个单元(Element),并在每个单元上指定有限个节点(Node),单元通过节点相 连构成整个有限元模型,用该模型代替实际结构进行结构分析。在对结构离散后,要求解的基本未知量就转变为各个节点位移(Ansys中称之为DOF(Degree Of Freedom),试想一下,节点的位移包括沿x,y,z轴的平动和转动,也就是节点的自由度),节点位移通过求解一系列代数方程组得到,在求得节点位移后,利用节点位移和应力、应变之间的关系矩阵就可以求出各个节点上的应力、应变,应用线性插值便可以获得单元内任意位置的位移、应力、应变等信息。 2、Ansys软件的发展近况 Ansys软件目前已发展到Ansys V12版本,从V10开始Ansys加入了一个新的工作环境Workbench,原先的Ansys被称为Ansys (classic),虽然操作界面不同,但两者的求解器是一样的。Ansys (classic)的前处理功能相对较弱(主要是建模方面),因而往往需要借助第三方软件,如CAD软件。也许是迫于另一个有限元分析软件ABQUS的竞争压力,Ansys推出了新的Workbench工作环境,Workbench在建模、划分网格、求解和后处理上都作了改进,尤其在建模和划分

ANSYS在玻璃钢-混凝土复合材料结构中的计算分析

ANSYS在玻璃钢-混凝土复合材料结构中的计算分析 殷波 (扬州大学水利与建筑工程学院土木工程系,扬州225009) 摘要:混凝土结构由于受荷载变化、材料、施工质量等因素影响,会造成结构的强度、刚度不足,玻璃钢-混凝土复合材料结构则改善其性能。本文通过ansys有限元软件计算,分析说明了玻璃钢-混凝土复合材料结构将有力的提高结构的强度和刚度。 关键词:ansys,混凝土、玻璃钢-混凝土、有限元、复合材料 ANSYS’s calculation in glass fibre reinforced plastics- concrete composite material structure YIN BO (Dept.of Civil Engin,Hydr and Civil Engin Coll, Yangzhou University, Yangzhou, 225009, China) Abstract: As the variety of load、material、constructional quality and so on ,concrete structure may be insufficient in intensity and rigidity. Glass fibre reinforced plastics- concrete composite material structure can improve its capacity. With the calculation of ansys, this paper indicate that glass fibre reinforced plastics-concrete composite material structure will raise the intensity and rigidity. Key words: ansys ,concrete, glass fibre reinforced plastics- concrete, finite element; composite material 1.前言 钢筋混凝土结构合理地利用钢筋和混凝土两种材料的力学性能,因而具有整体性、耐久性等优点。但同时也存在着自重大等缺点,在受弯构件中较显著。自重太大,对于设计大跨度结构、高层建筑以及抗震结构都很不利。针对普通钢筋混凝土梁以及经由玻璃钢加固的钢筋混凝土复合梁,通过ansys有限元软件计算分析说明了玻璃钢-混凝土复合材料结构将提高其强度和刚度。 2.用ANSYS进行有限元计算 有限元法是目前工程技术领域中实用性最强、应用最为广泛的数值计算方法。它的基本思想是将问题的求解域划分为一系列单元,单元之间靠节点连接。单元内部点的待求物理量可由单元节点物理量通过选定的函数关系插值求得。由于单元形状简单,易于由平衡关系或能量关系建立节点量之间的方程式,然后由单元方程再形成总体代数方程组,加入边界条件后即可对方程组求解。以ANSYS为代表的有限元分析软件具有以下优点:减少设计成本、缩短设计和分析的循环周期、增加产品和工程的可靠性、采用优化设计、降低材料的消耗和成本、在产品制造或工程施工前预先发现潜在的问题、可以进行模拟实验分析、进行机械事故分析,查找事故原因。ANSYS软件是集结构、热、流体、电磁、声学于一体的大型通用有限元分析软件,可广泛地应用于土木工程、

最新ansys_复合材料分析汇总

a n s y s_复合材料分析

第五章复合材料 5.1 复合材料的相关概念 复合材料作为结构应用已有相当长的历史。在现代,复合材料构件已被大量应用于飞行器结构、汽车、体育器材及许多消费产品中。 复合材料由一种以上具有不同结构性质的材料构成,它的主要优点是具有很高的比刚度(刚度与重量之比)。在工程应用中,典型复合材料有纤维和叠层型材料,如玻璃纤维、玻璃环氧树脂、石墨环氧树脂、硼环氧树脂等。 ANSYS程序中提供一种特殊单元--层单元来模拟复合材料。利用这些单元就可以作任意的结构分析了(包括非线性如大挠度和应力刚化等问题)。对于热、磁、电场分析,目前尚未提供层单元。 5.2 建立复合材料模型 与铁或钢等各向同性材料相比,建立复合材料的模型要复杂一些。由于各层材料性能为任意正交各向异性,材料性能与材料主轴取向有关,在定义各层材料的材料性能和方向时要特别注意。本节主要探讨如下问题: 选择合适的单元类型; 定义材料层; 确定失效准则; 应遵循的建模和后处理规则。 5.2.1 选择合适的单元类型 用于建立复合材料模型的单元类型有SHELL99、SHELL91、 SHELL181、SOLID46和SOLID191 五种单元。但 ANSYS/Professional 只能使用SHELL99 和 SHELL46 单元。具体应选择哪一类单元要根据具体应用和所需计算结果类型等来确定。所有的层单元允许失效准则计算。 1、SHELL99--线性层状结构壳单元 SHELL99 是一种八节点三维壳单元,每个节点有六个自由度。该单元主要适用于薄到中等厚度的板和壳结构,一般要求宽厚比应大于10。对于宽厚比小于10的结构,则应考虑选用 SOLID46 来建立模型。SHELL99 允许有多达 250 层的等厚材料层,或者 125 层厚度在单元面内呈现双线性变化的不等材料层。如果材料层大于 250,用户可通过输入自己的材料矩阵形式来建立模型。还可以通过一个选项将单元节点偏置到结构的表层或底层。 2、SHELL91--非线性层状结构壳单元 SHELL91 与 SHELL99 有些类似,只是它允许复合材料最多只有100 层,而且用户不能输入自己的材料性能矩阵。但是,SHELL91 支持塑性、

第二章 有限元法及其软件ANSYS在压电复合材料分析中的应用

第二章有限元法及其软件ANSYS在压电复合材料分析中的应用2.1 有限元法概述 在工程技术领域中有许多力学问题和场问题,例如固体力学中的应力应变场和位移场分析以及电磁学中的电磁场分析、振动模态分析等,都可以看作是在一定的边界条件下求解其基本微分方程的问题。虽然人们已经建立了它们的基本方程和边界条件,但只有少数简单的问题才能求出其解析解。。对于那些数学方程比较复杂,物力边界形状又不规则,承受任意载荷的问题,采用解析法求解在数学上往往会遇到难于克服的困难。通常对于这类问题,往往采用各种行之有效的数值计算方法来获得工程需要的数值解,比如有限元法。 有限元法的基本思想最早出现在20世纪40年代初 期,但是直到1960年,美国的Clough.R.W在一篇论文 中首次使用“有限元法”这个名词。在20世纪60年代 末70年代初,有限元法在理论上已基本成熟,并开始 出现商业化的有限元分析软件,如ANSYS等。 有限元发的基本思想是将连续的结构离散成有限 个单元,并在每一个单元中设定有限个节点,将连续体 看作是只在节点处相连接的一组单元的集合体;同时选 定场函数的节点值作为基本未知量,并在每一个单元中 假设一近似插值函数以表示单元中场函数的分布规律; 进而利用力学中某些变分原理去建立节点未知量的有 限元法方程,从而将一个连续域中的无限自由度问 题化为离散域中的有限自由度问题。一经求解就可以利 用解得的节点值和设定的插值函数确定单元上以致整 个集合体上的场函数。有限元求解程序的内部过程如图 2-1所示。 2.2 压电复合材料的有限元分析方法 传统上我们一般采用常用的细观力学方法分析压电复合材料的有效性能,主要目的是建立材料的宏观有效性能,包括弹性、压电和介电性能,与细观结构的定量关系,以指导材料的设计和制造。但是对于一般的细观力学方法来说,例如Dilute模型,自洽方法,Mori-Tanaka 模型和微分方法等,其建立起来的力学模型中涉及了大量复杂的积分和微分公式,用普通的解析法一般无法准确的求出正确解。例如Dunn和Taya使用常用的细观力学模型自洽方法,Mori-Tanaka模型和微分方法对压电复合材料的压电系数进行了预报,并与实验数据进行比较。结果表明:除了体积分数较小时,这些方法给出了比较接近的数值结果;但是当体积分数较大时,其数值结果就跟实验结果有很大的差别。 而压电复合材料作为两种和两种以上组成的宏观非均匀材料可以用合适的具有某种周期分布的微结构材料来表示,这样针对某一周期的非均匀材料单元,利用通常的有限元及边界元方法,可以数值上求得纤维、基体及界面处的应力分布,在此基础上可以预报复合材料的有效性能。这弥补了应用常规细观力学方法无法预报纤维或高体积分数及具有复杂微分结构材料等情况的不足。 2.2.1 基本耦合公式 压电线性理论的基本耦合公式(Tiersten,1969)如下:

ANSYS复合材料仿真分析

在ANSYS 中可以定义多种材料属性: 主菜单-> preprocesser -> Material Prop -> Material Models -> 打开Define Material Model Behavior 对话框 -> 顶部菜单中:Material -> New Model ... -> 弹出Define Material ID 对话框-> 定义更多的材料 ANSYS复合材料仿真分析 2009-05-23 23:31 复合材料,是由两种或两种以上性质不同的材料组成。主要组分是增强材料和基体材料。复合材料不仅保持了增强材料和基体材料本身的优点,而且通过各相组分性能的互补和关联,获得优异的性能。复合材料具有比强度大、比刚度高、抗疲劳性能好、各向异性、以及材料性能可设计的特点,应用于航空领域中,可以获得显著的减重效益,并改善结构性能。 目前,复合材料技术已成为影响飞机发展的关键技术之一,逐渐应用于飞机等结构的主承力构件中,西方先进战斗机上复合材料使用量已达结构总重量的25%以上。飞机结构中,复合材料最常见的结构形式有板壳、实体、夹层、杆梁等结构。板壳结构如机翼蒙皮,实体结构如结构连接件,夹层结构如某些薄翼型和楔型结构,杆梁结构如梁、肋、壁板。此外,采用缠绕工艺制造的筒身结构也可视为层合结构的一种形式。 一.复合材料设计分析与有限元方法 复合材料层合结构的设计,就是对铺层层数、铺层厚度及铺层角的设计。采用传统的等代设计(等刚度、等强度)、准网络设计等设计方法,复

合材料的优异性能难以充分发挥。在复合材料结构分析中,已经广泛采用有限元数值仿真分析,其基本原理在本质上与各向同性材料相同,只是离散方法和本构矩阵不同。复合材料有限元法中的离散化是双重的,包括了对结构的离散和每一铺层的离散。这样的离散可以使铺层的力学性能、铺层方向、铺层形式直接体现在刚度矩阵中。有限元分析软件,均把增强材料和基体复合在一起,讨论结构的宏观力学行为,因此可以忽略复合材料的多相性导致的微观力学行为,以每一铺层为分析单元。 二.ANSYS复合材料仿真技术及其在航空领域应用 复合材料具有各向异性、耦合效应、层间剪切等特殊性质,因此复合材料结构的精确仿真,已成为现代航空结构的迫切需求。 许多CAE程序都可以进行复合材料的分析,但是大多程序并没有提供完备的功能,使复合材料的精确仿真难以完成。如有些程序不提供非线性分析能力,有些不提供层间剪切应力的求解能力,有些不提供考虑材料失效破坏继续计算能力等等。ANSYS作为一款著名的商业化大型通用有限元软件,广泛应用于航空航天领域,为飞机结构中的复合材料层合结构分析提供了完整精确的解决方案。 1.复合材料的有限元模型建立针对飞机结构中的复合材料层合板、梁、实体以及加筋板等结构类型,ANSYS提供一种特殊的复合材料单元

Ansys复合材料结构分析操作指导书 - 副本

第四章复合材料计算实例 在有了前几章知识做铺垫,这一章我们来学习两个复合材料分析的例子,加深复合材料分析的理解,也希望读者能从中收获一些经验。在这里将第二章的流程图再次拿出来,进一步熟悉ANSYS有限元分析的基本过程。 图7 Ansys 结构分析流程图 4.1 层合板受压分析 4.1.1 问题描述 层合板指的是仅仅由FRP层叠而成的复合板材,中间不包含芯材,板材的性能不仅与纤维的弹性模量、剪切模量有关,还与纤维的铺层方向有着密切关系。本例中的板材有4层厚度为0.025m的单元板复合而成,单元板的铺层方向为0°、90°、90°、0°,见图13所示。单元板的材料属性见表4.1。 表 4.1 单元板材料属性 图13 复合材料板

4.1.2 求解步骤 根据问题描述,所要分析的问题为壳体结构的复合材料板,可以采用SOLID46单元建立3D有限元模型进行分析。结合图7的一般步骤进行分析。 步骤一:选取单元类型,设置单元实常数 ⑴、在开始一个新分析前,需要指定文件保存路径和文件名。 文件保存路径GUI:【Utility Menu】|【File】|【Change Directory】见图14 指定新的文件名GUI:【Utility Menu】|【File】|【Change Jobname】见图15所示 图14 指定文件保存路径 图15 修改文件名 ⑵、选取单元类型

1)选取单元类型的GUI操作:【Main Menu】|【Preprocessor】|【Element Type】|【Add/Edit/Delete】,执行后弹出Element Types对话框。2)在Element Types对话框点击Add定义新的单元类型,弹出“Library of Element Types”对话框,见图16所示,按图中所示选择,单元 类型参考号输入框中输入数字1。 图16 单元类型对话框 3)点击“OK”,回到“Element Types”对话框见图17所示,从图中可以看到,定义的单元类型参考号为1,单元类型对应为SOLID46。 图17 已经定义好的单元类型 4)接下来,还要对单元类型做一些选项设置,点击“Options”,弹出“SOLID46 element type options”对话框,在“Form of input”下拉 选择列表中有三个选项,分别是各材料层厚度相同、变厚度材料层、 自定义宏观材料本构矩阵,选择不同的选项,会导致后面需要输入 的材料参数不同。由于本例各层厚度相同,选择“Const thk layer” 项,点击“OK”,弹出“More SOLID46 element type options”,在 K8选项中选择“All layers”,然后单击OK,随后单击ElementTypes 对话框上的Close,关闭该对话框,完成单元类型选择,见图18。

ansys 复合材料分析

第五章复合材料 5.1 复合材料的相关概念 复合材料作为结构应用已有相当长的历史。在现代,复合材料构件已被大量应用于飞行器结构、汽车、体育器材及许多消费产品中。 复合材料由一种以上具有不同结构性质的材料构成,它的主要优点是具有很高的比刚度(刚度与重量之比)。在工程应用中,典型复合材料有纤维和叠层型材料,如玻璃纤维、玻璃环氧树脂、石墨环氧树脂、硼环氧树脂等。 ANSYS程序中提供一种特殊单元--层单元来模拟复合材料。利用这些单元就可以作任意的结构分析了(包括非线性如大挠度和应力刚化等问题)。对于热、磁、电场分析,目前尚未提供层单元。 5.2 建立复合材料模型 与铁或钢等各向同性材料相比,建立复合材料的模型要复杂一些。由于各层材料性能为任意正交各向异性,材料性能与材料主轴取向有关,在定义各层材料的材料性能和方向时要特别注意。本节主要探讨如下问题: 选择合适的单元类型; 定义材料层; 确定失效准则; 应遵循的建模和后处理规则。 5.2.1 选择合适的单元类型 用于建立复合材料模型的单元类型有SHELL99、SHELL91、SHELL181、SOLID46和SOLID191 五种单元。但 ANSYS/Professional 只能使用 SHELL99 和SHELL46 单元。具体应选择哪一类单元要根据具体应用和所需计算结果类型等来确定。所有的层单元允许失效准则计算。 1、SHELL99--线性层状结构壳单元 SHELL99 是一种八节点三维壳单元,每个节点有六个自由度。该单元主要适用于薄到中等厚度的板和壳结构,一般要求宽厚比应大于10。对于宽厚比小于10的结构,则应考虑选用 SOLID46 来建立模型。SHELL99 允许有多达 250 层的等厚材料层,或者 125 层厚度在单元面内呈现双线性变化的不等材料层。如果材料层大于 250,用户可通过输入自己的材料矩阵形式来建立模型。还可以通过一个选项将单元节点偏置到结构的表层或底层。 2、SHELL91--非线性层状结构壳单元 SHELL91 与 SHELL99 有些类似,只是它允许复合材料最多只有 100 层,而且用户不能输入自己的材料性能矩阵。但是,SHELL91 支持塑性、大应变行为

相关文档
最新文档