汽车方面的结构有限元分析.ppt

汽车结构有限元分析,主要参考资料,1王勖成，邵敏.有限单元法基本原理和数值方法M.清华大学出版社,20012周中坚，卢耀祖.机械与汽车结构的有限元分析M.同济大学出版社，19973朱伯芳.有限单元原理及其应用.中国水利水电出版社，19984钱伟长，叶开沅.弹性力学.科学出版社，1980.035杨橙，马力，王仲范.复合材料储能飞轮临界转速与极限转速研究J.中国机械工程，2003（18）6廖芳，马力，杨橙，何耀华.客车骨架强度失效分析J.客车技术与研究，2003（4）7唐少雄，马力，杨代华，柴苍修.内燃机凸轮机构摩擦学仿真设计建模研究J.内燃机工程，2002（6）8柴苍修，马力，王元良，李春光.农用运输车车架有限元强度分析J.拖拉机与农用运输车，2002（1）,主要讲述内容,有限元法概述弹性力学简介平面问题有限元法ANSYS的应用简介轴对称问题有限元法板壳问题有限元法实体问题有限元法组合结构有限元法有限元模态分析有限元疲劳分析汽车结构有限元分析实例,前言关于有限元法英文缩写FEM（FiniteElementMethod）应用中习惯称有限元分析是一种连续结构离散化数值计算方法上世纪五十年代由美国飞机工程师提出（1956年，Turner）FEM与CAECAE计算机辅助工程（ComputerAidedEngineering）CAE范围更广，还包含其它工程分析方法,概述,国内外发展应用简介方法提出后吸引大量数学力学家进行研究，至今理论相当成熟目前应用非常广泛国外发达国家有限元分析软件相当完善，商业化程度很高应用意愿和应用水平很高国内发展相对较晚研究应用水平发展很快没有公认的、广为流行的、高水平的、功能强大的商业化分析系统就全国范围来讲，汽车行业主动应用意愿与国外相比差距较大,车辆工程中有限元主要应用范围按学科分类弹性力学断裂力学塑性力学结构分析动力学流体力学热力学电磁学声学,几乎覆盖车辆工程中的各个学科，如果包含交叉学科和边缘学科，范围则更广,按汽车结构分析分类汽车结构强度和刚度分析静强度失效载荷产生超过屈服极限或强度极限的应力载荷可是真的静载，也可为动载峰值结构变形影响运动关系检查运动干涉变形过大失效问题世界著名案例：日本车制动油管,常见汽车零部件失效,静强度失效,疲劳失效,共振原因失效,汽车结构模态分析属动力学分析领域分析结构的固有特征固有频率固有振型模态阻尼模态刚度广范用于汽车结构动态特性设计广泛用于解决汽车结构振动噪声问题,系统固有特性,输入（激励）,输出（响应）,响应分析时间响应分析通常在时域内求系统的时间相应常用于疲劳分析、平顺性分析等频率响应和响应谱分析多用于随机振动,常参数线性系统脉冲响应函数频率响应函数传递函数,输入（激励）,输出（响应）,温度场计算包括对流、传导和辐射发动机温度场计算及散热问题保温车（冷藏车）隔热计算热应力分析噪声振动分析声场分析空腔共鸣制动噪声部件机械振动噪声汽车空气动力学分析降低空气阻尼、升力等通风换气改善车身表面流场发动机进气及排放,现在最高的汽车车速是多少？,你家新房空调放在什么地方最省电？,能否实现在车内后排两老板（或情侣）说话司机不易听见而司机说话前者容易听见,汽车结构碰撞、冲击分析碰撞安全性问题车身安全性设计碰撞事故模拟再现汽车零部件冲压成形成形部件设计模具设计分析开裂、起皱、回弹等问题（液）流场分析液化石油气罐车,乘车坐在什么地方（相对）最安全,汽车结构可靠性分析灵敏度分析和结构动态修改汽车结构优化设计按汽车结构分析的力学特征分类线性问题计算分析非线性问题的计算分析几何非线性材料非线性状态非线性接触问题,弹性力学基本概念基本假设假设物体是连续的假设物体是均质的假设物体是各向同性的假设物体是完全弹性的假设物体的位移和应变是微小的,第一章弹性力学简介,基本物理量外力分布力：体力、面力等集中力应力应变位移,x,z,y,几何方程几何方程：应变和位移的关系,x,y,回顾单向拉伸：,刚体位移：应变为零的位移,为积分常数,物理方程和弹性矩阵物理方程：描述应力和应变之间的关系弹性矩阵,对称,回顾单向拉伸虎克定理：,平衡方程边界条件力的边界条件位移边界条件,y,z,应力微元体平衡,*,弹性力学问题归纳（位移法）求解位移函数，它满足三个平衡方程六个几何方程六个物理方程力的边界条件和位移边界条件数学上是偏微分方程的边值问题全部点的位移集合反映结构的变形弹性结构上有无穷多个点，所以有无穷多个自由度由位移函数可求得应变由应变可求得应力,圣维南原理若把物体的一小部分边界上的面力变换为分布不同但静力等效的面力（主矢量相同，对于同一点的主矩也相同），那么近处的应力分布将有显著的改变，但远处所受的影响可忽略不计。,P,P,P,P,P,弹性体虚位移原理虚位移结构约束所允许的任何微小位移虚应变虚位移所产生的应变虚功原理外力在虚位移上所做的功等于弹性体内应力在相应的虚应变上所做的功,第二章有限元法基思想,引言弹性力学对结构分析的描述（求应力和变形）十五个偏微分方程满足边界条件连续体，无穷多个自由度机械结构复杂性几何形状复杂位移边界条件复杂所受载荷复杂结构材料复杂求不出解析解，须用数值计算方法,有限元基本方法结构离散化离散成有限个单元单元之间通过节点相连问题变成求节点位移有限个自由度网格越密精度越高计算时间越长所需计算机资源越大单元节点位移用矢量单元种类很多,结构,单元,节点,网格,单元分析假设位移模式得到假设的位移函数矩阵N任一点的位移用节点位移表示单元力学特性分析由几何方程得到应变矢量（B为应变矩阵）由物理方程得应力由虚位移原理得到单元刚度矩阵节点力和节点位移的关系,单元e,节点力,整体分析单刚迭加形成总刚K单元节点力矢量迭加形成结构节点力矢量形成结构线性方程组约束处理并求解方程处理位移约束条件求解线性方程组的全部节点位移根据所求节点位移计算应力分量,有限元应用实例（一）汽车安全气囊计算,有限元应用实例（二）动力响应1,有限元应用实例（三）动力响应2,有限元应用实例（四）接触问题,有限元应用实例（五）加工过程仿真,有限元应用实例（六）冲压成型1,有限元应用实例（七）冲压成型2,有限元应用实例（八）热轧,有限元应用实例（九）汽车碰撞1,有限元应用实例（十）汽车碰撞2,有限元应用实例（十一）失稳问题,有限元应用实例（十二）碰撞问题,有限元应用实例（十三）冲击1,有限元应用实例（十四）冲击2,有限元应用实例（十五）流体力学,有限元应用实例（十六）超弹性,有限元结构分类和计算步骤结构分类平面问题平面应力平面应变轴对称问题杆系结构桁架（平面、空间）刚架（平面、空间）板壳结构空间实体结构组合结构建立力学模型,计算过程前处理（建立计算模型）建立几何模型划分单元网格选择单元类型划分网格给定材料常数(E,u)给定单元实常数(厚度)施加载荷处理边界条件提交计算后处理观察分析位移结果观察分析各种应力其它结果,第三章平面问题的有限元法,平面应力问题网格划分单元分析整体分析举例及ANSYS应用初步平面应变问题,基本思想三维问题：十五个偏微分方程应力应变位移(先计算位移再计算应变和应力)二维问题平面应力问题：想法让与z有关的应力分量为零平面应变问题：想法让与z有关的应变分量为零位移分量,x,z,y,3.1平面应力问题,弹性力学平面应力问题方程平衡方程边界条件几何方程：三个偏微分方程物理方程：三个偏微分方程,平面应力问题特点（Planestress)结构特点载荷沿厚度方向不变化，其合力在中面内板的两面为自由面（无载荷作用）板关于中面对称板可以是变厚度的ta或tb,a,b,P,p,x,y,z,位移、应力和应变分量位移分量:应力分量应变分量,由虎克定律求得,平面应力问题的由来,物理方程和弹性矩阵,X,y,3.2单元网格划分,平面应力单元类型简介3节点三角形单元4节点4边形单元8节点4边形曲边单元节点位移分量每节点2个位移分量（自由度）x方向的位移u，y方向的位移v单元位移分量（4节点）,i,j,k,l,单元e,i,j,k,1,2,3,4,5,6,7,8,三角形单元,四边形单元,8节点单元,单元网格划分：生成单元节点信息应力梯度变化比较大的地方，网格应密一些有应力集中的地方，网格应密一些单元边界长度不要相差过大单元各边夹角不要太大集中载荷处要设置节点结构不同材料交界面处要设置节点并作为单元边界结构厚度突变处要设置节点并作为单元边界分布载荷突变处要设置节点施加位移约束处要设置节点注意单元间的连接,举例说明,材料A,材料B,界面,正确,这样不行,病态单元a-边长差别太大b-边长差别太大c-边夹角太大,a,b,c,单元节点信息节点信息单元拓扑信息,平板长2宽1,x,y,1,2,3,4,5,6,7,8,1,2,3,4,5,9,3.3三角形单元分析,单元的节点位移和节点力节点坐标:节点位移节点力目标：单元刚度矩阵K,单元位移模式单元位移模式的选取（线性模式）形函数和形函数矩阵形函数,代入,解得,整理,其中,(i=i,j,m),A为三角形面积,形函数矩阵形函数的性质在本节点值等于一,在它节点等于零权性,收敛性讨论单元内位移模式必须是连续的，公共边上位移必须协调位移模式必须反映单元的刚体位移位移模式必须反映单元的常应变可以证明三节点三角形单元是收敛的,i,j,m,p,i,j,m,p,单元刚度矩阵单元应变与节点位移的关系:应变矩阵,B为应变矩阵,单元应力与单元节点位移的关系：应力矩阵特点单元内各点的应力和应变是相等的三节点三角形单元为常应力（常应变）单元,S为应力矩阵,单元刚度矩阵及其特性单元刚度矩阵（单刚）,单刚的力学意义对角元素的力学意义为：使单元第i个自由度（位移分量）产生单位位移而其它位移分量均为零时需要在该自由度上所施加的力。非对角元素的力学意义为：使单元第j个自由度（位移分量）产生单位位移而其它位移分量均为零时第i个自由度上所产生的力。,(1,2),(3,4),(5,6),自由度编号,不动,1方向单位位移,这里（1方向）要加多大的力,(1,2),(3,4),(5,6),1方向单位位移,不动,这里（5方向）要加多大的力,单刚的特性单元刚度矩阵与位移模式有关单元刚度矩阵与单元形状、大小和方位有关单元刚度矩阵与单元的位置无关单元刚度矩阵是对称矩阵刚度矩阵是奇异矩阵,不存在逆矩阵,非节点载荷的移植（静力等效）集中力等效集中力最好在集中力处设置节点分布面力分布体力,i,j,m,x,y,P,i,j,m,q,单元迭加整体分析结构平衡方程：总体刚度矩阵（总刚）结构的节点力矢量:一维存储简介R存放力的矢量R（12）是第几号节点的哪个方向的力？第n号节点的Y方向的力是线性表R中的第几号元素？A存放节点坐标（三维）A（88）是第几号节点的哪个方向的坐标？第n号节点的坐标是线性表A中的第几号元素？总刚也是一维存储,单刚迭加形成总刚单刚为矩阵，子块迭加举例,1,2,3,4,5,6,1,2,3,4,5,12*12矩阵，6*6子块,总刚的特性矩阵元素的力学意义与单刚相同对角元数总是正值总体刚度矩阵为稀疏矩阵，非零元素均集中在对角元素附近减小相邻节点编号差值可减小带宽，节省存储单元总刚为对称矩阵总刚为奇异矩阵,带宽,对称,边界条件的处理划零置一法设已知边界条件为总刚的处理,对角元充1对应的行和列充0,充大数法已知边界条件为,M是很大的数，远大于其它元素如M=1.0E+30,应力计算及结果整理求解整体结构平衡方程的结构节点位移总刚用一维变带宽存储技术常采用波阵法求解根据节点位移求单元应力3节点三角形单元为常应力单元可以理解应力为单元中心处的应力结果的整理采用绕节点平均法求节点应力插值法求边界节点应力3点的应力可由2点和4点插值而得,1,2,3,4,1,3,2,4,目前国内常见的有限元系统专业有限元分析系统ANSYSADINANASTRANABQUSALGORSUPERSAPCAD软件挂带IDEAS软件中的有限元系统PRO/E软件中的有限元系统UG软件带的有限元系统,3.4平面应力问题举例及ANSYS应用初步,有限元软件评价要点单元库材料库算法库前后处理能力数据接口和数据转换标准与CAD软件UG、PRO/E、Parasolid、IGES有限元软件之间的数据接口IGES等数据转换标准操作方便性,ANSYS基本功能结构静力分析：弹性、塑性、蠕变、大变形、接触问题结构动力学分析：交变力、冲击或爆炸、随机力（地震）、其它瞬态力（如桥上的运动载荷）热分析：线性非线性热分析传导、对流、辐射电磁场分析：电感、电容、磁通量密度、涡流、电场分析、磁力分布、力、运动效应、电路和能量损失计算流体力学分析：瞬态或稳态声场分析压电分析：电子设备结构动态性能分析,ANSYS高级功能多物理场耦合分析：电磁热、压力结构等优化设计拓扑优化设计（外形优化）单元生死问题可扩展功能:连接用户自己的FORTRAN程序和子过程用户摩擦系数用户塑性屈服准则用户失效准则用户优化其它,ANSYS操作方式GUI方式：菜单、对话框操作命令方式：大约有1200多个命令程序方式：ADPL语言编程自动运行ANSYS几何建模点（keypoint）线面体,关键点Keypoint：几何建模用不参与有限元计算节点Node：有限元分析,解题过程及ANSYS应用有限元解题过程ANSYS初步应用,1000,500,P=500N,P=2.0e9Pa,E=2.1e11Pa泊松比0.3板厚30,计算举例（二）问题一方板,边长140mm,板厚10mm，板中心孔直径为20mm，两端受均匀拉伸分布力1.0E10Pa。材料弹性模量为we2.0E11Pa，泊松比为0.3。如图所示，计算结构应力和变形。,140,建立几何模型对称结构，只取其四分之一部分计算选择单元选平面问题4节点8自由度单元划分网格给定材料常数和单元厚度,施加载荷和约束提交计算,该线上各点X方向位移为零,该线上各点y方向位移为零,分布拉力,计算结果,VonMiss应力和结构变形原结构轮廓,局部应力放大VonMiss应力,平面应变结构特点(PlaneStrain)Z方向尺寸远大于x、y方向，横截面沿z轴不变化载荷平行于横截面，且沿z轴不变化任一横截面均可看成对称面（简化成平面问题）典型结构如大坝,x,y,z,y,x,3.5平面应变问题,位移、应力和应变分量位移分量:应变分量应力分量,由虎克定律求得,平面应变问题的由来,处理方法过程同平面应力计算时材料常数的处理理论公式上做如下变换程序应用中选择平面应力选项即可几何模型为结构的横截面,第三章轴对称问题的有限元法,结构特点几何结构绕轴线z对称（完整的旋转体）载荷绕结构对称约束绕轴线对称材料绕轴线对称力学特点Z轴横截面对称环向位移为零应力应变、位移只与r和z有关，与无关简化成平面问题,r,z,r,z,应力和应变分量应力分量：应变分量：单元节点位移单元类型举例（同平面应力问题）3节点、4节点、8节点平面问题单元等节点位移（同平面应力问题）每节点两个自由度，r和z方向位移载荷集中力分布面力体积力：重力、离心力,i,j,k,l,r,z,建模注意,机器不干,实心圆杆,空心圆筒,飞轮,考虑自重，是否轴对称问题,考虑自重，是否轴对称问题,计算举例问题一厚壁封闭容器，两端为半球形，中部为圆柱形，材料为普通碳素钢，其弹性模量为，泊松比为。已知圆柱段的长度为240mm，外径D100mm，内径d60mm。该容器以的转速绕其轴线旋转，容器内壁受的均匀内压。计算该容器的应力分布及变形。分析典型的轴对称问题可利用结构的对称性载荷包括内压和离心力,建立几何模型利用对称性，只取截面的12部分划分网格选择8节点单元给定单元材料常数质量密度查表得出,施加载荷和约束内压载荷离心力通过指定转速来施加,此线上各点r方向的位移为零,此线上各点z方向位移为零,内压,计算结果,VonMiss应力云图,变形图虚线为原结构,变形图网格线为变形后结构,第四章板壳问题有限元法,结构特点薄壁构件，t远小于结构边长受全方位载荷车身为典型的板壳结构有限元要点（通常情况）几何模型为板壳中面（t/2处）的形状几何模型无厚度单元和节点均在中面上最大应力发生在结构的上下表面,t,基本假设板壳中面法线在板壳变形后仍为直线，且垂直与变形后的中面板壳中面只有薄膜应力，弯曲应力为零板壳的上下表面上的应力为弯曲应力与中面薄膜应力之和弯曲应力沿截面线性变化,应力应变分量应力分量应变分量单元类型举例3节点、4节点单元8节点曲面单元单元自由度节点自由度（位移分量）单元自由度（8节点单元48个自由度）,1,5,2,6,3,7,4,8,载荷节点力矢量集中力、力矩分布面力分布体力边界条件简支固支已知位移或转角所需常数弹性模量泊松比板厚,固支,简支,板壳问题ANSYS应用举例一圆柱面顶盖薄壳，壳的两边支承在隔墙上，作为简支边。圆柱半径为1米，边长2米，圆心角90度。弹性模量为2.1e11,泊松比0.3，盖顶承受1e4Pa的法向均布压力。求盖的变形与应力。,边长两米简支,边长两米简支,板厚5mm盖顶受法向均布载荷,半径1米圆心角90度,计算举例（二）问题横截面为槽形的悬臂梁如图1所示。梁长1m，板厚5mm,其弹性模量为，泊松比，分布载荷的分布长度为500mm，集度为q=，集中载荷为p=2000N（作用点在上缘中点处），试计算结构的变形与应力。,建立几何模型当板壳板壳问题计算，几何模型为结构中面网格划分选择单元：8节点板壳原（shell）划分网格给定材料和单元常数,施加载荷与约束,集中力,分布载荷,固定端所有节点的六个自由度均被约束,计算结果,最大VonMiss应力处,VonMiss应力云图局部放大,结构变形图与VonMiss应力云图,原结构,变形后结构,第五章空间实体问题有限元法,结构特点3D实体有限元要点应力分量应变分量单元节点位移,x,y,z,单元举例8节点六面体单元20节点六面体等参元单元Tet10节点单元节点自由度为3X方向的位移uY方向的位移vZ方向的位移w20节点单元自由度为60所需材料常数弹性模量和泊松比载荷集中载荷分布面力分布体力,六面体单元,等参数单元等参数单元基本概念8节点平面问题等参数单元,1,2,3,4,5,6,7,8,x,y,1,2,3,4,5,6,7,8,1点坐标(-1,-1),4点坐标(-1,1),3点坐标(1,1),8点坐标(-1,0),变换,基本单元（母单元）,实际单元,设则可有位移变换和坐标变换其中形函数为,变换形式和参数完全一样,20节点空间等参单元,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,x,y,z,母单元,实际单元,设则可有位移变换和坐标变换其中形函数为,等参单元中的Jacobi矩阵和高斯积分Jacobi矩阵,J为jacobi矩阵,高斯积分,一维积分N个点,二维积分个点,三维积分类推：个积分点,还有4点积分(n=4)、5点积分(n=5)等,有限元网格自动生成的参数变换法方法原理简介几何结构的分块任何一个形状复杂的平面结构，都可以看成由若干个曲边四边形子块组成任何一个复杂的板壳结构，可以一些曲边四边形薄壳子块组成任何一个空间实体，也可以看作由一些曲面六面体子块组成,1,2,3,4,5,6,7,8,可以是真正具有四条边的四边形，也可以是人为定义的曲边四边形,每一个四边形（六面体）都可以通过前面介绍的参数变换方法由正方形（正六面体）的母单元变换得到对参数变化范围为-1,1的母单元划分网格，经过变换后即可得到子块的网格所有子块的网格划分完成后，即可得到整个结构的有限元网格,x,y,子块网格划分计算方法平面子块(3节点或4节点单元）用NX表示在边沿方向应划分的单元数用NY表示在边沿方向应划分的单元数用权因子控制方向所划分单元大小变化用权因子控制方向所划分单元大小变化子块内第K个节点的坐标可由下式计算其中初值,其中空间板壳子块,其它公式和平面子块一样,还可以生成8节点单元，只须增加单元边中间节点的计算即可,同理可扩充到3D六面体子块网格的划分，在此不再赘述,结构有限元网格的形成每个子块网格划分完成后，将其拼集在一起便可得到结构的有限元网格。由于每个子块都是单独划分的，在相邻子块的交线（面）上会出现一些公共节点，它们坐标一样但编号却不同，即产生一些虚假节点，必须将其消去，否则就无法用有限元程序进行正确的计算。设任意两个节点和，对给定的精度有：则认为节点和是同一个节点，必须将其中的一个节点消去。多余节点消除后，修改所有单元和节点编号即可得到结构的有限元网格。,子块A,子块B,实体问题计算举例问题图示U形夹左端固定，圆孔下半部受分布压力作用,，,图中长度单位为cm。用有限元法分析变形及应力。,建立几何模型由于对称性，只取结构的12分析即可网格划分单元选择选用Tet10节点单元划分网格给定材料常数,施加载荷和约束,对称面上各点垂直于对称面的位移为零施加约束,端面为固支，其上各点的所有自由度都必须约束住,分布载荷,计算结果,虚线为原结构,变形结构和VonMiss应力云图,原孔处的局部放大效果,第六章杆系结构有限元法,杆系结构桁架平面桁架空间桁架杆件与杆件间为铰接铰接点只传递力而不传递转矩每根杆件均为二力杆杆件不产生弯曲变形和弯曲应力有限元计算采用杆元（杆单元：bar）,桁架结构,刚架平面刚架空间刚架杆件与杆件间可理解为焊接连接点可传递力也可传递转矩刚架有限元分析采用梁元（beam）可当作刚架的常见结构高压线塔客车车身骨架管式摩托车车架自行车车架长江大桥,刚架结构,平面刚架的有限元法平面刚架单元（Beam）2节点单元局部坐标为单元坐标单元建立在杆中心线上单元为数学意义上的线单元节点位移单元节点力节点力与节点位移的关系,局部坐标系下的单元刚度矩阵可由材料力学直接推出也可假设位移模式和平面问题步骤一样推出,截面面积,弹性模量,对Z轴主惯性矩,单元长度,单元长度由机器根据两个节点的坐标自动计算,单元节点载荷节点集中力节点集中力矩单元分布力单元分布力矩重力,分布弯曲力矩,分布横向力,计算机可处理的线性分布可变成另外两种种具体形式1.均布分布力2.三角形分布,还可有分布轴力,单元刚度矩阵的坐标变换推导单元刚度矩阵时采用的是单元局部坐标系，它的坐标方向是由单元方向确定的，采用这样的坐标系，对所有单元可得到统一形式的单元刚度矩阵。但实际结构的每一杆件的方位都不一定相同，因此要将局部坐标系下的单元刚度矩阵和节点力变换到结构整体坐标系下才能迭加求解。设局部坐标系下的各量为：整体坐标系下对应的各量为：则有：,如果以为局部坐标与整体坐标之间的转换矩阵，则有：由此可得：即：其中,整体坐标系下的单元刚度矩阵,分别为两局部坐标轴在整体坐标系中的方向余旋,和,空间刚架的有限元法2节点空间梁元（Beam）2个节点i和j定义单元增加辅助节点k单元建立在杆中心线上单元为数学意义上的线单元节点位移单元节点力每节点6个自由度，为12自由度单元,单元刚度矩阵,对称,截面面积,弹性模量,剪切唐性模量,扭转惯性矩,对Z轴主惯性矩,对y轴主惯性矩,单元长度,单元长度由机器根据两个节点的坐标自动计算如果要考虑剪切的影响切，则还要给出对y和z轴的剪切影响系数以及沿y和z轴方向的有效抗剪面积共四个常数,单元节点载荷,单元节点载荷节点集中力节点集中力矩单元分布力单元分布力矩重力,分布弯曲力矩,分布横向力,计算机可处理的线性分布可变成另外两种种具体形式1.均布分布力2.三角形分布,还有分布力矩、,还可有分布轴力,单元矩阵的坐标变换,分别为三局部坐标轴在整体坐标系中的方向余旋,单元局部坐标系的确定辅助点k与节点i和j可确定局部坐标系节点i和j可确定轴方向K取在平面内由可得轴方向由和轴即可得到轴,结构整体坐标系,单元局部坐标系,辅助节点k的选取可以利用结构上已有的节点作为k点若k点不是结构上的点，其六个自由度须全部施加零约束k点不能在ij点的连线上,辅助点（又叫参考点）k点决定了杆件的放置方位,结构整体坐标系,如果给计算程序输入的对y轴和z轴主惯性矩大小（顺序）不变的话，下面k点选取的方位不同，则表示结构杆件放置的方位不一样。千万注意不要出错。,ANSYS中可处理的标准截面形状和参数（一）,ANSYS中可处理的标准截面形状和参数（二）,第七章组合问题的有限元法,热应力问题考虑热应力的物理方程其中为温度变化引起的应变式中为材料线膨胀系数，为温度的变化,考虑热应力时单元的虚功节点力与节点位移关系,相当于温度变化而产生的等效节点力，称为热载荷，在有限元计算中给定节点温度即可,预应力问题（初始应力）预应力钢板弹簧复合材料飞轮的预应力客车车身预应力蒙皮预应力钢箱桥梁有限元处理方法实质上与温度应力相同ANSYS可以方便处理预应力也可以用设置节点温度的方法获得预应力,杆件与块件的混合结构如图所示结构，由A和B两个实体用水平连杆相互连接而成，目的是要该结构的两部分在水平方向成为整体，具有较大刚度，而在垂直方向又相互独立，以消除不均匀沉陷所引起的相互影响。,杆单元节点三个自由度，实体单元节点也是三个自由度，问题可方便而自然地解决,支撑桩也可简化成杆元,连杆为杆元,实体单元,梁元和板壳元组合结构加筋薄壳结构问题结构以薄壳的刚度为主肋条只起加强筋的作用壳体与加强筋刚性连接加强筋截面某主方向与连接处处壳体中面法向一致i、j节点自由度为线性关系,i,j,薄壳采用板壳单元,加强筋采用梁元,梁元节点j当作板壳元节点i的从属节点,e,i点和j点的位移关系（约束方程）i点（主节点）的位移矢量j点（从节点）的位移矢量位移关系,变换矩阵,一般情况下的梁元和板壳元的组合梁单元节点为六个自由度板壳元节点也为六个自由度一般情况下可很协调地组合起来不同节点自由度的单元组合问题实体单元与板壳单元,实体，采用3D实体单元,板壳，采用板壳单元,地面上所有节点的全部自由度都被约束,按板壳组合结构算，如果不采取特殊的结构处理，计算伏明霞跳水，会发生什么情况呢？,x,y,z,实体单元、杆元与梁元,移动吊重,实体基础实体元,实体基础实体元,塔杆刚架梁元,绳索桁架杆元,一定要考虑地基、塔杆和绳索的应力和变形，结构该如何处理呢,塔杆刚架梁元,典型的特殊单元弹簧元两个节点每节点3个位移自由度有刚度和阻尼两个参数接触单元点接触单元面接触单元两物体相互靠近时接触时接触单元起作用（传递力）两物体相互脱离时单元不起作用（不传递力）绳索元两物体相互靠近时接触时单元不起作用（不传递力）两物体相互脱离时单元起作用（传递力）,自定义单元刚度矩阵读入单元刚度矩阵约束方程ANSYS节点间的约束方程简介ANSYS单元库ANSYS单元库简介,练习一边长为一米的方桌如图1(a)所示。四根桌腿为空心圆管，外径为50mm，壁厚为2mm；桌面为3mm厚的钢板（离地高度为1米），四边为矩形加强管，横截面高60mm，宽50mm,壁厚为2.5mm，其放置方位如图1(b)所示。桌面中心位置处有一钢块（长宽高：400-400-200mm），它由四个角点处的四个弹簧支撑在桌面上，钢块中心处作用有一垂直向下的集中载荷P=1000N。钢块的密度为7.8E1000kg/m3，弹簧的初始长度为60mm，刚度为100N/mm。材料弹性模量和泊松比分别为和。试计算结构的变形、内力与应力。,加强管,重钢块,桌面钢板,P,弹簧,钢架桌,第六章有限元动力学分析,有限元动力学方程和质量矩阵动力学方程M结构的质量矩阵C阻尼矩阵K刚度矩阵,质量矩阵结构质量矩阵由单元质量矩阵迭加而成一致质量矩阵对平面问题三节点单元,采用一致质量矩阵：1.可以得到更精确的振形2.计算所得频率值是结构真实频率的上界3.计算比较复杂，存储量较大,团聚质量矩阵根据精力学平行力分解原则求得等厚度平面三角形单元的团聚质量矩阵,1.计算比较简单2.存储量少3.有使结构频率计算值降低的趋势4.采用协调元时，由于协调单元有较高的刚度，会使结构计频率计算值偏高5.这两种相反因素相抵，有时可得到较好的固有频率计算值,模态分析（特征值问题）无阻尼自由振动解的形式广义特征值问题,广义特征值问题的解固有频率（n阶，n为结构自由度数）振型（特征矢量，n阶，n为结构自由度数）特征矢量的正交性一般情况下，工程上只对前几阶固有频率和振型感兴趣,特征值问题常见求解方法广义雅可比方法通过一系列变换将质量矩阵和刚度矩阵对角化后求解此法为求全部特征值的方法只求少数低阶特征值时，此法不合算逆迭代法通过一系列迭代求解特征值一般只用来求前三五阶特征对求高阶特征对时精度难以满足（精度较低）子空间迭代法求大型结构前几阶特征对很有效,计算举例问题一悬壁薄板长300mm，宽200mm，厚2mm，沿短边约束。，。求前四阶固有频率和振型。分析典型的板壳问题,几何建模划分网格采用8节点shell单元施加约束给定材料常数,计算结果:16.171Hz64.49Hz119.07Hz217.49Hz,第一阶16.17Hz,第二阶64.49Hz,第三阶119.07Hz,第四阶217.49,动力响应问题阻尼矩阵知道结构的各阶特征对和阻尼比后，阻尼矩阵为瑞雷阻尼矩阵,第i阶振型,第i阶固有频率,可由实验测得,只要两个阻尼比即可求得,响应问题常见求解方法振型迭加法此法适用于外力只激起较少振型的问题（如地震）或所需计算响应时间较长的情况逐步积分法适用于激发较多的振型（如冲击）适用于计算响应时间短促的情况ANSYS求时间响应问题谐振动瞬态响应,第七章疲劳问题有限元法,影响疲劳寿命的主要因素一般因素,表71影响疲劳强度的因素,经常碰到的四个影响因素,1）应力集中的影响应力集中对疲劳强度有显著影响，常用疲劳缺口系数来描述其影响；2）尺寸的影响试件尺寸越大疲劳强度就越低。通常标准试件比实际零部件的尺寸小，因此疲劳尺寸系数在疲劳分析中必须加以考虑。3）表面状态的影响疲劳裂纹通常产生于试件表面，因此，零部件的表面状态对其疲劳寿命有显著影响，这种影响程度用表面敏感系数来描述。4）载荷的影响绝大多数材料的疲劳强度是用标准试件在对称循环正弦波加载情况下得到的，而实际零部件所受载荷却十分复杂，载荷的种类、载荷的频率、平均应力、载荷波形以及载荷中间停歇和持续等因素对结构疲劳寿命均有影响。,确定疲劳寿命的方法,疲劳分析方法都包括三部分内容材料疲劳行为的描述；循环载荷下结构的响应；疲劳累积损伤法则,金属结构材料的SN曲线,断裂极限疲劳极限通常取值,SN曲线可用幂函数表示为,Miner线性累积损伤理论,目前有多种疲劳累积损伤法则，但对接触网零部件疲劳分析简单有效的是Miner线性累积损伤理论，这也是目前工程上广泛采用的理论，该理论在许多国际流行的专业有限元分析系统（如Ansys）中均可直接采用。根据Miner线性累积疲劳损伤理论，在大于疲劳极限而小于静拉伸对应的疲劳强度的应力的作用下，其循环次数达到时，零件的损伤程度为。设应力分为z级，则在各级应力的作用下的累积损伤度为：式中D表示疲劳损伤度。