液体晃动有限元模态分析及试验研究_李松

液体晃动有限元模态分析及试验研究_李松
液体晃动有限元模态分析及试验研究_李松

有限元填空选择题及答案

1有限元是近似求解_一般连续_场问题的数值方法 2有限元法将连续的求解域离散为若干个子域_,得到有限个单元,单元和单元之间用节点相连 3从选择未知量的角度来看,有限元法分为三类位移法. 力法混合法 4以_节点位移_为基本未知量的求解方法称为位移法. 5以_节点力_为基本未知量的求解方法称为力法. 6一部分以__节点位移__,另一部分以_节点力_为基本未知量的求解方法称为混合法. 7直梁在外力的作用下,横截面的内力有剪力_和_弯矩_两个. 8平面刚架结构在外力的作用下,横截面上的内力有轴力_ 、剪力_和弯矩. 9进行直梁有限元分析,平面刚架单元上每个节点的节点位移为挠度和转角 10平面刚架结构中,已知单元e的坐标变换矩阵[T e ]和在局 部坐标系x’O’y’下的单元刚度矩阵[K’]e ,则单元在真体坐标 系xOy下的单元刚度矩阵为_ [K]e = [T e ] T [K’] e [T e ] 13弹性力学问题的方程个数有15个,未知量的个数有15个. 14弹性力学平面问题的方程个数有8_个,未知量个数有8_个15几何方程是研究__应变___和_位移之间关系的方程 16物理方程是描述_应力_和_应变_关系的方程 17平衡方程反映了_应力__和_位移_之间关系的 18把经过物体内任意一点各个_ 截面上的应力状况叫做__该点_的应力状态 19形函数在单元上节点上的值,具有本点为_1_.它点为零的性质,并且在三角形单元的任一节点上,三个行函数之和为_1_ 20 形函数是_三角形_单元内部坐标的_线性位移_函数,他反映了单元的_位移_状态 21在进行节点编号时,要尽量使用同一单元的相邻节点的狭长的带状尽可能小,以使最大限度地缩小刚度矩阵的带宽,节省存储,提高计算效率. 22三角形单元的位移模式为_线性位移模式_- 23矩形单元的位移模式为__线性位移模式_ 24在选择多项式位移模式的阶次时,要求_所选的位移模式应该与局部坐标系的方位无关的性质为几何_各向同性 25单元刚度矩阵描述了_节点力_和_节点位移之间的关系 26在选择多项式作为单元的位移模式时,多项式阶次的确定,要考虑解答的收敛性,即要满足单元的_完备性和协调性要求27三节点三角形单元内的应力和应变是_常数,四节点矩形单元内的应力和应变是线性_变化的 28在矩形单元的边界上,位移是线性_变化的 29整体刚度是一个呈_ 狭长的带状_分布的稀疏矩阵 30整体刚度[K]是一个奇异阵,在排除刚体位移_后,它正义阵1从选择未知量的角度来看,有限元法可分为三类(力法,位移法,混合法) 2下列哪有限元特点的描述中,哪种说法是错误的(D需要使用于整个结构的插值函数) 3几何方程研究的是(A应变和位移)之间关系的方程式 4物理方程是描述(D应力和应变)关系的方程 5平衡方程研究的是(C应力和位移)之间关系的方程式 6在划分单元时,下列哪种说话是错误的(A一般首选矩形单元) 7下列哪种单元的单元刚度矩阵必须通过积分才能得到(D矩形单元) 8单元的刚度矩阵不取决于下列哪种因素(B单元位置) 9可以证明,在给定载荷的作用下,有限元计算模型的变形与实际结构变形之间的关系为(B前者小于后者) 10ANSYS按功能作用可分为若干个处理器,其中(B求解器)用于施加载荷和边界条件 11下列有关有限元分析法的描述中,哪种说话是错误的(B单元之间通过其边界连接成组合体) 12下列关于等参数单元的描述中,哪些说话是错误的(C将规则单元变换为不规则单元后,易于构造位移模式) 13从选择未知量的角度来看,有限元可以分为三类,混合法的未知量是(C节点力和节点位移) 14下列对有限元特点的描述中,哪种说话是错误的(B对有限元求解域问题没有较好的处理方法) 15在划分单元时,下列哪种说话错误(D自由端不能取为节点) 16对于平面问题,选择单元一般首选(D三角形单元或等参单元) 17下列哪种说法不是形函数的性质(C三角形单元任一条边上的形函数,与三角形的三个节点坐标都有关) 18下列四种假设中,哪种分析不属于分析弹性力学的基本假设(C大变形假设) 19下列四种假设中,哪种不属于分析弹性力学的基本假设(B 有限变形假设) 20下列关于三角形单元说法中哪种是错误的(C在单元的公共边上应力和应变的值是连续的) 21下列关于矩形单元的说法哪项是错误的(D其形函数是线形的) 22应用圣维南原理简化边界条件时,静力等效是指前后的力系的(D主矢量相同,对于同一点的主矩也相同) 24描述同一点的应力状态需要的应力分量是(C6个) 25在选择多项式作为单元的位移模式时.多项式阶次的确定,要考虑解答的收敛性,哪种说法不是单元必须满足的要求(D 对称性)

