首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >Python画弯矩图剪力图(一)

Python画弯矩图剪力图(一)

作者头像
fem178
发布2021-03-09 15:12:39
2K0
发布2021-03-09 15:12:39
举报

Python画弯矩图,都以分段直线近似代替曲线,分段越多,越接近真实的弯矩曲线,就像这样

import matplotlib.pyplot as plt
import numpy as np
fig, axs = plt.subplots(1, 1, constrained_layout=True)
axs.plot([0,1,5],[0,2,1], linewidth=2, color="dimgrey")
plt.show()

除了大致轮廓,还可以把细节完善一下。比如两端封闭起来,填充颜色等等。以受均布荷载的简支梁为例

import matplotlib.pyplot as plt
import numpy as np
plt.rcParams['font.sans-serif'] = ['SimHei']  # 正常显示中文标签
plt.rcParams['axes.unicode_minus'] = False  # 正常显示负号

l = 5
q = 20
x = np.linspace(0,l,20) #将梁20等分
M = q/2*(l*x-x**2)     #计算等分点弯矩M
V = q*(l/2-x)          #计算等分点剪力V

fig, axs = plt.subplots(2, 1, constrained_layout=True,figsize=(8,4) )
axs[0].invert_yaxis() 
axs[0].plot(x,M, label="M1", linewidth=3) #画弯矩图
axs[0].fill_between(x, 0, M, alpha=0.3, where=(M<=0) ) 
axs[0].fill_between(x, 0, M, alpha=0.3, where=(M>=0) )

axs[0].plot([0,5],[0,0], linewidth = 2, color = "dimgrey") #画梁
axs[0].set_ylabel('弯矩(kN.m)', fontsize = 16)
fig.suptitle('内力图', fontsize = 18)
x1 = np.array( [0] )
x2 = np.array( [l] )
x3 = np.hstack( (x1, x, x2) )
v1 = np.array( [ 0 ] )
v3 = np.hstack( (v1, V, v1) )
axs[1].plot(x3, v3, label="M1", linewidth=3)#画剪力图
axs[1].set_ylabel('剪力(kN)', fontsize = 16)
axs[1].fill_between(x, 0, V, alpha=0.3, where=(V<=0) ) #填充颜色
axs[1].fill_between(x, 0, V, alpha=0.3, where=(V>=0) )
axs[1].plot([0,5],[0,0], linewidth = 2, color = "dimgrey") #画梁

fig.savefig('./f2.png', dpi = 500) #保存图片 
plt.show()

再看一个多跨梁

弯矩图注意事项:

  • 正弯矩画在杆件的下方,负弯矩画在杆件的上方。
  • 使杆件下部受拉的弯矩为正,上部受拉的弯矩为负。
  • 弯矩图画在杆件纤维受拉的一侧。

剪力图注意事项:

  • 正剪力画在杆件的上方;
  • 负剪力画在杆件的下方;
  • 使杆件截面顺时针方向转动的剪力为正剪力;
  • 使杆件截面逆时针方向转动的剪力为负剪力;
  • 一般情况下,剪力与杆件所受外力的方向相反。

弯矩图是一条表示杆件不同截面弯矩的曲线。这里所说的曲线是广义的,它包括直线、折线和一般意义的曲线。弯矩图是对构件弯矩的图形表示,弯矩图画在受拉侧,无须标正负号。

本文参与 腾讯云自媒体分享计划,分享自微信公众号。
原始发表:2021-02-28,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 数值分析与有限元编程 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体分享计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档