%-- Unknown date --%
SS=find(xzp(:,1)>=A(1) & xzp(:,1)<=A(1))
xzp(SS,:)
SS=find(xzp(:,1)>=A(1) & xzp(:,2)<=A(2))
xzp(SS,:)
SS=find(xzp(:,1)>=A(1) & xzp(:,2)<=A(2))
telem
tnode
tnodepn
elemn
xzp
xzp(1:22)
xzpn(1:22)
elemn(1:22)
[xmin xmax ymin ymax]
A
elemn
find(xzp(:,1)==A(1) & xzp(:,2)==A(2))
elemn
xzp
xzpn
size(xzpn)
xzpn'
xzp'
elemn
xzpn'
xzpn(365,:)
xzpn(365:end,:)
elemn'
tnodep
tnodepn
telemp
telemp*4
j = 1
shape=feisoq48fdtnf(xzpn(j,1),xzpn(j,2))
xzpn(j,1)
xzpn(j,2)
j=2
shape=feisoq48fdtnf(xzpn(j,1),xzpn(j,2))
xzpn(j,1)
xzpn(j,2)
j=150
shape=feisoq48fdtnf(xzpn(j,1),xzpn(j,2))
fal reszleg(0.75)
userfilebak
clear
load embBREZA
clear
load embBREZA
REZA.sigmafi
REZA.sigma20ii*(REZA.KOi*2+1)/3
clear
load embBREZA
REZA.sigmafi = REZA.sigma20ii*(REZA.KOi*2+1)/3;
REZA.sigmafi
clear
load embBREZA
REZA.sigmafi = REZA.sigma20ii*(REZA.KOi*2+1)/3;
save embBREZAn
clear
load embBREZAn
who
REZA.sigmafi
batchembB
%-- 9/6/05 12:30 PM --%
%-- 9/6/05 12:37 PM --%
wex
cd('D:\programfinal\pfinal48\embankmentB\betachange2')
loadvarall tott/86400
load('D:\programfinal\pfinal48\embankmentB\betachange2\embB-coarsemesh-betavalue-4.mat')
embBwks
load('D:\programfinal\pfinal48\embankmentB\valoito\embankmentB-t2000days.mat')
embBwks
uiopen('D:\programfinal\pfinal48\embankmentB\valoito\pwp-bott.fig',1)
uiopen('D:\programfinal\pfinal48\embankmentB\valoito\pwp-bott1cf.fig',1)
load PWPEMbBnn
who
clear
load PWPEMbBnn
who
hold on
plot(atpt26.x,atpt122.y,'r)
plot(atpt26.x,atpt122.y,'r')
plot(atpt26.x,atpt93.y,'b')
clear all
load PWPEMbBnn
a=findobj(gca,'type','line')
Xd=get(a(1),'xdata')';
% % X2=get(a(2),'xdata');
Yd=get(a(1),'ydata')';
atpt122.x=Xd;
atpt122.y=Yd;
clear a Xd yd
who
clear a Xd Yd
who
save PWPEMbBnnf
loadr
embBwks
disarzn = [];
DD = unique(xz(:,2));
for i = 1:length(DD)
locadz(i) = find(xz(:,2) == DD(i),1);
end
depdz1 = xz(locadz,2);
% figure;plot(timefunc/86400,disarz(locadz,:)')
% logx;xlabel('Time (day)');ylabel('Vertical displacement (m)');xlim([1 2000])
% locadzh = find(xz(:,2) == xz(1,2));
% figure;plot(timefunc/86400,disarz(locadzh,:)')
% logx;xlabel('Time (day)');ylabel('Vertical displacement (m)');xlim([1 2000])
Ad = disarz(locadz(1),:);
for i = 1:length(locadz)
Bd = disarz(locadz(i),:);
for j = 1:length(Bd)
if Ad(j) == 0
Cd = 0;
else
Cd = Bd(j)/Ad(j);
end
fac1(i,j) = Cd;
end
end
for i = 1:size(fac1,1)
dzloca(i,:) = fac1(i,:)'.*Yd;
end
figure;plot(timefunc/86400,dzloca)
logx;xlabel('Time (day)');ylabel('Vertical displacement (m)');xlim([1 2000])
hold on;plot(timefunc/86400,disarz(locadz,:)','x')
elemn
elemnx
elemn
disarz(elemn(j),:)
xzele(elemn(j),:)
elemn(j)
elemn
j=132
xzele(elemn(j),:)
elemn(j)
elemn(66)
elemn(67)
elemn(68)
elemn(69)
elemn(369)
elemn(368)
j = 368;
xze = xzele(elemn(j),:)
elemn(j)
conn
disarz(nd,:)
B = disarz(nd,:)
shape
shape*B
shape.*B
shape*B
B*shape
B*shape'
disarzn
plot(timefunc/86400,disarz(1:22,:)','r');hold on
connpn
[0 0.45 0.9 1.45 1.9 2.85 3.8 4.5 5.3 6.1 6.9 8.5 10.5 12.75 15 17.2 19 21 25 29 34 40 46 53 62 72]
xzn(26,:)
xzpn(26,:)
plot(timefunc/86400,disarz(1:26,:)',':b');hold on
figure;plot(timefunc/86400,disarz(1:26,:)',':b');hold on
xzpn(1,:)
xzpn(2,:)
xzpn(3,:)
find(xz(:,2) == xz(1,2))
find(xzpn(:,2) == xzpn(1,2))
figure;plot(timefunc/86400,disarz(1:2,:)',':b');hold on
convert2elemdata
open gauss2elem8fdtnf
nodeselement
i = 1
A = nodeselement(i,connp)
i=100
A = nodeselement(i,connp)
tnodep
sum(porar(A,:))/nelp
for i = 1:tnodep
A = nodeselement(i,connp);
porarn(i,:) = sum(porar(A,:))/nelp;
end
sum(porar(A,:))/nelp
sum(porar(A,:))'/nelp
for i = 1:tnodep
A = nodeselement(i,connp);
porarn(i,:) = sum(porar(A,:))'/nelp;
end
clear porarn
for i = 1:tnodep
A = nodeselement(i,connp);
porarn(i,:) = sum(porar(A,:))'/nelp;
end
clear porarn
%porarf = zeros(telemp,size(porar,2));
for i = 1:tnodep
A = nodeselement(i,connp);
porarn(i,:) = sum(porar(A,:),1)/nelp;
end
clear porarn
%porarf = zeros(telemp,size(porar,2));
for i = 1:tnodep
A = nodeselement(i,connp);
porarn(i,:) = sum(porar(A,:),1)/length(A);%nelp;
end
porarn(i,:)
i
sum(porar(A,:),1)
porar(A,:)
telem
porar(147,:)
minn(porar(147,:))
max(porar(147,:))
max(porar(146,:))
min(porar(146,:))
min(porar(122,:))
min(porar(102,:))
min(porar(1,:))
min(porar(23,:))
min(porar(27,:))
min(porar(127,:))
min(porar(187,:))
min(porar(147,:))
void0telem
gammatelem
per0telem
voidtelemf
kappatelem
Cktelem
visplstr
R1ele0
RRele
stressele0
C31limitele
kappatelem
open backstepnnew8fdtnf
pertelem
Gtelem
xzelen
xzelen(1:telemn,1) = (xzpn(connpn(1:telemn,1),1)+xzpn(connpn(1:telemn,2),1)+xzpn(connpn(1:telemn,3),1)+xzpn(connpn(1:telemn,4),1))/4;
xzelen(1:telemn,2) = (xzpn(connpn(1:telemn,1),2)+xzpn(connpn(1:telemn,2),2)+xzpn(connpn(1:telemn,3),2)+xzpn(connpn(1:telemn,4),2))/4;
clear elemn2
xzelen(1:telemn,1) = (xzpn(connpn(1:telemn,1),1)+xzpn(connpn(1:telemn,2),1)+xzpn(connpn(1:telemn,3),1)+xzpn(connpn(1:telemn,4),1))/4;
xzelen(1:telemn,2) = (xzpn(connpn(1:telemn,1),2)+xzpn(connpn(1:telemn,2),2)+xzpn(connpn(1:telemn,3),2)+xzpn(connpn(1:telemn,4),2))/4;
clear elemn2
for j = 1:telemn
A2 = xzelen(j,:);elemn2(j) = -1;
DD2 = find(xzelen(:,1)==A2(1) & xzelen(:,2)==A2(2));
if isempty(DD2)
for i = 1:telem
C2 = connp(i,:);
x1 = xzp(C2,1);y1 = xzp(C2,2);
xmax = max(x1);xmin = min(x1);ymax = max(y1);ymin = min(y1);
%connpxzl(i,:) = [xmin xmax ymin ymax];
if A2(1) >= xmin & A2(1) <= xmax & A2(2) >= ymin & A2(2) <= ymax
elemn2(j) = i;
break
end
end
else
elemnx2(j) = DD2;
end
end
elemn2
for j = 1:telemn
A2 = xzelen(j,:);elemn2(j) = -1;
DD2 = find(xzpn(:,1)==A2(1) & xzpn(:,2)==A2(2));
if isempty(DD2)
for i = 1:telem
C2 = connp(i,:);
x1 = xzp(C2,1);y1 = xzp(C2,2);
xmax = max(x1);xmin = min(x1);ymax = max(y1);ymin = min(y1);
%connpxzl(i,:) = [xmin xmax ymin ymax];
if A2(1) >= xmin & A2(1) <= xmax & A2(2) >= ymin & A2(2) <= ymax
elemn2(j) = i;
break
end
end
else
elemnx2(j) = DD2;
end
end
clear elemn2 elemnx2
for j = 1:telemn
A2 = xzelen(j,:);elemn2(j) = -1;
DD2 = find(xzpn(:,1)==A2(1) & xzpn(:,2)==A2(2));
if isempty(DD2)
for i = 1:telem
C2 = connp(i,:);
x1 = xzp(C2,1);y1 = xzp(C2,2);
xmax = max(x1);xmin = min(x1);ymax = max(y1);ymin = min(y1);
%connpxzl(i,:) = [xmin xmax ymin ymax];
if A2(1) >= xmin & A2(1) <= xmax & A2(2) >= ymin & A2(2) <= ymax
elemn2(j) = i;
break
end
end
else
elemnx2(j) = DD2;
end
end
for j = 1:telemn
if elemn2(j) == -1
porarn(j,:) = porarN(elemnx2(j),:);
%disarxn(j,:) = disarx(elemnx(j),:);
else
xze = xzele(elemn2(j),:);
nd = connp(elemn2(j),:);
pwp = porarN(nd,:);
zs=-xze(2)+xzelen(j,2);
rs=-xze(1)+xzelen(j,1);
shape=feisoq48fdtnf(rs,zs);
porarn(j,:) = shape*pwp;
%disarxn(j,:) = shape*Bx;
end
end
figure;plot(timefunc/86400,porar(find(xzele(:,2) == xzele(1,2)),:));hold on
plot(timefunc/86400,porarn(find(xzelen(:,2) == xzelen(1,2)),:)','x');hold on
figure;plot(timefunc/86400,porar(find(xzele(:,1) == xzele(1,1)),:));hold on
plot(timefunc/86400,porarn(find(xzelen(:,1) == xzelen(1,1)),:)','x');hold on
for j = 1:telemn
if elemn2(j) == -1
porarn(j,:) = porarN(elemnx2(j),:);
voidtelemtn(j,:) = voidtelemtN(elemnx2(j),:);
visplstrtn(j,:) = visplstrtN(elemnx2(j),:);
gamvisplstrtn(j,:) = gamvisplstrtN(elemnx2(j),:);
esarelexn(j,:) = esarelexN(elemnx2(j),:);
esarelexzn(j,:) = esarelexzN(elemnx2(j),:);
esarelezn(j,:) = esarelezN(elemnx2(j),:);
esareleyn(j,:) = esareleyN(elemnx2(j),:);
strainelexn(j,:) = strainelexN(elemnx2(j),:);
strainelezn(j,:) = strainelezN(elemnx2(j),:);
strainelexzn(j,:) = strainelexzN(elemnx2(j),:);
else
xze = xzele(elemn2(j),:);
nd = connp(elemn2(j),:);
%pwp = porarN(nd,:);
zs=-xze(2)+xzelen(j,2);
rs=-xze(1)+xzelen(j,1);
shape=feisoq48fdtnf(rs,zs);
porarn(j,:) = shape*porarN(nd,:);
voidtelemtn(j,:) = shape*voidtelemtN(nd,:);
visplstrtn(j,:) = shape*visplstrtN(nd,:);
gamvisplstrtn(j,:) = shape*gamvisplstrtN(nd,:);
esarelexn(j,:) = shape*esarelexN(nd,:);
esarelexzn(j,:) = shape*esarelexzN(nd,:);
esarelezn(j,:) = shape*esarelezN(nd,:);
esareleyn(j,:) = shape*esareleyN(nd,:);
strainelexn(j,:) = shape*strainelexN(nd,:);
strainelezn(j,:) = shape*strainelezN(nd,:);
strainelexzn(j,:) = shape*strainelexzN(nd,:);
end
end
sum(Cktelem(A,:),2)/AA
sum(Cktelem(A,:)
sum(Cktelem(A,:))
Cktelem(A,:)
Cktelem
A
Cktelem(A)
sum(Cktelem(A),2)
Cktelem(A)
i=100
A = nodeselement(i,connp)
sum(Cktelem(A),2)/AA
Cktelem(A)
AA = length(A);
sum(Cktelem(A),2)/AA
sum(Gtelem(A),2)/AA
sum(Gtelem(A),1)/AA
for i = 1:tnodep
A = nodeselement(i,connp);
AA = length(A);%nelp
porarN(i,:) = sum(porar(A,:),1)/AA;
voidtelemtN(i,:) = sum(voidtelemt(A,:),1)/AA;
visplstrtN(i,:) = sum(visplstrt(A,:),1)/AA;
gamvisplstrtN(i,:) = sum(gamvisplstrt(A,:),1)/AA;
esarelexN(i,:) = sum(esarelex(A,:),1)/AA;
esarelezN(i,:) = sum(esarelez(A,:),1)/AA;
esarelexzN(i,:) = sum(esarelexz(A,:),1)/AA;
esareleyN(i,:) = sum(esareley(A,:),1)/AA;
strainelexN(i,:) = sum(strainelex(A,:),1)/AA;
strainelexzN(i,:) = sum(strainelexz(A,:),1)/AA;
strainelezN(i,:) = sum(strainelez(A,:),1)/AA;
strainelexN(i,:) = sum(strainelex(A,:),1)/AA;
CktelemN(i,:) = sum(Cktelem(A),2)/AA;
GtelemN(i,:) = sum(Gtelem(A),1)/AA;
pertelemN(i,:) = sum(pertelem(A),2)/AA;
meanstresseleN(i,:) = sum(meanstressele(A),1)/AA;
kappatelemN(i,:) = sum(kappatelem(A),2)/AA;
end
j=100
CktelemN(elemnx2(j),:)
shape*CktelemN(nd,:)
CktelemN
shape*GtelemN(nd,:)
shape*meanstresseleN(nd,:)
shape*kappatelemN(nd,:)
kappatelemN
xzo
open updatexzxzp8fdtnf
open backstepnnew8fdtnf
jilun
xzpfinal
isequal(xz8,xzfinal)
isequal(xz8,xzori)
isequal(xz,xzpori)
runandstruct
load embBREZAn
nargin
extravarifany
constime
wstmp
load('D:\programfinal\pfinal48\embankmentB\valoito\t1s.mat')
basevars
restrictjilun
r
matds = mainvsrm8fdtnf(REZA,'t1s','restrictjilun',1);
matds
matds.basevars
basevars = matds.basevars;
basevars.mpritelem
gamvisplstrt = gamvisplstrtn;kappatelem = kappatelemn;meanstressele = meanstresselen;porar = porarn;
strainelex = strainelexn;strainelexz = strainelexzn;strainelez = strainelezn;visplstrt = visplstrtn;voidtelemt = voidtelemtn;
clear CktelemN GtelemN disarxN disarzN esarelexN esarelexzN esareleyN esarelezN gamvisplstrtN kappatelemN lamdatelemN meanstresseleN porarN strainelexN strainelexzN strainelezN visplstrtN voidtelemtN
tnode8 = tnoden;
conn8 = connn;
tnode = tnodepn;
conn = connpn;
xzpori = xzpn;
xzori = xzn;
xzele = xzelen;
xz = xzpori;
xz8 = xzori;
displace(1:2:tnode*2,1) = disarx(:,end);
displace(2:2:tnode*2,1) = disarz(:,end);
[xzfinal,xzpfinal] = updatexzxzp8fdtnf(xz,xzp,conndis2pore,displace);
clear displace
clear j elemnx2 elemn2 elemnx elemn DD2 A2 C2 DD A C i telemn tnoden tnodepn connn connpn xzelen xmax xmin ymin ymax shape zs rs Bx nd xze Bz pwp x1 y1
clear jj layerbc ndx1 ndx2 multilayer ignore1x xzpn connpn telempn tnodepn conepn ndivlpn ndivwpn connp4n xzn connn telemn tnoden conen ndivln ndivwn connpfn
who
clear Cktelem Gtelem disarx disarz esarelex esarelexz esareley esarelez gamvisplstrt kappatelem lamdatelem meanstressele porar strainelex strainelexz strainelez visplstrt voidtelemt
Cktelem = Cktelemn;Gtelem = Gtelemn;disarx = disarxn;disarz = disarzn;esarelex = esarelexn;esarelexz = esarelexzn;esareley = esareleyn;esarelez = esarelezn;
gamvisplstrt = gamvisplstrtn;kappatelem = kappatelemn;meanstressele = meanstresselen;porar = porarn;
strainelex = strainelexn;strainelexz = strainelexzn;strainelez = strainelezn;visplstrt = visplstrtn;voidtelemt = voidtelemtn;
clear CktelemN GtelemN disarxN disarzN esarelexN esarelexzN esareleyN esarelezN gamvisplstrtN kappatelemN lamdatelemN meanstresseleN porarN strainelexN strainelexzN strainelezN visplstrtN voidtelemtN
tnode8 = tnoden;
conn8 = connn;
tnode = tnodepn;
conn = connpn;
xzpori = xzpn;
xzori = xzn;
xzele = xzelen;
xz = xzpori;
xz8 = xzori;
displace(1:2:tnode*2,1) = disarx(:,end);
displace(2:2:tnode*2,1) = disarz(:,end);
[xzfinal,xzpfinal] = updatexzxzp8fdtnf(xz,xzp,conndis2pore,displace);
matds
basevars
basevars.conndis2pore
[xzfinal,xzpfinal] = updatexzxzp8fdtnf(xzori,xzpori,basevars.conndis2pore,displace);
open mainvsrinitialization8fdtnf
disxp = disarx(:,end);
diszp = disarz(:,end);
% xzfinal(:,1) = xzori(:,1)+disx;
% xzfinal(:,2) = xzori(:,2)+disz;
xzpfinal(:,1) = xzpori(:,1)+disxp;
xzpfinal(:,2) = xzpori(:,2)+diszp;
xzpori(:,1)+disxp
xzpori(:,2)+diszp
xzpfinal(:,2) = xzpori(:,2)+diszp
clear xzpfinal
xzpfinal(:,1) = xzpori(:,1)+disxp;
xzpfinal(:,2) = xzpori(:,2)+diszp;
basevars.conndis2pore
tnode
tnode8
A = nodeselemnt(i,conn8)
A = nodeselement(i,conn8)
i=2
A = nodeselement(i,conn8)
conndis2poren = basevars.conndis2pore;
B = conn8(A,:)
displacen(conndis2poren) = displace
2*conndis2poren-1
gg=2*conndis2poren-1
gg'
displacen(2*conndis2poren-1) = disxp;
displacen(2*conndis2poren) = diszp;
displacen
2*conndis2poren
(displacen(B(1))+displacen(B(3)))/2
displacen(2*B(2)-1) = (displacen(2*B(1)-1)+displacen(2*B(3)-1))/2;
displacen(2*B(2)) = (displacen(2*B(1))+displacen(2*B(3)))/2;
for i = 1:tnode8
A = nodeselement(i,conn8);
B = conn8(A,:);
displacen(2*B(2)-1) = (displacen(2*B(1)-1)+displacen(2*B(3)-1))/2;
displacen(2*B(2)) = (displacen(2*B(1))+displacen(2*B(3)))/2;
displacen(2*B(4)-1) = (displacen(2*B(3)-1)+displacen(2*B(5)-1))/2;
displacen(2*B(4)) = (displacen(2*B(3))+displacen(2*B(5)))/2;
displacen(2*B(6)-1) = (displacen(2*B(5)-1)+displacen(2*B(7)-1))/2;
displacen(2*B(6)) = (displacen(2*B(5))+displacen(2*B(7)))/2;
displacen(2*B(8)-1) = (displacen(2*B(7)-1)+displacen(2*B(1)-1))/2;
displacen(2*B(8)) = (displacen(2*B(7))+displacen(2*B(1)))/2;
end
[xzfinal,xzpfinala] = updatexzxzp8fdtnf(xzori,xzpori,conndis2poren,displacen);
[xzfinal,xzpfinala] = updatexzxzp8fdtnf(xzori,xzpori,conndis2poren,displacen');
isequal(xzpfinal,xzpfinala)
tnode*2
tnode8*2
isnear(xzpfinal,xzpfinala)
xzpfinal(2,:)
xzpfinala(2,:)
clear xzpfinal displacen
xzpfinal(:,1) = xzpori(:,1)+disxp;
xzpfinal(:,2) = xzpori(:,2)+diszp;
conndis2poren = basevars.conndis2pore;
displacen(2*conndis2poren-1) = disxp;
displacen(2*conndis2poren) = diszp;
for i = 1:telem
A = nodeselement(i,conn8);
B = conn8(A,:);
displacen(2*B(2)-1) = (displacen(2*B(1)-1)+displacen(2*B(3)-1))/2;
displacen(2*B(2)) = (displacen(2*B(1))+displacen(2*B(3)))/2;
displacen(2*B(4)-1) = (displacen(2*B(3)-1)+displacen(2*B(5)-1))/2;
displacen(2*B(4)) = (displacen(2*B(3))+displacen(2*B(5)))/2;
displacen(2*B(6)-1) = (displacen(2*B(5)-1)+displacen(2*B(7)-1))/2;
displacen(2*B(6)) = (displacen(2*B(5))+displacen(2*B(7)))/2;
displacen(2*B(8)-1) = (displacen(2*B(7)-1)+displacen(2*B(1)-1))/2;
displacen(2*B(8)) = (displacen(2*B(7))+displacen(2*B(1)))/2;
end
[xzfinal,xzpfinala] = updatexzxzp8fdtnf(xzori,xzpori,conndis2poren,displacen');
isequal(xzpfinal,xzpfinala)
disxp = disarx(:,end);
diszp = disarz(:,end);
% xzfinal(:,1) = xzori(:,1)+disx;
% xzfinal(:,2) = xzori(:,2)+disz;
clear xzpfinal displacen
xzpfinal(:,1) = xzpori(:,1)+disxp;
xzpfinal(:,2) = xzpori(:,2)+diszp;
conndis2poren = basevars.conndis2pore;
displacen(2*conndis2poren-1) = disxp;
displacen(2*conndis2poren) = diszp;
for i = 1:telem
%A = nodeselement(i,conn8);
B = conn8(i,:);
displacen(2*B(2)-1) = (displacen(2*B(1)-1)+displacen(2*B(3)-1))/2;
displacen(2*B(2)) = (displacen(2*B(1))+displacen(2*B(3)))/2;
displacen(2*B(4)-1) = (displacen(2*B(3)-1)+displacen(2*B(5)-1))/2;
displacen(2*B(4)) = (displacen(2*B(3))+displacen(2*B(5)))/2;
displacen(2*B(6)-1) = (displacen(2*B(5)-1)+displacen(2*B(7)-1))/2;
displacen(2*B(6)) = (displacen(2*B(5))+displacen(2*B(7)))/2;
displacen(2*B(8)-1) = (displacen(2*B(7)-1)+displacen(2*B(1)-1))/2;
displacen(2*B(8)) = (displacen(2*B(7))+displacen(2*B(1)))/2;
end
[xzfinal,xzpfinala] = updatexzxzp8fdtnf(xzori,xzpori,conndis2poren,displacen');
isequal(xzpfinal,xzpfinala)
loadr
whos
whosr
whosr('size')
help whosr
whosr -b
whosr('-b')
whosr('b')
whos
open reducenargin
com.mathworks.mlservices.MLEditorServices.newDocument(whos)
loadr
whos
clear
clc
loadr
whos
telem
tnode
size(disarz)
tnode8
telem
telemn
struc2var(basevars)
whos
basevars.C1telem
mainvsrerr
who
load('D:\programfinal\pfinal48\embankmentB\valoito\embankmentB-t2000daysfinemesh.mat')
porar
wks2st
loadr
C1telem
telem
clear
clc
clear
loadr
C1telem
Cktelem
telem
telemp
tnode
tnodep
tnode8
whos
open mainvsrm8fdtnf
if removefemdata == 1
eval(['clear removefemdatalist acri ACi ACr ACRr removefemdata AC1 AC2 AC3 AC4 AC5 AC6 AC7 AC8 ' removefemdatalist]);
end
whos
restrictjilun
savf
(stresselez(:,end)+stresselez(:,end)+stresseley(:,end))/3
esarelez(:,end)
(esarelex(:,end)+esarelez(:,end)+esareley(:,end))/3
meanstresselea = (esarelex(:,end)+esarelez(:,end)+esareley(:,end))/3
[meanstressele meanstresselea]
whos
visplstr
[vranall,vranc] = vranchk4sstnf(telem,visplstrt(:,end),gamvisplstrt(:,end),visdepend,visplsrange);
[vranalla,vranc] = vranchk4sstnf(telem,visplstrt(:,end),gamvisplstrt(:,end),visdepend,visplsrange);
for i = 1:length(vranall)
vranall(4*i-3:4*i) = vranalla(i);
end
for i = 1:length(vranalla)
vranall(4*i-3:4*i) = vranalla(i);
end
[vranalla,vranc] = vranchk4sstnf(telem,visplstrt(:,end),gamvisplstrt(:,end),visdepend,visplsrange);
for i = 1:length(vranalla)
vranall(4*i-3:4*i) = vranalla(i);
end
who
loadr
whos
embBwks
size(porar)
size(timefunc)
deltall
maxtimestep
open anynoreal
find(hasz,pars)
ind(hasz==pars)
find(hasz==pars)
timefuncall
voidtelemt
ha = h{i}
i
i=12
ha = h{i}
hw = strfind(ha,'keep')
keepC31chkall
keepnewrezaall
newrezaall
embBwks
who
betavaluetelem
totalstep = length(timefunc)
clear totalstep
try,totalstep = evalin('base','totalstep');,catch,totalstep = length(timefunc);end
return
clear all
clc
embBwks
tota
clear all
clc
loadr
totalstep
loadr
who
set(gca,'xtick',[1,10,20,50,100,500,1000,2000])
logx
xlim([1 2000])
figforpublish(gcf,0,0,0,0,0,1,1)
reszleg(0.60)
open tmpc
legend boxoff
gtext('At Node 1')
gtext('At Element 131')
logx
xlim([1 2000])
set(gca,'xtick',[1,10,20,50,100,500,1000,2000])
figforpublish(gcf,0,0,0,0,0,1,1)
reszleg(0.60)
legend boxoff
logx
findtexthandles
a=findtexthandles
help findtexthandles
get(a)
edit gtools
gtext(x,y,lineno,'HorizontalAlignment','center','fontname','Proggy Clean','fontsize',8,'backgroundcolor','none','foregroundcolor','none')
gtext('HorizontalAlignment','center','fontname','Proggy Clean','fontsize',8,'backgroundcolor','none','foregroundcolor','none')
logx
xlim([1 2000])
set(gca,'xtick',[1,10,20,50,100,500,1000,2000])
figforpublish(gcf,0,0,0,0,0,1,1)
reszleg(0.60)
legend boxoff
logx
xlim([1 2000])
set(gca,'xtick',[1,10,20,50,100,500,1000,2000])
figforpublish(gcf,0,0,0,0,0,1,1)
reszleg(0.60)
legend boxoff
logx
xlim([1 2000])
set(gca,'xtick',[1,10,20,50,100,500,1000,2000])
figforpublish(gcf,0,0,0,0,0,1,1)
reszleg(0.60)
legend boxoff
logx
xlim([1 2000])
set(gca,'xtick',[1,10,20,50,100,500,1000,2000])
figforpublish(gcf,0,0,0,0,0,1,1)
reszleg(0.60)
legend boxoff
logx
xlim([1 2000])
set(gca,'xtick',[1,10,20,50,100,500,1000,2000])
figforpublish(gcf,0,0,0,0,0,1,1)
reszleg(0.60)
legend boxoff
userfilebak
size(betavalueall,1)
REZA.ndviw
REZA.ndivw
REZA.sigmaii
REZA.sigmafi
plotmesh(xz,conn)
sigmafitelem
sigmaftelem
combineandpcode batchembB
load('D:\pararesults\gfaceffect\Case3fdt-gfac2-gfacfac-0001.mat')
gfacfac
gfac
open delzer8fdtnf
%-- 9/7/05 11:55 AM --%
open batchembB
wex
open pcodemainvsr
%-- 9/7/05 1:10 PM --%
got D:\Temp\shortcutsb
cd D:\Temp\shortcutsb
clc
run('D:\Temp\shortcutsb\shortstool.m')
uigetdir(pwd,'reza')
uigetdir(pwd,'Select the directory where the history.mat will be saved')
wex(prefdir)
prefdir
%-- 9/7/05 1:31 PM --%
%-- 9/7/05 1:33 PM --%
%-- 9/7/05 1:35 PM --%
%-- 9/7/05 1:36 PM --%