有限元模态分析报告实例

ANSYS模态分析实例 5.2ANSYS建模 该课题研究的弹性联轴器造型如下图5.2: 在ANSYS中建立模型,先通过建立如5.2所式二分之一的剖面图,通过绕中轴线旋转建立模拟模型如下图5.3

5.3单元选择和网格划分 由于模型是三给实体模型,故考虑选择三维单元,模型中没有圆弧结构,用六面体单元划分网格不会产生不规则或者畸变的单元,使分析不能进行下去,所以采用六面体单元。经比较分析,决定采用六面体八结点单元SOLID185,用自由划分的方式划分模型实体。课题主要研究对象是联轴器中橡胶元件,在自由划分的时候,中间件2网格选择最小的网格,smart size设置为1,两端铁圈的smart size设置为6,网格划分后模型如图5.4。 5.4边界约束 建立柱坐标系R-θ-Z,如5-5所示,R为径间,Z为轴向

选择联轴器两个铁圈的端面,对其面上的节点进行坐标变换,变换到如图5.5所示的柱坐标系,约束节点R,Z方向的自由度,即节点只能绕Z轴线转 5.5联轴器模态分析 模态分析用于确定设计中的结构或者机器部件振动特性(固有频率和振型),也是瞬态变动力学分析和谐响应分析和谱分析的起点。 在模态分析中要注意:ANSYS模态分析是线性分析,任何非线性因素都会被忽略。因此在设置中间件2的材料属性时,选用elastic材料。 5.5.1联轴器材料的设置 材料参数设置如下表5-1: 表5.1材料参数设置 表5.1材料参数设置 铁圈1 中间件2 铁圈3 泊松比0.3 0.4997 0.3 弹性模量Mpa 2E5 1.274E3 2E5 密度kg/m 7900 1000 7900 5.5.2联轴器振动特性的有限元计算结果及说明 求解方法选择Damped方法,频率计算结果如表5-2,振型结果为图5.6: 表5.2固有频率 SET TEME/FREQ LOAO STEP SUBSTEP CUMULATIVE 1 40.199 1 1 1 1 73.63 2 1 2 2 3 132.42 1 3 3 4 197.34 1 4 4

有限元分析及应用大课后复习

有限元分析及应用作业报告

目录 有限元分析及应用作业报告....................................... I 目录 ........................................................ II 试题1 . (1) 一、问题描述 (1) 二、几何建模与分析 (2) 三、第1问的有限元建模及计算结果 (2) 四、第2问的有限元建模及计算结果 (7) 五、第3问的有限元建模及计算结果 (13) 六、总结和建议 (16) 试题5 (17) 一、问题的描述 (17) 二、几何建模与分析 (18) 三、有限元建模及计算结果分析 (18) 四、总结和建议 (26) 试题6 (27) 一、问题的描述 (27) 二、几何建模与分析 (27) 三、有限元建模及计算结果分析 (27) 五、总结和建议 (35)

试题1 一、问题描述 图示无限长刚性地基上的三角形大坝,受齐顶的水压力作用,试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析,并对以下几种计算方案进行比较: 1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算; 2)分别采用不同数量的三节点常应变单元计算; 3)当选常应变三角单元时,分别采用不同划分方案计算。 图1-1模型示意图及划分方案

二、几何建模与分析 图1-2力学模型 由于大坝长度>>横截面尺寸,且横截面沿长度方向保持不变,因此可将大坝看作无限长的实体模型,满足平面应变问题的几何条件;对截面进行受力分析,作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力,满足平面应变问题的载荷条件。因此该问题属于平面应变问题,大坝所受的载荷为面载荷,分布情况及方向如图1-2所示,建立几何模型,进行求解。 假设大坝的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比σ=0.3 三、第1问的有限元建模 本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。1)设置计算类型:两者因几何条件和载荷条件均满足平面应变问题,故均取Preferences为Structural 2)选择单元类型:三节点常应变单元选择的类型是PLANE42(Quad 4node42),该单元属于是四节点单元类型,在网格划分时可以对节点数目控制使其蜕化为三节点单元;六节点三角形单元选择的类型是PLANE183(Quad 8node183),该单元属于是八节点单元类型,在网格划分时可以对节点数目控制使其蜕化为六节点单元。因研究的问题为平面应变问题,故对Element behavior(K3)设置为plane strain。 3)定义材料参数:按以上假设大坝材料为钢,设定:ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 → OK 4)生成几何模型: a. 生成特征点:ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS→依次输入三个点的坐标:

