应变的计算方法

1、应变的计算方法本章介绍了几种网格应变的计算方法，通过分析网格变形的特点及规律，将 网格的变形分解为分别沿两个主应变的方向一次变形而得， 从而通过欧拉法推导 了有限应变解析的方网格应变计算方法， 并把三维空间网格的每个网格作为线性 孔斯曲面介绍了三维空间网格的应变计算方法。 此外还介绍了工程应变、 等效应 变和厚度的计算。4.2 基于欧拉法和有限应变理论解析的方网格计算方法根据有限应变的理论， 不同的应力加载可以获得相同的应变结果。 对于近似 于平面应力状态的板材成形来说， 每个单元体的应变主方向 (除去因为位移造成 的转动)在成形过程中保持不变。 这样就可以将应变分成不同的加载阶段， 利用 真

2、实应变的可叠加性，就可以推导出方网格变形的应变计算方法。连续体的有限变形有两种表述方法。 一种方法的相对位移计算是以变形前后 物体内一点作为参考点， 即以变形前的坐标作为自变量， 这种方法称为拉格朗日 法。另一种方法的相对位移计算是以变形后物体内一点作为参考点， 以及已变形 后的坐标作为自变量， 这种方法称为欧拉法 48 。这里给出基于欧拉法和有限应 变理论解析的方网格计算原理。4.2.1 方网格内部的变形 设任意方向正方形网格内接于圆网格， 将其变形过程分解为两个阶段， 如图4-5所示。第一个阶段沿着X方向变形，丫方向保持不变；第二个阶段沿着 丫方 向变形，X方向保持不变，即应变主方向与坐标

3、轴相平行。变形的结果使圆网格 变形为椭圆， 正方形网格变形为平行四边形 (假设单元网格内沿主应变方向的变 形是均匀的)(a) 初始网格 (b) 横向变形后的网格 (c) 纵向变形后的网格图 4-5 基于有限应变的网格分解变形过程4.2.2 应变主方向和真实应变的计算对于方网格中心的应变， 假设网格内部变形是均匀的， 所以变形前后四边形 对角线的交点就是网格中心， 对角线把方网格划分成四个三角形。 将变形后的网 格中心和变形前的网格中心重合，建立直角坐标系，如图 4-6 所示。图 4-6 以欧拉法建立的变形前后网格中心重合的坐标系统 根据欧拉方法，以变形之后的网格坐标来分析， 将主应变方向定为坐

4、标方向， 设X方向为主应变的方向，丫方向为主应变的方向，两个方向分别有拉形比： (4-20) 则两个方向的真实应变等于两次分别变形的叠加：（4-21） 设变形前方网格边长为，为所取初始三角形的直角边长，则有：取其中初始三角形，其变形后为，根据变形后的网格点坐标、，得到变形 后三角形边长为：（4-22） 沿两个主应变方向的拉形比为：（4-23） 已知：（4-24） 得：（4-25） 由此得到根据三角形计算出来的主应变的方向，进而可以求出主应变：（4-26） 根据四边形网格划分的三角形分别求出来的主应变的方向和大小， 就得到了 方网格中心0点的真实应变值。由于进行多次计算，四边形的网格都得到了利用

5、， 平均之后，计算的精度得以提高，减小了误差。4.2.3 网格点上的应变 以上应变的计算获得的都是方网格中心的应变值，对于网格点上的应变值， 则三角形的三个顶点都要取网格点才能计算出网格点上的应变值。图4-7中所示有9个网格点19构成四个网格四边形，A、B、C、D分别为 四个网格的中心。通过四个网格（ 1、 2、 4、 5），（2、 3、 5、 6），（4、 5、 7、 8）和（5、6、& 9）可以分别求出中心点 A B C D的应变值。那么网格点5 的应变值的获得有下面几种方法：一种是利用（ 1、3、7、9）四个点构成的四边形利用上小节所述的方法进行 计算，分别求出三角形（ 1、5、

