与有限元解法有关的问题

四、四、与有限元解法有关的问题

在讨论空间与平面问题时,为简化起见,坐标(x,y,z)改用(),u仍表示位移列矢量,其三分量用()代替(u,v,w);对体力f等荷载仍用()表示其分量;对积分,除列式演算外,一般只用一个积分号表示.

[变形能的正定性与刚性位移] 弹性系数矩阵D是正定的,从变形能V=0可推出,即对应变来说它是正定的,但是对位移来说,它却退化为非负的,即不能由推出,因为弹性体还存在无应变的位移即刚性位移.

1° 1° 一般空间问题设位移矢量场u的旋度rotu的分量为

它与位移分量以及应变分量之间的关系式为

因此,在小变形的假定下,应变表示整个弹性体除了作微小的位移外,还可作微小的转动,它的刚性位移可写成

式中表示平移的分量,表示绕三坐标轴的转角,它们都是任意微小的常量.总之,弹性体无应变的微小刚性位移共有六个自由度,设

则上式可简写成

式中就是第i坐标轴的单位矢量,就是绕第i坐标轴的旋臂矢量.于是

(i=1,2,3)

分别表示分布体力f对弹性体的合力与合力矩沿第i坐标轴方向的分量.对面力q等分布荷载也有类似的结果.

2° 2° 平面应力、应变问题由于,刚性位移简化为

其自由度等于3.分布荷载与,的内积也有类似的定义与含意.

3° 3° 薄板弯曲问题由于位移分量只有挠度,刚性位移可写成

