01-9结构有限元素法.docVIP

结构有限元素分析理论 绪言 有限元方法在工程技术中的广泛应用 静、动、固、流、宏、微工程理论 直接应用（刚度强度分析） 间接应用（基础工具，优化，可靠性） 关于有限元方法基础 i) 有限元方法是求解数学-物理方程（偏微分方程）系统的数值逼近方法。 e.g. a. 三维空间弹性方程系统 平衡方程： 记法： 张量求和记法： 几何方程： 物理方程： 边界条件： Dirichlet: 基本边界条件： Neumann: 自然边界条件： b. (任意截面的扭转问题) 二维Poisson方程： 基本边界条件：, , , ii) 有限元方法的分析特征 A）将一个偏微分方程系统的解可转换为一个泛函的变分问题。即： 求解泛函变分的驻值等价于有解微分方程。（数学前提，并不是有限元本身） e.g. a) …………………………………………………① b) …………………………② B) 有限元方法是直接求解泛函变分驻值的过程；泛函极小化序列的逼近过程。 过程思路：（Ritz 1908） ① 满足 ②将Un代入②式，得 ③求?J=0的解，获得CK C)计算大型稀疏代数方程组（计算机技术的贡献） 小结： {变分原理。（固体力学中的能量泛函变分原理）变分数值计算（数值积分等）胡海昌 {泛函空间、逼近理论（包括函数逼近理论） D)有限元方法的分类：位移法，力法，混合杂交法。 3．有限元方法及技术的发展 i)方法本身的完善： ii)数学理论基础的完善：{有限元空间的性质，逼近性，弱闭性，嵌入性，紧致性 iii)工程技术方面的完善： 4．本课程内容安排 i)变分法，（理论法，直接法） 9 ii)能量泛函变分原理 9 iii)单元设计（形状函数设计及逼近理论）及刚阵计算 15 iv) 大型方程组数值计算 6 5．课程要求 i) 主要参考书 ii) 强调自觉学习，作业按时完成 iii) 与“结构有限元程序设计”协调（工程结构分析，程序设计） iv) 闭卷考试 v) 上课时间安排 （2～8） 6小时/周 §1 变分原理与变分法 1.1 关于变分原理 i) 大自然总是以可能最好的方式安排一切（昼/夜，日/月，阴/阳，矛盾/统一，静止/运动），存在着各种安排原理。 对静止事物：最小能量原理，平衡原理，对称/相似原理（哲学意味） 对运动事物：能量守恒，动量守恒，熵增原理 变分原理是自然界静止状态中的一个普遍适应的数学定律，获称最小作用原理的数学规划。 例： ①光线是最短路径传播； ②光线入射角等于反射角，光线在反射中也是光最短路径（Heron） ③光线折射遵循时间最短的途径（Fermat） ④最小能量原理------------------------导致有限元法 ii) 变分法是变分原理的数学规划方法，是计算泛函驻值的数学方法 iii) 关于泛函及泛函驻值 a) 定义：空间（集合）元素与一个实数域间的映射关系（算子） {J:} b)举例：①矩阵范数：线性算子（矩阵）空间→数域 ‖A‖1 = ②函数的积分：函数空间→数域 注意：自变量是集合中的元素（定义域）；值域是实数域。 思考题： ① 判定下列那些是泛函： a.; b.; c.3x+5y=2; d. ② 试举另一泛函例。 c) 物理问题中的泛函举例 ①弹性地基梁的系统势能 i)梁的弯曲应变能： x=0,固支，x=l,自由（图） ii)弹性地基贮存的能量： iii)外力位能： iv)系统总的势能： 泛函的提法：有一个梁的挠度函数，就会有一个对应的系统势能。 泛函驻值提法：在满足边界条件的所有挠度函数中，找一个w(x),使系统势能泛函取最小值。 ②最速降线问题 问题：已知空间两点A,B,A高于B，要求在两点间连接一条曲线，使得有重物从A沿此曲线自由下滑时，从A到B所需时间最短（忽略摩擦力） 作法： i)通过A，B作一垂直于水平面的平面，取坐标系如图。B点坐标（a,b）,设曲线为y=y(x),并已知：x=0,y=0;x=a,y=b-----------------------------------------③ ii)建立泛函： a. 设P(x,y)是曲线上的点，P点的速度由能量守