6、3），（ 1、5、7），（ 3、5、9），（ 7、5、9） 的直角顶点 5 的应变值，然后再求平均值，从而获得主应变和主方向。如图 4-8 所示。图 4-7 网格点上应变的获得图 4-8 由斜侧方向点求网格点上应变图 4-9 由纵横方向点求网格点上应变另一种方法就是利用（ 2、4、8、6）四个点构成的四边形同样利用上小节所 述的方法进行计算，分别求出三角形（ 2、5、4），（ 2、5、6），（ 4、5、8） ， （ 8、5、6）的直角顶点 5 的应变值，然后再求平均值，从而获得主应变和主方 向。如图 4-9 所示。以上两种方法都仅利用了网格点 5周围的 4个网格点，而没有充分利用网格 点 5

7、四周的 8 八个网格点，因此从计算精度上来说，前一种方法在45°方向上精度比较高，而后一种方法在 90°方向上精度比较高。为了进一步提高计算精 度，可以将以上两种方法结合起来， 再一次进行平均， 由此获得由某个网格点连 同周围 8 个点计算出来的网格点应变值。 ?显然，这种方法所带来的问题就是需 要的计算时间相应地要增加很多。4.3 其它应变4.3.1 工程应变 工程应变虽然不具有叠加性质，但它比较直观，在工程实际中应用广泛。 如图 4-10 所示，是材料的原始长度，是材料的伸长量。则工程应变为伸长 量相对于原始长度的比值：(4-27)图 4-10 材料伸长变形 真实应变是

8、伸长比(拉形比)的自然对数，可以表示为：(4-28) 因此可以由真实应变获得工程应变：(4-29) 工程应变值如果为正，可以解释为伸长率，如果为负值，则为收缩率，而对 于厚度工程应变， 则可以解释为减薄率。 对于拉延成形， 材料变薄情况是最需要 关注的，用减薄率来反映材料的变形就非常直观。相应地，真实应变也可以表示为：(4-30)4.3.2 等效应变 等效应变也可以由真实最大主应变和最小主应变获得：(4-31)4.3.3 厚向应变 由体积不变条件可以得到真实厚向应变：(4-32) 材料的真实厚度可由真实应变获得：(4-33)4.3.4 增量应变 对于板材成形工序次数较多的零件，工序之间的应变变

9、化也是需要关注的。从上一次工序到下一道工序之间的变形的增量应变的计算方法和平面应变的计 算方法相似。但在这种情况下， 网格已经发生了变形， 原先变形的网格和再次变 形的网格如图 4-11 所示。图 4-11 增量变形前后三角形网格的形状文献 49 介绍一种通过坐标变换， 将变形前后的三角形的某个边重合而获得 的增量应变三角形节点计算方法：有四个张量 50,51 ：变形梯度张量：的转置：柯西- 格林变形张量：拉格朗日应变张量变形梯度张量 F 是从初始三角形的坐标到再次变形的三角形的线性变换， 有：(4-34) 为了简化，将初始三角形和变形后的三角形通过旋转是其中一个边和轴重 合，则有，另有变形梯

10、度张量 F 为：(4-35)柯西- 格林变形张量可以表达为：(4-36) 代入前式，有(4-37) 则拉格朗日应变就可以通过柯西 - 格林变形张量计算获得：其中是单位矩阵，拉格朗日应变可以以如下形式给出：(4-38) 有：(4-39)真实应变则可表达为：(4-40)则主应变的方位角为：(4-41)4.3.5 厚度和曲率修正 应变可以通过节点之间的距离计算出来。 初始的网格距离就是网格点之间的 距离，而在两个节点之间的新的距离在以上的计算过程中都是当作直线来计算的 的。实际上工件表面是有曲率的并且是变化的， 因此对于曲率形状变化剧烈的工 件，采用直线距离作为曲面上的距离所带的误差就比较大了。同样

11、，通过图像处理和三维重建的网格点都位于板料的表面， 计算的仅仅是 板料表面的应变值， 而板料是有一定的厚度的， 因此有必要进行厚度和曲率修正， 文献 49 介绍了一种进行曲率和半径修正的方法为： 当节点处曲率半径为， 初始 节点之间的距离为，厚度为，则修正的公式为：(4-42)在应变分析实际应用中， 为了获得比较高精度的测量结果或者为了准确地测 量应变变化梯度较大的区域， 如曲率半径比较小的区域， 采用方法使用比较小直 径的网格，能够获得更准确的结果。4.4 三维空间网格的应变计算 前面的描述都是基于平面应变的假设基础上， 实际上，在冲压后， 工件的形 状是三维立体的，印制在平面板料上网格经过

