应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab程序求助
81/1返回列表
查看: 1634  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

seaspringsun

银虫 (小有名气)

[求助] matlab程序求助 已有1人参与

下面这段程序，没有看懂，谁能帮忙看一下啊

%        APPENDIX   A:
%        THERMO-MECHANICAL   LOAD-TRANSFER   ANALYSIS   FOR   ENERGY   FOUNDATIONS
%
%   This   program   performs   T-z   analyses   for   an   energy   foundation
%   under   "Thermo-mechanical   loading"   based   on   the   Beta   method.
clear   all
format   long   e
%%----------------------------------------------------------------------
-%
%%---------------------------------INPUT   PARAMETERS
%%----------------------------------------------------------------------
-%
%%   SOIL   DESCRIPTION   %%
cop   =   0.3;   %c/p   ratio
gamsoil   =   18;%KN/M3
Phi   =   30;   %Drained   friction   angle
%%   FOUNDATION   DESCRIPTION   %%
L   =   10;%m
D   =   1;   %m
E   =   20000;%MPa
gampile   =        24;%KN/m3
alphat   =   10e-6;   %Coeffcient   of   thermal   expansion   for   R
%%   MECHANICAL   LOADING%%
P   =   500;%KN
%%   TEMPERATURE   CHANGE   %%
deltat   =   20;
%-----------------------------------------------------------------------
-%
%%   BETA   METHOD   PARAMETERS   %%
beta   =0.55;   %Emperical   reduction   factor
k   =   65;   %emperical   coefficient   representing   the   soil   resistance   to
expansion   of   the   foundation
%%   END   BEARING   FITTING   PARAMETERS
ab   =   0.02;
bb   =   0.9;
%%   SHEAR   RESISTANCE   FITTING   PARAMETERS
as   =   0.0035;
bs   =   0.9;
%%   HEAD-STRUCTURE   STIFFNESS
Kh=300000;%KN/m   (0.5   GPa/m   for   this   case)
%%   NUMBER   OF   ELEMENTS   %%
N   =   30;
%-----------------------------------------------------------------------
-%
%%
Ab   =   pi*D^2/4;   %   m2   -   Cross-sectional   area
Cs   =   pi*D;   %m   -   Circumference   of   pile
Wpile   =   gampile*Ab*L;   %KN   -   Weight   of   pile
Wdisp   =   gamsoil*Ab*L;   %kN   -   Weight   of   soil   displaced
Wp   =   Wpile   -   Wdisp;   %kN   -   Buoyant   unit   weight   of   pile
%-----------------------------------------------------------------------
-%
%-----------------------------------------------------------------------
-%
%%   MECHANICAL   TZ   ANALYIS   BASED   ON   "P"   %%

%   ULTIMATE   SIDE   SHEAR   RESISTANCE   %
K0=1-sin(Phi*pi/180);   %Coefficient   of   lateral   earth   pressure   at   rest
dL   =L/N;   %m   -   Length   of   each   soil   sublayer
h   =   0:dL:L;   %m   -   Depth   vector,   z=0   at   ground   surface
zmid   =   h(1:length(h)-1)+dL/2;   %m   -   Depth   of   each   soil   sublayer
sigv   =   gamsoil*zmid;   %m   -   Effecive   vertical   stress
fsM   =   beta*sigv*K0*tan(Phi*pi/180);   %KN/m2
As   =   Cs*dL*ones(length(zmid),1);%m2

69






QsM   =   fsM.*As';   %kN
QstotM   =   sum(QsM);   %kN

%   ULTIMATE   END   BEARING%
sc   =   1.2;
dc   =   1.5;
Nc   =   5;
cub   =   L*gamsoil*cop;   %KN/m2
Qp   =   cub*Ab*sc*dc*Nc;   %kN

%   ULTIMATE   CAPACITY
QuM   =   Qp   +   QstotM   -   Wp;   %kN

n   =   1:1:N;
Li   =   L/N;   %m
z   =   Li*(n-.5);   %m
Ki   =   Ab*E*1000/Li;   %KN/m
QbM   =   zeros(N,1);
QtM   =   zeros(N,1);
FsM   =   zeros(N,1);
FsmaxM   =   zeros(N,1);
pbM   =   zeros(N,1);
ptM   =   zeros(N,1);
psM   =   zeros(N,1);
QaveM   =   zeros(N,1);
deltaM   =   zeros(N,1);
sM   =   zeros(N,1);
sigmaM   =   zeros(N,1);
pb   =   0;
Dtot   =   100;
tol   =   1e-6;
PDif   =   500;   %   Previous   Difference
CDif   =   500;   %   Current   Difference
FsmaxM(1)=QsM(1);
for   k=2:1:N
FsmaxM(k)=FsmaxM(k-1)+QsM(k);
end
wxyz   =   0;
while   abs(CDif)<=   abs(PDif)
wxyz   =   wxyz   +   1;
for   i   =   N:-1:1
if   i   ==   N
deltaQt   =   10;   %Initialize   while   loop
QtM(i)   =   0;
pbM(i)   =   pb;
QbM(i)=Qp*pbM(i)/(ab+pbM(i)*bb);
%QbM(i)=0;
while   deltaQt   >   tol
QaveM(i)   =   (QtM(i)+QbM(i))/2;
deltaM(i)   =   QaveM(i)/Ki;
ptM(i)   =   pbM(i)   +   deltaM(i);
psM(i)   =   pbM(i)   +   deltaM(i)/2;   %   pt   +(pt-pb)/2;
FsM(i)   =   FsmaxM(i)*psM(i)/(as+psM(i)*bs);
QtnewM(i)   =   FsM(i)   +   QbM(i);
deltaQt   =   QtnewM(i)   -   QtM(i);
QtM(i)   =   QtnewM(i);
sM(i)   =   deltaM(i)*1e6/Li;%   micro   starin
sigmaM(i)   =   QaveM(i)/Ab;
end

