【原创】maxwell模型转为prony series的matlab代码

chunyu

下面的mtlab代码可以将图示的maxwell模型转换为Abaqus中viscoelastic材料模型（以prony series 表示）的输入参数。其中体积形变视为弹性，而剪切变形视为粘弹性，适合普通压力下的大多数高分子材料。
%By Zhang chunyu(g0306101@nus.edu.sg,chunyu79@hotmail.com)
%paramters of Maxell model
E0=2.82;%GPa
v0=0.41;
E1=13.05;%GPa
y1=137.32;%GPa.s
E2=8.48;%GPa
y2=893.43;%GPa.s
y0=6500.19;%GPa.s
  
  G0=E0/(2*(1+v0));
  G1=E1/(2*(1+v0));
  G2=E2/(2*(1+v0));
  k0=E0/(3*(1-2*v0));
  p1=(G0*G1*y2+G0*y0*G1+y0*G1*G2+G0*y1*G2+G0*y0*G2)/(G0*G1*G2);
  p2=(G0*y0*y1+y0*G1*y2+G0*y1*y2+G0*y0*y2+y0*y1*G2)/(G0*G1*G2);
  p3=(y0*y1*y2)/(G0*G1*G2);
  q1=2*y0;
  q2=(2*G0*y0*G1*y2+2*G0*y0*y1*G2)/(G0*G1*G2);
  q3=2*y0*y1*y2/(G1*G2);
   
  c=[p3 p2 p1 1];
  r=roots(c);
  t1=-r(1);
  t2=-r(2);
  t3=-r(3);
  h1=1/p3*(q1-q2*t1+q3*t1^2)/((t2-t1)*(t3-t1));
  h2=1/p3*(q1-q2*t2+q3*t2^2)/((t1-t2)*(t3-t2));
  h3=1/p3*(q1-q2*t3+q3*t3^2)/((t1-t3)*(t2-t3));
  % parameters of Prony Series
  GR0=h1+h2+h3
  g1=h1/GR0
  k1=0
  tao1=1/t1
  g2=h2/GR0
  k2=0
  tao2=1/t2
  g3=h3/GR0
  k3=0
  tao3=1/t3

capbull

好贴 学习 。。。。。

chenyuanjp

真太好了 非常感谢

capbull

LZ 配的图片是广义kelvin模型吧？？

ysfeng628

很精彩，但是能请问推导的原理是什么呢

afeibuaa

楼主强大！绝对好贴！

fzkjlt123

好帖，顶

蒲公英87625

怎样将PROE建的三维模型导入ANSOFT， 各位高手提供点帮助啊

nuaaer

啊，还有这样的好东西啊
以后多多交流啊
谢谢楼主

yyl2013121245

谢谢楼主

轩落子

谢谢楼主，感动

dingti

maxwell模型有这么多参数？  这好像是burgers模型吧

danielzlf

有问题请教大神，可否联系，谢谢！Q270296993

danielzlf

请问这个图和这个程序是相匹配的吗？同时图中显示的模型为Kelvin模型

小星星

谢谢楼主，感动

dc137937

楼主，程序和图片都是广义Kelvin模型的。广义Maxwell模型没有n0参数，是不是跟这个一样的原理？

Fwerrs

应该是转载过来的吧

ldisgod

好东西  学习了

ldisgod

这个，有点迷糊

yan_lee

谢谢楼主 学习了