maple可以计算砼的本构模型吗？
本帖最后由 xiaaiyin 于 2012-4-16 16:40 编辑核心混凝土的应力应变本构关系
把这些式子连在一起，最后得到应力应变曲线，请问各位高手，如何实现图片上的功能呢（从word里面粘贴时候没法复制公式，只能贴图片，请各位大侠见谅）
下列计算中忽略单位，需要的话，自行修改数值即可
restart;
unprotect(Zeta, alpha);
f := 405; f := 38.9; A := 2905.24; A := 37094.76; epsilon := 0.2e-2; Zeta := alpha*f/f;
alpha := A/A;
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^2; sigma := f*(1.194+(13/f)^.45*(-0.7485e-1*Zeta^2+.5789*Zeta));
piecewise(`and`(epsilon >= 0, epsilon <= epsilon), sigma*(A*epsilon/epsilon-B*(epsilon/epsilon)^2), epsilon >= epsilon, piecewise(Zeta0 >= 1.12, sigma*(1-q+q*(epsilon/epsilon)^(.1*Zeta0)), `and`(Zeta0 < 1.12, Zeta0 > 0), sigma*epsilon/(epsilon*(beta*(epsilon/epsilon-1)^2+epsilon/epsilon))));
sigma0 := unapply(%, Zeta0, epsilon);
plot(sigma0(Zeta, epsilon), epsilon = 0 .. 0.12e-1, gridlines = true);
这个曲线Maple可以画，
给出具体的fy,fck,alpha,As,Ac,我试试看 本帖最后由 xiaaiyin 于 2012-4-17 12:03 编辑
这个里面Fy=405MPA
Fck=38.9MPA（1MPA单位等于N除以MM的平方
As=2905.24平方毫米
Ac=37094.76平方毫米
ε0=0.002
εc最大取0.012
feiyuzhen 发表于 2012-4-17 13:47 static/image/common/back.gif
下列计算中忽略单位，需要的话，自行修改数值即可
谢谢feiyuzhen版主的热心解答，技术太牛了，简直是当代的欧拉先生~~
我在弱弱的问一句，怎么提取或输出这上面的一些数据，比如说，当应变εc是300E-6,500E--6，1000E-6,1000E-6，2000E-6等值对应输出应力。取得这些值后，输入有限元软件ansys才能够定义本构模型 feiyuzhen 版主你好，piecewise和unapply这两句是怎么用的，如果遇到不懂的语法有什么自带的语法帮助可以查询吗 第7行代码定义了函数sigma0,自变量是Zeta和epsilon,
直接将数值带入函数即可得到任意Zeta和epsilon对应的应力，
piecewise是定义分段函数的，unapply是定义函数的，
Maple自带的帮助文档详细说明了用法，直接拷贝粘贴练习一下就能理解了
页:
[1]