材料力学电算大作业.doc

材料力学电算大作业 题目： 班 号： 学 号： 姓 名： 指导教师： 完成时间： 评语：成绩（满分10分）： 签名： 时间：一、题目及算法（需详细说明题目中需要用到的计算公式） （1）题目具体要求 单向应力 输入 ：x方向的应力1.工字型截面对形心轴的惯性矩Iy、 Iz及截面面积A2. 相当长度系数、长度l3. 材料参数：E、p、s；线性经验公式系数a、b4.工作压力P及稳定安全系数n输出1. 绘制临界应力总图（图示）2. 临界压力Pcr或临界应力cr3. 稳定安全校核结果 （2）算法 1.柔度、判别柔度、可用直线公式最小柔度计算 其中，为杆件长度系数，l为杆件长度，I为截面惯性矩，A为截面面积，E是弹性模量，是比例极限，a是线性经验公式中的a，b是线性经验公式中的b，是屈服极限。2判断.之间的大小关系时，用欧拉公式进行计算；时，用直线经验公式进行计算；时，校核强度；3具体校核计算时，；将与工作压力P进行比较，若P.则校核合格，反之则不合格。时，用直线经验公式进行计算，将与工作压力P进行比较，若Plamda_p %欧拉方法计算 Fcr=pi2*E*I/(mu*l)2); sigma_cr=Fcr/A; if P(Fcr/n) set(handles.edit14,String,合格); else set(handles.edit14,String,不合格); endelseif lamda_slamdalamda_p %直线经验公式 sigma_cr=a-b*lamda; if type=1 Fcr=sigma_cr*(pi*dc2/4); elseif type=2 Fcr=sigma_cr*(pi*Dr2/4-pi*dr2/4); elseif type=3 Fcr=sigma_cr*(b*h); else Fcr=sigma_cr*A; end if P(Fcr/n) set(handles.edit14,String,合格);%稳定性校验结果：合格 else set(handles.edit14,String,不合格);%稳定性校验结果：不合格 endelseif 0lamdalamda_s %校核强度 if type=1 sigma_cr=p/(pi*dc2/4); elseif type=2 sigma_cr=p/(pi*Dr2/4-pi*dr2/4); elseif type=3 sigma_cr=p/(b*h); else sigma_cr=p/A; end if sigma_crsigma_s set(handles.edit15,String,合格); else set(handles.edit15,String,不合格); endendset(handles.edit17,String,sigma_cr/1000000);%绘制临界应力总图x=0:0.1:lamda_s;y=0*x+sigma_s;plot(x,y,LineWidth,2);hold on;x=lamda_s:0.1:lamda_p;y=a-b*x;plot(x,y,LineWidth,2);hold