【正文】 (550,i*22+10,str)。i19。 CString str。 CPen* OldPen。 //歐拉法結束 printf(% %\n,phi1*180/pi,omega1[i])。 fi=(Me[i]dJe[i]*omega1[i]*omega1[i]/2)/(Je[i]*omega1[i])。 dJe[i]=2*(J2*omega2*epsilon2+m2*(Vs2x*As2x+Vs2y*As2y)+m3*Vc*Ac)。 As2y=l1*sin(phi1)omega2*omega2*ls2*sin(phi2)+epsilon2*ls2*cos(phi2)。 Vs2y=l1*cos(phi1)+omega2*ls2*cos(phi2)。 Ac=l1*(cos(phi1phi2)/cos(phi2)+ lambda*cos(phi1)*cos(phi1)/(cos(phi2)*cos(phi2)*cos(phi2)))。//ω2*=λcosΦ1/cosΦ2 epsilon2=lambda*(sin(phi1)*cos(phi2)*cos(phi2)+lambda*sin(phi2)*cos(phi1)*cos(phi1)) /(cos(phi2)*cos(phi2)*cos(phi2))。//λ=L1/L2 phi2=asin(e/l2lambda*sin(phi1))。 double fi。 (50+phi1[i]*180/pi,somega1[i]*k)。 inertia(phi1[i],l1,l2,ls2,e,J01,J2,m2,m3)。i19。//縱坐標單位 (50,somega1[0]*k)。//畫出Y坐標軸 (50,200)。//畫出X坐標軸 (50,200)。 (50,200)。 OldPen=(amp。 CClientDC dc(this)。 CPen Pen。 omega1[0]=62。m2=5。J01=3。ls2=。//定義變量l1,l2,ls2,e,J01,J2,m2,m3 l1=。 }}// The system calls this to obtain the cursor to display while the user drags// the minimized window.HCURSOR CCrankshaft2Dlg::OnQueryDragIcon(){ return (HCURSOR) m_hIcon。 // Draw the icon (x, y, m_hIcon)。 int x = (() cxIcon + 1) / 2。 GetClientRect(amp。 int cyIcon = GetSystemMetrics(SM_CYICON)。 // device context for painting SendMessage(WM_ICONERASEBKGND, (WPARAM) (), 0)。 } else { CDialog::OnSysCommand(nID, lParam)。 0xFFF0) == IDM_ABOUTBOX) { CAboutDlg dlgAbout。 // Set small icon // TODO: Add extra initialization here return TRUE。s main window is not a dialog SetIcon(m_hIcon, TRUE)。 pSysMenuAppendMenu(MF_STRING, IDM_ABOUTBOX, strAboutMenu)。 (IDS_ABOUTBOX)。 CMenu* pSysMenu = GetSystemMenu(FALSE)。 0xFFF0) == IDM_ABOUTBOX)。 //{{AFX_DATA_MAP(CCrankshaft2Dlg) // NOTE: the ClassWizard will add DDX and DDV calls here //}}AFX_DATA_MAP}BEGIN_MESSAGE_MAP(CCrankshaft2Dlg, CDialog) //{{AFX_MSG_MAP(CCrankshaft2Dlg) ON_WM_SYSCOMMAND() ON_WM_PAINT() ON_WM_QUERYDRAGICON() ON_BN_CLICKED(IDC_BUTTON1, OnButton1) ON_BN_CLICKED(IDC_BUTTON2, OnButton2) //}}AFX_MSG_MAPEND_MESSAGE_MAP()/////////////////////////////////////////////////////////////////////////////// CCrankshaft2Dlg message handlersBOOL CCrankshaft2Dlg::OnInitDialog(){ CDialog::OnInitDialog()。 //{{AFX_DATA_MAP(CAboutDlg) //}}AFX_DATA_MAP}BEGIN_MESSAGE_MAP(CAboutDlg, CDialog) //{{AFX_MSG_MAP(CAboutDlg) // No message handlers //}}AFX_MSG_MAPEND_MESSAGE_MAP()/////////////////////////////////////////////////////////////////////////////// CCrankshaft2Dlg dialogCCrankshaft2Dlg::CCrankshaft2Dlg(CWnd* pParent /*=NULL*/) : CDialog(CCrankshaft2Dlg::IDD, pParent){ //{{AFX_DATA_INIT(CCrankshaft2Dlg) // NOTE: the ClassWizard will add member initialization here //}}AFX_DATA_INIT // Note that LoadIcon does not require a subsequent DestroyIcon in Win32 m_hIcon = AfxGetApp()LoadIcon(IDR_MAINFRAME)。 // DDX/DDV support //}}AFX_VIRTUAL// Implementationprotected: //{{AFX_MSG(CAboutDlg) //}}AFX_MSG DECLARE_MESSAGE_MAP()}。// Dialog Data //{{AFX_DATA(CAboutDlg) enum { IDD = IDD_ABOUTBOX }。附錄A1總的設計程序// : implementation file//include include //include include //include //ifdef _DEBUGdefine new DEBUG_NEWundef THIS_FILEstatic char THIS_FILE[] = __FILE__。 另外我還要感謝我將要面對的答辯老師，謝謝你們在百忙之中對我給予指教。在這里，我對我的老師馮老師表示深深的感謝。馮老師經常教導我們，這次課程設計重要，但是在這次設計過程中學會做事的技巧更重要。在建立了等效模型后，通過求解得出了曲柄在某個轉角時的速度，由于單自由度機構的運動只決定于一個坐標參數，所以由此確定了整個機構的運動。在曲柄機構動力學模型建立過程中用到了單自由度機械系統等效力學模型。執(zhí)行“”后，顯示如下界面 初始界面 顯示曲柄速度曲線界面點擊按鈕“曲柄速度曲線”后，顯示如下界面 曲柄速度曲線顯示 顯示曲柄在某轉角時速度界面先點擊按鈕“曲柄速度曲線”，然后點擊“速度數值顯示”后，顯示如下界面 曲柄在某轉角速度顯示5 結 論(1)基本結論本設計在機械系統動力學的思想和理論知識的基礎上，建立了曲柄機構動力學模型。 (550,i*22+10,str)。i19。 CString str。 CPen* OldPen。 (50+phi1[i]*180/pi,somega1[i]*k)。 inertia(phi1[i],l1,l2,ls2,e,J01,J2,m2,m3)。i19。//縱坐標單位 (50,somega1[0]*k)。//畫出Y坐標軸 (50,200)。//畫出X坐標軸 (50,200)。 (50,200)。 OldPen=(amp。 CClientDC dc(this)。 CPen Pen。 omega1[0]=62。m2=5。J01=3。ls2=。//定義變量l1,l2,ls2,e,J01,J2,m2,m3 l1=。 //歐拉法結束 printf(% %\n,phi1*180/pi,omega1[i])。 fi=(Me[i]dJe[i]*omega1[i]*omega1[i]/2)/(Je[i]*omega1[i])。 dJe[i]=2*(J2*omega2*epsilon2+m2*(Vs2x*As2x+Vs2y*As2y)+m3*Vc*Ac)。 As2y=l1*sin(phi1)omega2*omega2*ls2*sin(phi2)+epsilon2*ls2*cos(phi2)。 Vs2y=l1*cos(phi1)+omega2*ls2*cos(phi2)。 Ac=l1*(cos(phi1phi2)/cos(phi2)+lambda*cos(phi1)*cos(phi1)/(cos(phi2)*cos(phi2)*cos(phi2)))。//ω2*=λcosΦ1/cosΦ2 epsilon2=lambda*(sin(phi1)*cos(phi2)*cos(phi2)+lam