-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathHydraulic_hysteresis.m
104 lines (102 loc) · 2.15 KB
/
Hydraulic_hysteresis.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
close all;
clear all;
load('hyd_hys_1.txt');
z1=length(hyd_hys_1);
s1=(hyd_hys_1(1:z1,2));
smax1=max(s1);
sp1=s1/smax1*100;
T11=(hyd_hys_1(1:z1,3));
T12=(hyd_hys_1(1:z1,4));
t1=(.1:.1:z1/10)';
load('hyd_hys_2.txt');
z2=length(hyd_hys_2);
s2=(hyd_hys_2(1:z2,2));
smax2=max(s2);
sp2=s2/smax2*100;
T21=(hyd_hys_2(1:z2,3));
T22=(hyd_hys_2(1:z2,4));
t2=(.1:.1:z2/10)';
load('hyd_hys_3.txt');
z3=length(hyd_hys_3);
s3=(hyd_hys_3(1:z3,2));
smax3=max(s3);
sp3=s3/smax3*100;
T31=(hyd_hys_3(1:z3,3));
T32=(hyd_hys_3(1:z3,4));
t3=(.1:.1:z3/10)';
load('hyd_hys_4.txt');
z4=length(hyd_hys_4);
s4=(hyd_hys_4(1:z4,2));
smax4=max(s4);
sp4=s4/smax4*100;
T41=(hyd_hys_4(1:z4,3));
T42=(hyd_hys_4(1:z4,4));
t4=(.1:.1:z4/10)';
load('hyd_hys_5.txt');
z5=length(hyd_hys_5);
s5=(hyd_hys_5(1:z5,2));
smax5=max(s5);
sp5=s5/smax5*100;
T51=(hyd_hys_5(1:z5,3));
T52=(hyd_hys_5(1:z5,4));
t5=(.1:.1:z5/10)';
load('hyd_hys_6.txt');
z6=length(hyd_hys_6);
s6=(hyd_hys_6(1:z6,2));
smax6=max(s6);
sp6=s6/smax6*100;
T61=(hyd_hys_6(1:z6,3));
T62=(hyd_hys_6(1:z6,4));
t6=(.1:.1:z6/10)';
scrsz=get(0,'Screensize');
figure('Position',[1 scrsz(4) scrsz(3) scrsz(4)])
subplot(3,2,1);
plot(t1,sp1,t1,T11,t1,T12);
xlim([0 17500]);
ylim([-1 101]);
xlabel('Time / s');
ylabel('Data');
title('Hydraulic heating - hysteresis 1');
grid on;
subplot(3,2,2);
plot(t2,sp2,t2,T21,t2,T22);
xlim([0 15500]);
ylim([-1 101]);
xlabel('Time / s');
ylabel('Data');
title('Hydraulic heating - hysteresis 2');
grid on;
subplot(3,2,3);
plot(t3,sp3,t3,T31,t3,T32);
%xlim([0 52]);
ylim([-1 101]);
xlabel('Time / s');
ylabel('Data');
title('Hydraulic heating - hysteresis 3');
grid on;
subplot(3,2,4);
plot(t4,sp4,t4,T41,t4,T42);
xlim([0 18500]);
ylim([-1 101]);
xlabel('Time / s');
ylabel('Data');
title('Hydraulic heating - hysteresis 4');
grid on;
subplot(3,2,5);
plot(t5,sp5,t5,T51,t5,T52);
xlim([0 18500]);
ylim([-1 101]);
xlabel('Time / s');
ylabel('Data');
title('Hydraulic heating - hysteresis 5');
grid on;
subplot(3,2,6);
plot(t6,sp6,t6,T61,t6,T62);
xlim([0 17000]);
ylim([-1 101]);
xlabel('Time / s');
ylabel('Data');
title('Hydraulic heating - hysteresis 6');
grid on;
legend('Strain','Temperature 1','Temperature 2','Location','East');
legend('boxoff');