华科大有限元分析题及大作业题答案——船海专业(DOC)

姓名:学号:班级:

有限元分析及应用作业报告 一、问题描述 图示无限长刚性地基上的三角形大坝,受齐顶的水压力作用,试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析,并对以下几种计算方案进行比较: 1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算; 2)分别采用不同数量的三节点常应变单元计算; 3)当选常应变三角单元时,分别采用不同划分方案计算。

二、几何建模与分析 图1-2力学模型 由于大坝长度>>横截面尺寸,且横截面沿长度方向保持不变,因此可将大坝看作无限长的实体模型,满足平面应变问题的几何条件;对截面进行受力分析,作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力,满足平面应变问题的载荷条件。因此该问题属于平面应变问题,大坝所受的载荷为面载荷,分布情况及方向如图1-2所示,建立几何模型,进行求解。 假设大坝的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比σ=0.3 三、第1问的有限元建模 本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。 1)设置计算类型:两者因几何条件和载荷条件均满足平面应变问题,故均取Preferences为Structural 2)选择单元类型:三节点常应变单元选择的类型是PLANE42(Quad 4node42),该单元属于是四节点单元类型,在网格划分时可以对节点数目控制使其蜕化为三节点单元;六节点三角形单元选择的类型是PLANE183(Quad 8node183),该单元属于是八节点单元类型,在网格划分时可以对节点数目控制使其蜕化为六节点单元。因研究的问题为平面应变问题,故对Element behavior(K3)设置为plane strain。 3)定义材料参数 4)生成几何模 a. 生成特征点 b.生成坝体截面 5)网格化分:划分网格时,拾取所有线段设定input NDIV 为10,选择网格划分方式为Tri+Mapped,最后得到200个单元。 6)模型施加约束: 约束采用的是对底面BC全约束。 大坝所受载荷形式为Pressure,作用在AB面上,分析时施加在L AB上,方向水平向右,载荷大小沿L AB由小到大均匀分布(见图1-2)。以B为坐标原点,BA方向为纵轴y,则沿着y方向的受力大小可表示为: ρ(1) = gh P- =ρ g = - 10 {* } 98000 98000 (Y ) y

模态分析中的几个基本概念模态分析中的几个基本概念分析

模态分析中的几个基本概念 物体按照某一阶固有频率振动时,物体上各个点偏离平衡位置的位移是满足一定的比例关系的,可以用一个向量表示,这个就称之为模态。模态这个概念一般是在振动领域所用,你可以初步的理解为振动状态,我们都知道每个物体都具有自己的固有频率,在外力的激励作用下,物体会表现出不同的振动特性。一阶模态是外力的激励频率与物体固有频率相等的时候出现的,此时物体的振动形态叫做一阶振型或主振型;二阶模态是外力的激励频率是物体固有频率的两倍时候出现,此时的振动外形叫做二阶振型,以依次类推。一般来讲,外界激励的频率非常复杂,物体在这种复杂的外界激励下的振动反应是各阶振型的复合。模态是结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。有限元中模态分析的本质是求矩阵的特征值问题,所以“阶数”就是指特征值的个数。将特征值从小到大排列就是阶次。实际的分析对象是无限维的,所以其模态具有无穷阶。但是对于运动起主导作用的只是前面的几阶模态,所以计算时根据需要计算前几阶的。一个物体有很多个固有振动频率(理论上无穷多个),按照从小到大顺序,第一个就叫第一阶固有频率,依次类推。所以模态的阶数就是对应的固有频率的阶数。振型是指体系的一种固有的特性。它与固有频率相对应,即为对应固有频率体系自身振动的形态。每一阶固有频率都对应一种振型。振型与体系实际的振动形态不一定相同。振型对应于频率而言,一个固有频率对应于一个振型。按照频率从低到高的排列,来说第一振型,第二振型等等。此处的振型就是指在该固有频率下结构的振动形态,频率越高则振动周期越小。在实验中,我们就是通过用一定的频率对结构进行激振,观测相应点的位移状况,当观测点的位移达到最大时,此时频率即为固有频率。实际结构的振动形态并不是一个规则的形状,而是各阶振型相叠加的结果。 固有频率也称为自然频率( natural frequency)。物体做自由振动时,其位移随时间按正弦或余弦规律变化,振动的频率与初始条件无关,而仅与系统的固有特性有关(如质量、形状、材质等),称为固有频率,其对应周期称为固有周期。 物体做自由振动时,其位移随时间按正弦规律变化,又称为简谐振动。简谐振动的振幅及初相位与振动的初始条件有关,振动的周期或频率与初始条件无关,而与系统的固有特性有关,称为固有频率或者固有周期。 物体的频率与它的硬度、质量、外形尺寸有关,当其发生形变时,弹力使其恢复。弹力主要与尺寸和硬度有关,质量影响其加速度。同样外形时,硬度高的频率高,质量大的频率低。一个系统的质量分布,内部的弹性以及其他的力学性质决定 模态扩展是为了是结果在后处理器中观察而设置的,原因如下: 求解器的输出内容主要是固有频率,固有频率被写到输出文件Jobname.OUT 及振型文件Jobnmae.MODE 中,输出内容中也可以包含缩减的振型和参与因子表,这取决于对分析选项和输出控制的设置,由于振型现在还没有被写到数据库或结果文件中,因此不能对结果进行后处理,要进行后处理,必须对模态进行扩展。在模态分析中,我们用“扩展”这个词指将振型写入结果文件。也就是说,扩展模态不仅适用于Reduced 模态提取方法得到的缩减振型,而且也适用与其他模态提取方法得到的完整振型。因此,如果想在后处理器中观察振型,必须先扩展模态。谱分析中的模态合并是因为激励谱是其实是由一系列的激励组合成的一个谱,里面的频率不会是只有一个,而不同的激励频率对于结构产生的结果是不一样的,对于结果的贡献也是不一样的,所以要选择模态组合法对模态进行组合,得到最终的响应结果。

