有限元法基础与板料成形模拟要点.ppt

板料成形有限元法—TL法 单元插值关系 单元本构关系 为时间增量步?t内第二类Kirchhoff应力增量 为时间增量步?t内Green应变增量 为切线本构矩阵 板料成形有限元法—TL法 单元几何关系 t 时刻 t+?t 时刻 板料成形有限元法—TL法 线性化的TL法有限元平衡方程 KT为t时刻的单元切线刚度矩阵 Ks为初应力单元刚度矩阵 板料成形有限元法—UL法 UL法 UL法是取当前时刻构形作为参考构形，在所有的时间步长内的计算都参照当前时刻t构形来定义。则t+?t时刻的虚功方程可以表示为 板料成形有限元法—UL法 单元插值关系 单元本构关系 为时间增量步?t内第二类Kirchhoff应力增量 为时间增量步?t内Green应变增量 为切线本构矩阵 板料成形有限元法—UL法 单元几何关系 t 时刻 t+?t 时刻 板料成形有限元法—UL法 线性化的UL法有限元平衡方程 Ks为初应力单元刚度矩阵 Rs为初应力节点力向量 板料成形有限元法—虚功率增率法 虚功原理 由于Lagrange（第一类Kirchhoff）应力与速度对物质坐标偏导数是能量共轭的，可得时刻现实构形，以初始时刻为参考构形的弹塑性大变形虚功方程为 板料成形有限元法—虚功率增率法 单元插值关系 单元本构关系 单元内任意一点的速度v与单元节点速度向量ve的关系为 板料成形有限元法—虚功率增率法 单元几何关系 变形率张量dij为 为速度梯度张量 以当前时刻为参考构形时，变形率张量dij等于应变率张量 线性化的虚功率增率法有限元平衡方程 板料成形有限元法—虚功率增率法 板料成形有限元法—逆算法 虚功方程 在最终状态上建立虚功方程 经过单元离散化 板料成形有限元法—逆算法 单元插值关系 单元本构关系 Hencky全量形变理论可得本构方程 板料成形有限元法—逆算法 单元几何关系 是B-1张量的特征值,θ是第一主应变方向 板料成形有限元法—逆算法 逆算法有限元平衡方程 板料成形有限元法—动力显式算法 板料的运动学微分方程 系统的虚功方程 板料成形有限元法—动力显式算法 单元插值关系 单元几何关系 板料成形有限元法—动力显式算法 线性化有限元方程 M是质量矩阵 C是阻尼矩阵 P是节点外力向量 F是节点内力向量 通常，动力显式积分算法采用集中质量矩阵，既M是一对角矩阵，并取 板料成形有限元法—动力显式算法 则联立方程组变成（节点数节点*自由度数）个相互独立的方程 显式时间积分的中心差分算法 板料成形有限元法—动力显式算法 将节点速度和加速度用差分格式写成 显式时间积分的中心差分算法 板料成形有限元法—动力显式算法 得到节点位移和速度的显式计算格式，前提条件是已知前两步的位移和速度 临界时间步长的确定 板料成形有限元法—动力显式算法 由于中心差分算法是条件稳定的，为了保证系统计算的稳定性，对时间增量步长的大小必须加以限制 稳定性条件通常由系统的最高频率 决定 ?是最高模态中的临界阻尼 满足稳定性条件的时间增量步长可以由膨胀波沿网格中任意单元的最小穿越时间近似得到 临界时间步长的确定 板料成形有限元法—动力显式算法 c 为膨胀波在材料中的传播速度 为第n状态单元e的名义长度 板料成形有限元法—薄膜单元 局部坐标系以1节点为坐标原点，x轴与单元1、2边重合并指向2节点，z轴与单元法向量nz平行。单位法向量nz为 板料成形有限元法—薄膜单元 随体局部坐标系oxyz与空间整体坐标系OXYZ之间坐标转换矩阵 ?为 板料成形有限元法—薄膜单元 单元局部坐标自由度向量ue与整体坐标自由度向量Ue的变换关系为 板料成形有限元法—薄膜单元 四边形薄膜单元 板料成形有限元法—薄膜单元 局部坐标系以1节点为坐标原点，x轴与单元1、2边重合并指向2节点，z轴与单元法向量nz平行。单位法向量nz为 板料成形有限元法—薄膜单元 随体局部坐标系oxyz与空间整体坐标系OXYZ之间坐标转换矩阵 ?为 板料成形有限元法—薄膜单元 单元局部坐标自由度向量ue与整体坐标自由度向量Ue的变换关系为 板料成形有限元法—单元平衡方程 单元平衡方程 单元插值关系(局部坐标系) 单元几何关系(局部坐标系) 单元本构关系(局部坐标系) 最小势能原理(局部坐标系) 坐标变换关系 板料成形有限元法—单元平衡方程 代入局部坐标系的单元3个关系式, 得 代入坐标转换关系式, 得 板料成形有限元法—单元平衡方程 ke 为局部坐标系的单元刚度矩阵 f e 为局部坐标系的单元载荷向量 Ke 为整体坐标系的单元刚度矩阵 F e 为整体坐标系的单元载荷向量 板料成形有限