1. 背景介绍

应力-应变曲线能够有效地帮助我们分析材料的力学性能，应力指物体单位面积上所受到力的大小，应变指物体的相对变形量。正应变描述的是沿着某一方向上长度的变化量，而切应变描述的是角度的变化量。
根据力的大小可根据胡克定律理论上计算得到任意位置应变。对于现代的测量方式而言，主要有引伸计、应变片、光学测量等方式，光学测量因非接触、精度高的特点，在工业测量领域受到高度关注。

2 变形梯度张量

光学测量使用数字图像相关法(Digital Image Correlation, DIC)获得被测物体表面的变形情况，因此，计算的是表面应变(Surface Strain)。在连续介质力学中，描述物体运动的方式有两种，即Lagrange描述与Euler描述。前者以不同时刻下物体的位置为研究对象，而后者以固定空间内的变化为对象。本文中我们在Lagrange描述下进行分析。

上图中，参考位置PQ⃗\vec{PQ}PQ，在ttt刻为pq⃗\vec{pq}pq，P和Q的位移量分别为u(X)u(\bf{X})u(X),u(X+dX)u({\bf{X}}+d\bf{X})u(X+dX)，该过程的运动描述为x=χ(X,t){\bf{x}}=\chi({\bf{X}},t)x=χ(X,t)。根据应变的定义可知，要想得到应变，就需要建立dxd{\bf{x}}dx与dXd{\bf{X}}dX之间的关系，根据图示，很容易能够得到dx=dX+dud{\bf{x}}=d{\bf{X}}+d{\bf{u}}dx=dX+du，进一步可知

其中，dud{\bf{u}}du是相对位移，F\bf{F}F为变形梯度张量(Deformation Gradient Tensor)。根据运动描述x=χ(X,t){\bf{x}}=\chi({\bf{X}},t)x=χ(X,t)的全微分形式

由此前两式可知，变形梯度张量F\bf{F}F实际上是该物体运动描述函数的梯度。dxd{\bf{x}}dx如果理解为变形后pq⃗\vec{pq}pq的长度，那么变形前长度就是dXd{\bf{X}}dX，因此，只要能够得到F\bf{F}F，就能计算出应变大小。关于变形梯度张量F\bf{F}F的计算，若是已知运动方程，则只需要根据方程，求F=[∂x/∂X∂x/∂Y∂y/∂X∂y/∂Y]{\bf{F}}=\begin{bmatrix} \partial{x}/\partial{X} & \partial{x}/\partial{Y} \\ \partial{y}/\partial{X} & \partial{y}/\partial{Y} \\ \end{bmatrix}F=[∂x/∂X∂y/∂X∂x/∂Y∂y/∂Y]的偏导数矩阵即可，如下所示的计算例子。
前述内容是关于变形梯度张量的代数意义，而它在几何上的意义可以用公式表示为F=RU\bf{F=RU}F=RU它表示的是连续介质刚体旋转与变形的过程，如下图所示，立方体在变形矩阵U\bf{U}U作用下变形为长方体，其后做刚体旋转R\bf{R}R，最终为ttt时刻状态。或者也可以先旋转再变形，F=VR\bf{F=VR}F=VR。

3. 表面应变计方法一

实际情况运动方程是未知的，那应该如何求解F\bf{F}F。根据前面的式子可知，F=I+∇xu{\bf{F}}={\bf{I}}+\nabla_x{\bf{u}}F=I+∇xu如果能得到相对位移的梯度，那就能得到F\bf{F}F。使用DIC设备匹配散斑图像，重建能够得到匹配点的三维坐标，同时计算得到位移。但是，这些这些数据都是离散的，如何来求相对位移的梯度呢？常用的方法就是拟合，根据这些离散点的坐标为自变量，位移为因变量，拟合得到函数u(X,Y)u(X,Y)u(X,Y),v(X,Y)v(X,Y)v(X,Y)。然后，求函数梯度，即可得到F\bf{F}F。根据F\bf{F}F，计算得到格林-拉格朗日应变张量(Green-Lagrange tensor)，

3.1 局部坐标系建立

DIC结果得到的实际上是一个没有厚度的面片，Surface strain或者Plain strain就是说，应变计算是一个沿着表面的、二维的度量。利用3×33\times33×3，7×77\times77×7或者其他局部范围内的点来构造一个平面，以此建立局部坐标系，有两种方式：一是这些点的拟合平面；二是在局部中心点处的切平面。
采用构造切平面的方式，使用局部3×33\times33×3内点坐标，拟合二次曲面z=f(x,y)z=f(x,y)z=f(x,y)，法向量为

切向量有无数条，关于切向量的选取，只需要固定xxx或yyy为一常数，从而确定一个切向量。它的选取最好能与沿着物体变形的方向。

根据向量叉乘得到另一个坐标轴，由此，获得n⃗\vec{n}n, t⃗\vec{t}t, b⃗\vec{b}b单位向量构建的局部坐标系

3.2 位移函数拟合

值得注意的是，每个点的位移量值不是将点投影到局部坐标系，然后做差。这种做法只适合于小弯曲变形情况，在Andresen 1999年撰写的论文中就有提到在大弯曲变形情况下，这种做法准确性很低。Surface strain, 是沿着表面方向发生的变形量，准确的来讲，应该使用弧长而不是位移向量，如图所示。换句话说，相当于把面展开成平面。，准确的方法应该是在上一步建立的局部坐标下，将重建的三维点，根据拟合的二次曲面z=f(x,y)z=f(x,y)z=f(x,y)沿展到tobtobtob面内，最后在tobtobtob面内得到每一个点的“位移量”。