有限元法课后习题答案

1、有限元是近似求解一般连续场问题的数值方法 2、有限元法将连续的求解域离散为若干个子域,得到有限个单元,单元和单元之间用节点连接 3、直梁在外力的作用下,横截面的内力有剪力和弯矩两个. 4、平面刚架结构在外力的作用下,横截面上的内力有轴力、剪力、弯矩. 5、进行直梁有限元分析,平面刚架单元上每个节点的节点位移为挠度和转角 6、平面刚架有限元分析,节点位移有轴向位移、横向位移、转角。 7、在弹性和小变形下,节点力和节点位移关系是线性关系。 8、弹性力学问题的方程个数有15个,未知量个数有15个。 9、弹性力学平面问题方程个数有8,未知数8个。 10、几何方程是研究应变和位移之间关系的方程 11、物理方程是描述应力和应变关系的方程 12、平衡方程反映了应力和体力之间关系的 13、把经过物体内任意一点各个截面上的应力状况叫做一点的应力状态 14、9形函数在单元上节点上的值,具有本点为_1_.它点为零的性质,并且在三角形单元的任一节点上,三个行函数之和为_1_ 15、形函数是_三角形_单元内部坐标的_线性_函数,他反映了单元的_位移_状态 16、在进行节点编号时,同一单元的相邻节点的号码差尽量小. 17、三角形单元的位移模式为_线性位移模式_- 18、矩形单元的位移模式为__双线性位移模式_

19、在选择多项式位移模式的阶次时,要求_所选的位移模式应该与局部坐标系的方位无关的性质为几何_各向同性 20、单元刚度矩阵描述了_节点力_和_节点位移之间的关系 21、矩形单元边界上位移是连续变化的 1.诉述有限元法的定义 答: 有限元法是近似求解一般连续场问题的数值方法 2.有限元法的基本思想是什么 答: 首先,将表示结构的连续离散为若干个子域,单元之间通过其边界上的节点连接成组合体。其次,用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。 3.有限元法的分类和基本步骤有哪些 答: 分类: 位移法、力法、混合法;步骤: 结构的离散化,单元分析,单元集成,引入约束条件,求解线性方程组,得出节点位移。 4.有限元法有哪些优缺点 答: 优点:

有限元模态分析报告实例

