[工学]数值仿真技术基础1.ppt

过程装备数值模拟与仿真技术基础（工程硕士） 绪论 第一章 直接方法及应用 第二章 变分方法及应用 第三章 加权残值法及应用 第四章 能量方法及应用 绪 论 1、什么是数值模拟与仿真技术？ 数值模拟技术 数值仿真技术 2、数值模拟与仿真的工程意义 3、本课程想达到的目的 建立有限元仿真分析的理论及应用基础，建立有限元格式的方法（四种方法） 了解相关计算机软件；便于论文写作； 利用本课程的知识进行相关CAE编程； 4、解决问题的一般步骤（7步） 5、课程的要求及其说明 第一章：体现物理直观的直接方法 1、线性弹簧系统 拉压弹簧；拉压杆 2、输运系统 导热问题 流体输送和分配 电力网络系统 单元、节点受力分析(局部分析) 解之,得到:{p}={1.1,0.98,0.73,0.24}T Q1=0.3(m3/s)[正号表示流入系统] 5、一般平面杆件系统 平面刚架（梁单元）问题_桁架是刚架问题的特例。 平面刚架问题 1、局部坐标系弯剪单元刚度矩阵 2、局部拉压单元的平衡方程 3、整体坐标系中平面刚架单元的平衡方程 7、平面杆件问题举例 举例：二维平面弹簧问题 1、结构离散化，单元、节点编号，建立连通表。 2、建立整体坐标系中单元刚度矩阵 单元1，?=0°,c=1,s=0, 单元2， ?=-45°,c=s= 根据 注意:求解得到的，作为节点自由度的位移是沿着整体坐标方向的,而单元的节点位移是沿着局部坐标方向的.也就是说,节点自由度是在整体坐标的系统方程下求出的,当需要计算杆件的应力时,需要转化成局部坐标中的位移,利用刚度计算受力,从而计算杆件应力.整体坐标下的力不能直接计算构件的应力 局部位移 整体位移 坐标变换问题_局部坐标与整体坐标 局部坐标：便于单元列式定义的坐标系统；参数一般上加一横。 整体坐标：单元集成和系统求解采用的坐标系统； i、j为节点编号；m为单元编号 m可以为任意单元（弹簧、流管、杆件、导热、平面单元等等） 变换 关系 杆件系统问题_桁架、刚架（平面，空间问题）—平面桁架 局部坐标系中桁架单元仅受拉压载荷，可借用上述承受轴向力的单元刚度矩阵表达式： 按照局部坐标方向进行扩大，即增加局部y方向的分量后，得到： 简写为： 局部坐标杆端部位移与整体坐标杆端位移的关系： 简写为： 杆件系统问题—平面桁架_杆件应力计算问题 还得到局部坐标系中杆端节点力与总体坐标系中杆端节点力间关系： 综合上述结果，得到单元平衡方程的整体坐标表达式： 其中， 是整体坐标中的单元刚度阵。 下面给出总体坐标系中，桁架单元刚度矩阵的完整表达式： 利用本节公式重新求解前述桁架问题，从局部-整体-局部就可以对问题进行完整的求解。 桁架问题举例： 1、结构离散化、单元节点编号如图所示。 2、建立单元刚度矩阵 单元1：?=0o，c=1,s=0，i=1，j=2 单元2：?=45o，c=s=1/?2，i=2，j=3 单元3：?=-90o，c=0，s=-1，i=2，j=4 3、叠加整体刚度矩阵 4、建立系统方程： 5、约束和载荷处理： 6、求解位移： 7、求解约束反力： 8、单元局部力系和变形 单元1：?=0o，c=1,s=0，i=1，j=2 单元2：?=45o，c=s=1/?2，i=2，j=3 单元3：?=-90o，c=0，s=-1，i=2，j=4 9、内力与约束力图示 9、构件应力计算 由于： 根据： 得到： 由于力在这里的方向是沿着ij，因此，表示i节点压缩应力为正向，而应力习惯符号是拉为正，压缩为负；此外，局部坐标中的二力杆，受力大小相等，方向相反，因此，可以取其中与局部坐标一致的力进行应力计算。即可以反映拉正压负的符号规定 单元1： 单元2： 注意：为表达简便期间，某些局部坐标系参数上的一横已经去掉，应注意与整体坐标系的参数意义上面的区别。 局部坐标系中的梁单元刚度矩阵 前面已经得到 待求 根据载荷与变形的线性关系，可以得到： 当 时，得到单元刚度矩阵第一列元素： 这说明单元刚度阵第1列元素是在第一个自由度变形量取值为1，其余自由度的变形量为零的条件下，单元所有节点产生的力，换言之，kij就是单元在j自由度方向产生单位广义变形时，在各个i自由度方向上产生的广义节点力。由此，可直接从刚度阵元素物理意义导出其表达式。 展开得到： 广义力和广义位移 局部坐标系中： 3、局部坐标系中平面刚架单元的平衡方程 因为坐标变换矩阵[T]是正交矩阵，所以有[T]-1=[T]T，将[T]-1左乘上式，就得到： 其中， ，就是整体坐标系中的