简支梁(直圆杆)变形

简支梁(直圆杆)变形

ID:35196025

大小:232.25 KB

页数:8页

时间:2019-03-21

简支梁(直圆杆)变形_第1页
简支梁(直圆杆)变形_第2页
简支梁(直圆杆)变形_第3页
简支梁(直圆杆)变形_第4页
简支梁(直圆杆)变形_第5页
简支梁(直圆杆)变形_第6页
简支梁(直圆杆)变形_第7页
简支梁(直圆杆)变形_第8页
资源描述:

《简支梁(直圆杆)变形》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

1、简支梁(直圆杆)变形——材料力学计算机实验说明学院:专业:班级:姓名:学号:一、实验课题题目六梁的变形计算输入:1.支撑条件及梁的各区段长度输入(左、右固定端悬臂梁;简支梁;左、右、双外伸梁)2.各载荷大小、作用位置及方向(qi、ai、bi;pj、cj;mk、dk)3.截面参数输入(形状、尺寸及区段)即:圆截面-d,圆心位置;圆环截面-D,d,圆心位置;矩形截面h,b,形心位置;工字型截面-惯性矩IZ,抗弯截面模量Wz4.材料常数:弹性模量E输出:1.某截面挠度及转角2.最大挠度及最大转角3.梁挠曲线的大致形状(图示)二、程序适用条件简支梁为直圆杆,其上作用有集中

2、力、弯矩和均布力,其中均布力均匀分布在整个简支梁上。如下图1所示。b2a2b1a1qyxMF图1简支梁(直圆杆)载荷分布图参数说明:F:作用在简支梁上的集中力大小a1:集中力作用位置距离简支梁左端距离b1:集中力作用位置距离简支梁右端距离q:作用在整个简支梁上的均布力大小M:作用在简支梁上的弯矩大小a2:弯矩作用位置距离简支梁左端距离b2:弯矩作用位置距离简支梁右端距离一、程序代码注:MATLAB软件运行,其中E为材料弹性模量,I为惯性矩clc,clear;disp('请给定材料信息');%输入材料信息l=input('圆杆长度l(/m)=');d=input('

3、圆杆直径d(/m)=');E=input('弹性模量E(/GPa)=');I=double(d^4*3.14/32);disp('请给定各载荷大小及位置');%输入受力情况F=input('集中力大小F(/N)=');a1=input('集中力作用位置(/m)=');b1=l-a1;q=input('均布力大小q(/N/m)=');M=input('弯矩大小M(/N*m)=');a2=input('弯矩作用位置(/m)=');b2=l-a2;x1=0:0.01:a1;%F引入的挠度和转角vx1=(-F*b1*x1*l.^2+F*b1*x1.^3+F*b1.^3*x

4、1)*(1/(6*E*10^9*I*l));ix1=(-F*b1*l.^2+3*F*b1*x1.^2+F*b1.^3)*(1/(6*E*10^9*I*l));x2=a1:0.01:l;x2=x2(2:end);vx2=(-F*l*(x2-a1).^3-F*b1*(l.^2-b1.^2)*x2+F*b1*x2.^3)*(1/(6*E*10^9*I*l));ix2=(-3*F*l*(x2-a1).^2-F*b1*(l.^2-b1.^2)+3*F*b1*x2.^2)*(1/(6*E*10^9*I*l));v11=[vx1,vx2];i11=[ix1,ix2];x11=[

5、x1,x2];x3=0:0.01:l;%q引入的挠度和转角vx3=(-q*x3.*l^3+2*q*l*x3.^3-q*x3.^4)*(1/(24*E*10^9*I));ix3=(-q*l^3+6*q*l*x3.^2-4*q*x3.^3)*(1/(24*E*10^9*I));x4=0:0.01:a2;%M引入的挠度和转角vx4=(M*x4.*l^2-3*M*x4.*b2.^2-M*x4.^3)*(1/(6*E*10^9*I*l));ix4=(M*l^2-3*M*b2.^2-3*M*x4.^2)*(1/(6*E*10^9*I*l));x5=a2:0.01:l;x5=x

6、5(2:end);vx5=(-M*x5.^3+3*M*l*(x5-a2).^2+M*(l^2-3*b2^2)*x5)*(1/(6*E*10^9*I*l));ix5=(-3*M*x5.^2+6*M*l*(x5-a2)+M*(l^2-3*b2^2))*(1/(6*E*10^9*I*l));v22=[vx4,vx5];i22=[ix4,ix5];x22=[x4,x5];v33=v11+vx3+v22;%叠加i33=i11+ix3+i22;vm=max(abs(v33));%最大挠度temp1=find(abs(v33)==vm);str1=strcat(num2str(

7、v33(temp1)),'m');im=max(abs(i33));%最大转角temp2=find(abs(i33)==im);fprintf(sprintf('最大挠度(/m)=%f',v33(temp1)))fprintf(sprintf('最大转角(/rad)=%f',i33(temp2)))figure;plot(x11,v33),xlabel('x/m'),ylabel('v(x)/m')text(x3(temp1),v33(temp1),str1);holdon;plot(x3(temp1),v33(temp1),'.r');title('简支

8、梁(直圆杆

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。