ANSYS 模态分析实例 5.2ANSYS 建模 该课题研究的弹性联轴器造型如下图 5.2 : 图勺2弹性联轴器 1-联接柴油机大铁圈;茁橡胶膜片;3-联接电动机小铁圈 在ANSYS 中建立模型,先通过建立如 5.2所式二分之一的剖面图,通过绕中轴线 旋转建立模拟模型如下图 5.3资料个人收集整理,勿做商业用途 _.:q: 4 1(. 片三 _」」_止

5.3单元选择和网格划分 由于模型是三给实体模型,故考虑选择三维单元,模型中没有圆弧结构,用六面体单元划分网格不会产生不规则或者畸变的单元,使分析不能进行下去,所以采用六面体单元。经比较分析,决定采用六面体八结点单元SOLID185,用自由划分的方式划 分模型实体。课题主要研究对象是联轴器中橡胶元件,在自由划分的时候,中间件2 网格选择最小的网格,smart size设置为1,两端铁圈的smart size设置为6,网格划分 后模型如图5.4。资料个人收集整理,勿做商业用途 5.4边界约束 建立柱坐标系R- &Z,如5-5所示,R为径间,Z为轴向

选择联轴器两个铁圈的端面,对其面上的节点进行坐标变换,变换到如图5.5所示的柱坐标系,约束节点R,Z方向的自由度,即节点只能绕Z轴线转资料个人收集整理,勿做商业用途 5.5联轴器模态分析 模态分析用于确定设计中的结构或者机器部件振动特性(固有频率和振型),也是瞬态变动力学分析和谐响应分析和谱分析的起点。资料个人收集整理,勿做商业用途在模态分析中要注意:ANSYS模态分析是线性分析,任何非线性因素都会被忽略。因此在设置中间件2的材料属性时,选用elastic材料。资料个人收集整理,勿做商业用途 5.5.1联轴器材料的设置 材料参数设置如下表5-1 : 表5.1材料参数设置 表5.1材料参数设置 5.5.2联轴器振动特性的有限元计算结果及说明 求解方法选择Damped方法,频率计算结果如表5-2,振型结果为图5.6: 表5.2固有频率

模态分析有限元仿真分析学习心得

有限元仿真分析学习心得 1 有限元分析方法原理 有限元分析(FEA,Finite Element Analysis)利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。还利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。 有限元法是随着电子计算机发展而迅速发展起来的一种工程力学问题的数值求解方法。20世纪50年代初,它首先应用于连续体力学领域—飞机结构静、动态特性分析之中,用以求得结构的变形、应力、固有频率以及阵型。由于其方法的有效性,迅速被推广应用于机械结构分析中。随着电子计算机的发展,有限元法从固体力学领域扩展到流体力学、传热学、电磁学、生物工程学、声学等。 随着计算机科学与应用技术的发展,有限元理论日益完善,随之涌现了一大批通用和专业的有限元计算软件。其中,通用有限元软件以ANSYS,MSC公司旗下系列软件为杰出代表,专业软件以ABAQUS、LS-DYNA、Fluent、ADAMS 为代表。 ANSYS作为最著名通用和有效的商用有限元软件之一,集机构、传热、流体、电磁、碰撞爆破分析于一体,具有强大的前后处理及计算分析能力,能够进行多场耦合,结构-热、流体-结构、电-磁场的耦合处理求解等。 有限元分析一般由以下基本步骤组成: ①建立求解域,并将之离散化成有限个单元,即将问题分解成单元和节点; ②假定描述单元物理属性的形(shape)函数,即用一个近似的连续函数描述每个单元的解; ③建立单元刚度方程; ④组装单元,构造总刚度矩阵; ⑤应用边界条件和初值条件,施加载荷; ⑥求解线性或者非线性微分方程组得到节点值,如不同节点的位移; ⑦通过后处理获得最大应力、应变等信息。 结构的离散化是有限元的基础。所谓离散化就是将分析的结构分割成为有限

有限元复习题答案

1、何为有限元法?其基本思想是什么? 有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法,该方法以计算机为手段,采用分片近似,进而逼近整体的研究思想求解物理问题。 基本思想是化整为零集零为整。 2、为什么说有限元法是近似的方法,体现在哪里? 有两点:用离散单元的组合体来逼近原始结构,体现了几何上的近似;而用近似函数逼近未知变量在单元内的真实解,体现了数学上的近似。 3、单元、节点的概念? 节点:表达实际结构几何对象之间相互连接方式的概念 单元:网格划分中的每一个小部分称为单元,网格间相互联结点称为节点 4、有限元法分析过程可归纳为几个步骤? 结构离散化、单元分析、整体分析 5、有限元方法分几种?本课程讲授的是哪一种? 位移法、力法、混合法本课程讲授位移法 6、弹性力学的基本变量是什么?何为几何方程、物理方程及虚功方程?弹性矩阵的特点? 弹性力学变量:外力、应力、应变和位移。 描述弹性体应变分量与位移分量之间的方程称为几何方程;物理方程描述应力分量与应变分量之间的关系;弹性体上外力在虚位移发生过程中所做的虚功与储存在弹性体内的需应变能相等。 弹性矩阵由材料的弹性模量和泊松比确定,与坐标位置无关。 7、何为平面应力问题和平面应变问题? 平面应力问题:在结构上满足a几何条件:研究对象是等厚度薄板。b载荷条件:作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布,而在两板面无外力作用。 平面应变问题:满足a几何条件:长柱体,即长度方向的尺寸远远大于横截面的尺寸,且横截面沿长度方向不变。b载荷条件:作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力两条件的弹性力学问题。 1、何为结构的离散化?离散化的目的?何为有限元模型? ①离散化:把连续的结构看成由有限个单元组成的集合体。②目的:建立有限元计算模型③通常把由节点,单元及相应的节点载荷和节点约束构成的模型称为有限元模型2、结构离散化时,划分单元数目的多少以及疏密分布,将直接影响到什么?确定单元数量的原则?通常如何设置节点?

有限元分析大作业试题

有限元分析习题及大作业试题 要求:1)个人按上机指南步骤至少选择习题中3个习题独立完成,并将计算结果上交; 2)以小组为单位完成有限元分析计算; 3)以小组为单位编写计算分析报告; 4)计算分析报告应包括以下部分: A、问题描述及数学建模; B、有限元建模(单元选择、结点布置及规模、网格划分方 案、载荷及边界条件处理、求解控制) C、计算结果及结果分析(位移分析、应力分析、正确性分 析评判) D、多方案计算比较(结点规模增减对精度的影响分析、单 元改变对精度的影响分析、不同网格划分方案对结果的 影响分析等) E、建议与体会 4)11月1日前必须完成,并递交计算分析报告(报告要求打印)。

