哈工大材力上机实验报告

1、材料力学电算实验压杆的临界力计算 院 系： 机电工程学院 班 级： 设 计 者： 学 号： 指导教师： 张桂莲 软件要求： 设计时间： 一.概述：本程序使用Microsoft Visual Basic编写，可以对不同材料、不同约束类型、不同截面类型的压杆进行临界力的计算。杆件的参数可以输入，得出结果之后也可以清零。二、问题分析及相关公式：1、压杆稳定当短粗杆受压时(图1)，在压力F由小逐渐增大的过程中，杆件始终保持原有的直线平衡形式，直到压力F达到屈服强度载荷Fs (或抗压强度载荷Fb)，杆件发生强度破坏时为止。但是，如果用相同的材料，做一根与图1a所示的同样粗细而比较长的杆件(图1b)，当压

2、力F比较小时，这一较长的杆件尚能保持直线的平衡形式，而当压力F逐渐增大至某数值F1时，杆件将突然变弯，不再保持原有的直线平衡形式，因而丧失了承载能力。我们把受压直杆突然变弯的现象，称为丧失稳定或失稳。此时，F1可能远小于Fs (或Fb)。可见，细长杆在尚未产生强度破坏时，就因失稳而破坏。图1在研究压杆稳定时，我们用一微小横向干扰力使处于直线平衡状态的压杆偏离原有的位置，如图1所示。当轴向压力F由小变大的过程中，可以观察到：1)当压力值F1较小时，给其一横向干扰力，杆件偏离原来的平衡位置。若去掉横向干扰力后，压杆将在直线平衡位置左右摆动，最终将恢复到原来的直线平衡位置。2)当压力值F2超过其一限

3、度Fcr时，平衡状态的性质发生了质变。这时，只要有一轻微的横向干扰，压杆就会继续弯曲，不再恢复原状，。 3)界于前二者之间，存在着一种临界状态。当压力值正好等于Fcr时，一旦去掉横向干扰力，压杆将在微弯状态下达到新的平衡，既不恢复原状，也不再继续弯曲，。临界状态是杆件从稳定平衡向不稳定平衡转化的极限状态。压杆处于临界状态时的轴向压力称为临界力或临界载荷，用Fcr表示。2、 两端铰支细长压杆的临界力图2为一两端为球形铰支的细长压杆，其临界力公式为：图2 (1)式(1)又称为欧拉公式。3、不同杆端约束细长压杆的临界力（1）一端固定另一端自由细长压杆的临界力图3为端固定另一端自由的压杆。当压杆处于临

4、界状态时，它在曲线形式下保持平衡。将挠曲线AB对称于固定端A向下延长，如图中假想线所示。延长后挠曲线是一条半波正弦曲线，与本章第二节中两端铰支细长压杆的挠曲线一样。所以，对于端固定另一端自由且长为的压杆，其临界力等于两端铰支长为的压杆的临界力，即图3 图4 图5（2）两端固定细长压杆的临界力在这种杆端约束条件下，挠曲线如图4所示。该曲线的两个拐点C和D分别在距上、下端为处。居于中间的长度内，挠曲续是半波正弦曲线。所以，对于两端固定且长为的压杆，其临界力等于两端铰支长为的压杆的临界力，即（3） 一端固定另一端铰支细长压杆的临界力在这种杆端约束条件下，挠曲线形状如图5所示。在距铰支端B为处，该曲线

5、有一个拐点C。因此，在长度内，挠曲线是一条半波正弦曲线。所以，对于一端固定另一端铰支且长为的压杆，其临界力等于两端铰支长为的压杆的临界力，即综上所述，只要引入相当长度的概念，将压杆的实际长度转化为相当长度，便可将任何杆端约束条件的临界力统一写 (2)称为欧拉公式的一般形式。由式(2)可见，杆端约束对临界力的影响表现在系数上。称为长度系数，为压杆的相当长度，表示把长为的压杆折算成两端铰支压杆后的长度。几种常见约束情况下的长度系数列入表1中。表1 压杆的长度系数压杆的约束条件长度系数两端铰支一端固定，另一端自由两端固定一端固定，另一端铰支=1=2=1/20.7三、VB程序源代码/ 111Dlg.c

6、pp : implementation file/#include "stdafx.h"#include "111.h"#include "111Dlg.h"#include<math.h>#ifdef _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endif/ CAboutDlg dialog used for App Aboutclass CAboutDlg : public CDialogpublic:CAbout

