MATLAB 在实际问题中的分析与应用

《MATLAB 在实际问题中的分析与应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、MATLAB在实际问题中的分析与应用班级姓名号次目录1.实验目的2.叙述问题3.分析问题4.模型建立及求解4.1塑性区（Ri£r£Rc）4.2弹性区（coR£r£R）4.2.1应力分析4.2.2计算流程4.2.3程序编写4.2.4求解结果.4.2.5计算验证4.2.6应力求解4.3应力分布曲线的绘制4.3.1程序编写4.3.2图形输出5.结果分析6.实验总结摘要：Matlab因其强大的功能而在世界范围内得以流行，特别在工程领域的应用甚为广泛。本文从弹性力学中的一个问题出发，先简要叙述并分析了目前厚壁筒弹塑性应力分析中所遇到的问题及其Matlab在其中应用的思考；后建立了厚壁筒弹塑性区的模型

2、，进行了各应力分析，用Matlab中介绍的相关方法求解了弹塑性分界面半径并绘制了各应力关于半径的变化曲线，最后根据模型求解的结果进行了简要的分析。关键词：Matlab弹塑性应力厚壁筒前言自20世纪80年代以来，出现了多种科学计算语言，亦称数学软件，比较流行的有MATLAB、Mathematica、Maple等。因为他们具有功能强、效率高、简单易学等特点，在在许多领域等到广泛应用。MATLAB便是一种影响大、流行广的科学计算语言。MATLAB的语法规则简单，更加贴近人的思维方式。MATLAB是英文MATrixLABoratory(矩阵实验室)的缩写。自1984年由美国MathWorks公司推

3、向市场以来，得到了广泛的应用和发展。在欧美各高等院校MATLAB已经成为线性代数、自动控制理论、数字信号处理、时间序列分析、动态系统仿真、图像处理等诸多课程的基本教学工具，成为大学生、硕士生以及博士生必须掌握的基本技能。在设计研究单位和工业部门，MATLAB已被广泛的应用于研究和解决各种具体的工程问题。近年来，MATLAB在我国也开始流行，应用MATLAB的单位和个人急剧增加。可以预见，MATLAB将在我国科学研究和工程应用中发挥越来越大的作用。Matlab是当前数值计算方面应用地非常广泛的一种计算机软件，特别是在工程应用求解中发挥了重要作用。其所具有的浅显易懂的编程语言、强大的绘图功能、

4、大量的内部函数等都深深地吸引了我认真地去学习它。同时在上《过程装备力学基础》时，其中涉及有很多的问题是超越方程、微积分的问题，难以用普通的线性方法求解，而Matlab在此方面有强大的功能，特别是超越方程的精确求解以及图形的绘制方面。数学当中的绘制函数图象、绘制立体图形的交线（如绘制两个等直径圆柱体的交线）、求多项式的根等问题，这些问题如果依靠我们人工进行操作，则需要很多的时间和精力，当我们掌握了基本原理后，借助于MATLAB进行解决则会大大提高效率和精确度。1.实验目的1结合实际问题展现MATLAB在生活和学习方面的广泛应用2学会利用MATLAB编程并求解实际问题3学会并运用for循环和s

5、witch结构，以及MATLAB中已有函数如sum4了解单元数组cell和结构数组struct的作用，学会创建这些数组.2、叙述问题厚壁筒在承受内压载荷的作用下，随着压力的增加，筒壁应力不断增加。厚壁筒在承受逐渐增加压力的过程中，会经历弹性阶段、筒体部分屈服阶段、整体屈服阶段、材料硬化、筒体过度变形、直至爆破失效阶段。而在分析厚壁筒的弹塑性应力分布时，遇到两个问题：【1】弹塑性区分界面的半径的精确确定较为困难，因此半径是一个非线性方程；【2】弹塑性应力的分布曲线绘制难以精确，不能真实反映不同半径处的应力状况。3、分析问题针对厚壁筒弹塑性应力分析时遇到的两个问题，可利用Matlab软件中的相

6、关知识解决。对于弹塑性区分界面半径的确定，可利用方程求根的相关方法，从二分法、开方法、Newton法、Newton下山法以及弦截法中选取精度和收敛速度均较佳的方程求根数值方法进行求解；同时还可以利用Matlab中的Solve函数求其精确解。对于弹塑性应力分布曲线的绘制，可以根据求出的各应力分布的具体方程，用Plot函数进行绘制，同时要注意曲线后期的处理工作，以便更好地从图中得到各应力随半径分布的情况。4、模型建立及求解图2受内压厚壁圆筒图3弹性区图4塑性区为简化分析，假设厚壁圆筒为理想弹塑性体，不考虑材料在塑性变形过程中塑性强化，筒体仅受内压pi作用，筒体的内半径为Ri，外半径为Ro。初始

7、假设厚壁圆筒的内半径为Ri=152.5mm，外半径为Ro=254mm，内壁所受的压力为pi=340MPa，取其屈服极限为s=750MPa先分析筒体仅受内压Pi作用的情形。当内压pi大于弹性极限压力pe时，圆筒内壁的屈服区向外扩展，筒体沿壁可以分成塑性区和弹性区两个区域，其中内侧为塑性区，外侧为弹性区。假想两区域的交界圆面的半径为Rc，则塑性区的内外半径分别为Ri和Rc，承受的内外压力分别为Pi和Pc，见图4；弹性区的内外