习题及上机指南:(试题见上机指南) 例题1 坝体的有限元建模与受力分析 例题2 平板的有限元建模与变形分析 例题1:平板的有限元建模与变形分析 计算分析模型如图1-1 所示, 习题文件名: plane 0.5 m 0.5 m 0.5 m 0.5 m 板承受均布载荷:1.0e 5 P a 图1-1 受均布载荷作用的平板计算分析模型 1.1 进入ANSYS 程序 →ANSYSED 6.1 →Interactive →change the working directory into yours →input Initial jobname: plane →Run 1.2设置计算类型 ANSYS Main Menu : Preferences →select Structural → OK 1.3选择单元类型 ANSYS Main Menu : Preprocessor →Element T ype →Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK (back to Element T ypes window) → Options… →select K3: Plane stress w/thk →OK →Close (the Element T ype window) 1.4定义材料参数 ANSYS Main Menu : Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY :0.3 → OK 1.5定义实常数 ANSYS Main Menu: Preprocessor →Real Constant s… →Add … →select T ype 1→ OK →input THK:1 →OK →Close (the Real Constants Window)

有限元分析及应用大作业

有限元分析及应用大作业 作业要求: 1)个人按上机指南步骤至少选择习题中3个习题独立完成,并将计算结果上交; 也可根据自己科研工作给出计算实例。 2)以小组为单位完成有限元分析计算; 3)以小组为单位编写计算分析报告; 4)计算分析报告应包括以下部分: A、问题描述及数学建模; B、有限元建模(单元选择、结点布置及规模、网格划分方案、载荷及边界 条件处理、求解控制) C、计算结果及结果分析(位移分析、应力分析、正确性分析评判) D、多方案计算比较(结点规模增减对精度的影响分析、单元改变对精度的 影响分析、不同网格划分方案对结果的影响分析等) 题一:图示无限长刚性地基上的三角形大坝,受齐顶的水压力作用,试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析,并对以下几种计算方案进行比较: 1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算;(注意ANSYS中用四边形单元退化为三节点三角形单元) 2)分别采用不同数量的三节点常应变单元计算; 3)当选常应变三角单元时,分别采用不同划分方案计算。 解:1.建模: 由于大坝长度>>横截面尺寸,且横截面沿长度方向保持不变,因此可将大坝看作无限长的实体模型,满足平面应变问题的几何条件;对截面进行受力分析,作

用于大坝上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力,满足平面应变问题的载荷条件。因此该问题属于平面应变问题,大坝所受的载荷为面载荷,分布情况P=98000-9800*Y;建立几何模型,进行求解;假设大坝的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比σ=0.3; 2:有限元建模过程: 2.1 进入ANSYS : 程序→ANSYS APDL 15.0 2.2设置计算类型: ANSYS Main Menu: Preferences →select Structural →OK 2.3选择单元类型: ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 182(三节点常应变单元选择Solid Quad 4node 182,六节点三角形单元选择Solid Quad 8node 183)→OK (back to Element Types window) →Option →select K3: Plane Strain →OK→Close (the Element Type window) 2.4定义材料参数: ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK 2.5生成几何模型: 生成特征点: ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS →依次输入四个点的坐标:input:1(0,0),2(10,0),3(1,5),4(0.45,5) →OK 生成坝体截面: ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点,1(0,0),2(6,0),3(0,10) →OK 2.6 网格划分: ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条直角边:OK→input NDIV: 15 →Apply→依次拾取斜边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh:Areas, Shape: tri, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window) 2.7 模型施加约束: 给底边施加x和y方向的约束: ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On lines →pick the lines →OK →select Lab2:UX, UY →OK 给竖直边施加y方向的分布载荷: ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量;2) 在Result窗口中出现{X},写入所施加的载荷函数: 98000-9800*{Y};3) File>Save(文件扩展名:func) →返回:Parameters →Functions →Read from file:将需要的.func文件打开,参数名取meng,它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取竖直边;OK →在下拉列表框中,选择:Existing table →OK →选择需要的载荷为meng参数名→OK 2.8 分析计算: ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load