12、变形后变成了三维空间立体网格。 每个网格节点都具有三维坐标，网格应变的计算就必须考虑空间坐标的问题。平面正方形网格经过变形后变成了空间网格， 网格中心也从点变到了点。 如 图 4-12 所示(斜侧图)。图 4-12 平面正方形网格变形为空间网格 由于网格中心经过变形后不再是空间四边形的几何中心了， 因此在计算网格 应变时需要通过四个边长的变化确定网格中心点的坐标。 为了简化计算， 参照线 性孔斯曲面的定义我们可以获得四点构成的曲面片的方程， 从而求得变形后空间 网格中心点的坐标。4.4.1 线性 Coons 曲面1964年 S.A.Coons 提出了一种曲面分片、拼合造型的思想，他用四条边界

13、构造曲面片并通过叠加修正曲面片，产生满足用户需要的曲面 52 。线性 Coons 曲面，也称之为简单曲面， 是通过四条边界曲线构成的曲面。 若 给定四条边界曲线，且。在向进行线性插值，得到直纹面为：(4-43)在向进行线性插值，得到直纹面为：(4-44)如图4-13所示：图4-13 线性Coo ns曲面的生成则用四条边界曲线构造的曲面可用矩阵形式表示为：(4-45)4.4.2基于线性Coo ns曲面的四点空间网格应变的计算设：已知变形后空间网格四点的坐标为，则可将变形后的空间网格看 成四条边界曲线都是直线的线性 Coo ns曲面。如图4-14所示。图4-14 变形后的空间网格根据线性Coons

14、曲面的描述和网格区域内变形是均匀的假设前提，有变形后网格中心点坐标：(4-46)根据4.1.3所述的方法就可以计算出点的应变值。以三角形为例，有应变主方向矢量落在三角形所在的平面内，如图4-15所示：图4-15空间坐标三角形单元的计算(4-47)(4-48)有：(4-49)同样，其他三角形可以得到相类似的结果。由于四个三角形不共面，根据每个三角形计算出来的主应变矢量方向分别落在各个三角形所在的平面内，而一点的主应变方向及大小是唯一的，因此，需要对应变主矢进行平均求取网格中心点 的切向应变。但是由于考虑到计算精度和计算效率的问题，本文在实现时只计算一个三角形内部的应变。应力应变的切线模量和割线模

15、量1、初始切线模量应力值为零时的应力-应变曲线的正切2、切线模量一一某一应力级位处应力-应变曲线的斜率.3、 割线模量一一以某一应力值对应的曲线上的点同起始点相连的割线的斜率4、回弹模量一一应力卸除阶段应力-应变曲线的割线模量前三种模量取值时的应变值是包含残余应变和.回弹应变在内的总应变，而回弹 模量取值时已扣除残余应变后的回弹就变。 因此，将前三种模量笼统地称为土的 弹性模量显然是不合适的。而回弹模量能反映土所具有的那部分弹性性质， 所以, 在以弹性力学为理论基础的路面设计方法中， 往往将土的回弹模量视为土的弹性 模量，并且作为路面设计中的一项重要计算参数。取混凝土应力应变曲线在原点 0切线

