﻿功能:求力学里面应变(应力)矩阵的第二不变量

格式:MechI2(A)

说明:A是3*3的对称矩阵(本函数不检测你提供的矩阵是否对称).有力学基础的都知道,对于某点的应力(应变)矩阵,它对应着第一、第二、第三不变量.第一不变量,我们可以直接通过{Trace<矩阵运算\Trace>}函数就可以获取,对于第三不变量,我们直接使用函数{Det<矩阵运算\Det>}求方阵A的行列式即可.本函数就是求第二不变量的.

例子:

a =
[ 1.00000000000000     0.00000000000000     6.00000000000000
  0.00000000000000     5.80000000000000     3.00000000000000
  6.00000000000000     3.00000000000000     2.80000000000000 ]
mechi2(a)//执行命令,求a矩阵的第二不变量
ans =
[ -20.160000000000 ]

//By2012/03/01