游龙游龙 发表于 2014-7-7 17:15

麻烦大家帮忙看下mathematical编的程序，不知道哪里出错了

本帖最后由 游龙游龙 于 2014-7-8 14:45 编辑

麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
行列式求不出，画图找不到根
(*定义材料常数*)
Subscript[\, TL] = 0.25;
Subscript = 53.8*10^9;(*L方向的弹性模量*)
Subscript = 17.93*10^9;(*T方向的弹性模量*)
Subscript = 8.96*10^9;(*LT方向的剪切模量*)
Subscript[\, LT] =
  Subscript[\, TL]*Subscript/Subscript;
\ = 1 - Subscript[\, LT]*Subscript[\, TL];
Subscript = Subscript/\;
Subscript = Subscript[\, LT]*Subscript/\;
Subscript = Subscript/\;
Subscript = Subscript;
(*脱层几何参数*)
\ = 1600;(*脱层密度*)
c = 0.06;(*阻尼*)
t = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
l = 1;(*为梁总的长度*)
Subscript = 0.3;(*为梁1区的长度*)
Subscript = 0.3;(*为梁2区的长度*)
Subscript = 0.4;(*为梁3区的长度*)
Subscript = 0.4;(*为梁4区的长度*)
nn1 = h/t;(*nn1为总铺层厚度*)
nn2 = 2.0;(*为上子铺层层数*)
nn3 = h/t - nn2;(*下子铺层层数*)
Subscript = t*nn2;(*梁的上部子厚度*)
Subscript = h - Subscript;(*梁的下部子厚度*)
Subscript = h;(*梁1区的厚度*)
Subscript = h;(*梁2区的厚度*)
(*求各单层刚度*)
agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)


ck = Cos;(*列表求各铺层角度余弦*)
sk = Sin;(*列表求各铺层角度正弦*)
s11k = ck^4*Subscript + 2 ck^2*sk^2*Subscript +
             sk^4*Subscript +
   4 ck^2*sk^2*Subscript;(*求各区刚度*)