16、的斜率，作为混凝土的初始弹性模量， 简称弹性模量Ec,即：Ec = tg a 0Ec 初始弹性模量；a0原点切线的斜率夹角。当应力较大时，混凝土已进入弹塑性阶段，弹性模量已不能正确反映此时的应力 应变关系。比较精确的方法采用切线模量 Ec',即在应力应变曲线任一点处作一 切线。此切线的斜率即为该点的切线模量，其表达式为Ec' = tg a = d (T / d £割线模量是原点与某点连线即割线的斜率作为混凝土的割线模量，称为变形模量Ec，它的表达式为Ec" = tg a 1 =(T c /£ c它们之间的大小关系为：原点切线模量 割线弹性模量 切线

17、弹性模量。原 点弹性模量也就是我们平常用的混凝土弹性模量，值大约为 30000MPa左右， 与混凝土强度等级有关。土基的力学强度特性及其设计参数土的非线性特性土的非线性特性土基的模量1）初始切线模量一一应力值为零时的应力-应变曲线的正切，如图所示，代表加荷开始 时土的应力-应变关系。2）切线模量一一某一应力级位处应力-应变曲线的斜率，如图所示，反映土在该级位 应力-应变变化的精确关系。3）割线模量一一以某一应力值对应的曲线上的点同起始点相连的割线的斜率，如图所 示，反映在该应力级范围内的应力-应变关系的平均情况。4）回弹模量一一应力卸除阶段应力-应变曲线的割线模量，如图所示，反映土在回弹变 形

18、范围内的应力-应变关系的平均情况。土基的流变性质土的变形随时间变化的关系。土在荷载作用下的变形不仅与荷载大小有关，而且还与荷载作用的持续时间有关，是一种具有流变性质的材料。 土基回弹模量测试回弹模量能较好地反映土基所具有的部分弹性性质，可以用回弹模量表示土基在瞬时荷载作用下的可恢复变形性质。我国公路水泥混凝土路面、沥青路面设计方法都以回弹模 量E作为土基的刚度指标。测定时宜采用逐级加载一一卸载法（直径30.4cm的板）。每一级荷载经过加载和卸载， 取得稳定的回弹弯沉之后，再加下一级荷载，如此施加n级荷载后，即可点绘出荷载 -弯沉曲线。在多数情况下，试验曲线呈非线性。在确定模量时，可以根据土基实

19、际受的压力范围或 可能产生的弯沉范围在曲线上取值。路面设计中，按1mm线性归纳法来确定土基的回弹模量。a）柔性承载板 b）刚性承载板地基反应模量文克勒地基模型是原捷克斯洛伐克工程师文克勒 （Win kler）1876年提出的， 其基本假定是地基上任一点的弯沉仅与作用于该点的压力 p成正比，而与相邻点 处的压力无关。直径76cm的刚性板测定。当地基较软弱时，取l=0.127cm时相对应的压力 p计算地基反应模量；当地基较为坚硬时，取单位压力 p=0.07MPa时相对应的弯 沉值I计算地基反应模量。加州承载比CBR加州承载比CBR是美国加利福尼亚州提出的一种评定基层材料承载能力的试 验方法。承载能

20、力以材料抵抗局部荷载压入变形的能力表征，并采用标准碎石的承载能力为标准，以相对值的百分数表示CBR值。土基的设计参数的确定我国在测定土基回弹模量时，常采用直径 30.4cm的刚性承载板用加载一卸 载的试验方法。试验通常在不利时期进行，并取有 84.1 %概率的回弹模量值作 为土基回弹模量的计算值。规范给出我国土基回弹模量设计参数选用的建议值。土基回弹模量与CBR的关系一直是世界各国在路基土研究中比较关心的内 容。根据试验给出了国内外部分土基回弹模量与 CBR勺关系，设计是可以根据实 际参考选用。路面材料的几种强度（1）主要材料类型 松散颗粒型材料及块料 沥青结合类材料 无机结合料类材料 水泥混

21、凝土材料由于材料(整体性材料和非整体性材料)的基本性质和成型方式的不同， 各种路面结构具有不同的力学强度特性(即应力 -应变关系)，也使得路面具有 不同的使用品质和使用寿命。(2) 抗剪强度摩尔一库仑强度理论：匚口山屮其中c和V是表征路面材料抗剪强度的两项参数，c是材料的粘结力(kpa)，f' 是材料内摩阻角，对于土可以通过直剪试验得到；对于松散粒料无法做直接剪切 试验，可用三轴压缩试验测定。图2-11三轴试验确定c. ©(3) 抗压强度指试样在无侧向压力条件下，抵抗轴向压力的极限应力。材料经过标准成型 和养生后通过无侧限抗压试验测定的强度。一般有大试件、中试件和小试件。小试件指50mm*50m；m中试件指 100mm*100mm大试件指 150mm*150m>m(4) 抗拉强度气温变化会引起路面材料收缩，湿度变化能产生半刚性材料干缩，当收缩变 形受到约束，即在材料内产生拉应力，材料抗拉强度不足即可引起路面结构拉伸 断裂。路面材料的抗拉强度主要由混合料中的结合料粘结力提供，可采用直接拉伸或间接拉伸试验测定材料的抗拉强度。(5) 抗弯拉强度路面材料的实际工作状态是弯曲反复变化的， 结构层底首先容易出现局 部拉裂、产生弯曲断裂。材料的抗弯