弹塑性有限元法基本理论及模拟方法.pptVIP

4.1 非线性问题及分类 在分析线性弹性问题时，假定： 应力应变线性关系 结构位移很小（变形远小于物体的几何尺寸） 加载时边界条件的性质不变 如果不满足上述条件之一，就称为非线性问题 非线性结构的基本特征：变化的结构刚度 非线性问题可以分为三类： 材料非线性：体系的非线性由材料的应力应变关系的非线性引起。 如金属变形弹塑性行为、橡胶的超弹性行为等 几何非线性：结构的位移使体系的受力状态发生了显著的变化。 如板壳的大挠度问题 ——平衡方程必须建立于变形后的状态 接触非线性：接触状态的变化所引起。 如金属成形、跌落试验、多零件装配体等 接触非线性例子 随着有限元算法理论、计算机硬件和软件技术的进步及实际工业的需求，CAE技术的应用逐步由线性模拟为主向非线性模拟为主快速发展。 1969年，第一个商业非线性有限元程序——Marc诞生。 目前几乎所有的商业有限元软件都具备较强的非线性问题的分析求解能力。 非线性求解技术的先进性与稳健性已经成为衡量一个结构分析程序优劣的标准。 非线性问题的有限元求解方法 非线性方程(组)的求解方法 直接迭代法 Newton-Raphson迭代法 修正的Newton-Raphson迭代法 非线性问题通常采用增量法求解（追踪加载过程中应力和变形的演变历史。） 每个增量步采用Newton-Raphson迭代法 非线性方程的迭代求解方法 非线性方程组的迭代求解方法 非线性问题的增量法求解过程 (1) 将总的外力载荷分为一系列载荷段 (2) 在每一载荷段中进行迭代，直至收敛 (3) 所有载荷段循环，并将结果进行累加 (1) 将总的外力载荷分为一系列载荷段 (2) 在每一载荷段中进行迭代，直至收敛 4.2 材料非线性问题及分类 概念：由于材料的应力应变非线性关系引起的非线性。 分类： 不依赖时间的弹、塑性问题 非线性弹性——橡胶 弹塑性——冲压成形 依赖于时间的粘（弹、塑）性问题 蠕变——载荷不变，变形随时间继续变化 松弛——变形不变，应力随时间衰减 非线性弹性材料行为 弹塑性材料进入塑性的特征：载荷卸去后存在不可恢复的永久变形。 应力应变之间不是单值对应关系，与加载历史有关。 单轴应力状态下弹塑性材料行为 单轴（一维）应力状态下材料的应力应变行为可以从拉伸试验中获得。 单调加载 反向加载 一般应力状态下弹塑性材料行为 屈服准则（初始屈服条件） 硬化法则 （后继屈服函数、加载函数、加载曲面） 流动法则 加载、卸载准则 屈服准则（初始屈服条件） 在单向受力情况下，当应力达到材料的屈服强度时材料开始产生塑性变形。 对于一般复杂的应力状态，应力状态由六个应力分量决定时，显然不能根据某个单独应力分量的数值作为判断材料是否进入塑性变形的标准。为此，引入以应力分量为坐标的应力空间，根据代表不同应力路径的实验结果，可以定出从弹性阶段进入塑性阶段的各个界限，即屈服应力点。在应力空间中，这些屈服应力点形成一个区分弹性和塑性的分界面——屈服面。描述这个屈服面的数学表达式就是我们所要寻求的一般应力状态下的屈服准则。 硬化法则 塑性硬化法则规定了材料进入塑性变形后的后继屈服函数（又称加载函数或加载曲面） 各向同性硬化 运动硬化 混合硬化 流动法则 加载、卸载准则 4.3几何非线性问题及分类 概念：由于大位移、大转动而引起的非线性。 分类： 大位移、大转动、小应变问题 ——板壳的大挠度和后屈曲 大位移、大转动、大应变问题 ——薄板成形、弹性材料的受力 比较：线弹性 —几何非线性 线弹性：小变形假设——假定物体发生的位移远小于物体本身的几何尺寸，应变远小于1。建立平衡方程时不考虑物体位置和形状的变化。 几何非线性：物体发生有限变形——大位移、大转动的情况。建立平衡方程时必须考虑物体位置和形状的变化。 由上式可以看出，初载荷{R}不仅与加载前的应力有关，而且与该次加载引起的应变增量有关，即式： 的两边都含未知数。因此每此加载时，必须利用迭代法来求解。 具体迭代过程：可先取 ，则得 ， 求得初载矢量 ， 在载荷 下，解方程式： 此时方程可写成： (4-69) 同理可推得 的表达式 : 写成矩阵的形式，就是： 按照上述同样方法可得： 式中： 在塑性区： (4-39) 公式（4-39）的具体推导： 所以： 由普兰特尔—罗伊斯关系有： 写成矩阵形式为： 仿照前面，不难推得： 而： 其中： 对于平面应变问题，有 只需从式4-12 和式4-36中消去上述为零的分量，就可得到下列各式。 (4-40)