代码

import matplotlib.pyplot as plt
import mpl_toolkits.axisartist as AA
from mpl_toolkits.axes_grid1 import host_subplot%matplotlib inline
config = {"font.family": "serif","font.size": 14,"mathtext.fontset": "stix","font.serif": ["Times New Roman"],"xtick.direction": "in","ytick.direction": "in",
}
plt.rcParams.update(config)host = host_subplot(111, axes_class=AA.Axes)
# plt.subplots_adjust(right=0.75)
par1 = host.twinx()
par2 = host.twinx()
par3 = host.twinx()# set the position of axes
host.axis["top"].set_visible(True)offset = 0
new_fixed_axis = par1.get_grid_helper().new_fixed_axis
par1.axis["right"] = new_fixed_axis(loc="right", axes=par1, offset=(offset, 0))
par1.axis["right"].toggle(all=True)offset = 50
new_fixed_axis = par2.get_grid_helper().new_fixed_axis
par2.axis["right"] = new_fixed_axis(loc="right", axes=par2, offset=(offset, 0))
par2.axis["right"].toggle(all=True)
# par2.axis["right"].set_visible(False)offset = 110
new_fixed_axis = par3.get_grid_helper().new_fixed_axis
par3.axis["right"] = new_fixed_axis(loc="right", axes=par3, offset=(offset, 0))
par3.axis["right"].toggle(all=True)host.set_xlabel("Time (s)", fontsize=config["font.size"] + 2)
# host.set_ylabel("Shear stress, normal stress (MPa)")
host.text(-22, 0.35, 'Normal stress, ', color='c', va='center', rotation='vertical', fontsize=config["font.size"] + 2)
host.text(-22, 0.74, 'shear stress ', color='m', va='center', rotation='vertical', fontsize=config["font.size"] + 2)
host.text(-22, 1.01, '(MPa)', color='k', va='center', rotation='vertical', fontsize=config["font.size"] + 2)
par1.set_ylabel("Slip displacement (mm)", fontsize=config["font.size"] + 2)
par2.set_ylabel("Slip velocity (mm/s)", fontsize=config["font.size"] + 2)
par3.set_ylabel("Normal displacement (mm)", fontsize=config["font.size"] + 2)# plot
colors = ["c", "m", "b", "r", "k"]
p0 = host.plot(data2.index,data2["垂直応力(1)[MPa]"],"-",linewidth=1,color=colors[0],label="Normal stress",
)
p0 = host.plot(data2.index,data2["せん断応力[MPa]"],"-",linewidth=1,color=colors[1],label="Shear stress",
)
p1 = par1.plot(data2.index,data2["せん断変位[mm]"],linewidth=1,color=colors[2],label="shear displacement",
)
p2 = par2.plot(data2.index,data2["velocity[mm/s]"],linewidth=1,color=colors[3],label="shear velocity",
)
# par2.plot([0, 800], [0, 0], "k-", linewidth=0.8)
p3 = par3.plot(data2.index,data2["垂直全平均変位[mm]"],linewidth=1,color=colors[4],label="Normal displacement",
)# 轴名称颜色
host.axis["left"].label.set_color(p0[0].get_color())
par1.axis["right"].label.set_color(p1[0].get_color())
par2.axis["right"].label.set_color(p2[0].get_color())
par3.axis["right"].label.set_color(p3[0].get_color())# set limits of the axes
host.set_xlim(0, 200)
host.set_ylim(0, 1.2)
par1.set_ylim(0, 18)
par2.set_ylim(-0.05, 0.4)
par3.set_ylim(0, 0.6)par3.plot([45, 45], [0, 20], "k--", alpha=0.5)
# start and end time point for evaluating normal stiffness
par3.plot([120, 120], [0, 20], "k--", alpha=0.5, label='Start for computing normal stifness')
par3.plot([148, 148], [0, 20], "k--", alpha=0.5, label='End for computing normal stifness')
# start and end time point for shearing dilation
par3.plot([178, 178], [0, 20], "k--", alpha=0.5, label='Start for computing shear dilation')
par3.plot([187, 187], [0, 20], "k--", alpha=0.5, label='End for computing shear dilation')# plt.legend(ncol=1, bbox_to_anchor=(1.45, 0.7), loc=2, borderaxespad=0)
host.text(3, 1.12,'${\\rm JRC=3.21}, \\sigma_{\\rm n}^{\\rm 0}=1 \\ {\\rm MPa}$',fontsize=config["font.size"] + 2)
host.text(-25, 1.185, '(a)', fontsize=config["font.size"] + 2, fontweight='bold')plt.show()

成品图