A11 = \!\(
\*UnderoverscriptBox[\(\\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
A31 = (\!\(
\*UnderoverscriptBox[\(\\), \(k =
     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
A41 = (\!\(
\*UnderoverscriptBox[\(\\), \(k =
     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)


Print["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
D11 = \!\(
\*UnderoverscriptBox[\(\\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)
D21 = D11
D31 = \!\(
\*UnderoverscriptBox[\(\\), \(k =
    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
D41 = \!\(
\*UnderoverscriptBox[\(\\), \(k =
    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
k1 = Subscript/l; k2 = Subscript/l; k3 = Subscript[l, \
3]/l; k4 = Subscript/l; \ = h/l;
a1 = (\*k1)/(DD1*(\1)^3); a2 = (\*k2)/(
DD2*(\2)^3); a3 = (\*k3)/(
DD3*(\3)^3); a4 = (\*k4)/(DD4*(\4)^3);
\1 = h/Subscript; \2 = h/Subscript; \3 \
= h/Subscript; \4 = h/Subscript;
b1 = Subscript/Subscript; b2 = Subscript[l, \
1]/Subscript; b3 = Subscript/Subscript; b4 = \
Subscript/Subscript;
DD1 = D11/(A11*h^2); DD2 = D21/(A11*h^2); DD3 = D31/(A11*h^2);
DD4 = D41/(A11*h^2);
r11 = Power^2, (4)^-1]; r12 = Power[
a2*\^2, (4)^-1]; r13 = Power[
a3*\^2, (4)^-1]; r14 = Power^2, (4)^-1];
BB = {{1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0},
  {0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0},
  {0, Cosh, 0, 0, 0, Sinh, 0, 0, 0, Cos, 0, 0, 0,
   Sin, 0, 0},
  {0, Sinh, 0, 0, 0, Cosh, 0, 0, 0, -Sin, 0, 0, 0,
   Cos, 0, 0},
  {-(r11)^2*DD1*(\1)^2*Cosh,
   0, (r13)^2*DD3*(\3)^2, (r14)^2*DD4*(\4)^2,
   -(r11)^2*DD1*(\1)^2*Sinh, 0, 0,
   0, (r11)^2*DD1*(\1)^2*Cos,
   0, -(r13)^2*DD3*(\3)^2, -(r14)^2*DD4*(\4)^2, (r11)^2*
    DD1*(\1)^2*Sin, 0, 0, 0},
  {-(r11)^3*DD1*(\1)^3*Sinh, 0, 0,
   0, -(r11)^3*DD1*(\1)^3*Cosh,
   0, (r13)^3*DD3*(\3)^3, (r14)^3*DD4*(\4)^3, -(r11)^3*
    DD1*(\1)^3*Sin, 0, 0, 0,
   (r11)^3*DD1*(\1)^3*Cos,
   0, -(r13)^3*DD3*(\3)^3, -(r14)^3*DD4*(\4)^3},
  {0, -(r12)^2*DD2*(\2)^2, (r13)^2*DD3*(\3)^2*
    Cosh, (r14)^2*DD4*(\4)^2*Cosh,
   0, 0, (r13)^2*DD3*(\3)^2*Sinh, (r14)^2*DD4*(\4)^2*
    Sinh, 0, (r12)^2*DD2*(\2)^2,
   -(r13)^2*DD3*(\3)^2*Cos, -(r14)^2*DD4*(\4)^2*
    Cos, 0, 0,
   -(r13)^2*DD3*(\3)^2*Sin, -(r14)^2*DD4*(\4)^2*
    Sin},
  {0, 0, (r13)^3*DD3*(\3)^3*Sinh, (r14)^3*DD4*(\4)^3*
    Sinh, 0, -(r12)^3*DD2*(\2)^3,
   (r13)^3*DD3*(\3)^3*Cosh, (r14)^3*DD4*(\4)^3*
    Cosh, 0, 0, (r13)^3*DD3*(\3)^3*Sin,
   (r14)^3*DD4*(\4)^3*Sin,
   0, (r12)^3*DD2*(\2)^3, -(r13)^3*DD3*(\3)^3*Cos,
   -(r14)^3*DD4*(\4)^3*Cos},
  {Cosh, 0, -1, 0, Sinh, 0, 0, 0, Cos, 0, -1, 0,
   Sin, 0, 0, 0},
  {r11*Sinh, 0, 0, 0, r11*Cosh, 0, -r13*b1,
   0, -r11*Sin, 0, 0, 0, r11*Cos,
   0, -r13*b1},
  {Cosh, 0, 0, -1, Sinh, 0, 0, 0, Cos, 0, 0, -1,
   Sin, 0, 0, 0},
  {r11*Sinh, 0, 0, 0, r11*Cosh, 0,
   0, -r14*b2, -r11*Sin, 0, 0, 0, r11*Cos,
   0, 0, -r14*b2},
  {0, 1, -Cosh, 0, 0, 0, -Sinh, 0, 0, 1, -Cos, 0, 0, 0,
    0, -Sin},
  {0, 0, -r13*b3*Sinh, 0, 0, r12, -r13*b3*Cosh, 0, 0, 0,
   r13*b3*Sin,
   0, 0, r12, -r13*b3*Cos, 0},
  {0, 1, 0, -Cosh, 0, 0, 0, -Sinh, 0, 1, 0, -Cos, 0, 0,
    0, -Sin},
  {0, 0, 0, -r14*b4*Sinh, 0, r12, 0, -r14*b4*Cosh, 0, 0, 0,
   r14*b4*Sin,
   0, r12, 0, -r14*b4*Cos}}
aa1 = Det
Plot, -10, 10}]
ss = FindRoot, 2}][]

百年孤独 发表于 2014-7-8 09:30

这是什么语言编程？

游龙游龙 发表于 2014-7-8 14:44

百年孤独 发表于 2014-7-8 09:30 static/image/common/back.gif
这是什么语言编程？

mathematical编的
页: [1]
查看完整版本: 麻烦大家帮忙看下mathematical编的程序，不知道哪里出错了