else
deltaQt   =   10;   %Initialize   while   loop
pbM(i)   =   ptM(i+1);
QbM(i)   =   QtM(i+1);
70






QtM(i)   =   0;
while   deltaQt   >   tol
QaveM(i)   =   (QtM(i)+QbM(i))/2;
deltaM(i)   =   QaveM(i)/Ki;
ptM(i)   =   pbM(i)   +   deltaM(i);
psM(i)   =   pbM(i)   +   deltaM(i)/2;
FsM(i)=FsmaxM(i)*psM(i)/(as+psM(i)*bs);
QtnewM(i)   =   FsM(i)   +   QbM(i);
deltaQt   =   QtnewM(i)   -   QtM(i);
QtM(i)   =   QtnewM(i);
sM(i)   =   deltaM(i)*1e6/Li;
sigmaM(i)   =   QaveM(i)/Ab;
end

end

end
i=1:1:N;
FsMtot   =   sum(FsM(i));
PDif   =   CDif;
CDif   =   P-FsMtot-QbM(N);
disp('CDif');
disp(CDif);
disp('PDif');
disp(PDif);
pb=pb+1e-6;
end

pb=pb-1e-6;
disp('pb');
disp(pb);
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-04-11 21:05:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

hzj12399
2楼2013-04-11 22:52:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2013-04-13 13:34:11
ben_ladeng: 专家考核 2013-04-13 20:13:36
这么长段程序,如果不是你的专业的,基本不明白是在干啥
赞一下 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
3楼2013-04-12 16:26:25
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

seaspringsun

银虫 (小有名气)
引用回帖:
3楼: Originally posted by libralibra at 2013-04-12 16:26:25
这么长段程序,如果不是你的专业的,基本不明白是在干啥

说实话，这只是一小部分
赞一下 回复此楼
4楼2013-04-12 20:03:03
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

李加伟

荣誉版主 (文坛精英)

不过凡人

优秀版主优秀版主

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2013-04-13 13:34:21
总体的效果是什么？如果单段可以运行的话，将其他部分注释来看看效果会对你有帮助的。
这么长的程序，你根据他的注释能理解主要意思是在做能量形成的T-Z分析，每一段的主要作用也可以根据注释来看的。
想要知道具体每一步骤的作用就要对每个函数的作用进行分析，这个你可以查一下Matlab中的help函数看看
赞一下 回复此楼
人生海海，在于内心
5楼2013-04-13 12:40:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

zxhao168

捐助贵宾 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
理解程序切忌断章取义，函数的调用、参数的传递等的需要上下文的对应。
赞一下 回复此楼
6楼2013-04-13 20:42:32
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

quanxing

金虫 (小有名气)
第一：你给的程序段不全，需要详细的
第二：你的专业是什么？需要有点了解。
赞一下 回复此楼
好彩自然来
7楼2013-04-13 21:01:44
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xuanqi521

新虫 (著名写手)

【答案】应助回帖

师兄，你这个问题解决没有啊
赞一下 回复此楼
8楼2016-06-23 14:37:59
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 seaspringsun 的主题更新
81/1返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] Sci. Bull. 悲剧经验 +4 jyang1999 2026-05-16 5/250 2026-05-20 13:49 by jyang1999
[基金申请] 国自然评分 +3 无名者登山 2026-05-20 4/200 2026-05-20 13:29 by vito刘
[基金申请] 提交了我也来说说感想 +4 fummck 2026-05-20 4/200 2026-05-20 12:56 by Equinoxhua
[基金申请] 面上本子正文33页，违规吗？会被低分嘛？ +9 1234567wang 2026-05-17 11/550 2026-05-20 12:37 by zhuzg0628
[考博] 如果工作了想读博，可以边工作边读全日制嘛？ 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[有机交流] 如何实现卤原子转化 +3 BT20230424 2026-05-15 8/400 2026-05-20 09:32 by xtlilibin
[基金申请] 河北省自然科学基金 +3 Peterchao 2026-05-18 3/150 2026-05-20 08:57 by 霸_霸
[教师之家] 上海大学实验技术岗位非升即走 +12 嘻嘻哈哈乐呵呵 2026-05-15 13/650 2026-05-20 08:34 by xli1984
[文学芳草园] 献血感触 +6 呀呀好傻 2026-05-19 11/550 2026-05-19 22:26 by 而立得乐
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +3 光伏爱好者 2026-05-17 5/250 2026-05-19 19:13 by 光伏爱好者
[基金申请] 2026山东省优青 +3 luckinging 2026-05-14 4/200 2026-05-19 16:15 by 52wls
[考博] 找博士生导师 +6 小代想上岸 2026-05-15 7/350 2026-05-19 10:22 by free_fisher
[考博] 26/27申博自荐-锂/钠电池方向 5+3 狗头军师. 2026-05-15 4/200 2026-05-19 09:10 by moonboat
[基金申请] 同样的基金本子，换个专家直接从C变A！ (金币-10) +3 国自然国社科中 2026-05-19 3/150 2026-05-19 08:50 by Equinoxhua
[考博] 2026博士还有哪些学校有名额 +7 小王求读研 2026-05-15 8/400 2026-05-19 08:27 by zhyzzh
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[有机交流] 求助2，4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
信息提示
关闭
请填处理意见
关闭