各种模态分析方法总结及比较

各种模态分析方法总结与比较 一、模态分析 模态分析是计算或试验分析固有频率、阻尼比和模态振型这些模态参数的过程。 模态分析的理论经典定义:将线性定常系统振动微分方程组中的物理坐标变换为模态坐标,使方程组解耦,成为一组以模态坐标及模态参数描述的独立方程,以便求出系统的模态参数。坐标变换的变换矩阵为模态矩阵,其每列为模态振型。 模态分析是研究结构动力特性一种近代方法,是系统辨别方法在工程振动领域中的应用。模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。这些模态参数可以由计算或试验分析取得,这样一个计算或试验分析过程称为模态分析。这个分析过程如果是由有限元计算的方法取得的,则称为计算模记分析;如果通过试验将采集的系统输入与输出信号经过参数识别获得模态参数,称为试验模态分析。通常,模态分析都是指试验模态分析。振动模态是弹性结构的固有的、整体的特性。如果通过模态分析方法搞清楚了结构物在某一易受影响的频率范围内各阶主要模态的特性,就可能预言结构在此频段内在外部或内部各种振源作用下实际振动响应。因此,模态分析是结构动态设计及设备的故障诊断的重要方法。 模态分析最终目标是在识别出系统的模态参数,为结构系统的振动特性分析、振动故障诊断和预报以及结构动力特性的优化设计提供依据。二、各模态分析方法的总结

