有限元基本概念

有限元基本概念第1页，课件共140页，创作于2023年2月有限元方法的分类有限元位移法的基本概念依据求解问题的路径不同，有限元方法大致可分为：位移法：以位移为基本未知量力法：应力为基本未知量混合法：部分以位移；部分以应力为基本未知量①将连续的弹性体离散将无限自由度问题简化为有限的自由度问题②选择表示一个单元内连续位移的试函数，并根据插值理论将单元内的位移用节点位移表征，也可以说通过一定的分析，将试函数中的待定参数变为单元节点位移，显然如果求出节点位移就可以得到单元内部各点的位移进而求出单元内各点的应力及应变。第2页，课件共140页，创作于2023年2月③利用李兹法建立基于节点位移试函数下的结构势能并利用泛函求极值的方法，即最小势能原理建立关于节点力与节点位移间的关系。④利用节点力平衡的条件将所有的单元集合成总的港督矩阵，形成代数方程组。实际上就是建立单元刚度矩阵：总刚度矩阵：总节点位移列阵：总节点载荷向量：第3页，课件共140页，创作于2023年2月⑦将求得的节点位移带入前面选择的单元位移函数中，求得单元内各点的位移、应变及应力。⑤根据边界条件进行约束处理⑥求解代数方程组得到节点位移。在有限元中：不同的问题节点位移自由度不同，对应节点力不同，单元的刚度矩阵不同；同样问题位移函数不同，单元刚度矩阵也不相同。第4页，课件共140页，创作于2023年2月常见的单元类型①一维单元单元的形状是一条直线，每个单元有两个节点（在两端），或者内部设置节点的单元。一般对应杆系结构，或柔索结构。杆系结构：单一拉压杆件、两个方向的弯曲杆件、弯曲及扭转杆件等；柔索结构：考虑轴向拉伸具有较小抗弯及抗扭刚度的结构。第5页，课件共140页，创作于2023年2月②二维单元二维单元的形状及种类有：三角形、矩形及任意四边形。每种形状的单元节点数可以是各边交点形成的节点，也可以在各边中或单元内部增加节点，以便提高计算的精度。三节点六节点七节点可用于：平面问题、薄板弯曲问题等第6页，课件共140页，创作于2023年2月矩形单元任意四边形六节点曲边三角形等参元八节点曲边四边形等参元第7页，课件共140页，创作于2023年2月③三维单元一般用于空间实体问题常见最常见的空间实体单元有：四节点四面体单元、八节点六面体单元同样为了提高精度，可以在每个棱边上或单元内增设节点，成为精度更高的单元。四节点四面体单元八节点长方体单元第8页，课件共140页，创作于2023年2月曲面四面体单元二十节点曲面六面体单元三角形截面四边形截面第9页，课件共140页，创作于2023年2月总结：位移有限元理论求解的基本思路①假设单元的位移函数；将单元内部位移函数用节点位移表示；将单元内形变用节点位移表示；将单元内应力用节点位移表示；③应用李兹法求出单元刚度矩阵②单元内部各物理量用节点位移表示第10页，课件共140页，创作于2023年2月④利用节点力的平衡方程式得总刚度矩阵⑤利用虚功原理将单元内及边界上的分布力等效至节点上。⑥求解第11页，课件共140页，创作于2023年2月设定单元位移函数时常使用的插值多项式①插值函数使用的原因用多项式在一些点上的已知函数值，插值近似表示该多项的方法，称之为插值函数方法。而常用代数多项式作为近似连续函数，称此多项式为插值多项式。②多项式插值方法的分类分两类：a.只要求插值多项式本身在插值点上取已知的函数值；b.不仅要求插值多项式本身在插值点上取已知的函数值；而且要求插值多项式在插值点上取已知的函数值导数值。拉格朗日插值埃尔米特插值第12页，课件共140页，创作于2023年2月③拉格朗日的插值方法为了避免求解n+1阶线性代数方程组的麻烦，采用以下插值方法：构造一个特殊的插值多项式：其中：为n次多项式显然：此时满足条件：且的线性组合仍然是n次多项式，可见就是要求的第13页，课件共140页，创作于2023年2月为保证条件成立，得：可令：根据上述条件拉格朗日插值多项式第14页，课件共140页，创作于2023年2月拉格朗日插值函数可表示为：第15页，课件共140页，创作于2023年2月④艾尔米特插值拉格朗日插值的缺点：插值函数在基点的导数要通过插值函数的表达式来进行计算，而一阶艾尔米特插值则把基点的导数直接取为插值的条件。设函数在n个互不相同点上的取值为：导数取值为：寻找一个插值多项式使它在所有的插值点上均满足：第16页，课件共140页，创作于2023年2月艾尔米特插值构造原理由于插值点上给定共2n个值，因此可以确定2n-1多项式与拉格朗日多项式求解方法一样，可以用函数在n个基点上的值及导数值与2n个多项式线性组合而成。其中：第17页，课件共140页，创作于2023年2月经推导：其中：第18页，课件共140页，创作于2023年2月单元位移的一般表达形式：选定单元的类型和位移函数之后，我们可将单元位移表示为以单元节点位移为插值点的插值函数，并用矩阵表示之。单元中任意点的位移：其中：——形函数——单元节点位移弹性体位移有限元的基本表达关键问题有两个：位移函数的选取、单元刚度矩阵的求解第19页，课件共140页，创作于2023年2月其中：几何矩阵的一般表达形式：第20页，课件共140页，创作于2023年2月根据物理方程——应力与应变的关系：其中：——初应力向量；——材料的弹性矩阵；第21页，课件共140页，创作于2023年2月刚度方程的基本表达形式——推导思路①写出用矩阵形式表达的弹性体单位体积总势能；②写出整个单元的总势能，进而得出整个结构总势能；③根据变分原理，总势能一阶变分=0；得出单元刚度矩阵单元刚度矩阵求解的基本形式第22页，课件共140页，创作于2023年2月体积力面积力直接作用于节点集中力以及作用于单元上的集中力外部载荷向量的分析：一般结构受到的力：分别表示为：为此可采用下式表达力函数：第23页，课件共140页，创作于2023年2月体积应变能一般表达：单元的总势能一般表达：第24页，课件共140页，创作于2023年2月M个单元的总势能：考虑：第25页，课件共140页，创作于2023年2月根据变分原理，一阶变分等于零结构总势能：结构总的力函数可表示为：第26页，课件共140页，创作于2023年2月得到下列线性代数方程组：单元的总刚度矩阵由单元刚度矩阵组成：等效节点力：单元初始应力对应的等效节点力体积力引起的等效节点力面积力引起的等效节点力集中引起的等效节点力第27页，课件共140页，创作于2023年2月上述结构一阶变分=0得到的线性代数方程组又可表示为：其中等效节点力：对于集中作用外力我们可详细的将其分成两种：