7、Dlg();/ Dialog Data/AFX_DATA(CAboutDlg)enum IDD = IDD_ABOUTBOX ;/AFX_DATA/ ClassWizard generated virtual function overrides/AFX_VIRTUAL(CAboutDlg)protected:virtual void DoDataExchange(CDataExchange* pDX); / DDX/DDV support/AFX_VIRTUAL/ Implementationprotected:/AFX_MSG(CAboutDlg)/AFX_MSGDECLARE_MESSA

8、GE_MAP();CAboutDlg:CAboutDlg() : CDialog(CAboutDlg:IDD)/AFX_DATA_INIT(CAboutDlg)/AFX_DATA_INITvoid CAboutDlg:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);/AFX_DATA_MAP(CAboutDlg)/AFX_DATA_MAPBEGIN_MESSAGE_MAP(CAboutDlg, CDialog)/AFX_MSG_MAP(CAboutDlg)/ No message handlers/AFX_MSG_MA

9、PEND_MESSAGE_MAP()/ CMy111Dlg dialogCMy111Dlg:CMy111Dlg(CWnd* pParent /*=NULL*/): CDialog(CMy111Dlg:IDD, pParent)/AFX_DATA_INIT(CMy111Dlg)m_A1 = FALSE;m_A2 = FALSE;m_A3 = FALSE;m_A4 = FALSE;m_A5 = FALSE;m_B1 = FALSE;m_B2 = FALSE;m_B3 = FALSE;m_C1 = FALSE;m_C2 = FALSE;m_C3 = FALSE;m_b = 0.0f;m_h = 0.

10、0f;m_DD = 0.0f;m_D = 0.0f;m_d = 0.0f;m_l = 0.0f;m_r = 0.0;m_Fcr = 0.0;/AFX_DATA_INIT/ Note that LoadIcon does not require a subsequent DestroyIcon in Win32m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);void CMy111Dlg:DoDataExchange(CDataExchange* pDX)CDialog:DoDataExchange(pDX);/AFX_DATA_MAP(CMy1

11、11Dlg)DDX_Check(pDX, IDC_CHECK1, m_A1);DDX_Check(pDX, IDC_CHECK2, m_A2);DDX_Check(pDX, IDC_CHECK3, m_A3);DDX_Check(pDX, IDC_CHECK4, m_A4);DDX_Check(pDX, IDC_CHECK5, m_A5);DDX_Check(pDX, IDC_CHECK6, m_B1);DDX_Check(pDX, IDC_CHECK7, m_B2);DDX_Check(pDX, IDC_CHECK8, m_B3);DDX_Check(pDX, IDC_CHECK9, m_C

12、1);DDX_Check(pDX, IDC_CHECK10, m_C2);DDX_Check(pDX, IDC_CHECK11, m_C3);DDX_Text(pDX, IDC_EDIT1, m_b);DDX_Text(pDX, IDC_EDIT2, m_h);DDX_Text(pDX, IDC_EDIT3, m_DD);DDX_Text(pDX, IDC_EDIT4, m_D);DDX_Text(pDX, IDC_EDIT5, m_d);DDX_Text(pDX, IDC_EDIT6, m_l);DDX_Text(pDX, IDC_EDIT7, m_r);DDX_Text(pDX, IDC_

13、EDIT8, m_Fcr);/AFX_DATA_MAPBEGIN_MESSAGE_MAP(CMy111Dlg, CDialog)/AFX_MSG_MAP(CMy111Dlg)ON_WM_SYSCOMMAND()ON_WM_PAINT()ON_WM_QUERYDRAGICON()ON_BN_CLICKED(IDC_BUTTON1, OnButton1)ON_BN_CLICKED(IDC_BUTTON2, OnButton2)ON_BN_CLICKED(IDC_BUTTON3, OnButton3)ON_BN_CLICKED(IDC_CHECK1, OnCheck1)ON_BN_CLICKED(I

14、DC_CHECK10, OnCheck10)ON_BN_CLICKED(IDC_CHECK11, OnCheck11)ON_BN_CLICKED(IDC_CHECK2, OnCheck2)ON_BN_CLICKED(IDC_CHECK3, OnCheck3)ON_BN_CLICKED(IDC_CHECK4, OnCheck4)ON_BN_CLICKED(IDC_CHECK5, OnCheck5)ON_BN_CLICKED(IDC_CHECK6, OnCheck6)ON_BN_CLICKED(IDC_CHECK7, OnCheck7)ON_BN_CLICKED(IDC_CHECK8, OnChe

15、ck8)ON_BN_CLICKED(IDC_CHECK9, OnCheck9)/AFX_MSG_MAPEND_MESSAGE_MAP()/ CMy111Dlg message handlersBOOL CMy111Dlg:OnInitDialog()CDialog:OnInitDialog();/ Add "About." menu item to system menu./ IDM_ABOUTBOX must be in the system command range.ASSERT(IDM_ABOUTBOX & 0xFFF0) = IDM_ABOUTBOX);A