(一)单自由度法 一般来说,一个系统的动态响应是它的若干阶模态振型的叠加。但是如果假定在给定的频带内只有一个模态是重要的,那么该模态的参数可以单独确定。以这个假定为根据的模态参数识别方法叫做单自由度(SDOF)法n1。在给定的频带范围内,结构的动态特性的时域表达表示近似为: ()[]}{}{T R R t r Q e t h r ψψλ= 2-1 而频域表示则近似为: ()[]}}{ {()[]2ωλωψψωLR UR j Q j h r t r r r -+-= 2-2 单自由度系统是一种很快速的方法,几乎不需要什么计算时间和计算机内存。 这种单自由度的假定只有当系统的各阶模态能够很好解耦时才是正确的。然而实际情况通常并不是这样的,所以就需要用包含若干模态的模型对测得的数据进行近似,同时识别这些参数的模态,就是所谓的多自由度(MDOF)法。 单自由度算法运算速度很快,几乎不需要什么计算和计算机内存,因此在当前小型二通道或四通道傅立叶分析仪中,都把这种方法做成内置选项。然而随着计算机的发展,内存不断扩大,计算速度越来越快,在大多数实际应用中,单自由度方法已经让位给更加复杂的多自由度方法。 1、峰值检测 峰值检测是一种单自由度方法,它是频域中的模态模型为根据对系统极点进行局部估计(固有频率和阻尼)。峰值检测方法基于这样的事实:在固有频率附近,频响函数通过自己的极值,此时其实部为零(同相部分最

有限元分析及应用例子FEM14

第9章受内外压筒体的有限元建模与应力变形分析(Project 2) 计算分析模型如图9-1 所示, 习题文件名: cylinder。 X (a) σO=100N/mm2 σI =200N/mm2 γ =7.85g/cm3 μ =0.3 E =210000N/mm2 (b) 图9-1 计算分析模型 9.1进入ANSYS 程序→ANSYSED 6.1ed →Interactive →change the working directory into yours→input Initial jobname: cylinder→Run 9.2 设置计算类型 ANSYS Main Menu: Preferences…→select Structural →OK 9.3 选择单元类型 ANSYS Main Menu: Preprocessor → Element Type →Add/Edit/Delete… → Add… →select Solid Quad 4node 42 →Apply →select Solid Brick 8node 45 → OK → Close (the Element

Types window) 9.4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Materials Models →Structural→Lineal →Elastic→Isotropic…→input EX:2.1e5, PRXY:0.3→ OK 关闭材料定义窗口 9.5构造筒体模型 ?生成模型截平面 ANSYS Main Menu: Preprocessor →Modeling→Create →Keypoints →In Active CS… →按次序输入横截平面的十个特征点和旋转对称轴上两点坐标(十个特征点:(300,0,0), (480,0,0), (480,100,0), (400,100,0), (400,700,0), (480,700,0), (480,800,0), (300,800,0), (300,650,0), (300,150,0),对称轴上两点:(0,0,0), (0,800,0))(每次输入完毕,用Apply结束,0可以不输入) →Cancel (back to Create window) →-Areas- Arbitrary → Through KPs →依次连接截面边线上的十个特征点(注意在选完第10点后结束,不要再选第1点)→ OK ?对平面进行网格划分 ANSYS Main Menu: Preprocessor →Meshing→Mesh Tool →(Size Controls) Globl: Set →input SIZE (element edge length): 50 →OK (back to MeshTool window)→Mesh → Pick All (in Picking Menu) → Close( the MeshTool window) ?用旋转法生成筒体模型 ANSYS Main Menu: Preprocessor →Modeling→Operate →Extrude→Elem Ext Opts→select TYPE:SOLID 45→Element sizing options for extrusion No. Elem divs: 1→OK (back to Extrude window)→Areas →About Axis →Pick All(in Picking Menu)→OK→Pick the two keypoints (11,12) of the Symmetrical Axis → OK→input ARC: 90; NSEG: 3→ OK 9.6 模型加位移约束 ANSYS Main Menu: Solution→Define Loads →Apply→Structural→Displacement ?两截面分别加Z, X方向的约束 ANSYS Utility Menu: Select → Entities…→Nodes → By Location →select X coordinates →input 0→ OK (back to Displacement window)→On Nodes → Pick All(in Picking Menu) → select Lab2:UX →OK →ANSYS Utility Menu: Select → Everything ANSYS Utility Menu: Select → Entities…→ Nodes → By Location →select Z coordinates →input 0→ OK (back to Displacement window)→On Nodes →Pick All(in Picking Menu) → select Lab2:UZ →OK →ANSYS Utility Menu: Select →Everything ?底面加Y方向的约束 ANSYS Utility Menu: Select → Entities… → Nodes → By Location →select Y coordinates →input 0→ OK (back to Displacement window)→On Nodes →Pick All(in Picking Menu) →

有限元分析及其应用思考题附答案2012

有限元分析及其应用-2010 思考题: 1、有限元法的基本思想是什么?有限元法的基本步骤有那些?其中“离散”的含义是什 么?是如何将无限自由度问题转化为有限自由度问题的? 答:基本思想:几何离散和分片插值。 基本步骤:结构离散、单元分析和整体分析。 离散的含义:用假想的线或面将连续物体分割成由有限个单元组成的集合,且单元之间仅在节点处连接,单元之间的作用仅由节点传递。当单元趋近无限小,节点无限多,则这种离散结构将趋近于实际的连续结构。 2、有限元法与经典的差分法、里兹法有何区别? 区别:差分法:均匀离散求解域,差分代替微分,要求规则边界,几何形状复杂精度较低; 里兹法:根据描述问题的微分方程和相应的定解构造等价的泛函表达式,求得近似解; 有限元:基于变分法,采用分片近似进而逼近总体的求解微分方程的数值计算方法。 3、一根单位长度重量为q的悬挂直杆,上端固定,下端受垂直向下的外力P,试 1)建立其受拉伸的微分方程及边界条件; 2)构造其泛函形式; 3)基于有限元基本思想和泛函求极值构造其有限元的计算格式(即最小势能原理)。4、以简单实例为对象,分别按虚功原理和变分原理导出有限元法的基本格式(单元刚度矩 阵)。 5、什么是节点力和节点载荷?两者有何区别? 答:节点力:单元与单元之间通过节点相互作用 节点载荷:作用于节点上的外载 6、单元刚度矩阵和整体刚度矩阵各有何特点?其中每个矩阵元素的物理意义是什么(按自 由度和节点解释)? 答:单元刚度矩阵:对称性、奇异性、主对角线恒为正 整体刚度矩阵:对称性、奇异性、主对角线恒为正、稀疏性、带状性。 Kij,表示j节点产生单位位移、其他节点位移为零时作用i节点的力,节点力等于节点位移与单元刚度元素乘积之和。 7、单元的形函数具有什么特点?有哪些性质? 答:形函数的特点:Ni为x,y的坐标函数,与位移函数有相同的阶次。 形函数Ni在i节点的值为1,而在其他节点上的值为0; 单元内任一点的形函数之和恒等于1; 形函数的值在0~1间变化。 8、描述弹性体的基本变量是什么?基本方程有哪些组成? 答:基本变量:外力、应力、应变、位移 基本方程:平衡方程、几何方程、物理方程、几何条件 9、何谓应力、应变、位移的概念?应力与强度是什么关系? 答:应力:lim△Q/△A=S △A→0 应变:物体形状的改变 位移:弹性体内质点位置的变化 10、问题的微分方程提法、等效积分提法和泛函变分提法之间有何关系?何谓“强形 式”?何谓“弱形式”,两者有何区别?建立弱形式的关键步骤是什么?

相关文档
最新文档