①直接作用于节点上的外力,可在方程中加入该力

②作用于单元内部的外力——此时我们也要求出其相应的等效节点力，该等效力具有单元的意义。第28页，课件共140页，创作于2023年2月值得注意的是：

节点力的等效处理过程中：我们使用了形函数，由于形函数是我们假设的近似函数，所以等效的过程是一个近似计算的过程。第29页，课件共140页，创作于2023年2月边界条件的处理边界条件处理的必要性：总刚度矩阵为奇异矩阵，要将边界的条件代入消除奇异性。从位移边界而言一般有几种类型的边界条件：①指定某些节点具有不变的确定性位移；②某些点上的位移=零③某节点为弹性边界——节点的反力与节点位移具有线性或非线性的关系。第30页，课件共140页，创作于2023年2月数学处理上一般采用的数学手段：当某一位移等于零时，如第31页，课件共140页，创作于2023年2月当某一位移为指定位移时，如可采用乘大数方法第32页，课件共140页，创作于2023年2月当某一点弹性支持时，可采用直接修正对应节点所在行的对角线上的刚度系数方法。如在i节点上有弹性约束：对应支座反力与位移的关系为：第33页，课件共140页，创作于2023年2月空间梁单元的节点位移、节点力、位移函数、及刚度矩阵。（1）空间杆单元单元位移：u表征杆件应变的物理量为节点位移个数：每个节点6个位移个数，节点力的个数：每个节点6个节点力。两节点单元共有12个自由度第34页，课件共140页，创作于2023年2月第35页，课件共140页，创作于2023年2月（2）空间梁单元空间梁单元的位移：位移函数的选取：且12个待定参数，利用节点位移进行插值求解第36页，课件共140页，创作于2023年2月表征其变形的几何方程特性的物理量：第37页，课件共140页，创作于2023年2月单元刚度矩阵：第38页，课件共140页，创作于2023年2月不考虑横向剪切影响：第39页，课件共140页，创作于2023年2月第40页，课件共140页，创作于2023年2月第41页，课件共140页，创作于2023年2月（3）平面问题的有限元节点位移与节点力oyxmijvmumvjujviui第42页，课件共140页，创作于2023年2月位移函数式中：将三节点i,j,m坐标代入（1）式：将单元内部位移用节点位移表示之—广义坐标系第43页，课件共140页，创作于2023年2月

