2012年计算力学第三课.ppt

哈尔滨工程大学航天与建筑工程学院 王滨生 E-mail: wang_bs@126.com; 本次课的主要问题： 有限元的基本概念 什么是位移函数 位移函数的选取原则、收敛性讨论;单元和节点 自然单元：桁架 分割单元：自然分割单元：刚架、阶梯形梁等 任意分割单元：板。; 单元划分原则： ⑴划分单元的个数，兼顾精度与速度。 ⑵划分单元的大小，可根据部位不同有所不同。 ⑶划分单元的形状，但要注意。必须节点与节点相连。 ⑷单元各边的长度不要相差太大，否则将影响求解精度。 ⑸尽量把集中力或集中力偶的作用点取为节点。 ⑹尽量利用对称性，以减少计算量。; 节点?：分主外节点、副主外节点和内节点。 节点力和节点荷载 节点力：相邻单元之间的相互作用是通过节点来实现的，这种通过节点的相互作用力叫做节点力，即节点内力。 如Xij，X表示水平方向，第一个下标表示节点编号，第二个下标表示单元编号。 节点荷载：作用在节点上的外荷载称为节点荷载。包括 ⑴作用在节点上的外力； ⑵单元的分布力(按静力等效原则)等效到节点上的外力。 静力等效原则 : ⑴刚体：等效前后具有相同的主矢和主矩。 ⑵弹性体：等效前后两者在虚位移上的的虚功相等。;将单元的重力G向各节点移置 在线性位移模式下，该单元的重心 C的虚位移为 根据虚功相等，有 同理; 单元重心C的相应虚位移为 根据虚功相等，有 同理;位移函数 定义：用以表示单元内的位移或位移场的近似函数，称为位移函数。 选取原则：位移函数应满足 ⑴位移函数在单元内部是连续的 ⑵两相邻单元在交界处的位移是连续的; 选取方法?? ⑴广义坐标法(待定系数法) 如一维单元 其中 如二维单元 其中 {F}按Pascal(贾宪——杨辉?——帕斯卡)三角形选取 ; ⑵插值函数法(试凑法)? 将位移函数用节点位移乘上一个已知插值函数之积的和来表示。 如一维单元 如二维单元 *前者用途较广，但求广义坐标，须计算一些矩阵的逆阵，如果矩阵奇异，则逆阵不存在。 ;平衡问题 一般来说 ⑴单元内部的平衡条件往往不能满足? ⑵单元之间的平衡条件常常得不到满足? ⑶节点力的平衡条件是满足的 *如果单元的组合数量越来越多，任何破坏连续性和平衡条件的情况都会趋于缓和。;收敛性讨论 问题 ⑴由于受计算机字长的限制的误差? ⑵进行相差悬殊的数值作加减运算等会产生误差 ⑶边界上用直线代替曲线导致离散化模型与实际物体的误差? ⑷单元位移函数与实际单元的位移场一致的误差? *在单元内部和单元边界上的连续条件都得到满足的情况下，所取单元尺寸逐步减小时，有限元解就可收敛于问题的精确解。; 收敛条件 ⑴充分条件：单元内部连续+单元边界连续 ⑵必要条件：单元内部连续 收敛准则 将两个连续性条件用以下三个判别准则： ⑴单元的位移函数在单元内部连续，在两相邻单元的公共边界上必须协调。(后者表示既不重叠也不分离) ⑵单元的位移函数应包括常应变状态。