假设点的位移模式为如下二次形式

现在已知每一个点在tobtobtob面内的坐标(x1,x2)(x_1,x_2)(x1,x2)，以及该点的在x1x_1x1方向的“位移”u(x1,x2)u(x_1,x_2)u(x1,x2)，如果选取的是3×33\times33×3的点，9个方程，6个未知系数，利用最小二乘法求解，得到函数的表达式。

3.3 位移函数偏导

由前述得到位移函数的表达式之后，也就能计算得到∂u/∂x1=β2\partial u/\partial x_1 = \beta_2∂u/∂x1=β2∂u/∂x2=β3\partial u/\partial x_2 = \beta_3∂u/∂x2=β3同理得到另一个方向的偏导数，带入上述变形梯度张量F\bf{F}F的计算中，最终求解格林应变张量。根据3.1~3.3的步骤，遍历DIC重建得到的每一个点，由此可以获得物体表面变形的应变分布。

4. 表面应变计算方法二

通过第3节中的计算过程，其实可以看出来，连续变形过程的应变计算本质上讲就是位移的变化率，用公式来说就是strains=∇displacementsstrains = \nabla displacementsstrains=∇displacements这里通过搜集的资料，给大家分享第二种表面应变计算方法，其本质没有变，只是策略是利用了有限元理论中，三角形单元的应变计算。
既然DIC设备重建得到的是点的信息，那就可以将这些点连接为一个个的三角形，三角形内任意点的位移可以根据如下位移模式得到

方程中含有6个未知数，三角形的每一个定点坐标与位移量大小都是DIC获取的，每一点能够建立两个方程，那么3个点正好6个方程，可以求解a1,a2,...,a6a_1,a_2,...,a_6a1,a2,...,a6。同样地，就可以得到位移的梯度
∂u/∂x=a2,∂u/∂y=a3\partial u/\partial x = a_2, \partial u/\partial y = a_3∂u/∂x=a2,∂u/∂y=a3∂v/∂x=a5,∂v/∂y=a6\partial v/\partial x = a_5, \partial v/\partial y = a_6∂v/∂x=a5,∂v/∂y=a6代入变形梯度张量公式与格林应变公式，即可得到结果。三角形单元是一次形式的单元，单元内任意一点的应变是一个常数，那么在四个三角形中心点处的应变，需要求平均值。
同样地，遍历所有点，构造周围四个三角形，得到该点格林应变。由上述过程可以看出，该方法的前提就是获取点之间的间距必须足够小，否则对于大变形的情况精度会很差，而方法一的误差则主要来自于拟合误差。

表面应变计算（Surface strain calculation）相关推荐

GEO-BYB型振弦式表面应变计的主要作用
振弦式表面应变计是适用于长期布设在水工结构物或其它结构物的表面,测量结构物表面的应变量,并可同步测量布设点的温度的机器. GEO-BYB型振弦式应变计加装配套附件可组成表面应变计.钢板计.无应力计等多 ...
matlab分析应力应变图,应变分析 (strain analysis).ppt
应变分析 (strain analysis) 第三章应变分析 (strain analysis) 3.1 几何方程(geometry equation) 第三章应变分析 (strain analy ...
牛顿迭代法进行低周疲劳应力应变计算 C++编程函数
参考网上的牛顿迭代法C++程序写了两个计算函数输入初始迭代值得到结果计算公式可以根据实际情况更改 #include<stdio.h> #include<math.h> ...
晶圆测试软件,半导体晶圆接触角测量仪 Wafer水滴测试
SDC-600半导体晶圆接触角测量仪 Wafer水滴测试的详细资料: 半导体晶圆接触角测量仪 Wafer水滴测试一·概述(Summary) 接触角(contact angle)是指在气.液.固三相交 ...
给力干货！一文教你如何通过Jade计算晶粒尺寸、微观应变和结晶度
当晶粒尺寸小于100nm或者样品中存在微观应变时都会引起衍射峰的宽化.样品衍射峰的宽化(FW(S)).衍射峰的并高宽(FWHM)和衍射峰宽度(FW(I))之间存在如下的关系: 式中,D称为反卷积参数, ...
fluent瞬态计算终止条件在哪里设置_fluent模拟表面喷涂
幸运草喷涂是通过喷枪或碟式雾化器,借助于压力或离心力,分散成均匀而微细的雾滴,施涂于被涂物表面的涂装方法. 这一篇文章我们一起学习一下fluent对喷涂过程的仿真计算. 模型是我在网上下载的一个手机 ...
matlab求三维图形表面积,基于MATLAB的不规则平面立体表面积计算
·12· 设计与研究机械 2012 年第 2 期总第 39 卷 --------------- 收稿日期:2011-06-10 作者简介:和蕴锋(1988-),女,河南漯河人,硕士研究生,主要研究 ...
计算声学·物理基础篇1
计算声学 @rz_wen 第一部分·物理基础流体力学符号说明 c i j k i , j , k = 1 , 2 , 3 c_{ijk} \qquad i,j,k=1,2,3 cijki,j,k ...
Science评论：量子计算目前最大的挑战，在0和1之间
2020-09-18 14:20 导语:前路漫漫~ 作者 | 青暮编辑 | 丛末 2019年10月,谷歌的一项关于量子计算的研究登上了Nature封面.谷歌声称用53个量子比特的量子计算机 ...

表面应变计算（Surface strain calculation）