maple可以计算砼的本构模型吗？

xiaaiyin

核心混凝土的应力应变本构关系

把这些式子连在一起，最后得到应力应变曲线，请问各位高手，如何实现图片上的功能呢（从word里面粘贴时候没法复制公式，只能贴图片，请各位大侠见谅）

feiyuzhen

下列计算中忽略单位，需要的话，自行修改数值即可

1. 
   
2. restart;
   
3. unprotect(Zeta, alpha);
   
4. f[y] := 405; f[ck] := 38.9; A[s] := 2905.24; A[c] := 37094.76; epsilon[0] := 0.2e-2; Zeta := alpha*f[y]/f[ck];
   
5. alpha := A[s]/A[c];
   
6. k := .1*Zeta^.745; A := 2-k; B := 1-k; q := k/(.2+.1*Zeta); beta := 0.5e-3*0.236e-4^(.25+(Zeta-.5)^7)*f[ck]^2; sigma[0] := f[ck]*(1.194+(13/f[ck])^.45*(-0.7485e-1*Zeta^2+.5789*Zeta));
   
7. piecewise(`and`(epsilon[c] >= 0, epsilon[c] <= epsilon[0]), sigma[0]*(A*epsilon[c]/epsilon[0]-B*(epsilon[c]/epsilon[0])^2), epsilon[c] >= epsilon[0], piecewise(Zeta0 >= 1.12, sigma[0]*(1-q+q*(epsilon[c]/epsilon[0])^(.1*Zeta0)), `and`(Zeta0 < 1.12, Zeta0 > 0), sigma[0]*epsilon[c]/(epsilon[0]*(beta*(epsilon[c]/epsilon[0]-1)^2+epsilon[c]/epsilon[0]))));
   
8. sigma0 := unapply(%, Zeta0, epsilon[c]);
   
9. plot(sigma0(Zeta, epsilon[c]), epsilon[c] = 0 .. 0.12e-1, gridlines = true);

复制代码

feiyuzhen

这个曲线Maple可以画，
给出具体的fy,fck,alpha,As,Ac,我试试看

xiaaiyin

这个里面Fy=405MPA
Fck=38.9MPA（1MPA单位等于N除以MM的平方
As=2905.24平方毫米
Ac=37094.76平方毫米
ε0=0.002
εc最大取0.012

xiaaiyin

feiyuzhen 发表于 2012-4-17 13:47
下列计算中忽略单位，需要的话，自行修改数值即可

谢谢feiyuzhen版主的热心解答，技术太牛了，简直是当代的欧拉先生~~
我在弱弱的问一句，怎么提取或输出这上面的一些数据，比如说，当应变εc是300E-6,500E--6，1000E-6,1000E-6，2000E-6等值对应输出应力。取得这些值后，输入有限元软件ansys才能够定义本构模型

xiaaiyin

feiyuzhen 版主你好，piecewise和unapply这两句是怎么用的，如果遇到不懂的语法有什么自带的语法帮助可以查询吗

feiyuzhen

第7行代码定义了函数sigma0,自变量是Zeta和epsilon[c],
直接将数值带入函数即可得到任意Zeta和epsilon[c]对应的应力，
piecewise是定义分段函数的，unapply是定义函数的，
Maple自带的帮助文档详细说明了用法，直接拷贝粘贴练习一下就能理解了