材料力学剪力图弯矩图绘制有详细的程序Word文件下载.docx
- 文档编号:22402969
- 上传时间:2023-02-03
- 格式:DOCX
- 页数:24
- 大小:18.06KB
材料力学剪力图弯矩图绘制有详细的程序Word文件下载.docx
《材料力学剪力图弯矩图绘制有详细的程序Word文件下载.docx》由会员分享,可在线阅读,更多相关《材料力学剪力图弯矩图绘制有详细的程序Word文件下载.docx(24页珍藏版)》请在冰豆网上搜索。
2.4左端外伸梁QMDWZ
2.5右端外伸梁QMDWY
2.6两端外伸梁QMDWL
functionMM=QMM(n,x1,a,M,MM)
forj=1:
n
ifx1(j)==a
n1=j;
end
end
MM(n1:
n)=MM(n1:
n)+M;
function[QQ,MM]=QMP(n,x1,b,P,QQ,MM)
ifx1(j)==b;
QQ(n1:
n)=QQ(n1:
n)-P;
n)-P*(x1(n1:
n)-b);
function[QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM)
ifx1(j)>
c
QQ(j)=QQ(j)-q*(x1(j)-c);
MM(j)=MM(j)-0.5*q*(x1(j)-c)^2;
d
QQ(j)=QQ(j)+q*(x1(j)-d);
MM(j)=MM(j)+0.5*q*(x1(j)-d)^2;
function[QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM)
XQM=[x1'
QQ'
MM'
];
[Qmax,i]=max(QQ);
Q1=[Qmax,x1(i)];
[Qmin,i]=min(QQ);
Q2=[Qmin,x1(i)];
[Mmax,i]=max(MM);
M1=[Mmax,x1(i)];
[Mmin,i]=min(MM);
M2=[Mmin,x1(i)];
disp('
剪力极值及位置'
)
QDX=[Q1;
Q2]
弯矩极值及位置'
MDX=[M1;
M2]
t1=findobj(0,'
Tag'
'
text31'
);
str=num2str(Q1
(1));
set(t1,'
String'
str);
t2=findobj(0,'
text39'
str=num2str(Q1
(2));
set(t2,'
t3=findobj(0,'
text32'
str=num2str(Q2
(1));
set(t3,'
t4=findobj(0,'
text40'
str=num2str(Q2
(2));
set(t4,'
m1=findobj(0,'
text33'
str=num2str(M1
(1));
set(m1,'
m2=findobj(0,'
text41'
str=num2str(M1
(2));
set(m2,'
m3=findobj(0,'
text34'
str=num2str(M2
(1));
set(m3,'
m4=findobj(0,'
text42'
str=num2str(M2
(2));
set(m4,'
functionTU_QM(x1,QQ,MM)
h1=findobj(0,'
axes1'
axes(h1);
plot(x1,QQ);
grid;
title('
剪力图'
h2=findobj(0,'
axes2'
axes(h2);
plot(x1,MM);
弯矩图'
functionXQM=QMDJ(x,MPQ)
[n,m]=size(x);
L=x(m);
x1=[];
fori=1:
m-1
x1=[x1,linspace(x(i),x(i+1),50)];
MM=zeros(size(x1));
QQ=zeros(size(x1));
[m,t]=size(MPQ);
[t,n]=size(x1);
m
switchMPQ(i,1)
case1
M=MPQ(i,2);
a=MPQ(i,3);
RA=-M/L;
QQ=QQ+RA;
MM=MM+RA*x1;
ifa>
0&
a<
L
MM=QMM(n,x1,a,M,MM);
ifa==0
MM=MM+M;
case2
P=MPQ(i,2);
b=MPQ(i,3);
RA=(L-b)*P/L;
ifb>
b<
[QQ,MM]=QMP(n,x1,b,P,QQ,MM);
case3
q=MPQ(i,2);
c=MPQ(i,3);
d=MPQ(i,4);
RA=(L-0.5*(c+d))*q*(d-c)/L;
MM=MM+RA*x1+MA;
[QQ,MM]=QMQ(n,x1,c,d,q,QQ,MM);
[QDX,MDX,XQM]=MAX_MIN(x1,QQ,MM);
TU_QM(x1,QQ,MM);
disp('
梁的有限元分析结果'
位置-----------剪力----------弯矩'
functionXQM=QMDXZ(x,MPQ)
[PN,t]=size(MPQ);
PN
MM=MM-M;
ifa==L
RA=P;
MA=-P*b;
RA=q*(d-c);
MA=-0.5*q*(d-c)*(d+c);
functionXQM=QMDXY(x,MPQ)
ifb==0
QQ=QQ-P
MM=MM-P*x1;
functionXQM=QMDWZ(x,L1,MPQ)
RA=P*(L-b)/(L-L1);
[QQ,MM]=QMP(n,x1,L1,-RA,QQ,MM);
QQ=QQ-P;
b=(c+d)*0.5;
P=(d-c)*q;
functionXQM=QMDWY(x,L1,MPQ)
RA=-M/L1;
RB=-RA;
MM=MM+x1*RA;
RA=P*(L1-b)/L1;
RB=P*b/L1;
[QQ,MM]=QMP(n,x1,L1,-RB,QQ,MM);
functionXQM=QMDWL(x,L1,L2,MPQ)
RA=-M/(L2-L1);
LL=L2-L1;
bb=b-L1;
RA=P*(LL-bb)/LL;
RB=P*bb/LL;
[QQ,MM]=QMP(n,x1,L2,-RB,QQ,MM);
untitled.m
functionvarargout=untitled(varargin)
%UNTITLEDM-fileforuntitled.fig
%UNTITLED,byitself,createsanewUNTITLEDorraisestheexisting
%singleton*.
%
%H=UNTITLEDreturnsthehandletoanewUNTITLEDorthehandleto
%theexistingsingleton*.
%UNTITLED('
CALLBACK'
hObject,eventData,handles,...)callsthelocal
%functionnamedCALLBACKinUNTITLED.Mwiththegiveninputarguments.
Property'
Value'
...)createsanewUNTITLEDorraisesthe
%existingsingleton*.Startingfromtheleft,propertyvaluepairsare
%appliedtotheGUIbeforeuntitled_OpeningFunctiongetscalled.An
%unrecognizedpropertynameorinvalidvaluemakespropertyapplication
%stop.Allinputsarepassedtountitled_OpeningFcnviavarargin.
%*SeeGUIOptionsonGUIDE'
sToolsmenu.Choose"
GUIallowsonlyone
%instancetorun(singleton)"
.
%Seealso:
GUIDE,GUIDATA,GUIHANDLES
%Edittheabovetexttomodifytheresponsetohelpuntitled
%LastModifiedbyGUIDEv2.503-Jun-200823:
12:
06
%Begininitializationcode-DONOTEDIT
gui_Singleton=1;
gui_State=struct('
gui_Name'
mfilename,...
'
gui_Singleton'
gui_Singleton,...
gui_OpeningFcn'
@untitled_OpeningFcn,...
gui_OutputFcn'
@untitled_OutputFcn,...
gui_LayoutFcn'
[],...
gui_Callback'
[]);
ifnargin&
&
ischar(varargin{1})
gui_State.gui_Callback=str2func(varargin{1});
ifnargout
[varargout{1:
nargout}]=gui_mainfcn(gui_State,varargin{:
});
else
gui_mainfcn(gui_State,varargin{:
%Endinitializationcode-DONOTEDIT
%---Executesjustbeforeuntitledismadevisible.
functionuntitled_OpeningFcn(hObject,eventdata,handles,varargin)
%Thisfunctionhasnooutputargs,seeOutputFcn.
%hObjecthandletofigure
%eventdatareserved-tobedefinedinafutureversionofMATLAB
%handlesstructurewithhandlesanduserdata(seeGUIDATA)
%varargincommandlineargumentstountitled(seeVARARGIN)
%Choosedefaultcommandlineoutputforuntitled
handles.output=hObject;
%Updatehandlesstructure
guidata(hObject,handles);
%UIWAITmakesuntitledwaitforuserresponse(seeUIRESUME)
%uiwait(handles.figure1);
%---Outputsfromthisfunctionarereturnedtothecommandline.
functionvarargout=untitled_OutputFcn(hObject,eventdata,handles)
%varargoutcellarrayforreturningoutputargs(seeVARARGOUT);
%Getdefaultcommandlineoutputfromhandlesstructure
varargout{1}=handles.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 材料力学 剪力 弯矩 绘制 详细 程序