数学模型和数值分析方法.pptVIP

2有限元法的基本概念-直接刚度法12.2有限元法通过一个例子来介绍直接刚度法，说明有限元法求解的一般步骤。例：考虑一个变截面杆，如图示。杆的一端固定，另一端承受P=1000N的载荷，杆的顶部宽w1=2cm，杆的底部宽w2=1cm，杆的厚度t=0.125cm，长度L=10cm，杆的弹性模量E=10.4×105MPa。试分析该杆沿长度方向不同位置的变形情况，质量忽略不计。1.前处理过程(1)求解区域离散化将求解的问题分解为结点和单元第29页，共36页，星期日，2025年，2月5日3.2数学建模软件简介2.2有限元法(2)建立单元位移方程先考虑一截均一长度为l、截面为A的固体单元受外力F的变形情况。以弹簧的变形来模拟固体单元的变形，等价刚度为：keq=AE/l前面问题可近似将杆模型化为不同截面的等截面杆的串联。杆在顶端固定，静态平衡时，可获得如下方程：10000-k1k1+k2-k2000-k2k2+k3-k3000-k3k3+k4-k4000-k4k4u1u2u3u4u50000P=第30页，共36页，星期日，2025年，2月5日2.2有限元法(3)建立单元刚度方程每个单元有两个结点，每个结点可建立两个方程。方程包含结点位移和单元的刚度。结点i和i+1处可有方程：(4)单元集成将上式用于所有单元，进行集成，得到总体刚度矩阵：keq-keq-keqkeqfifi+1uiui+1=k1-k1000-k1k1+k2-k2000-k2k2+k3-k3000-k3k3+k4-k4000-k4k4=K(G)(5)施加边界条件和载荷2.求解阶段3.后处理阶段第31页，共36页，星期日，2025年，2月5日(一)加权余量法工程或物理学中的许多问题，通常是以未知场函数应满足的微分方程和边界条件的形式提出来的，可以一般地表示为未知函数u应满足微分方程组：2.2有限元法2有限元法的基本理论2域Ω可以是体积域，面积域等，如图示，同时未知函数u还应满足边界条件：Υ是域Ω的边界(2-1)(2-2)第32页，共36页，星期日，2025年，2月5日2.2有限元法要求解得未知函数u可以是标量场（如温度场），也可是几个变量组成的向量场（如位移、应变、应力等）。A、B是对于独立变量的微分算子。微分方程组是（2-1）在域Ω中每一点都必须为零，因此：(2-3)其中V是函数向量，它是和微分方程个数相等的任意函数。同理，假如边界条件式（2-2）也同时在边界上的每一点都得到满足，对于一组任意函数V应有：(2-4)因此若积分形式对于所有的V和V都成立，式（2-5）称为微分方程的等效积分形式。(2-5)第33页，共36页，星期日，2025年，2月5日2.2有限元法对于微分方程式（2-1）和边界条件式（2-2）所表达的物理问题，未知函数u可以用近似函数来表示。近似函数是一族带有待定参数的已知函数，一般形式是：(2-6)式中，ai是待定参数，Ni是称为试探函数的已知函数。在式（2-5）中，用n个规定函数来代替任意函数V和V，即：(2-7)可以得到近似的等效积分形式：权函数采用使余量的加权积分为零的方法来求微分方程的近似解的方法称为加权余量法。(2-8)第34页，共36页，星期日，2025年，2月5日2.2有限元法按照对权函数的不同选择，加权余量法分为以下几种：1.配点法若域Ω是独立坐标x的函数，δ(x-xj)则有如下性质：当x≠xj时，Wj=0，有这种方法相当于简单地强迫余量在域内n个点上等于零。2.子域法在n个子域Ωj内Wj=I，在子域Ωj以外Wj=0。此方法的实质是强迫余量在n个子域Ωj的积分为零。3.最小二乘法近似解取则权函数此方法是使