有限元数学原理.pptVIP

几何方程：ab的变形为：因此正应变为：其中ρ为曲率半径。物理方程：由广义虎克定律有整理得边界条件。弯矩x方向平衡y方向平衡。其中c1,c2,c3是待定系数。最后可得求解方程得讨论应变能外力功势能由于单元有四个位移分量，可设梁单元的位移模式v(x)为包含4个待定常数的三次多项式：有限元分析步骤-单元分析式中根据梁的平面假定可知梁单元的轴向应变为：这里利用平面假设（变形后横截面仍保持平面，与纵线正交）如图：从而可以由单向虎克定律得出单元的轴向应力：由虚功原理可以推得2弹性力学问题近似求解的加权残值法1直接针对原始三大类方程边界条件下求解三大类变量往往是非常困难的，尤其是当几何形状和边界条件比较复杂时，一般求不出相应的解析解。如果事先假定满足一定边界条件的试函数，再在此基础上进行近似求解，则可以大大降低求解难度。这种试函数方法可以使得求解过程比较规范和简单，并有一定的适应性，但是求解的精度有所降低。2试函数方法的基本原理：先假定满足一定边界条件的试函数，然后将其带入需要求解的方程中（控制方程），通过使与原来的方程的误差残值最小来确定试函数中的待定系数。为了提高解或逼近精度，可以采用较多项数的试函数来进行计算，这种方法叫做加权残值法。加权残值法WRM:Galerkin加权残值残值最小二乘法12345设满足以下方程和边界条件的位移场为公式中的L为微分算子。由于平面弯曲梁的平衡方程为故对于更一般的情形，设有一组满足所有边界条件的试函数，将其线性组合为新的试函数其中c1,c2,c3…cn为待定系数。假设能找到事先满足式中的边界条件的一个试函数，将其带入到控制方程，则一定存在残差，记为将试函数代入原始方程组，则必有残差，真实的c1,c2,c3…cn使得残值的积分为零，即其中w1,w2,w3…wn为权函数。以上为关于c1,c2,c3…cn的方程组，由上式可以求出他们，最后由线性组合形式的试函数得到真实。如果将权函数w1,w2,w3…wn取为?1,?2,?3…?n，则该方法称作伽辽金法。代入控制方程得残差由Galerkin加权残值方程分析可得，试函数取为求解上式可得代入控制方程得残差由Galerkin加权残值方程分析可得取另一个试函数为求解上式可得同时满足面力边界条件，根据Galerkin法分析可得仅仅取1项试函数时，由伽辽金加权残数法得到的结果与精确解得相对误差为0.3861%。02几种函数结果比较01由以上比较可以看出，此解法精度还是比较高的。04仅仅取2项试函数时，由伽辽金加权残数法得到的结果与精确解得相对误差为-0.027%。03同样，设有满足所有边界条件的试函数，若将试函数代入原始方程中，则必有残差值，真实地待定常数使得残差值平方的加权积分取极小值，即其中w为权函数，一般可以取1。将上式进一步具体化，即对上式取极值，有这是一组关于c1,c2,c3…cn的方程组，由上式可以求出它们，从而得到真实的试函数。受均布外载荷简支梁的残值最小二乘法求解将其带入到控制方程，则一定存在残差，取权函数为1，则残差平方的积分为由最小二乘法有由上式可以解出与用伽辽金加权残值法相同的结果。试函数取为平面问题有：三大应力分量，三个应变分量，以及两个位移分量。三个物理方程，三个几何方程，两个平衡方程，三类边界条件基于三大物理量之间的关系（见三大方程），可以得到基于位移表达的平衡方程（控制方程）如下：加上边界条件，就构成了基于位移求解的平面应力问题的所有方程。设有位移的试函数为式中的cui和cvi是待定常数，ψui和ψvi是满足所有边界条件的基底函数。将试函数代入控制方程，使得其残值在加权积分下为零，即可得到Garlerkin加权残值方程，得到一组关于待定系数的方程组。同样，也可以得到空间问题的加权残值法，包括伽辽金加权残值法和最小二乘法。可以看出，弹性问题的加权残值法的特点为：1、试函数要满足所有边界条件，即位移边界条件和力边界条件2、积分中试函数的最高阶导数较高（对梁的弯曲问题，导数为4阶，而对于其他一般弹性问题为2阶），因此试函数对连续性要求较高3、整个方法为计算一个全场的积分4、有求取积分问题的最小值，将原方程得求解化为线性方程组的求解5、整个方法的处理流程比较规范6、难点在于如何在全场范围内寻找同时满足所有边界条件的具有较高连续性的试函数？3.3弹性问题近似求解最小势能原理及其变分原理最小势能原