按关键词阅读: 构造 函数 单元
6.5.1 单元载荷的移置根据虚位移原理 , 等效结点力所做的功与作用在单元上的集中力、表面力和体积力在任何虚位移上所做的功相等 , 由此确定等效结点力的大小 。
对于平面三角形单元 , 有(5.132)式中 , 单元结点虚位移列阵 , 单元内任一点的虚位移列阵;等号左边表示单元的等效结点力Re 所做的虚功;等号右边第一项是集中力G所做的虚功 , 等号右边第二项是面力q所做的虚功 , 积分沿着单元的边界进行;等号右边第三项表示体积力p所做的虚功 , 积分遍及整个单元;t为单元的厚度 , 假定为常量 。
用形函数矩阵表示的单元位移模式方程为(5.133)代入式(3-132) 。
20、 , 注意到结点虚位移列阵d*e可以提到积分号的外面 , 于是有(5.134)注意到(d * e ) T 的任意性 , 上式化简为R e = F e +Q e +P e (5.135)其中F e = N T G (5.136)(5.137)(5.138)式(5.134)右端括号中的第一项与结点虚位移相乘等于集中力所做的虚功 , 它是单元上的集中力移置到结点上所得到的等效结点力 , 它是一个61阶的列阵 , 记为F e 。
同理 , 式(5.134)右端括号中的第二项是单元上的表面力移置到结点上所得到的等效结点力 , 记为Qe;第三项是单元上的体积力移置到结点上所得到的等效结点力 , 记为P e 。
6.5.2 结构整体载荷列阵的形成结构 。
21、载荷列阵由所有单元的等效结点载荷列阵叠加得到 。
注意到叠加过程中相互联接的单元之间存在大小相等方向相反的作用力和反作用力 , 它们之间相互抵消 , 因此 , 结构载荷列阵中只有与外载荷有关的结点有值 。
下面逐项进行讨论 。
(1)集中力的等效载荷列阵逐点合成各单元的等效结点力 , 并按结点号码的顺序进行排列 , 组成结构的集中力等效载荷列阵 , 即(5.139)上式中 , 单元e的集中力的等效结点力为(记单元结点局部编号为i,j,m)(5.140)式中(i, j, m) (5.141)式中 , (Ni )c 、(Nj )c 、(Nm )c 为形函数在集中力作用点处的值 。
(2)表面力的等效载荷列阵把作用在单元边界上的表面力移置到结点上 。
22、 , 得到各单元的表面力的等效结点力 。
按照结点号码的顺序进行排列 , 逐个结点叠加合成后 , 组成结构表面力的等效载荷列阵 , 即(5.142)式中 , (5.143)由于作用在单元边界上的内力在合成过程中已相互抵消 , 上式中的结点力只由作用在结构边界上的表面力所引起 。
(3)体积力的等效载荷列阵与表面力类似 , 体积力的等效载荷列阵也是由单元体积力的等效结点力按结点号码顺序排列 , 在各结点处合成得到(5.144)式中 , 单元e的体积力的等效结点力为(5.145)6.5.3载荷移置与静力等效关系上述基于形函数的载荷等效所得到的结果与按照静力学的平行力分解原理得到的结果完全一致 。
例如 , 如图3-19所示的单元e , 在ij边上作用有 。
23、表面力 。
假设ij边的长度为l , 其上任一点P距结点i的距离为s 。
根据面积坐标的概念 , 有 ,, (a)代入式(5.137) , 求得单元表面力的等效结点力(b)可见 , 求得的结果与按照静力等效原理将表面力q向结点i及j分解所得到的分力完全相同 。
图3-19 表面力等效示意再如 , 从图3-20所示的单元e的A点处取体积微元tdxdy , 作用在其上的体积力为ptdxdy , 为便于分析 , 认为力的作用方向与单元平面垂直 。
根据平行力分解原理 , 对jm边取力矩 , 求得结点i处的分力为(c)整个单元e的体积力在结点i处的分力为(d)类似地 , 分别对im及ij边取力矩 , 可得到结点j和结点m处的分力(e)(f)图3-20 体积力等效示意 。
24、因此 , 对于平面三角形单元 , 按照静力学中平行力的分解原理所得到的结点力与按照虚功原理求得的结点力完全一致 , 在实际计算等效结点力时 , 可以直接应用静力学中有关平行力分解的结果 。
例如 , 对均质等厚度的三角形单元所受的重力 , 只要把1/3的重量直接加到每个结点上 , 对于作用在长度为l的ij边上强度为q的均布表面力 , 可以直接把 (qtl) /2 移置到结点i和j上 。
习题3-1 解释基本概念:位移插值函数、位移模式、有限元解的收敛准则、位移解的下限性质 。
3-2 简答下列问题:什么是有限元解的收敛性?什么是解的收敛准则?什么是形函数?它有什么性质?如何建立有限元的形函数?3-7 推导基于形状函数和结点的一维线性插值 。
【单元|单元形函数的构造】25、格式 , 将结果表示成矩阵形式 。
3-8 横截面面积为常数的弹性杆两端固定 , 杆长为 , 弹性杆各处受相同的体积力作用 , 试采用3个长为的线性单元 , 用形函数(不用插值多项式)给出Rayleigh-Ritz解的表达式 。
来源:(未知)
【学习资料】网址:/a/2021/0318/0021711081.html
标题:单元|单元形函数的构造( 四 )