16、SSERT(IDM_ABOUTBOX < 0xF000);CMenu* pSysMenu = GetSystemMenu(FALSE);if (pSysMenu != NULL)CString strAboutMenu;strAboutMenu.LoadString(IDS_ABOUTBOX);if (!strAboutMenu.IsEmpty()pSysMenu->AppendMenu(MF_SEPARATOR);pSysMenu->AppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu);/ Set the icon for th

17、is dialog. The framework does this automatically/ when the application's main window is not a dialogSetIcon(m_hIcon, TRUE);/ Set big iconSetIcon(m_hIcon, FALSE);/ Set small icon/ TODO: Add extra initialization herereturn TRUE; / return TRUE unless you set the focus to a controlvoid CMy111Dlg:OnS

18、ysCommand(UINT nID, LPARAM lParam)if (nID & 0xFFF0) = IDM_ABOUTBOX)CAboutDlg dlgAbout;dlgAbout.DoModal();elseCDialog:OnSysCommand(nID, lParam);/ If you add a minimize button to your dialog, you will need the code below/ to draw the icon. For MFC applications using the document/view model,/ this

19、is automatically done for you by the framework.void CMy111Dlg:OnPaint() if (IsIconic()CPaintDC dc(this); / device context for paintingSendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);/ Center icon in client rectangleint cxIcon = GetSystemMetrics(SM_CXICON);int cyIcon = GetSystemMetrics(SM

20、_CYICON);CRect rect;GetClientRect(&rect);int x = (rect.Width() - cxIcon + 1) / 2;int y = (rect.Height() - cyIcon + 1) / 2;/ Draw the icondc.DrawIcon(x, y, m_hIcon);elseCDialog:OnPaint();/ The system calls this to obtain the cursor to display while the user drags/ the minimized window.HCURSOR CMy

21、111Dlg:OnQueryDragIcon()return (HCURSOR) m_hIcon;void CMy111Dlg:OnButton1() /计算柔度UpdateData();if(m_B1)/两端固定，长度因数=0.5if(m_C1)/矩形截面m_r=1.732*m_l/m_h;/柔度=（l）/i,其中惯性矩i=h/(2*1.732)if(m_C2)/实心圆柱截面m_r=2*m_l/m_DD;/惯性矩i=D/4if(m_C3)/空心圆柱截面m_r=2*m_l/sqrt(m_D*m_D+m_d*m_d);/惯性矩i=sqrt(D*D+d*d)/4if(m_B2)/一端固定，一端铰支

22、，长度因数=0.7if(m_C1)m_r=1.732*1.4*m_l/m_h;if(m_C2)m_r=2.8*m_l/m_DD;if(m_C3)m_r=2.8*m_l/sqrt(m_D*m_D+m_d*m_d);if(m_B3)/两端铰支，长度因数=1if(m_C1)m_r=1.732*2*m_l/m_h;if(m_C2)m_r=4*m_l/m_DD;if(m_C3)m_r=4*m_l/sqrt(m_D*m_D+m_d*m_d);UpdateData(FALSE);void CMy111Dlg:OnButton2()/计算临界力UpdateData();if(m_A1)/材料为硅钢，p=100

23、，s=60，E=150GPa，a=578MPa，b=3.744MPaif (m_r>=100)/>=p，为细长杆，用欧拉公式if(m_C1)/矩形截面m_Fcr=9.87*150000*m_b*m_h/(m_r*m_r);/欧拉公式cr=*E/(*)，Fcr=cr*A，其中A为截面面积m_Fcr=7.75*150000*m_DD*m_DD/(m_r*m_r);if(m_C3)/空心圆柱截面m_Fcr=7.75*150000*(m_D*m_D-m_d*m_d)/(m_r*m_r);if(m_r<100)&&(m_r>60)/S<<p,为中长杆，