(即）

其自由度等于3.内积

分别表示分布荷载q沿中面的合力以及绕轴的合力矩.

4° 4° 圆柱扭转问题根据假定w=0与广义剪应变,即扭转率可得无应变的刚性位移为.相应的旋臂矢量为.这时,位移也可写成

(即 )

其自由度等于1.设在端面给定分布荷载切力Q(其分布,严格地说,应与剪应力相同),则内积就是沿端面的合力矩M.

[变分问题的可解性] 这里只讨论一般空间问题、薄板弯曲问题与圆柱扭转问题,其他问题类推.

1°1° 一般空间问题先考虑边界纯荷载支承的情况.从§1可以看出,求解位移函数u使弹性体的总势能达到极小值,即

的变分问题可以归结为求解变分方程

(对一切）

的问题,其中位移u的变分力学上称为虚位移,而积分

称为虚功泛函,它是位移u和虚位移的双线性泛函.

上面已经指出,变形能只能得出是任意刚性的虚位移.反之,对一切非刚性的虚位移..因此利用等式

可以证明,变分问题有解的充分必要条件是

(对一切刚性的虚位移）

由于刚性位移的自由度为6,又可分写为6个条件

这表明,要使弹性体在纯荷载支承下达到平衡状态的充分必要条件就是外力对整个弹性体所产生的合力与合力矩都等于零.

反过来,设变分问题有一解,则位移与任意微小的刚性位移之和+也是它的解,这表明弹性体的平衡状态可以相差一个微小的刚性位移,具有6个自由度.换句话说,解的一般形式可写成

如果边界上还有一部分受几何约束或弹性支承,则有解的充分必要条件的个数或自由度都会相应地减少,只要这部分边界条件使弹性体不可能有微小的刚性位移,则问题的解存在而且唯一.例如当u在部分边界上给定位移值时,变分方程改为

（对于)

只要包含不共线的三点,条件就排除一切刚性的虚位移,解存在的充分必要条件的个数等于零了.如果在部分边界上施加弹性支承,则从变形能

与二积分的正定性可得

而C的正定性又推出.因此,只要包含不共线的三点,那末变形能对位移u就是正定的,问题的解也就存在而且唯一.

2° 2° 薄板弯曲问题考虑纯荷载支承的情况,其变形能与外力势能可用广义应变表示如下

变分问题

可归结为求解变分方程

(对一切虚位移)

同样[,]=0只能得出广义应变,刚性虚位移可表示为

其自由度等于3,而变分方程有解的充分必要条件是

(对一切刚性虚位移)

由于上式可分写为三个条件:

前一条件是沿板厚度方向的外力要取得平衡使板不产生移动;后二条件分别表示外力对板所产生的绕x,y轴方向的力矩要取得平衡,使板不产生转动.在这些条件下,板的平衡状态也不唯一而具有3个自由度,即问题的解可写成

如果边界上还有一部分受几何约束或弹性支承,则可解的充分必要条件的个数或自由度都会减少或等于零.例如,只要有不共线的三点给定挠度,或给定一点的挠度与两个绕不同方向的转角就使板不可能有微小的刚性位移,这时问题的解存在而且唯一.对部分边界施加弹性反力或反力矩也会提高变形能的正定性.例如在部分边界上给定弹性反力-cw与反力矩,则从可得Bwº0 与,即使只是一个点也减少两个自由度,如果有两个不同的外法线方向,则自由度为零.这时[w, w]对w是正定的,问题的解也就存在而且唯一.

3^o 圆柱扭转问题考虑纯荷载支承的情况,其变形能与外力势能可用广义位移(即扭转角)ω表示如下

变分问题

可归结为求解变分方程

(对一切)

同样, [,]=0只能得出扭转率.而刚性(广义)虚位移可表示为,其自由度等于1.变分方程有解的充分必要条件是

(对一切绕z轴的扭转角)

这表示沿整个柱面的扭矩要取得平衡.在这个条件下,圆柱的平衡状态可相差一微小的扭转角,即其解可写成

只要有一端受几何约束,自由度就消失;如果有一端(例如z=l)给予弹性支承,则相应的变形能(例如)对就是正定的了.对于这些情况,问题的解存在而且唯一

[变分问题与微分方程定解问题的等价性] 弹性体内部力的平衡方程是在平衡状态下对其微元进行均值分析而推导出来的;同样,其边界条件也是边界微元上内力与外力达到平衡的体现.因此,弹性力学中各种微分方程定解问题就是要找出弹性体因连续变形而产生的内、外力平衡状态.变分问题则是以整个弹性体的能量的形式来表达这种平衡状态,不论是用变形能泛函或虚功泛函,其实质是一样的.这里还要指出它与微分方程定解问题的等价性.下面仍以一般空间问题、薄板弯曲问题与圆柱扭转问题为范例作简单介绍.

1° 1° 一般空间问题力的平衡方程

可以看作一组关于位移分量()的二阶微分方程,其边界条件仍按三种支承方式分别表达如下:

它们包含位移分量及其一阶偏导数.

现在从虚功泛函的形式出发,考虑如下的变分方程

（对)

对应用高斯公式进行分部积分可得

从和变分方程就得到在W内力的平衡方程和,上的边界条件.至于在上的几何约束条件由于在P(u)中没有反映出来,需要强加上去,即在解变分方程时,除了对变分加上条件外,还得要求其满足.反过来,以满足的变分乘上弹性体W内部力的平衡方程与上的边界条件并且作相应的积分,应用高斯公式显然也得出上述的变分方程.

2°2°薄板弯曲问题以中面的挠度w(x,y)为待定函数,弹性体内部力的平衡方程在推导过程已自行满足,而以上板面的边界条件作为w所满足的微分方程

对应的边界条件是由板的边缘截面内外力平衡来确定.由于支承方式的复杂性,下面先对纯弹性支承的情况的虚功泛函作分部积分得出一般形式,然后再来考虑其他支承对各项积分的影响.设

代以B,D及内力等表达式并分部积分,可以证明

式中为边界上沿s正向依序排成的拐点.

由于荷载支承可看作弹性支承的特例(即c或的情况),只要考虑在部分边界上分别给定挠度和转角:

(另一条件只能是弹性反力矩)

(另一条件只能是弹性反力)

(可以重迭)并在拐点有集中力(可以取零,特别当)的情况,这时相应的外力势能可写成

同样在变形能中应加上条件:

对照上述二式,可知变分问题:满足强加的条件

且使

(对)

的解w,是下面微分方程的定解问题的解:

反过来,这个问题的解w显然也是上述变分问题的解,因此它们是等价的.

3°3° 圆柱扭转问题以柱体横截面的扭转角w(z)为待定函数,其广义力(即扭矩)的平衡方程为沿柱长的微分方程:

先考虑纯弹性支承的情况,其边界条件可写成

相应的虚功泛函与外力势能为

对前一式作分部积分,再对照后一式不难看出求解变分方程

（对一切）

的问题同上述二阶常微分方程的定解问题是等价的.端点受荷载支承的情况可看作其特例(即的情况).对于端点受几何约束的情况,只要把变分方程中与该端有关的各项去掉并对ω及改为强加条件,两个问题仍然等价.例如求满足及变分方程

的函数ω(z)的问题与下面微分方程的定解问题是等价的:

[热效应与热当量荷载] 在弹性问题中,弹性体由于受热引起各微元的胀缩,从而调整弹性体内部应力或应变的分布,这种作用称为热效应.设为零应变状态的温度分布;T为小变形后的温度分布,由于温升,弹性体将向周围作均匀的膨胀(或收缩,当).在任何直角坐标系中,其位移分量只与有关,即只有线应变而无剪应变(i¹j),而且在小变形假定下,它们与有线性关系: