Matlab定积分为什么会有虚数,Mathematica求定积分，结果有虚数怎么处理？ - 程序语言 - 小木虫 - 学术 科研 互动社区...

Mathematica求积分结果有虚数怎么处理？是错了吗？怎么求解结果才没有虚数！(求大神解释怎么处理)

程序如下：

(*定义材料常数*)(*脱层深度为2，长度为0 .2*)

Subscript[\[Nu], TL] = 0.25;

Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)

Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)

Subscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)

Subscript[\[Nu], LT] =

Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];

\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];

Subscript[C, L] = Subscript[EE, L]/\[Mu];

Subscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];

Subscript[C, T] = Subscript[EE, T]/\[Mu];

Subscript[C, LT] = Subscript[G, LT];

(*脱层几何参数*)

\[Rho] = 1600;(*脱层密度*)

c = 0.1;(*阻尼*)

t = 0.01;(*t每层厚度*)

h = 0.05;(*h为梁总厚度*)

l = 1;(*为梁总的长度*)

Subscript[l, 1] = 0.35;(*为梁1区的长度*)

Subscript[l, 2] = 0.35;(*为梁2区的长度*)

Subscript[l, 3] = 0.3;(*为梁3区的长度*)

Subscript[l, 4] = 0.3;(*为梁4区的长度*)

nn1 = h/t;(*nn1为总铺层厚度*)

nn2 = 2.0;(*为上子铺层层数*)

nn3 = h/t - nn2;(*下子铺层层数*)

Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)

Subscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)

Subscript[h, 1] = h;(*梁1区的厚度*)

Subscript[h, 2] = h;(*梁2区的厚度*)

(*求各单层刚度*)

agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180,

0 Pi/180};(*列表各铺层角度*)

ck = Cos[agk];(*列表求各铺层角度余弦*)

sk = Sin[agk];(*列表求各铺层角度正弦*)

s11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +

sk^4*Subscript[C, T] +

4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)

A11 = \!\(

\*UnderoverscriptBox[\(\[Sum]\), \(k =

1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\);(*一区拉伸刚度*)

A21 = A11;

A31 = (\!\(

\*UnderoverscriptBox[\(\[Sum]\), \(k =

1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\));(*三区拉伸刚度*)

A41 = (\!\(

\*UnderoverscriptBox[\(\[Sum]\), \(k =

nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*

t\)\));(*四区拉伸刚度*)

Print[&