《有限元方法及软件应用有限元平面问题3》由会员分享，可在线阅读，更多相关《有限元方法及软件应用有限元平面问题3（18页珍藏版）》请在金锄头文库上搜索。
1、Fnie EemenMeho3.面力的移置设三角形单元某边界s 上受面力q 作用,分量为，,则取d 则由一般公式： 积分在边界s上以上三种载荷的等效节点荷载由公式e导出一般我们称： 为荷载移量的一般公式:几点阐明:1 虚功等效静力等效。 唯一性2 一般3 更多节点的单元公式形式不变，但不同4 虽然公式e导出但对于面力和体力的计算都是很麻烦和困难的 为,y 的函数,若p, q再为 x， y 的函数则更难,且单移分限不好定。因此，我们将来还要进一步把这个问题解决好。四 三角形单元的面积坐标（自然坐标,局部坐标)1 面积坐标的定义:图示三角形单元I ,j,k中任意一点m ,其位置可由xoy坐标系中两
2、个坐标来拟定,即m(x,y)若我们连接，,，则形成了3 个小三角形jm， ikm,jk.则有:若m(x，y)拟定im, ikm, j.面积拟定。反之,ijm, im,m.面积拟定(x,）拟定(用同底等高的概念解释！!）因此，三角形单元内任一点可以我们如何用三角形面积来描述m点的位置呢?定义:节点I对边为底的三角形面积为；节点j对边为底的三角形面积为;节点k对边为底的三角形面积为；设三角形单元的面积为A令 (-3）则三个比值，,称为三角形单元中m点的面积坐标.三角形面积坐标的性质：1 面积坐标为三角形单元的局部坐标,与三角形的形状及位置无关。其定义域为 ;2 三个面积坐标之和:+.即只有两个面积
3、坐标是独立的。(2-38)证明：+=+=(+)= （亦可几何解释）。3 三角形单元内与jk边平行的直线上各点相似(轮换）。(同底等高三角形=)4 形心处的面积坐标为： =13 （2-39)5 三角形单元节点的面积坐标为: (2-）证:节点I: =A. =0.3三角形面积坐标与直角坐标及形函数的关系下面我们来推导面积坐标与直角坐标的关系：设m点的坐标为m（x,y)，m为任一点则：= (）（）+()+(）显然： , ，=(） （2-41) 与体现式比较可知:三节点三角形单元的面积坐标就是其形函数。（对于一般的状况：面积坐标永远是线性坐标而形函数可以是非线性的,后来我们可以把形函数用面积坐标表达)即
4、=,， （24)具有的所有性质 式(241)还可写成矩阵的形式： 直面 （2-44)这就是直角坐标与面积坐标的转换关系。下面的成果留给人们自己证明: 面直 （2-45）4 面积坐标函数的运算我们可以不加证明得地给出面积坐标函数的微积运算成果。(证明复杂麻烦用函数等）.偏导 设z=f(，) g(，y） (I= I , ,)则: (24）. 面积分 (47）其中,为正整数; 0!1,A：三角形面积ex： (I ,j ,k)3.线积分： (s为直线长) (2-48）以上公式要会用 注意表达的边五 三角形单元的荷载移置有了面积坐标与形函数的关系，我们即可对荷载移置进行计算了。1 集中力的移置 设m点作
5、用有集中力m点的形函数为:（I= I ,j ,k)等效节点荷载为：这就是三角形单元内点作用有的等效节点荷载。只要计算出(I= ，j ，k)即可。作为特例，考虑三角形单元形心处重力的移置。形心坐标：=0 =-R故：重力作用于形心时各节点均担。2. 体积力的移置设单元作用有体力 则等效节点荷载为： 若为x，y的函数，则把用面积坐标表达(转换)在常体力的作用下有:=即：常体力作用下，总体力均分三节点。2 面力的移置。 设三角形单元I ,j边上作用有梯形分布的面力q 由面力移置公式得:（可分别由节点合力表达及用节点分力表达) = （为合力，非分力）则 = q为x, y 的函数，把q表达面积坐标的函数有
6、q= ,在门边上是线性坐标，可运用两点式方程写出。 则：=同理：=注意到在s上=0=0故:=或:此法:1 避免复杂的分离。 2 便于编程计算。特例:若分布荷载为三角形分布。令(或）则有： (近端为2 ，远端为1)阐明:用以上积分的措施求等效节点荷载合用于任意节点的三角形单元，形函数也未必是线性的。六 三角形单元节点荷载的形成 通过荷载解决后，我们已把非节点的荷载转化为常点荷载。实际计算的荷载为：计算荷载=原节点荷载+等效节点荷载即： (2-4）等效节点荷载要注意：1 同步奉献的问题2 用哪个单元计算的问题。七计算成果的整顿：有限元计算提供的成果一般为：1。节点位移 。单元应力1 节点位移的解决
7、:一般把节点位移按比例标出，提供出构造常点位移分布规律1 连成折线（线性位移函数）2 连成光滑曲线（实际变形）2 单元应力的解决：输出的单元应力一般为,，（形心处)（三节点单元为常应力元,无所谓）1 变换为单元的主应力。, (材力）2 变换为节点应力的主应力。 (） （平均法) Ex. 节点5的应力为：即:= （x,， xy )(,)(,)然后标出应力变化曲线。计算成果的工作量随构造的单元,节点划分增长面增大。要关注的是：)位移的变化规律 2）应力的最大值及发生地点。小概念：位移最大的地方，应力未必最大。八：有限元计算小结:1 基本原理:持续法有限个节点连接,有限大小的单元的组合法。建立的节点
8、位移为未知数,总刚为系的阶线性代数方程组。2 研究措施(拟定)节点位移单元位移单元应变单元应力 单刚总刚计算荷载（等效节点荷载） 约束解决求解方程整顿成果.解答特点:1.假定单元内的位移分布规律，近似离散的数值解。2.误差重要来源于:构造离散(持续离散),假定位移分布。3收敛性:单元缩小划分细密收敛于精确解。Cap3. 平面问题较精密单元的分析(矩形,高阶单元，等参单元). 问题的提出在三角形单元中,我们假定位移函数是线性的。即:单元内的位移按线性规律变化。这是最简朴,最基本的一种有限单元。而实际构造中在外载荷的作用下位移分布常非按线性变化。设单元位移曲线为图示的f(x)显然，用线性插值解的精
9、度较差。提示解的精度的措施:1 增长单元，节点数（工作量大，费用高)2 提高插值阶数因此，提出了用高阶插值,高阶单元的问题我们人们都懂得,位移是一种持续函数（持续体),而任意的持续函数都是可以展成幂级数,用幂级数来表达的。因此，一种单元内的位移分布为f(x)时，我们就可以取级数的前几项来表达它。用二次三次函数来插值,以改善计算成果。至于等参数单元(等参单元）是一种为清除曲边误差而出的一种单元。如果实际构造为曲线边界,则无论如何提高位移函数（插值)的阶数也不能使解得到多大的改善。有限个直边替代曲边,究竟是替代，而不会是相等。为理解决曲边问题,人们提出了等参单元的概念。有平面等参单元,空间等参单元
10、等。我们只向人们简介平面等参单元，以供理解。3-2 四节点矩形单元的有限元分析。 矩形单元常用于规则边界的有限元分析,它也是常用的一种有限单元。一 单元的位移函数。设单元e为矩形单元，边长为2a，2b;其节点为I，j,k,m；为研究以便我们取局部坐标系x-y如图(原点在形心)。1 单元的自由度及位移函数：个节点,每个节点2个自由度（位移)u，v，则单元的总自由度为8个。为保证单元的收敛性准则，位移函数必须保证有常数项，线性项。设位移函数为:(对称性与坐标选择无关） （U 有关坐标对称) （31)其中xy项是根据pasca三角形及y的对称性选用（选二次项尚有协调问题,选，项不行）这样,已知(i,
11、2,3，4)这八个节点位移（i1,8）2.形函数的推导：同理： 我们不难从前4个方程中解出,，,具体做法： +：v+: -: - : v- v =u=令: （3-2）则: (3-3)式（3-2）称为节点矩形单元的形函数;式（2)尚可写为: （II， j, k， ) （3-4）式(3)为节点位移表达的单元位移函数。式（3-3）还可以写成矩阵的形式 = （3-5)式中: (3-6)二 形函数的性质:1形函数(I =I, j，,m)在节点 上=,在其他节点上=0 (轮换）即在 节点I ：1 节点j： （)= = 节点k: 节点m： 证明: （=I, j, k，)在节点ix,y=时:=1, k, m各
12、节点至少有一种节点坐标x-或y=故=0 (i, k， m)同理可得到所有成果。24个形函数之和：+证明：写出形函数(，的符号与该点的，相似) += =1因此:4个形函数只有3 个是独立的。在 j边上,0（节点除外），=0 （轮换）(一条边上的四个)证明:在I 边上,=, （x-)在I, 边上,= (x-） I， j边上 y=-=0同理：I, j边上：=0。 证毕.在4条边界上的性质(节点除外) 在涉及节点I 的边界上,0，否则=0。（轮换）(四条边上的一种）证明:性质3 的另一种表述 在涉及节点 的边界上：, 或 y=显然有: ()或 (节点j, m除外）在不涉及节点I 的边界上:X= 或y=- 显然：=0得证由以上的性质，我们可以描述的几何图形。三 位移函数的性质1 位移函数是双线性的位移函数： 显然, u, v涉及坐标x y的二次项 ，但当x=n 或y=cost时,V都是一种线性函数。 即：在单元内任一点，无论沿x方向变化（此时y=const）或沿y方向变化(
