有限元分析第三章

合集下载

有限元的分析过程(第三章)22

有限元的分析过程(第三章)22
0 0 0 L 1 1 0 L
c ( x [ x , x ])
n i i i 0 L
2-1
0 L 1 2 3
其中 ( x [ x , x ]) 为所采用的基底函数,它定义在全域 [ x , x ]上, c , c , c ...
i 0 L
为展开的系数。 第二种是基于子域 [ x , x ] 上的分段展开形式,若采用线性函数,有
上式中的
为节点位移,
为节点力,可以看出,图2-9
分别就单元①、②、③写出了各自的节点力,如对于节点2,即写出了单元①中节 点力 ,又给出了单元②中节点2的节点力 可以看出,在单元组装后,实
3有限元分析的基本流程:
际上只需要合成后的节点力; 因此,今后只需要对各个单元的刚度系数按对应节点位移的位置进行组装,
【典型例题】2.1(1) 一个一维函数的两种展开方式的比较
设有一个一维函数 f ( x), x [ x , x ] ,分析它的展开与逼近形式。
0 L
解答:首先考虑基于全域的展开形式,如采用傅立叶级数展开,则有:
f ( x) c ( x [ x , x ]) c ( x [ x , x ]) ...
4有限元分析的特点:
有限元分析的最大特点就是标准化和规范化,这种特点使
得大规模分析和计算成为可能,当采用了现代化的计算机以及 所编制的软件作为实现平台时,则复杂工程问题的大规模分析 就变为了现实。 实现有限元分析标准化和规范化的载体就是单元,这就需 要我们构建起各种各样的具有代表性的单元,一旦有了这些单 元,就好像建筑施工中有了一些标准的预制构件(如梁、楼板
综合分段函数描述的优势和问题,只要采用功能完善 的软件以及能够进行高速处理的计算机,就可以完全发挥 “化繁为简”策略的优势,有限元分析的概念就在于此。

有限元法基本原理及应用第3章重庆大学龙雪峰

有限元法基本原理及应用第3章重庆大学龙雪峰

有限元原理及应用
第三章 弹性力学有限元法
• 3.单元分析 • 单元分析包括位移模式选择,单元力学分析两个内容。 • 位移模式也称位移函数或插值函数,在有限元位移法中是 以节点位移为基本未知量,再由这些节点位移插值得到单 元内任意一点的位移值。单元的位移模式一般采用多项式, 因为多项式计算简便,并且随着项数的增加,可以逼近任 何一段光滑的函数曲线。 • 单元力学分析 根据所选单元的节点数和单元材料性质, 应用弹性力学几何方程和物理方程得到单元刚度矩阵。由 于连续体离散化后假定力是通过节点在单元间传递的,因 此要利用插值函数把作用在单元上的体积力、面积力和集 中力按静力等效原则移到节点上。
Hale Waihona Puke 有限元原理及应用第三章 弹性力学有限元法
• 5.结果后处理和分析 • 求解线性方程组得到位移矢量后,由几何和物理关系可以 得到应变和应力。 • 由于应变(应力)来自位移的微分可能导致单元间应力不 连续,这会使应力计算误差较大,要在节点附近进行平均 化处理。 • 通过后处理还可得到位移、应变和应力的最大最小值及其 所在位臵以及主应力、主应变或其它定义的等效应力。 • 结果的输出可以应用图表、动画等各种方式。最后还要对 这些结果进行分析以指导工程设计、产品开发等等。
有限元原理及应用第三章弹性力学有限元法?如果挠度与板厚相比不再为小量如金属板当挠度如果挠度与板厚相比不再为小量如金属板当挠度ww与板厚tt的关系在范围内板的中面应变就不能忽略如图的关系在范围内板的中面应变就不能忽略如图35所示面内的两个自由度也要一并考虑所示面内的两个自由度也要一并考虑导致单元的每个节点上a四边形弯曲单元b三角形弯曲单元图34薄板弯曲单元导致单元的每个节点上就要有五个自由度此类单元一般称为薄板单元
有限元原理及应用

结构分析的有限元法-第三章

结构分析的有限元法-第三章

式中
H 1 u B A yH v
(3.32)

H 0 u H 0 v 0 0 0 0 1 0 0 2 0 6x
(3.33)
单元刚度矩阵
再次应用式(2.70),并进行一系列的积分运算,可以得出单元刚度矩阵的显式如下:
l
K
e
E d A B B d x
0 1 l
Av
1
2 l
0 0 1 l 2 1 l
(3.21)
MATLAB不仅可以进行数值运算,也能进行符号运算。如式(3.20)中的矩 阵Au和Av的求逆运算,我们可以在MATLAB的命令窗口下输入 >> syms L >> Au = [ 1 0 1 L ] ; >> Av = [ 1 0 0 0 0 1 0 0 1 L L^2 L^3 0 1 2*L 3*L^2] ; 第一句是定义符号变量L,后面定义两个矩阵Au和Av。然后我们再输入下 面求逆的命令 >> inv(Au) ans = 0 1 1 [ 1, 0] Au [ -1/L, 1/L] 1 l 1 l >> inv(Av) ans = 0 0 1 [ 1, 0, 0, 0] 0 1 0 1 [ 0, 1, 0, 0] A v 2 2 3 l 2 l 3 l [ -3/L^2, -2/L, 3/L^2, -1/L] 3 2 3 1 l 2 l [ 2/L^3, 1/L^2, -2/L^3, 1/L^2] 2 l
根据材料力学的有关知识,我们可以立刻写出杆单元的结点位移与结点力 之间的关系为
FNi EA l (u i u j ) FNj EA l (u j u i )

4-第3章 有限元分析的力学基础

4-第3章  有限元分析的力学基础

第3章有限元分析的力学基础由固体材料组成的具有一定形状的物体在一定约束边界下(外力、温度、位移约束等)将产生变形(deformation),该物体中任意一个位置的材料都处于复杂的受力状态之中,本章将定义用于刻画任意形状弹性变形体的力学变量和表达这些变量之间的关系。

具体地,将在五个简化条件下,定义有关位移、变形、力的三大类变量,推导这些变量之间的三大类方程,给出典型的两类边界条件,本章的主要内容就是弹性力学中的基础部分。

3.1 变形体的描述、变量定义、分量表达3.1.1 变形体在外力的作用下,若物体内任意两点之间发生相对移动,这样的物体叫做变形体(deformed body),它与材料的物理性质密切相关。

如果从几何形状的复杂程度来考虑,变形体又可分为简单形状变形体和任意形状变形体。

简单变形体如杆、梁、柱等,材料力学和结构力学研究的主要对象就是简单变形体,而弹性力学则处理任意形状变形体。

有限元方法所处理的对象为任意形状变形体,因而,弹性力学中有关变量和方程的描述将是有限元方法的重要基础。

3.1.2 基本变量当一个变形体受到外界的作用(如外力或约束的作用)时,如何来描述它?首先,我们可以观察到物体在受力后产生了内部和外部位置的变化,因此,物体各点的位移应该是最直接的变量,它将受到物体的形状、组成物体的材质以及外力的影响,变形体的完整描述如图3.1所示。

图3.1 变形体的描述描述位移是最直接的,因为可以直接观测,描述力和材料特性是间接的,需要我们去定义新的变量,如图3.2所示,可以看出应包括位移、变形程度、受力状态这三个方面的变量,当然,还应有材料参数来描述物体的材料特性。

图3.2 变形体的描述及所需要的变量总之,在材料确定的情况下,基本的力学变量应该有:·位移(displacement) (描述物体变形后的位置)·应变(strain) (描述物体的变形程度)·应力(stress) (描述物体的受力状态)对于任意形状的变形体,我们希望建立的方程具有普遍性和通用性,因此,采用微小体元(representative volume) dxdydz的分析方法来定义位移、应变、应力这三类变量。

第三章MATLAB有限元分析与应用

第三章MATLAB有限元分析与应用

第三章MATLAB有限元分析与应用有限元分析(Finite Element Analysis, FEA)是一种工程计算方法,用于解决结构力学和流体力学等问题。

它将一个复杂的结构分割成多个简单的离散单元,通过建立数学模型和求解方程组,得到结构的力学、热力学和流体力学等性能参数。

MATLAB是一种功能强大的数学计算软件,具有直观的用户界面和丰富的工具箱,可以方便地进行有限元分析。

本章将介绍在MATLAB中进行有限元分析的基本步骤和方法,以及一些常见的应用例子。

首先,进行有限元分析需要将结构进行离散化。

常用的离散化方法有节点法和单元法。

节点法是将结构的几何形状划分为小的节点,并在节点上进行计算。

单元法是将结构划分为多个小的单元,并在每个单元内进行计算。

在MATLAB中,可以通过创建节点和单元的矩阵来描述结构和单元的关系。

例如,创建一个2D结构形式的节点矩阵:nodes = [0 0; 1 0; 0 1; 1 1];然后,通过创建描述节点连接关系的矩阵,来定义结构的单元:elements = [1 2 3; 2 4 3];这里的每一行代表一个单元,数字表示节点的编号。

接下来,需要定义材料的力学参数和边界条件。

材料的力学参数包括弹性模量、泊松比等。

边界条件包括支座约束和加载条件。

在MATLAB中,可以通过定义力学参数和边界条件的向量来描述。

例如,定义弹性模量和泊松比的向量:E=[200e9200e9];%弹性模量nu = [0.3 0.3]; % 泊松比定义支座约束的向量(1表示固定,0表示自由):constraints = [1 1; 0 0; 0 1; 0 1];定义加载条件的向量(包括点力和面力):最后,通过求解方程组得到结构的应力和位移等结果。

在MATLAB中,可以利用有限元分析工具箱中的函数进行计算。

例如,可以使用“assem”函数将节点和单元的信息组装成方程组,并使用“solveq”函数求解方程组。

有限元分析法第3章 杆单元

有限元分析法第3章 杆单元

提示: 1)本例中单元应力的计算采用了材料力学中的方法,与采 用有限元单元应力公式 E EBd 的结果相同。 2)对锥形杆,单元截面积可用平均值。 3)求应力之前需要求出节点位移——有限元位移法。
第三章
杆单元
§ 3 –1
习题2:
一维等截面杆单元
已知:
求:杆两端的支反力

第三章 杆单元
u2
v2 u3 v3
1 1 1 1 EA 1 1 1 1 2 L 1 1 1 1 1 1 1 1
第三章 杆单元
§ 3 –2
二维空间中的杆单元
将单元1,2的刚度方程扩张到系统规模(6阶), 相加后引入节点平衡条件:
第三章
杆单元
§ 3 –2
0 1 0 0 0 1 1 0 1 1 1 0
u1 v1 u2
v2
1 1 1 1 EA 1 1 1 1 2 L 1 1 1 1 1 1 1 1
第三章 杆单元
§ 3 –2
单元2:2-3
135,l
按公式计算杆应力:
二维空间中的杆单元
得:
0 E 2 L 0 1 1 1 1 2 ( P1 P2 ) 1 L 2 EA P 2A 1 P2
P 1 E 2 L P2 1 1 1 1 2 ( P1 P2 ) 2 L 2 EA 0 2 A 0
第三章
杆单元
§ 3 –2
二维空间中的杆单元
节点位移向量的坐标变换:
~ d i Tdi
第三章
杆单元
§ 3 –2
二维空间中的杆单元

第3章 有限元分析的数学求解原理-三大步骤

第3章 有限元分析的数学求解原理-三大步骤

U x x y y z z xy xy yz yz zx zx dV
X u Y v Z w dV X u Y v Z w d W
V V
用 * 表示;引起的虚 应变分量用 * 表示
j Vj
Ui
i Vi


0 X
y
¼ 1-9 Í

ui* * vi wi* * * u j , v* j w*j

x* * y * z * * xy *yz * 18 zx
19
7.间接解法:最小势能原理
20
最小势能原理
W U 0
最小势能原理就是说当一个体系的势能最小时,系统会处于稳定 平衡状态。或者说在所有几何可能位移中,真实位移使得总势能取最小值
0 表明在满足位移边界条件的所有可能位移 最小势能原理: 中,实际发生的位移使弹性体的势能最小。即对于稳定平衡状态,实 际发生的位移使弹性体总势能取极小值。显然,最小势能原理与虚功 原理完全等价。 n m
虚功原理的矩阵表示
在虚位移发生时,外力在虚位移上的虚功是:
* 式中

U i u i* V i v i* W i w i* U j u *j V j v *j W j w *j
* 是 的转置矩阵。
T

*
F
T
同样,在虚位移发生时,在弹性体单位体积内,应力在虚应变上的虚 功是: * * * * * * * T x x y y z z xy xy yz yz zx zx
27
⑴解析法

第三章+有限元法的分析基础

第三章+有限元法的分析基础

第三章有限元方法的分析基础如前所述,有限元分析主要包含以下步骤:(1)求解域的离散;(2)选择场函数;(3)单元刚度矩阵的建立;(4)单刚集成结构总体刚度矩阵;(5整体结构平衡方程求解;(6) 按单元计算场内的各种物理量。

本章将概述有限元分析各阶段的基本要求。

第一节结构的离散化将求解域离散为子域(有限元)是有限元分析的第一步,这也意味着用一个具有有限自由度数目的系统来代替具有无限自由度数目的系统。

离散的实施主要取决于对物理问题的认知和物理模型的建立,在选择单元的形状、尺寸、数量和排列时必须以能正确反映待解决问题为前提,尽可能以更好的精度和更高的计算速率精确地模拟原问题。

3.1.1基本单元形状对任一个给定的物体,必须靠工程实际或研究判断力来选择适当的单元进行离散化。

在大多数情况下,单元类型的选择取决于物体的几何形状以及描述系统所需要的独立的空间坐标数。

图3.1.1、图3.1.2、图3.1.3分别示出了某些常用的一维、二维和三维单元.图3.1.1 一维单元当几何形状、材料性质和其他参数(如应力、位移)仅需用一个空间坐标描述时,可以采用如图3.1.1所示的一维单元。

虽然这种单元有横截面面积,但一般在示意图中都用线段来表示。

在某些问题中,单元的横截面面积可沿长度变化。

当问题的几何形状和其他参数可以用二个独立的空间变量结点来描述时,可以采用图3.1.2所示的二维单元。

二维分析中常用的基本单元是三角形单元,虽然四边形(或其特殊形式矩形或平行四边形)单元可以用二个或四个三角形单元集合而成(如图3.1.3所示),但在某些情况下用四边形(或矩形,平行四边形)单元仍然是有利的。

图3.1.2二维单元149150图3.1.3 由二个或四个三角形单元集合成的四边形单元如果物体的几何形状,材料性质和其他参数可以用三个独立的空间坐标来描述,就可以采用图3.1.4所示的三维单元来离散物体。

与二维问题中的三角形单元类似,基本三维单元是四面体单元,但在某些情况下用六面体单元会更有利。

第三章+matlab有限元分析与应用

第三章+matlab有限元分析与应用
优化设计
在满足一定约束条件下,寻找使某个或多个设 计指标达到最优的设计方案的过程。
目标函数
用于衡量设计方案优劣的数学表达式,通常是 最小化或最大化的某个性能指标。
约束条件
限制设计方案选择的条件,包括设计变量的上下界、设计变量的关系等。
基于Matlab的有限元优化设计方法
MATLAB优化工具箱
提供了一系列用于求解各种优化问题的函数和算法,包括线性规划、非线性规划、混合 整数规划等。
有限元模型
由一组离散化的元素组成,每个 元素代表系统的一部分,并具有 特定的属性和行为。
节点
元素之间的连接点,用于传递力 和位移。
有限元分析的基本步骤
前处理
01
建立有限元模型,包括定义元素类型、几何形状、材料属性、
边界条件和载荷等。
Байду номын сангаас求解
02
应用数学方程求解有限元模型的节点位移和应力分布。
后处理
03
对于一些复杂模型,如具有非线性、大变形、多 材料等特性,建模难度大,需要发展更高级的建 模方法和技术。
数据安全与隐私保护
在进行有限元分析时,需要处理大量的数据,如 何保证数据的安全和隐私保护是一个重要的问题 。需要采取有效的数据加密和保护措施来确保数 据的安全性和隐私性。
未来发展方向与展望
跨学科融合
结果后处理
显示结果
使用Matlab的图形功能,如`plot`、`mesh`等,绘制 结果的可视化图像。
分析结果
对结果进行详细的分析,如查看位移分布、应力分布 等。
结果优化
根据分析结果,对模型进行优化设计,以提高性能或 降低成本。
03
有限元分析实例
Chapter

东南大学 有限元分析课程 第三章 轴对称问题和空间问题有限元法

东南大学 有限元分析课程 第三章 轴对称问题和空间问题有限元法

B = Bi
Bj
Bm

0 0 ( s = i , j , m) cs bs
K e = 2π rc B T DB
单元刚度矩阵的分块矩阵近似表达式为:
K sp = 2π rc BsT DB p bs bp + f s f p + A1 (bs f p + f s bp ) + A2 cs c p = A1 (cs bp + cs f p ) + A2bs c p A1 (bs c p + f s c p ) + A2 cs bp cs c p + A2bs bp
式中: 式中:
A1 =
µ
1− µ
A2 =
1 − 2µ 2(1 − µ )
由于几何矩阵中的元素不是常量,单元刚度矩阵需要通过积分得到, 为简化计算可以用三角形单元形心位置的坐标 rc , z c 代替 B 矩阵中的变 量,将单元中的r和z近似地当作常量,并且分别等于 rc , z c 。
1 r ≈ rc = ( ri + rj + rm ) 3
K e = 2π ∫∫ B T DBrdrdz
单元刚度矩阵的分块矩阵为, 单元刚度矩阵的分块矩阵为,
K sp = 2π ∫∫ BsT DB p rdrdz ( s, p = i, j , m)
由于几何矩阵中的元素不是常量,单元刚度矩阵需要通过积分得到, 为简化计算可以用三角形单元形心位置的坐标 rc , z c 代替 B 矩阵中的变 量 r, z 。
1 rc = ( ri + r j + rm ) 3
1 z c = ( zi + z j + z m ) 3

MATLAB有限元分析与应用可编辑全文

MATLAB有限元分析与应用可编辑全文

%
modulus of elasticity E, cross-sectional
%
area A, and length L. The size of the
%
element stiffness matrix is 2 x 2.
y = [E*A/L -E*A/L ; -E*A/L E*A/L];
2019/11/28
y = k * u/A;
2019/11/28
18
§3-2 线性杆元
3、实例计算分析应用
如图所示二线性杆元结构,假定E=210MPa,A=0.003m^2,P=10kN, 节点3的右位移为0.002m。
求:系统的整体刚度矩阵; 节点2的位移; 节点1、3的支反力; 每个杆件的应力
解:
步骤1:离散化域
%LinearBarElementStresses This function returns the element nodal
%
stress vector given the element stiffness
%
matrix k, the element nodal displacement
%
vector u, and the cross-sectional area A.
? ?
?
? ?
?
10??
630000 ????0.002?? ?? F3 ??
线性杆元也是总体和局部坐标一致的一维有限单元,用线性函数描述
每个线性杆元有两个节点(node)
? EA
单刚矩阵为:k
?
? ?
L
???
EA L
?
EA L
? ? ?

第三章 有限元和刚性体

第三章 有限元和刚性体

第三章有限元和刚性体有限元和刚性体是ABAQUS模型的基本元素。

有限元是可变形的。

刚性体不可变形,只作刚体运动。

用户对有限元多少有一些了解,而对有限元软件中刚性体概念可能感到有些陌生。

ABAQUS软件引入刚性体是为了提高计算效率。

在ABAQUS中任何物体或物体的局部都可以定义为刚性体,大多数单元类型可为刚性体定义所用(参考ABAQUS/Stardard用户手册24-1节所列的表)刚性体与变形体结合的优点是对刚性体运动的完全描述可只用一个参考点,所以不超过六个自由度即可定位一个刚性体。

与此相反,变形单元需要许多自由度,并要求大量的单元计算来确定变形。

当某一局部变形可忽略或对它的变形不感兴趣时,则可在建模时构造一个刚性体部件,这样做就会极大地节省计算时间而并不影响整体结果。

3.1 有限元ABAQUS有各种各样的单元,其庞大的单元库提供了一套强大的工具来解决许多不同类型的问题,本节介绍影响单元特性的五个方面问题。

3.1.1 单元的表征每一个单元都由下面几个特性来表征:·单元族·自由度(和单元族直接相关)·节点数·数学描述,即单元列式·积分ABAQUS中每一种单元都有自己特有的名字,例如T2D2,S4R和C3D81。

单元的名字标志着一种单元的五个特性。

单元命名的规则将在本章里说明。

单元族图3-1给出了应力分析中最常用的单元族。

单元族之间一个明显的区别是每一个单元族所假定的几何类型不同。

图3-1 常用单元族在这本指南里将用到的单元族有实体单元、壳单元、梁单元、桁架和刚性体单元,这些将在以后的各章里详细讨论。

其它的单元族在这本指南中没有讲到;如果对应用它们感兴趣,请查阅ABAQUS/Standard用户手册的Part V。

单元名字里开始的字母标志着这种单元属于哪一个单元族。

例如,S4R中的S表示它是壳单元,C3D81中的C表示它是实体单元。

自由度自由度(dof)是分析中计算的基本变量。

有限元分析基础

有限元分析基础
13
第二章 结构几何构造分析
(a) 结构本身可变 (b) 缺少必要的约束条件 (c) 约束汇交于一点 图2-1 几何可变结构
14
第二章 结构几何构造分析
2.2 结构计算基本知识
2.2.1 结构计算简图
实际结构总是很复杂的,完全按照结构的实际情况 进行力学分析是不可能的,也是不必要的,因此在对实 际结构进行力学计算之前,必须将其作合理的简化,使 之成为既反映实际结构的受力状态与特点,又便于计算 的几何图形。这种被抽象化了的简单的理想图形称之为 结构的计算简图,有时也称为结构的力学模型。 结构计算所常用的结点和支座的简化形式:
对称结构在正对称载荷下,对称轴截面上只能产生 正对称的位移,反对称的位移为零;对称结构在反对称 载荷下,对称轴截面上只有反对称的位移,正对称的位 移为零。 (1) 具有奇数跨的刚架
① 正对称载荷作用
(a) 对称刚架
(b) 变形状态分析
(c) 对称性利用
图2-22对称性利用示意图
19
第二章 结构几何构造分析
单元结点位移条件
当 x0 时
v vi,
v x
i
当 xl
时 v vj,
v x
j
1 vi
2 i
3
3 l2
vi v j
1 l
2i
j
4
2 l3
vi v j
1 l2
i j
34
第三章 杆系结构静力分析的有限单元法
a. 杆件的转折点、汇交点、自由端、集中载荷作用 点、支承点以及沿杆长截面突变处等均可设置成结点。 这些结点都是根据结构本身特点来确定的。
b. 结构中两个结点间的每一个等截面直杆可以设置 为一个单元。 变换为作用在结点上的等效结点载荷。

patran培训教材(有限元分析)

patran培训教材(有限元分析)

目录第一章Patran基础知识 (2)第二章悬臂梁的有限元建模与变形分析 (17)第三章受热载荷作用的薄板的有限元建模与温度场求解 (31)第四章带孔平板的受力分析(平面) (36)第五章厚壁圆筒的受内压作用时的应力分析 (44)第六章受压力载荷作用时板的受力分析 (51)第七章板的模态分析 (57)第八章板的瞬态响应分析 (62)第九章板的频率响应分析 (67)第十章提取车架中性面的模态分析 (72)第一章Patran基础知识一.Patran的用户界面介绍Patran具有良好的用户界面,清晰、简单、易于使用且方便记忆,其用户界面如图1-1所示。

图1-1 patran界面按照各部分的功能,可将Patran界面划分为四个区域:菜单和工具栏区、操作面板区、图形编辑区、信息显示和命令行输入区。

下面,就分别对这几个区域进行介绍。

1.菜单和工具栏区如图1-2所示,patran 的界面上有一行菜单,两行工具栏。

图1-2 菜单工具栏Patran 的菜单是该软件的重要组成部分,使用菜单项,可以完成多设置和操作。

本来,菜单与各种工具是配合使用的,两者是不能独立区分的。

这里对菜单栏进行简单的介绍,一般情况下,Patran 有九个主菜单项,如图1-2所示,文件管理(File )菜单主要用于Patran 数据库文件的打开/关闭,同时也用来从其他CAD 系统输入模型;组(Group )菜单主要用于组的操作,作用类似CAD 系统中的“层”;视窗管理(Viewport )菜单用于视窗设置;视图操作(Viewing )菜单用于图形显示设置,包括了工具栏中一些工具的功能;元素显示管理(Display )菜单用于设置各种元素的显示方式;参数设置(Preferences )菜单用于选择求解器,定制用户自己的环境等操作;工具选项(Tools )菜单中提供了许多非常有用的工具;在线帮助(Help )菜单为使用者提供在线帮助。

工具栏各工具功能见表一:表一 Patran 工具栏各工具功能列表 菜单栏应用菜单按钮工具栏2.操作面板区,图形编辑区和信息显示和命令输入区由工具按钮和菜单选项打开的各种面板一般都显示在patran界面的右侧,即操作面板区。

第三章 杆梁结构的有限元分析原理

第三章 杆梁结构的有限元分析原理

function y = LinearBarElementForces(k,u) %LinearBarElementForces This function returns the element nodal % force vector given the element stiffness % matrix k and the element nodal displacement % vector u. y = k * u;
function y = PlaneTrussElementLength(x1,y1,x2,y2) %PlaneTrussElementLength This function returns the length of the % plane truss element whose first node has % coordinates (x1,y1) and second node has % coordinates (x2,y2). y = sqrt((x2-x1)*(x2-x1) + (y2-y1)*(y2-y1));
function y = PlaneTrussElementStiffness(E,A,L, theta) %PlaneTrussElementStiffness This function returns the element % stiffness matrix for a plane truss % element with modulus of elasticity E, % cross-sectional area A, length L, and % angle theta (in degrees). % The size of the element stiffness % matrix is 4 x 4. x = theta*pi/180; C = cos(x); S = sin(x); y = E*A/L*[C*C C*S -C*C -C*S ; C*S S*S -C*S -S*S ; -C*C -C*S C*C C*S ; -C*S -S*S C*S S*S];

有限元ansys电磁场分析详解第三章-1

有限元ansys电磁场分析详解第三章-1

•电流分布 选择导体
Postproc>elec&mag calc>current
用 ANSYS/Emag进行电磁场分析
3.1-40
• 利单元表数据JT (实数解)看电流等值图 Postproc>plot results>elem table
用 ANSYS/Emag进行电磁场分析
3.1-41
• 检察与外加电流相差90度相位的场量(虚数解) Postproc>by load step>
用 ANSYS/Emag进行电磁场分析
用 ANSYS/Emag进行电磁场分析
3.1-3
二维轴对称有限元模型
导电杆 通量平行条件
绞线圈
电流密度: 1E6 A/m2 频率: 100Hz
通量垂直条件
用 ANSYS/Emag进行电磁场分析
3.1-4
• 二种求解结果: – 实数解: 线圈励磁相位0度 – 虚数解:相位差90度
实数解
虚数解
用 ANSYS/Emag进行电磁场分析
3.1-5
• 利用这两种求解结果,任何时间处的场量都能用迭加的方法来生成
执行动画文件:acaz.avi观察场动画
用 ANSYS/Emag进行电磁场分析
3.1-6
• 根据Faradays 定律,线圈中的时变电流会在导体中感生电流
执行动画文件acjt.avi ,观察电流动画
用 ANSYS/Emag进行电磁场分析
打开材料号显示
用 ANSYS/Emag进行电磁场分析
3.1-31
• 模拟端部条件需要耦合电压(VOLT)自由度 – 选择导体节点
用 ANSYS/Emag进行电磁场分析
3.1-32

汽车驱动桥壳的有限元建模与分析_第三章驱动桥壳有限元模型的建模_36_63

汽车驱动桥壳的有限元建模与分析_第三章驱动桥壳有限元模型的建模_36_63

29第三章 驱动桥壳有限元模型的建模 作为MSC.NASTRAN 的前后处理器,MSC.PATRAN 是工业领域最著名的并行框架式有限元前后处理和分析系统。

在驱动桥壳几何模型的基础上,本章将探讨应用MSC.PATRAN 建立驱动桥壳有限元模型的问题。

3.1 导入驱动桥壳几何模型到MSC.PATRAN中 3.1.1 驱动桥壳几何模型的存储 前一章已经采用CAD 软件Pro/E 建立了所研究驱动桥壳的几何模型,为将几何模型导入到MSC.PATRAN 中,需要将在Pro/E 中建立的几何模型存成一定格式的数据。

STEP 格式是国际标准化组织(ISO )于1984 年提出的关于产品数据的交换标准,全称是“产品模型数据交换标准(Standard for Exchange of Product Model Data )”。

与IGES 数据格式相比,STEP 数据格式模型的数据不易丢失,导入速度较快,因此,将在Pro/E 中建立的几何模型存成STEP 数据格式。

图3-1 New Database对话框 图3-2 New Model Preference菜单 3.1.2 MSC.PATRAN模型数据库文件的建立 (1)启动MSC.PATRAN ,选择“File ”菜单中的“New ”命令,或直接在工具栏上单击按钮 ,出现图3-1所示对话框;30(2)在文件名输入框中输入:CA141_Housing.db ,单击“OK ”按钮确认,即建立新的PATRAN 模型数据库文件,如图3-1所示;(3)建立新的数据库文件后,会出现New Model Preference 菜单,使菜单的内容与图3-2所示一致,单击“OK ”按钮确认。

3.1.3 驱动桥壳几何模型的导入 (1)由MSC.PATRAN 菜单File/Import 打开输入模型对话框,在“Object ”中选择“Model ”,在“Source”中选择“STEP”,即确定模型导入的数据格式是STEP 格式,如图3-3所示;(2)在“File Type ”中选择AP203类型;(3)选择要输入的文件,单击“Apply ”按钮,输入几何模型;(4)MSC.PATRAN 弹出一个模型输入统计报告,导入完成。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。


由此可知,两个主应力就是正应力的最大和最小 值。
过一点任意斜面的应力极值
再求切应力的极值。 将sx=s1,sy=s2 ,txy=0代入切应力公式(2-5), 并利用两个方向余弦的平方和为1,得

由此可知,当 l2=0.5 ,s1≥s2 时,切应力的最大和 最小值如下,其作用平面的法线方向与x轴和y轴成 45°角:
3、研究的方法:有较大的区别。 虽然都从静力学、几何学与物理学三方面进行研究,但 是在建立这三方面条件时,采用了不同的分析方法。
材料力学是对构件的整个截面来建立这些条件的,因而 要常常引用一些截面的变形状况或应力情况的假设。这样虽 然大大简化了数学推演,但是得出的结果往往是近似的,而 不是精确的。 而弹性力学是对构件的无限小单元体来建立这些条件的, 因而无须引用那些假设,分析的方法比较严密,得出的结论 也比较精确。所以,可以用弹性力学的解答来估计材料力学 解答的精确程度,并确定它们的适用范围。
n x y
平面AB上的切应力tn即为
上面所求的全应力P向切线方 向的投影: t n mpx lp y 或
2 2 2 t n px py s n
过一点任意斜面的主应力与主方向 问题3:若经过该点的某一斜面上的切应力为0,求
此斜面上的主应力s和应力主方向a ? 设如图所示的斜面上切应 力为0,则该面上的全应力等 于正应力,也等于主应力, 于是有 p x s nl sl
t xy l (s y s )m 0
2 xy
s (s x s y )s (s xs y t ) 0
展开得平面问题的主应力特征方程:
s I1s I 2 0
2
I1 s x s y I 2 s xs y t
2 xy
由求根公式有:
2 1
I1 I 4 I 2 s x s y s x s y 2 2 s 1,2 ( ) t xy 2 2 2
1、平面应力问题

构件几何特征:很薄的等厚度薄板。 厚度为h远远小于结构另外两个方向的 尺度。薄板的中面为平面。
外力与约束:其所受体力、面力和
约束均平行于中面Oxy面内,并沿厚 度方向Oz不变。而且薄板的两个表面 不受外力作用。因此应力沿厚度方向 不变。


表面面力边界条件:表面不受外力作
s z z t 0 t zx z t 0 t zy z t 0
2、由通过中心C点并平行于z 轴的直线为转轴,列出力矩的 平衡条件,并利用小变形假设, 可推导出“切应力互等定理”, 即 txy=tyx 3、由x轴和y轴两个方向的平 面力系的平衡条件,可推导出 “平衡微分方程”,即 s x t yx
平衡微分方程:注意事项

列平衡条件时,应力和体力应分别乘以其作用 面积和体积,才能得到合力; 应用了两个基本假设:连续性假设(不同面间 应力分量采用泰勒级数展开)和小变形假设(受力 变形前后微分体尺寸不变),这也是其适用的条件。
两个应力主方向是相互垂直 的
过一点任意斜面的应力极值 问题4、已知任一点处两个主应力s1和s2,及其应力
主方向,可求得经过该点正应力、切应力的最大和最 小值。 为了分析简便,选取x轴和y轴分别与两个应力主方 向一致,则该点的应力分量为 sx=s1, sy=s2 , txy=0

先求正应力的极值。 上式代入正应力公式(2-4),并利用两个方向余 弦平方和为1,得 sn=(s1-s2)l2+ s2

平面应力问题和平面应变问题的平衡微分方程 相同(平面应变问题中的正应力sz不影响方程的推 导)
3.3.2 几何方程——应变与位移的关系
3.3.3 物理方程——应力应变关系
1、平面应力问题
平面应力问题条件:
很薄的等厚度薄板,厚 度为h远远小于结构另外两 个方向的尺度。其所受体力、 面力和约束均平行于板面, 即只是Oxy面内的量,并沿 厚度方向不变。薄板的两个 表面不受任何外力和约束的 作用。
一点应力状态分析就是求解上述有关应力
分量,具体为:已知任一点处坐标面上的应 力分量sx,sy和txy,求解如下四个问题:
1:求经过该点、平行于z轴而斜交于x轴和y轴的任何 斜面上的应力p? 2:求经过该点、平行于z轴而斜交于x轴和y轴的任何 斜面上的正应力sn和切应力tn ? 3:若经过该点的某一斜面上的切应力为0,求此斜面 上的主应力s和应力主方向a ?

平衡微分方程中各个量的量纲都相同,其中第 一式的各项为x方向的力,第二项为y方向的力;
平衡微分方程:注意事项

平面问题的平衡微分方程有2个方程,但包含 有3个未知函数,只根据静力学条件无法定解,即 是超静定的。要想定解,还必须考虑几何学和物 理学方面的条件。

平衡微分方程表示了平面区域内任意点的微分 单元体的平衡条件,必然保证任一有限大部分和 整个区域是满足平衡条件的,因而所考虑的静力 学条件是严格和精确的;
3:主应力s和应力主方向a :
s x s y s x s y 2 2 s 1,2 ( ) t xy 2 2 t xy s1 s x tan a1 , tan a 2 t xy s 2 s y
4:经过该点的正应力sn和切应力tn 的最大和最小 值: (s ) s , s
n 极值 1 2
(t n )极值
(s 1 s 2 ) (s 1 s 2 ) , 2 2
3.3 弹性力学的平衡方程、几何方程及物理方程

平面问题的平衡微分方程是考虑平面问题的静力 学条件,根据弹性体内微分单元的静力平衡条件来推 导出应力分量与体力分量之间的关系。
如图,在弹性体内任一点
2 2 2
1、平面应力问题

应力分量分布特点:由于板很薄,外力沿厚度均匀分 布,同时应力沿厚度还是连续分布的,因此应力分量也沿 厚度均匀分布,所以板中各点均有:
sz 0
s x f x, y
t xz 0
t yz 0
t xy hx, y
因此只剩下Oxy面内的三个应力分量,且只是坐标x, y的函 数,沿厚度方向Oz不变,即
s y gx, y
p y s n m sm
又由于有
px s xl t xy m p y t xy l s y m
过一点任意斜面的主应力与主方向
从而有关于方向余弦l,m的线性方程组: 有
2
(s x s )l t xy m 0
t xy m s s x l t xy s s y
平面问题中一点应力状态分析
4:求经过该点的正应力sn和切应力tn 的最大和最小 值?
过一点任意斜面的全应力 问题1:已知任一点处坐标面上的应力分量sx,sy和
txy,求经过该点、平行于z轴而斜交于x轴和y轴的任何
斜面上的应力p?
取如图所示的微分三角板或
三棱柱PAB,当平面AB无限接 近于P点时,该平面上的应力即 为所求。 根据该微分单元的力系平衡 条件,在x和y轴方向上合力为0, 从而有:
过一点任意斜面的主应力与主方向
下面求应力主方向。
将所求主应力s1代入第一个方程:
m1 s 1 s x tan a1 l1 t xy
2
将所求主应力s2代入第二个方 程: t m
tan a 2
l2

xy
s 2 s y
(s x s )l t xy m 0
t xy l (s y s )m 0
(3) 物体是均匀的,也就是说整个物体是由同一种材料组成的。 这样,整个物体的所有各部分才具有相同的物理性质,因而 物体的弹性常数(弹性模量和波桑系数)才不随位置座标而变。
(4) 物体是各向同性的,也就是说物体内每一点各个不同方 向的物理性质和机械性质都是相同的。
(5) 物体的变形是微小的,亦即当物体受力以后,整个物体 所有各点的位移都远小于物体的原有尺寸,因而应变和转角 都远小于1,这样,在考虑物体变形以后的平衡状态时,可以 用变形前的尺寸来代替变形后的尺寸,而不致有显著的误差; 并且,在考虑物体的变形时,应变和转角的平方项或乘积项 都可以略去不计,这就使得弹性力学中的微分方程都成为线 性方程。
取一微小的正平行六面体, 其x、y方向的尺寸分别为dx、 dy,为计算方便,设它在z方 向的尺寸为单位长度1。
平面问题的平衡微分方程
由于六面体是微小的,各面上的应力可
认为是均匀分布,且作用于对应面的中心。
一般而论,应力分量是变量x 和y的函数,作用于左右两对 面或上下两对面的应力分量 不完全相同,具有微小的差 量。
第三章 弹性力学基础知识
弹性力学—区别与联系—材料力学 1、研究的内容:基本上没有什么区别。 弹性力学也是研究弹性体在外力作用下的平衡和运动, 以及由此产生的应力和变形。 2、研究的对象:变形体但也有区别。 材料力学基本上只研究杆、梁、柱、轴等杆状构件, 即长度远大于宽度和厚度的构件。弹性力学虽然也研究杆 状构件,但还研究材料力学无法研究的板与壳及其它实体 结构,即两个尺寸远大于第三个尺寸,或三个尺寸相当的 构件。
Fx 0 p x s xl t xy m Fy 0 p y t xy l s y m
过一点任意斜面的正应力与切应力 问题2:求经过该点、平行于z轴而斜交于x轴和y轴
的任何斜面上的正应力和切应力?
平面AB上的正应力sn即为
上面所求的全应力p向法线方 向n的投影: s lp mp
总之,弹性力学与材料力学既有联系又有区别。它们都 同属于固体力学领域,但弹性力学比材料力学,研究的对象 更普遍,分析的方法更严密,研究的结果更精确,因而应用 的范围更广泛。 但是,弹性力学也有其固有的弱点。由于研究对象的变 形状态较复杂,处理的方法又较严谨,因而解算问题时,往 往需要冗长的数学运算。因此为了简化计算,便于数学处理, 它仍然保留了材料力学中关于材料性质的假定:
相关文档
最新文档