24、用经验公式if(m_C1)/矩形截面m_Fcr=(578-3.744*m_r)*m_b*m_h;/经验公式cr=a-b*，Fcr=cr*A，其中A为截面面积if(m_C2)/实心圆柱截面m_Fcr=(578-3.744*m_r)*0.785*m_DD*m_DD;if(m_C3)/空心圆柱截面m_Fcr=(578-3.744*m_r)*0.785*(m_D*m_D-m_d*m_d);if(m_r<=60)/<=s，为短粗杆，用强度理论if(m_C1)/矩形截面m_Fcr=353*m_b*m_h;/cr=S,其中S为屈服极限，Fcr=cr*A，其中A为截面面积if(m_C2)/实心圆柱

25、截面m_Fcr=353*0.785*m_DD*m_DD;if(m_C3)/空心圆柱截面m_Fcr=353*0.785*(m_D*m_D-m_d*m_d);if(m_A2)/材料为Q235钢，p=100，s=57，E=200GPa，a=304MPa，b=1.12MPaif (m_r>=100)if(m_C1)m_Fcr=9.87*200000*m_b*m_h/(m_r*m_r);if(m_C2)m_Fcr=7.75*200000*m_DD*m_DD/(m_r*m_r);if(m_C3)m_Fcr=7.75*200000*(m_D*m_D-m_d*m_d)/(m_r*m_r);if(m_r&

26、lt;100)&&(m_r>57)if(m_C1)m_Fcr=(304-1.12*m_r)*m_b*m_h;if(m_C2)m_Fcr=(304-1.12*m_r)*0.785*m_DD*m_DD;if(m_C3)m_Fcr=(304-1.12*m_r)*0.785*(m_D*m_D-m_d*m_d);if(m_r<=57)if(m_C1)m_Fcr=240.16*m_b*m_h;if(m_C2)m_Fcr=240.16*0.785*m_DD*m_DD;if(m_C3)m_Fcr=240.16*0.785*(m_D*m_D-m_d*m_d);if(m_A3)/材料为

27、优质碳钢，p=86，s=60，E=200GPa，a=461MPa，b=2.568MPaif (m_r>=86)if(m_C1)m_Fcr=9.87*200000*m_b*m_h/(m_r*m_r);if(m_C2)m_Fcr=7.75*200000*m_DD*m_DD/(m_r*m_r);if(m_C3)m_Fcr=7.75*200000*(m_D*m_D-m_d*m_d)/(m_r*m_r);if(m_r<86)&&(m_r>60)if(m_C1)m_Fcr=(461-2.568*m_r)*m_b*m_h;if(m_C2)m_Fcr=(461-2.568*m

28、_r)*0.785*m_DD*m_DD;if(m_C3)m_Fcr=(461-2.568*m_r)*0.785*(m_D*m_D-m_d*m_d);if(m_r<=60)if(m_C1)m_Fcr=306*m_b*m_h;if(m_C2)m_Fcr=306*0.785*m_DD*m_DD;if(m_C3)m_Fcr=306*0.785*(m_D*m_D-m_d*m_d);if(m_A4)/材料为铝合金，p=50，s=0，E=70GPa，a=372MPa，b=2.14MPaif (m_r>=50)if(m_C1)m_Fcr=9.87*70000*m_b*m_h/(m_r*m_r);i

29、f(m_C2)m_Fcr=7.75*70000*m_DD*m_DD/(m_r*m_r);if(m_C3)m_Fcr=7.75*70000*(m_D*m_D-m_d*m_d)/(m_r*m_r);if(m_r<50)/因铝合金较软，即s=0，故不会被压断，只会产生失稳变形。if(m_C1)m_Fcr=(372-2.14*m_r)*m_b*m_h;if(m_C2)m_Fcr=(372-2.14*m_r)*0.785*m_DD*m_DD;if(m_C3)m_Fcr=(372-2.14*m_r)*0.785*(m_D*m_D-m_d*m_d);if(m_A5)/材料为铸铁，E=150GPa，a=578MPa，b=3.744MPa，强度极限S=400MPa；因铸铁为脆性材料，故只能被压断，不会产生失稳变形。if(m_C1)m_Fcr=400*m_b*m_h;/cr=S=400MPa，Fcr=cr*A，其中A为截面面积if(m_C2)m_Fcr=400*0.785*m_DD*m_DD;if(m_C3)m_Fcr=400*0.785*(m_D*m_D-m_d*m_d);UpdateData(FALSE);v