2013-计算力学习题讲解.docVIP

幻灯片1 计算力学习题讲解 东南大学土木工程学院 幻灯片2 1，试述弹性力学中按位移求解与有限元单元法中按位移求解之间的异同点。 幻灯片3 1，试述弹性力学中按位移求解与有限元单元法中按位移求解之间的异同点。 弹性力学 有限单元法 物理模型 基本方程 解法 解答形式 解答精度 幻灯片4 1，试述弹性力学中按位移求解与有限元单元法中按位移求解之间的异同点。 弹性力学 有限单元法 物理模型 连续体 基本方程 几何方程 物理方程 平衡微分方程 解法 解微分方程 解答形式 用函数表示 解答精度 精确解 幻灯片5 1，试述弹性力学中按位移求解与有限元单元法中按位移求解之间的异同点。 弹性力学 有限单元法 物理模型 连续体 离散化结构 基本方程 几何方程 物理方程 平衡微分方程 几何方程 物理方程 结点平衡方程 解法 解微分方程 解代数方程 解答形式 用函数表示 用数值表示 解答精度 精确解 近似解 幻灯片6 2，以位移为基本未知量的有限元法其解具有下限性质，试证明之。 解：系统总位能的离散形式 求解的方程 在平衡情况下，系统总位能等于负的应变能。在有限元解中，由于假定的近似位移模式一般来说总是与精确解有差别的。 幻灯片7 2，以位移为基本未知量的有限元法其解具有下限性质，试证明之。 根据最小势能原理，得到的系统的总位能总会比真正的总位能要大，故 幻灯片8 2，以位移为基本未知量的有限元法其解具有下限性质，试证明之。 幻灯片9 3，某平面结构采用四结点矩形单元和三结点三角形单元建立有限元计算模型，其如图所示。试求结点2的等效荷载列阵{R2}。 平衡所用到的的量均要属于节点的量，如单元位移、单元力。载荷亦应如此，必须将体积力、表面力转化到节点上去，成为等效荷载（载荷）。 幻灯片10 把总等效荷载? R ?e分解成体积力、表面力和集中力的等效荷载之和，有 ?FV?——单元上体积力的等效荷载 ?FS?——单元上表面力的等效荷载 ?pC?——单元上节点上的集中力 1、体积力的等效荷载 2、表面力的等效荷载 幻灯片11 3，某平面结构采用四结点矩形单元和三结点三角形单元建立有限元计算模型，其如图所示。试求结点2的等效荷载列阵{R2}。 ② ① 幻灯片12 单元①荷载作用在1-2边上，故等效结点力只与1、2号结点有关， ② ① 形函数 线性分布面力 幻灯片13 ② ① 幻灯片14 单元②形函数 故结点2的等效荷载列阵为 幻灯片15 4，图中两个三角形单元组成平行四边形，已知单元①按局部编码i , j , m 的单元刚度矩阵和应力矩阵是 按图中所示单元②的局部编码写出单元刚度矩阵和应力矩阵。 幻灯片16 平面问题的单元刚度矩阵[k]不随单元（或坐标轴）的平行移动或作n?角度（n为整数）的转动而改变。 幻灯片17 5，如图所示8结点矩形单元（每边中点为结点），3点为坐标原点，a=b=2，单元厚度为t。 （1）求该单元的位移函数和形函数，并检验其是否满足收敛性条件。 （2）求在2-6-3边作用均布水平荷载q时的等效结点荷载。 幻灯片18 Pascal三角形 幻灯片19 幻灯片20 课本内容 (1) 构造角节点的插值函数：双线性乘积，暂时保证在本节点处为1,在其余角节点处为0 幻灯片21 课本内容 (2) 构造边中节点的插值函数：某个方向的一次项与另一个方向的Lagrange多项式的乘积。保证在本节点处为1,在其余所有节点处为0，但在单元内部不为0。若没有中心点(单元内），由此构造的插值函数为最终的结果。如有单元中心点，则需对边中节点修正 幻灯片22 课本内容 (3) 使角节点的插值函数在边中点处等于0 幻灯片23 局部坐标 幻灯片24 （2）荷载作用在2-3边上，故等效结点力只与2、6、3号结点有关 幻灯片25