有限元线性静态分析
有限元分析FEA
广州有道计算机科技有限公司有限元分析FEA有限元法(FEA,Finite Element Analysis)的基本概念是用较简单的问题代替复杂问题后再求解。
它将求解域看成是由许多称为有限元的小的互连子域组成,对每一单元假定一个合适的(较简单的)近似解,然后推导求解这个域总的满足条件(如结构的平衡条件),从而得到问题的解。
这个解不是准确解,而是近似解,因为实际问题被较简单的问题所代替。
由于大多数实际问题难以得到准确解,而有限元不仅计算精度高,而且能适应各种复杂形状,因而成为行之有效的工程分析手段。
有限元分析(FEA,Finite Element Analysis)利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。
还利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。
大型通用有限元商业软件:如ANSYS可以分析多学科的问题,例如:机械、电磁、热力学等;电机有限元分析软件NASTRAN等。
还有三维结构设计方面的UG、CATIA、Proe等都是比较强大的。
国产有限元软件:FEPG、SciFEA、,JiFEX、KMAS等有限元法:把求解区域看作由许多小的在节点处相互连接的单元(子域)所构成,其模型给出基本方程的分片(子域)近似解,由于单元(子域)可以被分割成各种形状和大小不同的尺寸,所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件。
有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
有限元法的收敛性是指:当网格逐渐加密时,有限元解答的序列收敛到精确解;或者当单元尺寸固定时,每个单元的自由度数越多,有限元的解答就越趋近于精确解。
ANSYS Workbench 17·0有限元分析:第5章-线性静态结构分析
第5章 线性静态结构分析 在工程应用中,经常会遇到计算在固定不变的载荷作用下的结构效应,主要有平面应力、平面应变、轴对称、梁及桁架分析、壳分析、接触分析等问题的求解,这些问题均是线性静态结构问题,线性静态结构分析是有限元(★ 掌握线性静态结构分析的基本过程。
5.1 线性静态结构分析概述线性静态结构分析(Lines Static Structural Analysis )用于计算在固定不变的载荷作用下结构的效应,它不考虑惯性和阻尼的影响,如结构随时间变化载荷等情况。
静力分析可以计算固定不变的惯性载荷对结构的影响(如重力和离心力),以及可以近似为等价静力作用的随时间变化载荷(如通常在许多建筑规范中所定义的等价静力风载和地震载荷)。
在经典力学理论中,物体的动力学通用方程为:[]()[]()[]{}(){}M x C x K x F t ++=&&&其中[]M 为质量矩阵,[]C 为阻尼矩阵,[]K 为刚度系数矩阵,{}x 为位移矢量,{}F 为力矢量。
在线性静态结构分析中,力与时间无关,因此位移{}x 可以由下面的矩阵方程解出:[]{}{}K x F =在线性静态结构分析中,假设[]K 为一常量矩阵且必须是连续的,材料必须满足线弹性、小变形理论,边界条件允许包含非线性的边界条件,{}F 为静态加载到模型上的力,该力不随时间变化,不包括惯性影响因素(质量、阻尼等)。
静力分析用于计算由不包括惯性和阻尼效应的载荷作用于结构或部件上引起的位移、应力、应变和力等。
假定载荷和响应是固定不变的,即假定载荷和结构的响应随时间的变化而缓慢变化。
静力分析所施加的载荷包括:ANSYS Workbench 17.0有限元分析从入门到精通外部施加的作用力和压力。
稳态的惯性力(如中力和离心力)。
位移载荷。
温度载荷。
5.2 线性静态结构的分析流程在ANSYS Workbench 左侧工具箱中Analysis Systems 下的Static Structural 上按住鼠标左键拖动到项目管理区,或双击Static Structural 选项,即可创建静态结构分析项目,如图5-1所示。
IDESA有限元分析_第6篇第33章 装配体有限元模型
第33章MasterFEM 教程:装配体的有限元分析本章介绍怎样定义装配体的有限元模型进行分析,用户将学会:∙定义一个装配体的有限元模型。
∙施加约束和载荷。
∙定义并产生网格。
∙线性静态分析。
∙评估分析结果。
设置设置工作单位为mm (milli newton).导入Archive文件,把Frame assembly拿上工作台。
FileImport注意事项:Archive文件为压缩文件,放在教学光盘中。
路径:光盘根目录下/frame.arc文件内容为已做好的装配体模型文件。
·1·创建有限元模型转变到Boundary Conditions任务中创建此装配体的有限元模型文件:选择Create FE Model图标在FE Model Create 菜单里选择get图标·2··3·在Select Parent Part 菜单里选择Frame (assembly) 选择OK命名有限元模型文件:FEM打开Geometry Based Analysis Only 选项 选择OK命名新的几何模型文件:Part 选择OK设置放入装配体到抽屉中: 选择Manage Bins 图标在Manage Bins 菜单里选择Frame (assembly) 选择Put Away 命令·4·选择Dismiss 命令更改零件显示过滤器: 选择Display Filter 图标在Display Filter 菜单里选择Parts 关闭Coordinate Systems 开关 打开Local Origin 开关 每次都选择OK施加边界条件全约束leg 底边在Top piece 顶部创建轴对称约束: X translation free Y translation free Z translation constant X rotation constant Y rotation constantZ rotation free·5·在Side piece 顶部创建轴对称约束:X translation constant Y translation free Z translation free X rotation feeY rotation constant Z rotation constant在表面上定义均布载荷1000 mN/mm^2在同样的载荷集中,创建10G 的加速度。
有限元分析
工程分析模块概述 三维模型相关处理 有限元分析过程
分析模块概述
工程分析
Catia有限元分析包括静态分析(Static Analysis)和动态分析。
动态分析
Catia动态分析又分为自由固有频率分析和约 束固有频率分析。
常用动态分析
专业有限元软件如ANSYS、ABAQUS、 NASTRAN的动态分析包括模态分析、谐分析、谱 分析和瞬态分析等。
生成分析报告 可将有限元模型、多种工况及其分析结果 汇总,产生分析报告。
分析过程3-后处理1
变形云图
分析过程3-后处理2
等效应力云图
分析过程3-后处理3
分析过程3-后处理4
粒子云图
总结与练习
Summary Practice1 如右图示模型 进行静力分析 约束模态分析 Practice2 对实例分别 更改约束条件 和载荷类型数据 进行分析仿真 The end! Thank you!
分析过程3-前处理1
单元类型及尺寸
分析过程3-前处理2
网格预览
分析过程3-前处理3
约束/载荷 Case定义
分析过程3-求解及后处理
工况求解与结果预览
针对同一个有限元网格模型,可以定义不同 的边界条件和载荷类型模拟多种工况,并进工况 合并。分析结果以静态云图和动画的形式,反 出模型的最大最小量值。
3D模型与有限元模型
3D模型的简化 工程数据准备 工况仿真模拟
有限元分析过程1
确定工程分析类型 导入几何模型(已赋予材质) 设置单元类型进行网格划分 模拟施加约束边界条件 模拟施加不同类型的载荷 工况求解与结果预览 生成分析报告
材料力学中的有限元方法分析
材料力学中的有限元方法分析材料力学是研究物质初始状态至最终破坏状态之间的力学行为及其规律的科学。
有限元分析是一种数值计算方法,可以求解各种工程问题的数学模型。
有限元方法在材料力学研究中有着重要的应用,本文将从有限元方法的基本原理、材料力学中的有限元分析、有限元模拟在材料力学中的应用等方面进行分析。
一、有限元方法的基本原理有限元方法是一种通过建立复杂结构的有限元模型,将一个复杂的连续问题转化为离散问题来求解的方法。
其基本思想是将一个连续物体分割成很多小的单元,使用一些简单的解析方法求解每个小单元内的力学问题,然后将所有小单元的解组合在一起来求解整体力学问题。
有限元方法求解的过程分为以下基本步骤:1.建立有限元模型2.离散化3.施加约束4.建立刚度矩阵和荷载向量5.求解未知量二、材料力学中的有限元分析材料力学中的有限元分析是指通过有限元方法对材料力学问题进行分析、计算和评估的方法。
材料力学问题中的目标是通过施加荷载或外界力,来得到物体内部的应力和应变状态,以及其随时间和载荷变化的规律。
在建立材料力学有限元模型时,需要考虑以下因素:1.应力集中和应变集中的位置和程度2.物理边界和几何结构3.材料的力学性质和力学参数材料力学中的有限元分析包含以下几个方面:1.静态分析:研究物体在静态等效荷载下的应力状态,计算物体的静态变形。
2.动态分析:研究物体在动态载荷下的应力和应变状态,计算物体的动力响应。
3.疲劳分析:研究物体在周期性载荷下的损伤状态、损伤机理和寿命预估。
4.热力耦合分析:研究物体在温度场和应力场的共同作用下的应力和应变状态。
5.多物理场分析:研究物体在电、磁、声、液、气、红外、光、辐射等多个物理场的共同作用下的应力和应变状态。
三、有限元模拟在材料力学中的应用有限元模拟在材料力学中的应用范围非常广泛,包括了以下几个方面:1.材料的结构设计和分析2.材料的性质和参数的测试和评估3.材料的制造和加工工艺的模拟4.材料的破坏和损伤机理的研究5.材料的寿命评估和振动疲劳分析最终,有限元分析的结果可以在材料设计、材料优化和制造流程等方面提供准确的数据支持,帮助人们更好地理解材料的力学行为和性质,促进材料科学的发展。
有限元静力分析
10
二、静力分析中的载荷
11
பைடு நூலகம்
静力分析中的载荷
Definition
静力分析用于计算由那些不包括惯性和阻尼效应的 载荷作用于结构或部件上引起的位移,应力,应变和 力。固定不变的载荷和响应是一种假定;即假定载荷 和结构的响应随时间的变化非常缓慢。静力分析所 施加的载荷包括:
外部施加的作用力和压力
稳态的惯性力(如重力和离心力)
模态分析---用于计算结构的固有频率和模态。 谐波分析---用于确定结构在随时间正弦变化的载荷作用下的响应。 瞬态动力分析---用于计算结构在随时间任意变化的载荷作用下的响应,
并且可计及上述提到的静力分析中所有的非线性性质。
谱分析---是模态分析的应用拓广,用于计算由于响应谱或PSD输入
(随机振动)引起的应力和应变。
23
状态变化(包括接触)
Lesson Objectives
许多普通结构表现出一种与状态相关的非线性行为,例如,一 根只能拉伸的电缆可能是松散的,也可能是绷紧的。轴承 套可能是接触的,也可能是不接触的。冻土可能是冻结的, 也可能是融化的。这些系统的刚度由于系统状态的改变 在不同的值之间突然变化。状态改变也许和载荷直接有 关(如在电缆情况中), 也可能由某种外部原因引起 (如在冻土中的紊乱热力学条件)。ANSYS程序中单元 的激活与杀死选项用来给这种状态的变化建模。
位移载荷
温度载荷
12
三、线性静力分析和非线性静 力分析
13
线性静力分析和非线性静力分析
Definition
静力分析既可以是线性的也可以是非线性的。非 线性静力分析包括所有的非线性类型:大变形,塑性, 蠕变,应力刚化,接触(间隙)单元,超弹性单元 等。本节主要讨论线性静力分析,非线性静力分析 在下一节中介绍.
有限元方法分类
有限元方法分类
有限元方法是一种强大的数值分析工具,广泛应用于工程计算、物理模拟等领域。
按照不同的分类方式,有限元方法可以划分为多个类别:
1. 按求解问题类型划分:结构力学有限元、热传导有限元、电磁场有限元、流体力学有限元、声学有限元等,分别对应于解决固体结构应力变形、热量传递、电磁场分布、流体流动以及声音传播等问题。
2. 按单元性质划分:线性有限元和非线性有限元。
线性有限元处理的是线性问题,如弹性力学中的小变形问题;非线性有限元则是针对材料非线性、几何非线性等问题。
3. 按时间因素考虑划分:静态有限元分析和动态有限元分析。
静态分析处理稳态问题,不考虑随时间变化的影响;动态分析则考虑了随时间演变的效应,如瞬态动力响应。
4. 按离散形式划分:等参有限元、非等参有限元。
等参有限元在单元内部采用一致的坐标变换,非等参有限元则根据实际情况灵活选择节点和形状函数。
5. 按求解流程划分:直接法有限元和迭代法有限元。
直接法直接求解全局刚度矩阵,而迭代法则通过多次迭代逐步逼近解。
总之,有限元方法因其灵活性和普适性,能够处理各类复杂的物理问题,已成为现代工程与科学研究中不可或缺的分析手段。
SolidWorks Simulation 培训——线性静态分析
什么是有限元
• 广义的有限元分析:在运用有限元分析软件,对实际产品进行相应 的力学分析基础之上,通过产品实验的结合,并对实验数据及仿真 数据的归纳总结,得到企业产品的一系列性能数据。
• 对工程师的要求:专业的分析工程师,拥有深厚的力学专业功底和 实验数据处理能力
• 给企业带来的效益:创新,技术自主,企业知识库 • 企业遇到的问题:时间周期过长
部分行业实际案例
线性动力学分析案例:钱江电器抗震分析
钱江电器于2013年初向智利出口某型号变压器,由于智利方需要 一份依据IEEE Std 693-2005国际标准的有限元抗震分析报告,因此钱江 电器请求我们华睿提供技术支持,在指定时间内成功完成分析并提交 报告,促成钱江电器完成订单。
高压管套 变压器器身
案例:1992年,法国在进行210次核试验之后,宣布其核武器更新将依 靠数值仿真计划来实现,该项目15年总投资210亿欧元,这标志着发达 国家在复杂武器工程分析方面已经进入了大规模并行计算时代。
油箱
散热器 底座
底座
橡胶件非线性静态分析案例:皮尔轴承橡胶密封
上海皮尔轴承有限公司隶属于美国皮尔轴承公司的全资子公司,于
2008年被全球轴承业巨头SKF收购。
2013年中旬和我们合作,解决橡胶件件密封及磨损相关问题解决
方案。
轴
承
防
尘
盖
装
配
分
析 数
实验
据
全
局
规
仿真
划
温度 磨损量 摩擦力 接触压强 变形过程
不锈钢的化学性能:不容易在高温发生化学反应 手柄温度安全
锅底受热不均
什么是有限元
铜的化学性能:容易在高温发生化学反应 手柄温度过高
悬臂梁—有限元ABAQUS线性静力学分析实例分解
线性静力学分析实例——以悬臂梁为例线性静力学问题是简单且常见的有限元分析类型,不涉及任何非线性(材料非线性、几何非线性、接触等),也不考虑惯性及时间相关的材料属性。
在ABAQUS 中,该类问题通常采用静态通用(Static,General)分析步或静态线性摄动(Static,Linear perturbation)分析步进行分析。
线性静力学问题很容易求解,往往用户更关系的是计算效率和求解效率,希望在获得较高精度的前提下尽量缩短计算时间,特别是大型模型。
这主要取决于网格的划分,包括种子的设置、网格控制和单元类型的选取。
在一般的分析中,应尽量选用精度和效率都较高的二次四边形/六面体单元,在主要的分析部位设置较密的种子;若主要分析部位的网格没有大的扭曲,使用非协调单元(如CPS4I、C3D8I)的性价比很高。
对于复杂模型,可以采用分割模型的方法划分二次四边形/六面体单元;有时分割过程过于繁琐,用户可以采用精度较高的二次三角形/四面体单元进行网格划分。
悬臂梁的线性静力学分析1.1 问题的描述一悬臂梁左端受固定约束,右端自由,结构尺寸如图1-1 所示,求梁受载后的Mises 应力、位移分布。
材料性质:弹性模量 E 2e3 ,泊松比0.3均布载荷:F=103N图1-1 悬臂梁受均布载荷图1.2 启动ABAQUS 启动ABAQUS 有两种方法,用户可以任选一种。
(1)在Windows 操作系统中单击“开始” --“程序” --ABAQUS 6.10 --ABAQUS/CAE 。
(2)在操作系统的DOS 窗口中输入命令:abaqus cae。
启动ABAQUS/CAE 后,在出现的Start Section(开始任务)对话框中选择Create Model Database。
1.3 创建部件在ABAQUS/CAE 顶部的环境栏中,可以看到模块列表:Module :Part,这表示当前处在Par(t 部件)模块,在这个模块中可以定义模型各部分的几何形体。
有限元分析
有限元分析什么是有限兀方法?一个多用途的有限元法计算机设计程序,用于解决简单的线性静态分 析问题,也可用于解决复杂的非线性动态分析问题。
为什么使用有限元分析方法:由于常见的实验测试方法(如应变计、激光全息和光栅云纹等)对试 样制备要求高,且提供的测量值只能是平均值或断面值,而且小尺寸 封装的焊点应力应变复杂,所以,实验测试方法在电子封装可靠性这 方面的研究进展地比较慢,而通过有限元软件即可实现的数值模拟方 法发展迅速。
过程:建立BGA 单个焊点的模型 ------ 加载温度循环载荷 ---- 设定边界条件 -- 有限元分析 ---- 应力应变云图 --- 失效位置分析Mold Com poundDieDie Auac hSolder心州闻 MaskLami nate图:BGA 模型示意图林处丼单元计鼻程单凡•门出rnr世救对算求駢瓒「I泄汀汁厂广图:有限元的计算流程在利用有限元分析时借助于ANSYS软件,利用ANSYS模拟的流程如下图所示前块图:ANSYS模拟的一般流程为什么进行有限元分析?热循环实验能够直接得到焊点裂纹的产生位置。
但是,相当费时,费钱和费力。
因此,评价SMT焊点可靠性的方法更多是基于焊点的寿命预测模型和有限元模拟来预测焊点的寿命。
对于钎料合金,最根本的寿命预测方程为Man so n-Coffi n方程相关参数失效循环数;△ £ —循环应变范围;B, C —经验常数什么是纳米压痕仪?根据载荷-位移曲线和压头的面积函数计算瞬时接触面积纳米压痕测试的结果仅为一组载荷-位移曲线,并不能直接得到材料的力学性能参数优点:(1)避免了直接成像带来的误差(2)对压头载荷随压入深度连续变化进行精确测量(3)能够在相对有限的材料体积内产生很高的应变、应力。
纳米压痕仪的组成:特殊几何形状的压头、加载的施力马达、精确测量位移的传感器压头形式:金字塔形状的Berkovich压头、直角立方压头和球形压头纳米压痕测试所采用的压头为Berkovich 形压头Permanentmag 口已柱Depth sensorIndenterMitajieiic shielding ' Suppoiring springSNfReference lingBerkovich 形压头:相关参数:P是载荷,A是该载荷下投影接触面积,S是弹性接触刚度,B = 1.034, E r是弹性模量纳米压痕法采用对载荷-位移曲线的直接分析得到压痕的接触面积•次完幣的加栈和科載的压痕过程的载荷一位移曲线的示总图压痕法测量塑性变形的机理以通过压痕P-H 曲线来反推得到被压材料的弹塑性性能,Strain* c图:金属应力一应变的弹塑性行为w p从P ―― H 图中我们可以知道以下参数: 是C 、一^、H max 、P max 、刚度S 、弹性模量、w t 屈服应力、和加工硬化指数。
有限元第五讲 结构线性静力分析
4.1结构静力分析过程与步骤
一般包括建立模型、施加载荷并求解和检查结果3个步骤 4.1.1 建立模型 主要包括定义单元类型、单元实常数、材料属性和几何模型等。 建立模型注意事项: (1)单元类型必须指定为线性或非线性结构单元类型。 (2)材料属性可为线性或非线性、各向同性或正交各向异性、常量或与温 度相关的量等。 (3)必须定义杨氏模量和泊松比。 (4)对于诸如重力等惯性载荷,必须定义能计算出质量的参数,如密度等。 (5)对热载荷,必须要定义热膨胀系数。 (6)对应力、应变感兴趣的区域,网格划分比仅对位移感兴趣的区域要密。 (7)如果分析中包含非线性因素,网格应划分到能捕捉非线性因素影响的 程度。
4.1结构静力分析过程与步骤
结构静力分析主要用来分析由于稳态外载荷所引起的系统或 零部件的位移、应力、应变和作用力,很适合求解惯性及阻 尼的时间相关作用对结构响应的影响并不显著的问题,其中 稳态载荷主要包括外部施加的力和压力、稳态的惯性力,如 重力和旋转速度、施加位移、温度和热量等。静力分析可分 为线性静力分析和非线性静力分析。
4.1.2 施加载荷并求解
2.在模型上施加载荷 用户能够将载荷施加在几何模型或有限元模型上。结构静力分析的载荷 类型主要包括位移、力或力矩、压力、温度、流通量、重力和旋转角速 度等。 GUI路径为:Main Menu>Solution>DefineLoads>Apply。 指定载荷步选项主要包括普通和非线性选项,其中普通选项包括对载荷 步终止时间(Time)、对热应变计算的参考温度(Reference temperature)和 用于轴对称单元的摸态数 (Mode number)等;非线性选项包括对下面选 项的设置:时间子步数、时间步长、渐变加载还是阶跃加载、是否采用 自动时间步跟踪、平衡迭代的最大数、收敛精度、矫正预测、线搜索、 蠕变准则、求解终止选项、数据和结果文件的输入输出,以及结果外插 法等。
有限元-结构静力学分析
轧制磷青铜 轧制锰黄铜
铸铝青铜 硬铝合金 冷拔黄铜 轧制纯铜
轧制锌 轧制铝
铅 钢 铝 铸铁 不锈钢 镁 镍 玻璃 黄铜 铜 右墨 钛 钨 木材
弹性模量E GPa
115~160 151~160 200~220
210 175 115 110 105 71 91~99 110 84 69 17 207 71.7 100 190 44.8 207 46.2 106 119 36.5 102.04 344.7 11
现在有限元静定、超静定全部都可以方便计算了。
杆件的结构静力分析分类
杆系结构还可分为平面结构和空间结构。当结构的全部杆 件、支座及作用力均位于同一平面时,称结构为平面结构; 否则即为空间结构。工程中的绝大多数结构都是空间结构。 但在许多情况下往往可以引入一些适当的假定,把它们简化 为平面结构,从而避免复杂的计算并取得精度符合工程要求 的结果。在计算机发展后,习惯上常简化为平面结构的桁架 和刚架(见框架)等,已逐步转向按空间结构计算。
0.42 0.29 0.33 0.211 0.305 0.35 0.291 0.245 0.324 0.326 0.425 0.3 0.28 0.33
第二部分 杆件的结构分析
杆件的结构静力分析分类
杆件分析主要见于大型钢结构中的分析,如果都使用 实体模型的话,模型将非常大。
杆系结构分为静定结构和超静定结构。凡是仅用静力平衡原理即可 求出结构的全部内力和反力时,称结构为静定结构;否则为超静定结构。 超静定结构可用力法、位移法或混合法等求解。在求得内力后,静定结 构和超静定结构均可用位移计算公式或其他方法求得结构中任意指定点 的位移。较复杂的超静定结构,由于其计算工作量很大,在20世纪30~ 50年代期间,曾发展了许多近似法、渐近法及实用的简化方法。这些方 法在当时曾解决过许多工程结构的计算问题,也推动了结构力学的发展。 但随着电子计算机的发展和普及,适合于计算机的矩阵力法、矩阵位移 法及有限元法等已成为分析复杂问题的主要方法。
有限元分析法概述
第十一章 有限元分析方法概述1、基本概念有限元分析方法是随着电子计算机的发展而迅速发展起来的一种现代没计计算方法。
它是20世纪50年代首先在连续体力学领域—飞机结构静、动态特性分析中应用的一种有效的数值分析方法,随后很快就广泛地应用于求解热传导、电磁场、流体力学等连续性问题。
在工程分析和科学研究中,常常会遇到大量的由常微分方程、偏微分方程及相应的边界条件描述的场问题,如位移场、应力场和温度场等问题。
求解这类场问题的方法主要有两种:用解析法求得精确解;用数值解法求其近似解。
应该指出,能用解析法求出精确解的只是方程性质比较简单且几何边界相当规则的少数问题。
而对于绝大多数问题,则很少能得出解析解。
这就需要研究它的数值解法,以求出近似解。
目前工程中实用的数值解法主要有三种:有限差分法、有限元法和边界元法。
其中,以有限元法通用性最好,解题效率高,目前在工程中的应用最为广泛。
下面通过一个具体例子,分别采用解析法和数值解法进行求解,从而体会一下有限元分析方法的含义及其相关的一些基本概念。
如下图所示为一变横截面杆,杆的一端固定,另一端承受负荷P ,试求杆沿长度方向任一截面的变形大小。
其中,杆的上边宽度为1w ,下边宽度为2w ,厚度为t ,长度为L ,杆的材料弹性模量为E 。
已知P =4450N ,1w =50mm ,2w =25mm ,t =3mm ,L =250mm ,E =72GPa 。
① 采用解析法精确求解假设杆任一横截面面积为)(y A ,其上平均应力为σ,应变为ε。
根据静力平衡条件有:0)(=-y A P σ根据虎克定律有:εσE =而任一横截面面积为:t y L w w w y A )()(121-+= 任一横截面产生的应变为:dydu=ε将上述方程代入静力平衡条件,进行变换后有:dy y EA Pdu )(=沿杆的长度方向对上式两边进行积分,可得:⎰⎰⎰-+==y yudy y Lw w w Et P dy y EA P du 01210)()(将)(y A 表达式代入上式,并对两边进行积分,得杆沿长度方向任一横截面的变形量:]ln )[ln()()(112112w y Lw w w w w Et PL y u --+-=当y 分别取0、62.5、125、187.5、250值时,变截面杆相应横截面处的沿杆长方向的变形量分别为:m u m u m u m u m u 6564636211080.142 ;1083.96 ;1027.59 ;1051.27 ;0----⨯=⨯=⨯=⨯==② 采用数值解法近似求解将变横截面杆沿长度方向分成独立的4小段,每一小段采用等截面直杆近似,等截面直杆的横截面面积为相应的变截面杆横截面面积的平均面积表示,每一小段称为一个单元,小段之间通过节点连接起来。
详解Catia静态有限元分析
Catia静态有限元分析指南注意:在进行有限元分析之前,必须赋予零件材质属性。
切换到GPS模块时出现的对话框说明如下:缺省情况下,CATIA会自动计算并为每个零件赋予网格特性。
网格特征可以删除和添加。
一、模型管理创建四面体网格,用于3D体单元网格划分。
创建2D面网格,用于面和板壳单元网格划分。
创建1D网格,用于线和梁单元网格划分。
修改局部网格大小,达到网格划分不同密度的需要。
修改网格类型,分为线性和非线性两种。
创建局部网格塌陷。
创建实体特性,缺省情况下,CATIA自动为part赋予实体特性。
创建壳单元特性。
创建梁单元特性,分为以下几种:圆柱,参数R。
管状,参数R i和R o。
矩形,参数H和L。
匣形,参数L i、L e、H i和H e。
U形梁,参数H、L和T。
I形梁,参数H、L、T l和T h。
T形梁,参数H、L、T h和T l。
X形梁,参数H、L、T h和T l。
用户自定义的梁。
输入梁的参数数值。
创建导入的梁特性。
检查模型,可以检查特性、连接和网格等方面,建议在进行计算之前进行模型的检查。
二、网格规范创建适应性框,来修改网格规格。
三、群组群组功能可以使你生成一组点、线、面和体的映像,方便操作。
群组点。
群组线。
群组面。
群组体。
四、连接特性创建滑动连接,在共同的接触面上,垂线方向上两个体扣紧,切线方向上可以相互滑动。
创建接触连接,防止体在彼此共同接触面上分离。
创建扣紧连接,使体在共同面上扣紧。
创建压力装配连接,防止体在彼此共同接触面上分离。
创建螺钉固定连接,防止体在彼此共同接触面上分离。
创建刚性连接,在体之间的共有边界上创建硬性的紧扣连接,表现就好像共有面见具有无穷的刚性。
创建柔性连接,在体之间的共有边界上创建紧扣连接,表现好像它们之间是柔软的。
创建虚拟刚性螺钉连接,只考虑使用螺钉装配式的拉紧压力,而不包括螺钉。
创建虚拟柔性螺钉连接,在一装配系统中指定边界作用。
自定义间隔连接,在一定的距离之内,指定单元的类型和关联特性。
连杆有限元分析
研究与试验连 杆 有 限 元 分 析包头职业技术学院 文占科 连杆是柴油机的主要零件之一。
它在柴油机中,把作用于活塞的膨胀气体压力传给曲轴,又受曲轴的驱动而带动活塞压缩气缸中的气体。
连杆在工作中承受着急剧变化的动载荷。
在连杆设计过程中,为了减小其尺寸,保证安全,故对它进行有限元分析,得出位移和应力分布结果,以便对设计方案的刚度和强度有一个较为准确的估计,从而合理地改进和优化设计方案,提高设计的效率和可靠性。
1 三维实体模型的建立及软硬件条件本文研究的连杆实体模型采用笛卡儿坐标系,全部采用基于特征的参数化建模技术进行三维实体建模。
连杆本身结构形状比较简单,为了能最大限度地使分析结果接近实际情况,在实体建模过程中,对连杆所作的简化很少,仅对小油孔等作必要的简化。
由于采用参数化特征建模,后续工作将能更方便而快捷地进行。
因为设计者既可对箱体连杆总体模型进行多层次的变结构、变参数改型设计,也可对各零部件模型进行低层次(也是更大自由度)的变结构、变参数改型设计,以便充分发挥设计者的创造性。
有限元分析全部内容在大型CAD/CAE/ CAM集成化软件I-DEAS Master Series5上进行,在微型计算机上完成。
其基本配置为: CPU主频为600MHz,内存为128MB,硬盘20GB。
2 基于实体施加边界条件为使分析结果最大限度地接近真实情况,所有边界条件的施加,全部采用通过实体施加的方式,由程序自动完成实体边界条件与等效有限元边界条件之间的转化。
基本参数:缸径D=110mm,转速n= 2600rpm,最大爆发压力P2=130bar,连杆两孔中心距L=210mm,连杆大端重量W3=1. 03kg,连杆小端重量W2=0.83kg,活塞组重量W1=2.25kg,曲轴半径R=66mm。
根据连杆的工作形式,本文选择最大受拉和最大受压两种工况进行计算。
连杆最大受拉工况出现在排气冲程的上止点处,此时连杆小端孔表面受活塞组W1的最大往复惯性力: P1=(1+R/L)×W1×R×(3.1416×N/ 30)2=14481N连杆大端孔所受的惯性力为:P M=P1+P2+P3=[(1+R/L)×(W1+ W2)+W3]×R×(3.1416×N/30)2=24858N 连杆最大受压工况以最高转速时的爆发压力来计算。
Catia静态有限元分析指南
Catia静态有限元分析指南注意:在进行有限元分析之前,必须赋予零件材质属性。
切换到GPS模块时出现的对话框说明如下:缺省情况下,CATIA会自动计算并为每个零件赋予网格特性。
网格特征可以删除和添加。
一、模型管理创建四面体网格,用于3D体单元网格划分。
创建2D面网格,用于面和板壳单元网格划分。
创建1D网格,用于线和梁单元网格划分。
修改局部网格大小,达到网格划分不同密度的需要。
修改网格类型,分为线性和非线性两种。
创建局部网格塌陷。
创建实体特性,缺省情况下,CATIA自动为part赋予实体特性。
创建壳单元特性。
创建梁单元特性,分为以下几种:圆柱,参数R。
管状,参数R i和R o。
矩形,参数H和L。
匣形,参数L i、L e、H i和H e。
U形梁,参数H、L和T。
I形梁,参数H、L、T l和T h。
T形梁,参数H、L、T h和T l。
X形梁,参数H、L、T h和T l。
用户自定义的梁。
输入梁的参数数值。
创建导入的梁特性。
检查模型,可以检查特性、连接和网格等方面,建议在进行计算之前进行模型的检查。
二、网格规范创建适应性框,来修改网格规格。
三、群组群组功能可以使你生成一组点、线、面和体的映像,方便操作。
群组点。
群组线。
群组面。
群组体。
四、连接特性创建滑动连接,在共同的接触面上,垂线方向上两个体扣紧,切线方向上可以相互滑动。
创建接触连接,防止体在彼此共同接触面上分离。
创建扣紧连接,使体在共同面上扣紧。
创建压力装配连接,防止体在彼此共同接触面上分离。
创建螺钉固定连接,防止体在彼此共同接触面上分离。
创建刚性连接,在体之间的共有边界上创建硬性的紧扣连接,表现就好像共有面见具有无穷的刚性。
创建柔性连接,在体之间的共有边界上创建紧扣连接,表现好像它们之间是柔软的。
创建虚拟刚性螺钉连接,只考虑使用螺钉装配式的拉紧压力,而不包括螺钉。
创建虚拟柔性螺钉连接,在一装配系统中指定边界作用。
自定义间隔连接,在一定的距离之内,指定单元的类型和关联特性。
基于有限元的汽车车架静态分析
J n.2 0 u 07
文 章 编 号 :10 —4 4 2 0 )2—0 1 —0 0 9 4 X(0 7 0 12 5
基 于 有 限元 的汽 车 车架 静 态分 析
刘新 田 ,黄 虎 刘 长虹 , ,郭 辉 ,范 平 清
( 上海 工 程 技 术 大 学 汽 车 工 程 学 院 , 海 2 1 2 ; . 一 一 研 究 所 , 海 2 0 9 ) 1 上 060 2七 上 00 0
相应部分转 角和角刚度 。通过分析 计算 , 可知该车 架的静应 力和扭 转刚度符合设 计要 求 。
关键 词 : 应 力 ; 转 刚度 ;有 限元 分析 静 扭
中图 分类号 : 4 3 3 U 6 .2 文 献标志 码 : A
St tc An lss o t m o r m e Ba e n t e FEA a i ay i fAu o bi F a s d o h l e
S c n l h r n e l we e f r e n e rwh es we e f e To s n si n s ft e fo t mi d e a d e o dy t efo twh es r o c d a d r a e l r i d. r i t f e s o h r n , d l n x o f
LI X i i n U nta , HU A NG u H ,LI Ch ng o g U a h n ,GU O u , FAN ng ng H i Pi qi
( . ol e f uo b eE gne i ,S a g a Unvri f n ier gS ine S a g a 2 1 2 C i ; 1 C l g tmo i n i r g h n h i i s yo gn e n c c , h n h i 0 6 0, hn e oA l e n e t E i e a
有限元静力学及动力学分析
1. 基本方程和术语
通用运动方程: 假定为自由振动并忽略阻尼:
M u Cu Ku F(t
M u Ku 0
假定为谐运动u = u0cos( t) :
(K 2M u0 0
1. 这个方程的根是 i, 即特征值, i 的范围从1到自由度的数目,特征
值的平方根是i , 它是结构的自然圆周频率(弧度/秒),并可得出
i j
(t) (t)
Fi (t)e [K ]e{ (t)}e
(3)整体分析
(4)通用运动方程
1.结构离散 与静力分析相同,选用适当的单元类型将连
续的弹性体离散成有限多个单元和节点。 2.单元分析
从离散的弹性体中任 意取出一个单元。利用给 定的位移插值方式表示单 元 内 任 一 点 的 位 移 {δ(t)}e, 进而确定节点的速度和加 速度。 3.整体分析 利用各节点处的变形协调条件和动力平衡条件即达朗贝尔原理,建立整体刚 度方程;
1.建模 2.加载求解 3.检查分析结果
2.1动力学有限元分析原理
1) 动力学分析的原因 2) 动力学有限元分析引例 3) 动力学分析的定义和目的 4) 动力学分析类型
1) 动力学分析的原因
静力分析也许能确保一个结
构可以承受稳定载荷的条件,但 这些还远远不够,尤其在载荷随 时间变化时更是如此。
无需 添加
选模态 设阶次 算结果
不添约束时,前6阶振型为刚体位移,固有频率均为0.
第三节 谐响应分析
1.定义和目的 2.术语和概念 3.谐响应分析求解方法 4.谐响应分析步骤 5. 实例-弹簧质量系统谐响应分析
1、定义与目的
1、定义与目的
1、引例-货车有阻尼受迫振动
在o< η <1之间, η大,即运行速度高,则受迫振动 振幅大;在η >1时, η 大,即运行速度高,则受迫振 动振幅小。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
这是最常用的分析类型。所有航空,汽车,海洋和土木工程行业会用线性静态分析。
HyperMesh线性静态分析举例
理论结果 在机械学里,静态状况意思就是平衡力和扭矩(V=0)。静态载荷应该没有变化。要是载荷变化很慢,结
构反应可能考虑用静态分析确定。但是载荷变化很快的话(相对于结构的反应能力)就需要用动态分析确定。 计算静态问题,所有有限元分析求解器会算下面的方程: K: 整体刚度矩阵 x: 需要确定反应的位移矢量 f: 施加到结构的外力矢量
静态 静态分析有两个条件: 1. 力是静态的,也就是说,力是常数,静负载
2. 平衡条件:∑ Force = 0, ∑ Moments = 0
∑Fx = 0
∑x = 0
∑Fy = 0
∑My = 0
∑Fz = 0
∑Mz = 0
有限元模性在每个节点应该满足这个条件。在每个位置,外部力的总和应该等于反作用力总和。
第四步- 定义约束载荷集,设置模型约束
1. 在 Model Browser 里右键点击然后选 Create > LoadCollector 打开 Create LoadCollector 窗口。 2. Name 输入 SPC,Card Image 选成 none,然后点 Create。
创建这个载荷集以后,在 Model Browser 可以看得到约束 SPC 在 LoadCollector 里。
5. 点击 Create。 6. 重复这个流程,这次节点 2 和 3 只选 dof2 和 dof3。模型现在是这样的:
现在模型只有两个自由度,Ux2 和 Ux3。 第五步 – 设定力载荷集,施加载荷 1. 创建载荷集 Force 选 no card image。
创建这个载荷集以后,可以在 Model Browser 里看到 Load Collector 组里的 Force。
2. 现在应该在节点 3 创建 X 负方向的力。用 BCs > Create > Forces 下面的 Forces 面板。 3. Nodes 选 node3。magnitude 输入-10,然后选 x-axis,具体如下:
4. 点击 Create 创建力。
第六步 – 定义载荷步。 1. 选 Setup > Create > LoadSteps 创建载荷步。 这样就会进入 LoadSteps 面板,这里可以选约束、载荷和分析类型。
8. 重复这个过程,创建第二部件 A2,选 PROD 属性,面积等于 78.54。也指定 steel 材料。记得把部件 选成另外颜色。
创建这两个部件以后,应该可以在 Model Browser 里看到两个叫 Component 和 Assembly Hierarchy 的新 组,两个里面都有 A1 和 A2。
11.2 HyperMesh设定线性静态分析
用 HyperMesh设定RADIOSS BulkData做静态分析基本过程如下: 1. 从 Model Browser,右键点击 Create > Material。
2. 然后输入材料名称,选 card image (mat 1) 再点击 Create。
3. 在 HyperView 的工具栏,点击 Contour, 后点击 Apply。这就会生成云图如下:
,选 Displacement 和 Mag 作为 Result Type, 然
4. 可以用 Measure 面板, ,比较这个结果和第一部分里面的结果。点击 Add 然后把 Measure Type 变成 Nodal Contour,然后选节点 2 和 3。
3. 然后需要输入材料杨氏模量 [E]和泊松比 [Nu],再点击 return.
4. 要是只有一维单元,泊松比可以不输入。 5. 现在应该可以看到在 Model Browser 里的 Material 组下面的 steel 材料。
第二步 创建新的部件和属性,然后指定材料 1. 从 Model Browser,右键点击选 Create > Component。 2. 然后输入部件名称:如 A1,注意目前没有指定的属性或者材料。
可以看得到,Measure 里的数值跟第一部分里的是一样的。
很多复杂问题可以简化为一个静态系统。也可以用静态分析先查一下非线性特性,是不是需要做更复杂的 分析类型。材料密度是不需要的,除非要看部件质量影响。做静态分析中需要指定材料的弹性模量和泊松比。 评估静态分析结果时,需要注意线性假设有效性。要考虑一下:
两个圆:Ri=50 mm; Ra=250 mm。
两圆中间创建填充表面。 然后沿表面法向拉 5mm,形成封闭的 3D 几何。
网格划分 用 Solid Map Mesh 面板创建实体单元。
确定圆盘中心和柱坐标系 模型里需要指定柱坐标系才能计算径向应力和切向应力(Sxx 和 Syy)。
可以用 Distance 面板确定圆盘中心。注意用三节点子面板。 在内圈选三个节点,点击 circle center. 重复上一步,确定另外表面中心。 用这两临时节点确定柱坐标系 z 方向,选 Geometry > Create > Systems > Node Reference。
11.4 项目:旋转圆盘分析
下面演示如何设定旋转圆盘线性静态分析。
图形:3D FEM 旋转圆盘模型。Ri=50 mm, Ra=250 mm, 厚度 T= 5 mm, 材料属性: E= 210.000 N/mm2, Rho=7.9e-9 t/mm3, 平均单元大小10 mm
几何 首先需要建圆盘模型,一个办法是创建两个圆:
相对于结构大小,变形和旋转够不够小? 超过材料屈服应力了没有? 部件与部件,部件与边界有没有穿透,这将改变系统载荷的路径。
11.3 线性静态分析教程和交互式视频
首先需要用邮箱地址在 HyperWorks 用户中心网站注册才能看下面的视频和教程。 推荐教程: •RD-1000: Linear Static Analysis of a Plate with a Hole 推荐视频 (10-15 分钟,不需要安装 HyperWorks): •Loads •Loadstep/Subcase •Input File •Log file, Error Messages •Direct Frequency Response Analysis of a Flat Plate •Direct Transient Dynamic Analysis of a Bracket •Convergence of finite elements – by Prof. J. Chessa, Texas 网络研讨会: •Linear Solutions with RADIOSS •NVH Solutions with RADIOSS
第七步 - 提交分析
1. 到 Analysis 页里,选 Radioss 面板。 2. 点击 save as… 然后给 Radioss 文件制定存储路径和后缀为.fem 的文件名。 3. 确认设置跟下面一样,然后选 Radioss 提交分析。
第八步 – 结果后处理
1. 计算完成以后,看结果最容易办法就是点击 Hyperview 绿色按钮。 2. 这就会打开 HyperView 并且加载结果.h3d 文件。
XI 线性静态分析
11.1 线性静态分析
应力
软件根据这个路劲进行线性 静态计算
实际应力-应变曲线
应变
线性 线性表示材料线性弹性行为。也就是说,应力-应变曲线的线性部分,一个直线按照胡克定律:σ=εE。这
条直线方程就是 y=mx跟原点有交叉。曲线的斜率就是弹性模量,E,也是常数。在真实情况下,过了屈服应 力,材料会按照非线性曲线走,可是求解器还是会按照线性曲线走。超过极限应力,部件应该破裂,可是线性 静态分析不会显示失败。它只会显示没有破坏的部件,在失败位置显示有很高的应力。也可能会看到不现实的 大变形。工程师需要比较最大应力和屈服应力或者极限应力,然后决定部件是安全还是失败的。所以工程师必 须考虑材料和载荷在线性静态情况是否可行。
3. 在我们例子里需要创建 CROD 单元。用 Element Types 面板 (Mesh > Assign > Element Type)将杆 单元设置成 CROD。
4. Create rod 面板可以从 Mesh > Create > 1D Elements > Rods 下拉菜单选择。 5. 在 Create rod 面板里,从节点 1 到节点 2 创建一个杆单元,放在 A1 部件里,指定 A1 属性(用 Model Browser 选 A1 部件当作 Current Component)。注意要先选属性,再选节点,确保指定合适的属性。 6. 重复第五步,让 A2 部件当 Current Component,然后从节点 2 到节点 3 创建杆单元,指定 A2 属性。
3. 选 BCs > Create > Constraints 创建单点约束。这些约束会固定不允许动的自由度。 4. 在 RADIOSS 里,杆单元是空间单元,有三个自由度 (Ux, Uy, 和 Uz). 为了复制第一部分简单模型,
我们需要消除所有 Uy,Uz 自由度和节点 1 上的 Ux 自由度。在 Constraints 面板里,选节点 1,然后只选 dof1, dof2,dof3。
我们现在考虑一个简单的静态例子来帮助我们理解静态分析。例子如下:
这是二级柱,两个均等部分,有不一样的横截面面积。 下面说的办法可以扩展到任何问题。先需要用节点和单元建模结构。明显这个问题至少要两个单元,每一 部分一个单元。 我们先来考虑下面的有限元模型解决这个线性静态问题。在这个例子中,我们将使用杆单元,如下图片:
第三步 创建有限元网格和合适属性 这一步取决于模型。可以从导入的 CAD 几何开始,或者导入另外 FEA 模型,或者这两办法并用。在我们