| 查看: 1505 | 回复: 7 | ||
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); |
回复此楼» 猜你喜欢
291求调剂
已经有7人回复
-
300求调剂
已经有8人回复
-
求调剂推荐
已经有7人回复
-
289 分105500药学专硕求调剂(找B区学校)
已经有4人回复
-
0854求调剂
已经有13人回复
-
初试324 中药学 一志愿天中医 求调剂
已经有3人回复
-
药学求调剂
已经有14人回复
-
327求调剂
已经有27人回复
-
急需调剂
已经有5人回复
-
271求调剂
已经有33人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助MATLAB做三维图问题
已经有4人回复
- 求助一个matlab关于分类讨论的问题
已经有9人回复
- matlab画图求助
已经有11人回复
- 求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟
已经有10人回复
- 求助matlab软件压缩包
已经有4人回复
- 求助matlab编程
已经有13人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】matlab中如何将程序的结果输出到一个文件中
已经有11人回复
- 【求助】请教把java程序转化成matlab程序
已经有9人回复
- 【求助】初学matlab
已经有11人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
- matlab程序求助!
已经有3人回复
 |
2楼2013-04-11 22:52:33
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
3楼2013-04-12 16:26:25
seaspringsun
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 412.3
- 帖子: 252
- 在线: 1070.1小时
- 虫号: 2222038
- 注册: 2013-01-04
- 专业: 工程地质学
4楼2013-04-12 20:03:03
李加伟
荣誉版主 (文坛精英)
不过凡人
- 应助: 258 (大学生)
- 贵宾: 1.257
- 金币: 24684.9
- 散金: 8829
- 红花: 146
- 沙发: 13
- 帖子: 12823
- 在线: 3507.4小时
- 虫号: 2359263
- 注册: 2013-03-19
- 性别: GG
- 专业: 可再生与替代能源利用中的
- 管辖: 论文道贺祈福
5楼2013-04-13 12:40:15
6楼2013-04-13 20:42:32
7楼2013-04-13 21:01:44
8楼2016-06-23 14:37:59
5