张量分析论文

张量分析课程报告1知识总结1.1指标符号例如,三维空间任意一点p在笛卡儿坐标系（321,,xxx），若是再推广到比三维更高的空间时不好描述了。因此，发展了另一种记法指标记法。在三维空间力里，矢量有三个分量，采用一般的指标将它们用一个简单的分量进行缩写。因此在指标记法里边用指标符号表示为（ix，i=1,2,3）。一个n维空间的矢量（nxxxx,,,,321）也可用分量表示为（nixi,,2,1,）。其中i—指标（取值范围为小于或等于n的所有正整数）n—维数1.1.1求和约定和哑指标求和约定是指标记法的补充。若在一项中，只要一个下标在同一式子中重复出现，则表示要对这个指标从1,2,3......n求和。要表示求和nnxaxaxaS2211，可表示为njjjniiixaxaS11，约定：jjiixaxaS，（用拉丁字母表示3维，希腊字母表2维）。其中求和指标与所用的字母无关指标重复只能一次。对于双重求和，3131ijjiijyxA，其中，333323321331322322221221311321121111yxAyxAyxAyxAyxAyxAyxAyxAyxAyxAjiij可表示为kjiijkzyxA，代表27项的和式。1.1.2自由指标333323213123232221211313212111bxAxAxAbxAxAxAbxAxAxA可以简写为ijijbxA，其中j——哑指标i——自由指标，在每一项中只出现一次，一个公式中必须相同1.1.3Kronecker-符号和置换符号(Ricci符号)（1）Kronecker-符号定义张量分析课程报告首先是标积，从物理学知道，一个力矢量f与一个位移矢量s,可以确定一个标量，即功W，cosWfsfsgg其中记作fs．所以又称点积。用指标符号，则当用基矢分别表示fs,时，它们的点积记为由于1e，2e，3e是相互垂直的单位矢量，由点积的定义，知当i=j时，ij的分量是1；当ij时，ij的分量为0。即013,2,1,01133132232112332211时，有当当当jijijiijji克朗内克（Kronecker）符号ij可看作是一个单位矩阵的缩写形式，即100010001333231232221131211ij当将1、2、3赋值给i时，这一点很容易被验证，于是得到的分量分别为1v，2v，3v所以ijjivv，可见，最终的结果是由于在数值变换上用i代替j。所以，显而易见，将ij应用于vj只是将vj中的j用i置换；因此ij符号通常称为置换算子。（2）置换符号(Ricci符号)等若有两个或三个指标相若若2,3,1,3,1,2,1,2,3,,2,1,3,1,3,2,3,2,1,,011kjikjieijk交错张量eijk还为缩写提供了另一种方法。例如，叉积可以写为ijkjkievwe。注意求和约定。三个矢量U，V，W的点积和叉积可以得到几种有意义的乘积形式：张量分析课程报告下面的关系式成立并且有用：以UVW,为边的平行六面体的体积或者该体积的负值，这要根据UVW,是不是构成右手坐标系而定。1.2矢量的基本运算在三维空间中,任意矢量都可以表示为三个基矢量的线性组合321,,eee，iieaeaeaeaa332211其中ai为矢量a在基矢量ei下的分解系数,也称矢量的分量（1）矢量点积可表示为ijjiee，jjiiijjijjiibababaebeaba(2)矢量叉积可表示为kijkjieeee其中kikiee，kjkjeekijktijttjsirrstjjjiiijieeeeeeeeeee321321321故可得：jiijkkkjiijkkijkjijijijjiibaeccebaeeebaeebaebeaba(3)矢量的混合积可表示为rrkjiijkecebaecba=krrjiijkcbae=kjiijkcbae其中有krijrkjieeeeee=ijkrkijree（ijke表示Ricci符号）1.3坐标变换张量分析课程报告1.4梯度、散度1.4.1标量场的梯度假定在空间某区域定义一个标量，那么可以得到分别对三个坐标的导数，即321ezeyexgrad1.4.2矢量的散度zuyuxudivzyx=ujjiijjeueu,ug是一个标量，在空间任一点，它只有一个值，不像矢量那样有三个分量。1.5笛卡尔张量1.5.1张量的概念与表示方法矢量是比标量更复杂的一种物理量或几何量。自然界还有比矢量更复杂的量，如弹性体中一点的应力状态，就有正应力x，y，z和剪应力xy，xz，yz，zx，yx，zy共九个分量值。这样一种量，叫做张量。先从矢量v的表示方法考虑，从有向线段的图示法出发，应用平行四边形合张量分析课程报告张量分析课程报告1.5.2张量的代数运算（1）加(减)法jijijijijieeTeeBABAT)(（2）矢量与张量的点积(点乘)张量的乘法，又叫张量的外积或直积。任何阶的几个张量都可施行乘法运算。其意义是第一个张量的每一个分量乘以第二个张量的每一个分量，不难证明它们组成的集合仍是一个张量，叫做原两个张量的积张量。积张量的阶数等于两相乘张量的阶数之和。矢量与张量点乘的结果仍为张量,新张量b比原张量T的阶数降低一阶左点乘beTaeeTeaTakijjkikjjkii)()(右点乘jkikijkkjiijeaTeaeeTaT)()(ceaTijijaTTa（只有对称张量两者才相等）（3）矢量与张量的叉积矢量与张量叉乘的结果仍为张量,新张量与原张量同阶左叉乘krijrjkikjjkiieeeTaeeTeaTa)()(AeeTaekrjkiijr右叉乘rjkrikijkkjiijeeeaTeaeeTaT)()(BeeaTerikijjkr（4）两个张量的点积两个张量点积的结果仍为张量。新张量的阶数是原两个张量的阶数之和减2SeeeeBAeeeeBAeeeBeeeABAtsjitkskijtskrjitrskijtsrtrskjikij)()(两个二阶张量点积的结果为一个新的二阶张量,这相当于矩阵相乘（5）张量的缩并在张量的不变性记法中,将某两个基矢量点乘,其结果是一个较原张量低二阶的新张量,这种运算称为缩并，张量的缩并是张量特有的又一个代数运算。对阶张量进行缩并，就是对其中两相同的指标按求和约定求和。不难证明，缩张量分析课程报告并之后仍是张量，其阶数比原张量低某个偶数，这要看它是对几对指标缩并而定。二阶张量的缩并，是一个标量。低于二阶的张量（如矢量）不能进行缩并运算。jiijeeAA332211AAAAAeeAAiiijijjiij（6）指标置换这是张量所特有的代数运算之一，也是最简单的张量代数运算，如kjiijkeeeAA若对该张量的分量中任意两个指标交换次序,得到一个与原张量同阶的新张量kjiijkkjijikeeeBeeeAkjiijkkjijikkijijkeeeBeeeAeeeA如果一个张量只是对某一对特定指标对称（或者斜对称），则称之为对这对指标对称的（或者斜对称）张量。如果在一个坐标系中，一个张量对某一对指标对称（或者斜对称），那么在所有的坐标系中，它对该对指标都对称（或者斜对称）。张量分析课程报告2知识应用————井筒岩体裂隙等效渗透系数张量数值法研究2.1井筒岩体淋水现状受地质、水文、施工等多种复杂因素的影响，煤矿在井筒建设之初就有淋水现象。随时间推移，淋水点越来越多，淋水量不断增大。副井井筒主要出水层位是在进入到主要含水层后开始，明显出水点集中在井壁破裂处，局部裂隙、明显出水点和井壁破裂严重的位置基本对应，大多出现在井筒的东北和西南方向，裂隙井筒岩体渗透性及其随着应力、温度的影响受到广泛关注。图2.1副井井筒内罐道梁和支护体系在淋水后锈蚀沿整个副井井深，在井壁的东北和西南出现明显的压剪破坏裂缝，东北方位裂缝多数右上到左下方向，而西北方位裂缝多数是左上到右下，呈现出与地层侏罗系地层X型压剪共轭裂缝相对应的特征。随着井筒深度增加，裂缝倾角（初始近70）有逐步变缓趋势。井壁裂缝沿井深的分布广，井深50m以下直到马头门上方均布。图2.2井深80~200m副井井壁混凝土破坏段主要位于马头门向上100m左右范围，从地质柱状图可知，这一范围恰好处于1煤至5煤的含煤地层。从2012年4月检查情况看，局部区域的井壁厚度甚至不足15cm，井壁破坏严重，井壁后煤层清晰可见。由于煤层的自身承载力能力弱，且具有遇水膨胀的特性，井壁的淋水渗透到煤层引起向井筒内的膨胀变形进一步加剧了井壁破坏。为研究岩体温度与裂隙渗透性变化关系，进行为期3个月的静水压力及温张量分析课程报告度对含单裂隙石英岩渗透性实验，并分析其变化机制开展了不同温度及应力作用下人工裂隙渗流试验，提出了岩样裂隙结构面温度-应力-水力耦合本构关系式。2.2井筒岩体裂隙渗透特性试验研究（1）试验装置及试验方案首先开展了大理岩人工裂隙渗透率随应力及温度变化的试验研究，获得加卸载过程中裂隙渗透率的演化规律；其次通过数值方法研究了某裂隙岩体特定张开度条件下其等效渗透系数的尺寸效应及各向异性，获得了该裂隙岩体的等效渗透系数REV及渗透张量。试样取井筒中部大理岩，将试样加工成直径49.07mm，高80mm圆柱体后，采用巴西劈裂法制作人工裂隙，制成后的裂隙试样见图2.3，其中裂隙试样为一个整体，裂隙闭合，人力难以将裂隙两侧岩块分开。图2.3试样巴西劈裂试验形成的裂缝（2）试验结果分析表2.1不同静水压力试样D-C10的渗透试验结果记整体坐标系下的渗透系数矩阵为K，局部坐标系下的渗透系数矩阵为K＇，OXY逆时针旋转角度为oxy，x轴与X，Y轴的方向余弦分别记为1cosl，1sinm，y轴与X，Y轴的方向余弦分别记为2sinl，2cosm。张量分析课程报告K和K的分量之间满足坐标转换关系：（1）其中1122ijlmlm为整体坐标系和局部坐标系的转换系数。如果以主渗透系数1k，2k的方向为整体坐标系的X轴和Y轴的方向，则逆时针旋转度的局部坐标系x'轴方向的渗透系数为（2）根据椭圆方程的极坐标形式，有（3）对比式（2），令121ka，221kb，21kl。则可以看出渗透系数张量的分量的某种形式可以用椭圆来描述，该椭圆称为渗透椭圆（见图2.4），渗透椭圆的主轴分别为11k，21k，沿着角度方位径长为1k。图2.4渗透椭圆示意图图2.5为根据渗透椭圆方程（3）拟合出来的椭圆，其中沿着角度方位径张量分析课程报告长为1k，根据椭圆拟合结果得到渗透张量主值和主方向，最大渗最大渗透系数为6.1×10-10m/s，方向沿着x轴逆时针旋转15.4o，最小渗透系数为了3.22×10-10m/s方向沿着y轴逆时针旋转15.4o。利用式（1）求得当前坐标系下渗透张量的各个分量，等效渗透系数张量为105.90.74100.743.42K（m/s）最大渗最大渗透系数为6.1×10-10m/s，方向沿着x轴逆时针旋转15.4o，最小渗透系数为了3.22×10-10m/s方向沿着y轴逆时针旋转15.4o。图2.6为渗透系数张量分量xxk的数值计算结果和渗透系数张量计算结果对比。图2.5渗透椭圆拟合曲线图2.6分量xxk各个方向拟合值和数值计算值比较张量分析课程报告（3）结论首先进行了大理岩