材料力学上机报告VB.docVIP

Harbin Institute of Technology 材料力学上机大作业 题目名称： 梁的剪力、弯矩与挠度计算与绘图 班 号： 材科一班（1219001） 学 号： 1121900131 姓 名： 刘 豪 指导教师： 曲焱喆 时 间： 2014年6月11日 一、 一、 问题的提出 以VB6.0为基础，利用叠加原理计算梁在多种载荷下受力的弯曲问题，并作图。 1. 可选梁类型：简支梁与悬臂梁； 2. 可选截面：矩形截面，圆形截面，圆环截面； 3. 可选载荷：集中力，集中力偶，均布力； 4. 输出量：（n个集中力，h个集中力偶，m段均布共同作用下的任意截面的）结构图、剪力图、弯矩图、挠度图；最大剪力值、最大切应力值、最大正应力值、最大弯矩值、最大挠度值及其所在位置；任意一点的剪力值、弯矩值、挠度值（支持滚动查看）； 5. 可对常见输入错误经行判断并提醒； 二、 数学模型 将悬臂梁与简支梁杆件受力简化，逐步计算每一个载荷单独作用所产生的剪力、弯矩与挠度，再用叠加原理求出全部载荷共同作用时产生的剪力、弯矩与挠度，然后作图。 挠度计算使用近似微分方程求解所得公式计算，只适用于线弹性、小变形杆件。 三、 程序流程图 主要变量： 1.全局变量 1）动态数组：F0()、M0()、V0()叠加后的剪力值、弯矩值与挠度值；Fs()、Mz()、V()载荷单独作用时剪力值、弯矩值与挠度值；F()、q() 、M()输入的集中力、力偶与分布力大小；a1()集中力作用位置，a2()力偶作用位置，a3()均布力初始位置，b3()均布力结束位置。 2）普通变量：h 、b矩形宽度与高度；D0圆直径；D1、d2圆环外径与内径；E杨氏模量；L杆件长度；n1、n2 、n3集中力、力偶与均布力的数量；a任一点位置。 2.局部变量 Fm 、Mm 、Vm剪力、弯矩与挠度最值； c1 c2 c3 最值位置。 四、 典型算例 1.简支梁；圆环截面（外径200mm，内径100mm；长度4000mm弹性模量200Gpa；集中力1个（F=400N，位置500mm）；集中力偶2个（M1=300N·m,位置1000mm；M2=-300N·m,位置2600mm）；均布力1个（q=-200N/m,起始位置1500mm，结束位置3500mm）；任一点位置2025mm 2.悬臂梁；矩形截面（宽度100mm，高度400mm；长度5000mm弹性模量120Gpa；集中力2个（F1=400N，位置800mm；F2=-600N,位置3000mm）；集中力偶2个（M1=500N·m,位置1800mm；M2=-300N·m,位置3700mm）；均布力1个（q=-400N/m,起始位置1500mm，结束位置4500mm）；任一点位置3456mm 五、 心得体会 将近两个月时间里，我从VB6.0的基础开始学起，逐步掌握编程所需的全部知识。在后来开始编程的最初一段时间里，遇到许多困难，错误不断。因而再次利用程序课本学习，不断改正错误。在编程期间多次更改想法进行尝试，不断完善提高，一步步做到现在这个效果，总的来说还是想当满意的。下面谈谈程序的优缺点。 优点： 1.能够指导用户按照步骤进行数据输入，避免错误； 2.能够对用户输入的数据进行初步判断，并提示错误，暂停程序，让用户能够修改错误； 3.能够计算最大切应力与最大正应力，方便用户进行强度判断； 4.能够利用滚动条查询任一点的剪力、弯矩与挠度; 5.计算之前先做出结构图，让用户进一步确认输入的数据是否正确。 缺点 1.长度精确到毫米，因为当杆件长度小于50mm时，图像与最值点误差较大。 2.没有对输入非数字字符情况进行限制。 六、 附录（程序主体部分） Public F0() As Single, M0() As Single, V0() As Single, Fs() As Single, Mz() As Single, V() As Single Public F() As Single, q() As Single, M() As Single, a1() As Single, a2() As Single, a3() As Single, b3() As Single Public h!, b!, D0!, D1!, d2!, E!, n1%, n2%, n3%, L!, y%, a! Private Sub Command1_Click() Picture1.Cls Command2.Enabled = True Command2.