悬臂梁 弹性力学.doc

PAGE / NUMPAGES 《弹性理论及其工程应用》 课程三级项目说明书 学生姓名：　 李 志 鹏 专业班级： 10级工设一班 指导教师： 周 庆 田 　　 得 　 分: 　 　　 　 一、设计任务 使用matlaｂ软件对端部受集中载荷的悬臂梁进行数值分析 具体内容 对悬臂梁进行应力及位移分析,并以云图形式给出结果。 由图形结果确定梁最易折断部分. １。首先讨论梁内应力分布。 　　 其边界条件为： 　 　（σ）=0; 　 （τ)=０; 　　 （σ)＝0; 　 F＝ - σ== 　 (a) （为应力函数) 双调和方程为:＋ +＝0　(b） 通过对(a）、(b)两式积分可得: 　 (c) （d） 　2。系数的确定 　 由上述边界条件及（c）、（d）可得： 　 ; ; 　 （　为截面对中性轴的截面二次矩【惯性矩】） 　至此,所有常数均已求出,于是可得应力场: ３.然后讨论梁内位移分布 (1）应用应变位移关系及胡克定律,由应力场方程可得出: 　　　 通过对上式积分得到位移表达式： （p） 常数由阻止梁在Ｏｘy面内作刚体运动所必需的三个约束条件连确定,如在固定端（处） 有： 　 　　｝ 　 （q） 代入式（p）求得： 于是可得梁的位移场方程: 在固定端，由位移场方程可得: 　 梁轴的竖向位移为： 而端部挠度为： 上式等号右边的第一项，是我们熟悉的材料力学中所得到的悬臂梁端部挠度的结果.第二项显然是切力对挠度的影响.而这部分与弯曲的影响之比为 所以当2h〈＜时,梁的挠度主要由于弯曲所引起的.由此可见,材料力学中所得到的结果，对于细而长的梁是足够精确的。 二、设计过程（流程，代码） 1。流程 （1)确定悬臂梁各项参数 Ｅ＝2。01e５ 　弹性模量,MＰa v＝0。3 　 　泊松比 deltａ=1０　 梁厚度, mm h=15 　　 　　梁半高， ｍm leｎgth=500 　梁的跨度， mm F＝30　 　 集中载荷，N I=2/3*dｅltａ＊ｈ^3； ％截面惯性矩 （2）根据弹性力学公式编写代码 2．代码 　　 clc ｃleaｒ disｐ(梁弹性平面弯曲场变量可视化) Ｅ=2．01e5 　　%弹性模量,ＭPa ｖ=0。3 　 　 %泊松比 ｄeｌta＝10 %梁厚度，　ｍm h= 　0。５ 　 %梁半高， dm ｌengｔh=10　　 ％梁的跨度,　dm F=3０ ％集中载荷,N I=２/３＊dｅｌta＊h^３; ％截面惯性矩 Ｎx=10；Ny=20; %横纵方向离散数量 ｆoｒ ｉ=1：Ｎx+１ x=0+length／Ｎx＊(i-1)；　 ％横向离散 　ｆｏr j=１：Ny+1　 ｙ＝—h+2＊h/Ny＊（ｊ－1）； 　％纵向离散 　 Ｔ(i，j)=-F*（ｈ^2－y^２）/（２*Ｉ）； %计算切应力 　 　 fx（i，ｊ）＝-F＊x＊y／Ｉ； %计算ｘ方向应力 　 uｘ(ｉ，ｊ)=F/（E＊Ｉ)＊（leｎgth^2-x＾2）+（２+v)＊F＊y^３/（6＊E＊I）； ％计算位移ux 　 ｕy(ｉ，ｊ)=F/(E＊I)＊（x^2／６+lenｇth＾3/3+x/2＊（v*y^2-length＾2）＋ｈ^2＊（1+ｖ)*（ｌｅnｇth-x)); %计算位移uy end end x=0：length/Nx:length； ｙ=—h：2＊h/Ｎy:ｈ； [X,Ｙ］=meｓhgrid（ｘ’,y)； ｓuｂplot(2２1) suｒｆ（Ｘ,Ｙ，ux） colｏrbar ｔｉｔle（＇x方向位移） ｓｕｂploｔ（222) surf（X，Y，uy＇) coloｒｂar ｔiｔle（’y方向位移’） subploｔ（223） sｕrf（X，Y，fx’） coｌｏrbａr title（’轴线方向正应力’） sｕbｐｌｏt(２２4） surｆ（X，Y，T) titｌe（’切应力’) colｏrbａr 三、设计心得体会 采用ＭＡＴLＡB数值分析软件的三维绘图功能对简单的结构进行应力、位移分析,让计算结果直观明了，便于理性理解。计算机技术在机械设计的应用越来越广泛,分析软件就是很具有代表性的一部分,人力的解析计算耗费人力和时间，还不能保证足够的准确性，分析软件节省时间而且分析更全面，结果更直观。所以，多学习和了解各个相关方面的软件不仅有助于以后知识的学习，更重要的是帮助我们与工厂的实际接轨，增加自己的知识筹码。 书山有路勤为径,学海无涯苦作舟。我们有理由相信,当我们没有停下追逐梦想的脚步