简记为：其中：（3）第44页，课件共140页，创作于2023年2月位移矩阵为三角形i,j,m面积（4）得节点位移与插值函数的乘积之和第45页，课件共140页，创作于2023年2月——形函数矩阵或位移矩阵第46页，课件共140页，创作于2023年2月单元应变(几何矩阵）用弹性理论平面应力问题的几何方程式，可得单元应变：或：式中：几何矩阵用节点位移表达的单元应变第47页，课件共140页，创作于2023年2月单元应变（用节点位移表示的单元应力）根据虎克定律的矩阵表示式应力矩阵第48页，课件共140页，创作于2023年2月单元刚度矩阵（表示节点位移与节点力关系矩阵）求解单元刚度矩阵的方法：虚功原理基本公式：给节点虚位移：真实应变：即：节点力在虚位移上所作的虚功=虚位移引起单元内部的虚应变能单元上真实的节点位移：真实应力：相应的虚应变：节点力在虚位移上所作的虚功：虚位移引起单元内部的虚应变能第49页，课件共140页，创作于2023年2月进行比较得：式中单元刚度矩阵：将[B]、[D]代入上式：分割子矩阵对称、奇异、稀疏矩阵第50页，课件共140页，创作于2023年2月§8-5结构刚度矩阵

本节通过单元节点的力的平衡关系来建立结构的平衡式，包括结构刚度矩阵的建立。取出节点i，列出x,y方向力的平衡方程式：(1)(2)(1)(2)immjjmnnii第51页，课件共140页，创作于2023年2月该结构共有两个单元，外力只作用于i节点之上。对于其它节点同样可列出相应的方程式。将这些方程式合并一齐用矩阵表达，形成整个结构的节点力平衡方程。其形式如下：外力列阵，每一个节点有2行。应包括：直接作用在节点上的外力、支座反力及等效节点力。各单元因节点发生可能位移而产生的节点力之合。可由各单元刚度矩阵依对号入座方式形成其中：第52页，课件共140页，创作于2023年2月式中[Kij]为单元刚度阵的子矩阵，上式可简记为：[K]—结构总刚度矩阵具有n个节点的结构,总节点力平衡方程式为：注意：总刚度矩阵具有与上章所述相同的性质。对称性、稀疏性、奇异性。行列个数2n第53页，课件共140页，创作于2023年2月

§8-6外载荷处理单元中分布力的移置单元上受到的体积力与面积力等效到节点上力的计算（1）单元上体积力的等效计算计算原则虚功原理。当两个力系在同一个可能发生的虚位移上所做的虚功相等时，这两个力系静力等效。设：三角形平面单元受均匀分布力，其合力Q作用于形心处第54页，课件共140页，创作于2023年2月外力作用下：给节点i单位的虚位移：有：第55页，课件共140页，创作于2023年2月单元位移函数及位移函数的选取①单元位移函数坐标系两种坐标系：自然坐标系、广义坐标系自然坐标系：指无因次的，坐标变化范围在±1之间的坐标系。又称面积坐标广义坐标系：一般笛卡尔坐标系、极坐标系柱坐标系。自然坐标系与广义坐标系的关系：可从广义坐标系下的位移推得自然坐标系下的位移。第56页，课件共140页，创作于2023年2月以平面问题为例加以说明：广义坐标系下，二维单元的位移函数可表示为：为了保持几何不变性，选取多项式时要保持其完整性。参照帕斯卡三角形选取多项式。第57页，课件共140页，创作于2023年2月帕斯卡三角形第58页，课件共140页，创作于2023年2月多项式的项数——广义坐标的数目，取决于单元的自由度，一般是相等的。自然坐标下，二维单元的位移函数可直接表示为：节点位移与插值函数的乘积之和。形函数：由于插值函数反映了单元位移的形态所以又被称为形函数。第59页，课件共140页，创作于2023年2月——单元的节点位移向量；——形函数矩阵；n——单元节点个数；第60页，课件共140页，创作于2023年2月单元位移函数收敛必须满足的三个条件①连续性相邻单元无分离、无重叠②完整性位移函数含有：常位移项及常应变项③几何各向同性单元位移函数不因坐标变化而发生改变也就是保持对坐标的对称性。协调元与完备元的概念第61页，课件共140页，创作于2023年2月三角形单元是什么类型的单元？(1)(2)ijmnijmp矩形单元的位移函数第62页，课件共140页，创作于2023年2月将[B]、[D]代入上式：第63页，课件共140页，创作于2023年2月ijmp单元类型的讨论局部坐标转换问题与协调性问题第64页，课件共140页，创作于2023年2月（3）薄板小挠度弯曲问题薄板弯曲问题的位移：与空间问题不同，表征单元内部一点位移的物理量与另外位移的一阶导数相关。薄板小挠度弯曲问题的基本假定条件：①弯曲后板的中面保持无应力状态，即无面内应力。②弯曲板满足直法线假定条件。③挤压应力引起的变形可略去不计。3.1弹性的基本理论第65页，课件共140页，创作于2023年2月形变为：一般我们可以看到表征弯曲板形变特征的矢量应力：第66页，课件共140页，创作于2023年2月如图:建立坐标系oxztab

ijkm第67页，课件共140页，创作于2023年2月应力在薄板的横截面上合力应力可以通过先求内力矩得到称其为内力矩矢量第68页，课件共140页，创作于2023年2月薄板有限元的基本概念①单元位移函数小挠度弯曲薄板的应力、应变、以及位移均与挠度函数的二阶导数相关。如果在构造位移函数时，要求函数不仅在单元内部并且在相邻单元的边界上均要满足挠度位移函数及其一阶导数连续的条件，难度是非常大的，所以对于本问题的单元位移函数具有以下特点：仅满足位移函数的完备性、几何同向性而不满足在相邻边界上位移的连续性。对应的单元为非协调元位移函数的特点：第69页，课件共140页，创作于2023年2月③三角形板单元：9个矩形板单元：12个当单元边长选取足够小时，这两种单元可以用于求解薄壳问题。②每个节点3个常用的薄板单元形式：三角形板单元、矩形板单元节点自由度：第70页，课件共140页，创作于2023年2月一个三角形单元的节点位移矩阵与节点力矩阵：：一个矩形单元的节点位移矩阵与节点力矩阵：第71页，课件共140页，创作于2023年2月③矩形弯曲单元的位移函数：讨论：所选位移函数的特点④位移矩阵——形函数的求解非协调元第72页，课件共140页，创作于2023年2月将节点坐标即对应节点位移代入可出系数可整理得到w的插值函数以及形函数单元内部挠度位移：第73页，课件共140页，创作于2023年2月⑤矩形板的几何矩阵⑥矩形板的内力矩阵内力矩阵第74页，课件共140页，创作于2023年2月⑥应力矩阵或内力矩阵：第75页，课件共140页，创作于2023年2月第76页，课件共140页，创作于2023年2月第77页，课件共140页，创作于2023年2月⑦单元刚度矩阵基本公式：给节点虚位移：真实应变：即：节点力在虚位移上所作的虚功=虚位移引起单元内部的虚应变能单元上真实的节点位移：真实应力：相应的虚应变：节点力在虚位移上所作的虚功：虚位移引起单元内部的虚应变能：推导：第78页，课件共140页，创作于2023年2月单元刚度矩阵第79页，课件共140页，创作于2023年2月（4）空间四面体单元一边描述空间实体位移为三个量：在x、y、z方向上的线位移。单元的节点自由度：除了将位移u、v、w外。有时还可以将位移对坐标的偏导数也取为节点的自由度，这样可以增加节点自由度，提高插值位移函数的精度。所以对于空间四面体单元：可以有自由度数为12、48空间六面体单元：可以有自由度数为24、96每个节点自由度3个或12个。第80页，课件共140页，创作于2023年2月四节点四面体单元1234第81页，课件共140页，创作于2023年2月12自由度四面体的节点位移：位移函数：其中:形变为常形变第82页，课件共140页，创作于2023年2月几何矩阵：其中:形变为常形变第83页，课件共140页，创作于2023年2月物理方程：对称应力：第84页，课件共140页，创作于2023年2月单元刚度矩阵分割自矩阵形式第85页，课件共140页，创作于2023年2月第86页，课件共140页，创作于2023年2月有限元在组合结构——船体的应用1、船舶结构有限元计算的主要结构形式①单纯的杆系结构②局部板架结构——平面板架结构甲板板架船底板架舷侧板架结构船底板架结构③刚架结构空间梁结构④空间梁结构⑤空间板、梁结构⑥全船板、梁结构第87页，课件共140页，创作于2023年2月2、一般结构的单元选择与网格划分的原则依据结构受力的特点以及考虑问题的出发点不同来进行选取。网格划分依据应力分布的突变程度进行划分。建立的有限元模型必须是稳定的结构。3、边界条件处理满足结构不能出现刚体位移；局部结构计算是选取模型简化的区域足够大，以至边界的约束对结构计算结果不会产生过大的影响。对称性的使用；能够较真实地反映实际约束状态第88页，课件共140页，创作于2023年2月4、作用的外力简化问题应尽量与实际情况接近，并依据所选单元的特点加以简化第89页，课件共140页，创作于2023年2月第90页，课件共140页，创作于2023年2月不稳定的结构稳定的结构第91页，课件共140页，创作于2023年2月三维梁单元及壳体单元第92页，课件共140页，创作于2023年2月平面单元梁单元第93页，课件共140页，创作于2023年2月第94页，课件共140页，创作于2023年2月单元划分的疏密程度是否合理将影响计算精度与计算时间第95页，课件共140页，创作于2023年2月第96页，课件共140页，创作于2023年2月第97页，课件共140页，创作于2023年2月第98页，课件共140页，创作于2023年2月第99页，课件共140页，创作于2023年2月第100页，课件共140页，创作于2023年2月一个强框架的典型网格第101页，课件共140页，创作于2023年2月第102页，课件共140页，创作于2023年2月第103页，课件共140页，创作于2023年2月第104页，课件共140页，创作于2023年2月第105页，课件共140页，创作于2023年2月含有极细网格的整体舱段模型部分第106页，课件共140页，创作于2023年2月第107页，课件共140页，创作于2023年2月有限元模型的纵向范围第108页，课件共140页，创作于2023年2月有限元模型的示例第109页，课件共140页，创作于2023年2月有限元模型的示例第110页，课件共140页，创作于2023年2月定义热点：热点系指可能产生疲劳裂纹的部位。名义应力：名义应力系指考虑宏观几何影响但不考虑由于结构不连续性和焊缝存在而引起的应力集中的构件的应力。名义应力可用粗网格有限元分析方法求得，也可用简化方法求得。第111页，课件共140页，创作于2023年2月热点应力：热点应力系指在热点处的局部应力。热点应力考虑到连接处的几何形状引起的结构不连续性所产生应力集中的影响，但不考虑焊缝的影响热点应力可通过规定的细化网格有限元分析，也可以名义应力乘以应力集中系数求得。第112页，课件共140页，创作于2023年2月切口应力：切口应力应通过将热点应力乘以疲劳切口因子求得。切口应力系指焊趾处的峰值应力，考虑了结构几何形状和焊