基于灵敏度分析的有限元模型修正

回归模型结果分析

回归模型结果分析为了提高回归模型的准确性，上文中我们分别按月份、颜色比、退偏振比三种情况进行回归建模，从以上的分析结果看来，按月份划分建立的回归模型反演效果较好。为了更好地对不同情况下得到的回归模型及反演结果进行对比，我们把相同情况下得到的所有反演结果表示在一张图上，并与相应的太阳光度计观测值进行对比分析。（a）

（b）（c）

图4.1 图4.1中（a）、（b）、（c）三幅图为分别按月份、颜色比和退偏振比建立回归模型后得出的所有颗粒物体积浓度的反演结果与相应太阳光度计观测值的对比分析图。图（a）数据的样本容量为250，图（b）和图（c）的样本容量为150，虽然图（a）样本容量多，但是与图（b）和图（c）相比，图（a）中数据更为集中，大部分数据的反演结果与太阳光度计观测值接近，出现误差的数据少且误差小，图（c）的反演结果略优于图（b），总体来说按月份建立的颗粒物体积浓度的回归模型最准确，而按颜色比建立的回归模型准确性较差。（a）

（b）（c）图4.2

图4.2中(a)、（b）、（c）三幅图为分别按月份、颜色比和退偏振比建立回归模型后得出的所有有效粒子半径的反演结果与相应太阳光度计观测值的对比分析图。图（a）样本容量较多且数据比较集中，但有一部分数据反演结果明显偏小，严重影响了回归模型的准确性，图（b）数据较离散，部分数据误差大，线性相关系数较小，图（c）个别数据误差大，虽然数据集中程度没有图（a）好。但是数据横纵坐标的差异比其他两幅图小。在确定最优样本容量时，我们发现随着样本容量的增加，线性相关系数减小，所以在无法统一样本容量且线性相关系数差异不大的情况下无法确定在哪种情况下建立的回归模型最准确。所以在建立有效粒子半径的回归模型时，我们可以按月份建立回归模型，也可以按退偏振比建立回归模型。

回归模型分析

新疆财经大学实验报告课程名称：统计学实验项目名称：回归模型分析姓名： lili 学号： 20000000 班级：工商2011-2班指导教师： 2014 年5 月

新疆财经大学实验报告

附:实验数据。

1、作散点图,加趋势线， 2、建立回归模型（用公式编辑器写），对模型进行统计检验。解释模型意义SUMMARY OUTPUT 回归统计 Multiple R 0.974111881 R Square 0.948893956 Adjusted R Square 0.947131679 标准误差527.4648386 观测值31 方差分析 df SS MS F Significance F 回归分析 1 149806425.5 149806426 538.4476 2.82E-20 残差29 8068355.522 278219.156 总计30 157874781.1 Coefficients 标准误差t Stat P-value Lower 95% Upper 95% Intercept 121.5246471 365.0193913 0.33292655 0.741585 -625.024 X Variable 1 1.270433698 0.054749518 23.2044728 2.82E-20 1.158458

RESIDUAL OUTPUT 观测值预测 Y 残差标准残差 1 14252.56 -369.959 -0.71338 2 10116.66 196.2382 0.378401 3 7032.43 206.6701 0.398516 4 6607.597 412.4032 0.795225 5 7006.005 6.895144 0.013296 6 7843.094 -602.494 -1.16177 7 7098.874 -93.6736 -0.18063 8 6493.004 185.8963 0.358458 9 14147.49 720.0062 1.388367 10 8644.356 618.1438 1.191949 11 12461.12 717.8799 1.384267 12 6555.382 244.618 0.47169 13 9467.216 532.2839 1.026388 14 6365.198 536.2019 1.033943 15 7832.295 567.6051 1.094497 16 6399.5 526.5002 1.015235 17 7697.502 -375.502 -0.72407 18 7871.17 -171.17 -0.33006 19 12363.8 16.59511 0.032 20 7443.669 341.3307 0.658178 21 7111.959 147.341 0.284113 22 9164.599 -1070.9 -2.06498 23 7490.04 -448.14 -0.86414 24 6408.901 160.099 0.308714 25 7774.109 -130.509 -0.25166 26 10342.54 -1577.04 -3.04097 27 7362.997 -462.997 -0.89278 28 6852.282 -195.082 -0.37617 29 6982.121 -236.821 -0.45665 30 6893.317 -362.817 -0.69961 31 7260.6 -39.5998 -0.07636 y=β0+β1x y=121.225+1.27X 3、求相关系数与方向说明数意根据以上的结果，0《r≤1,这表明x与y之间正线性相关，因为r=0.9741可视为高度相关；

有限元分析报告样本

《有限元分析》报告基本要求： 1. 以个人为单位完成有限元分析计算，并将计算结果上交；（不允许出现相同的分析模型，如相同两人均为不及格） 2. 以个人为单位撰写计算分析报告； 3. 按下列模板格式完成分析报告； 4. 计算结果要求提交电子版，报告要求提交电子版和纸质版。（以上文字在报告中可删除）《有限元分析》报告一、问题描述（要求：应结合图对问题进行详细描述，同时应清楚阐述所研究问题的受力状况和约束情况。图应清楚、明晰，且有必要的尺寸数据。）一个平面刚架右端固定，在左端施加一个y 方向的-3000N 的力P1，中间施加一个Y 方向的-1000N 的力P2，试以静力来分析，求解各接点的位移。已知组成刚架的各梁除梁长外，其余的几何特性相同。横截面积：A=0.0072 m2 横截高度：H=0.42m 惯性矩：I=0.0021028m4ｘ弹性模量： E=2.06ｘ10n/ m2/ 泊松比：u=0.3 二、数学模型（要求：针对问题描述给出相应的数学模型，应包含示意图，示意图中应有必要的尺寸数据；如进行了简化等处理，此处还应给出文字说明。）（此图仅为例题）

三、有限元建模（具体步骤以自己实际分析过程为主，需截图操作过程）用ANSYS 分析平面刚架 1．设定分析模块选择菜单路径：MainMenu—preference 弹出“PRreferences for GUI Filtering”对话框，如图示，在对话框中选取：Structural”,单击[OK]按钮，完成选择。 2．选择单元类型并定义单元的实常数（1）新建单元类型并定（2）定义单元的实常数在”Real Constants for BEAM3”对话框的AREA中输入“0。0072”在IZZ 中输入“0。0002108”，在HEIGHT中输入“0.42”。其他的3个常数不定义。单击[OK]按钮，完成选择 3．定义材料属性在”Define Material Model Behavier”对话框的”Material Models Available”中，依次双击“Structural→Linear→Elastic→Isotropic”如图

多元线性回归模型的案例分析

1. 表1列出了某地区家庭人均鸡肉年消费量Y 与家庭月平均收入X ，鸡肉价格P 1，猪肉价格P 2与牛肉价格P 3的相关数据。年份 Y/千克 X/ 元 P 1/(元/千克) P 2/(元/千克) P 3/(元/千克) 年份 Y/千克 X/元 P 1/(元/ 千克) P 2/(元/ 千克) P 3/(元/千克) 1980 2.78 397 4.22 5.07 7.83 1992 4.18 911 3.97 7.91 11.40 1981 2.99 413 3.81 5.20 7.92 1993 4.04 931 5.21 9.54 12.41 1982 2.98 439 4.03 5.40 7.92 1994 4.07 1021 4.89 9.42 12.76 1983 3.08 459 3.95 5.53 7.92 1995 4.01 1165 5.83 12.35 14.29 1984 3.12 492 3.73 5.47 7.74 1996 4.27 1349 5.79 12.99 14.36 1985 3.33 528 3.81 6.37 8.02 1997 4.41 1449 5.67 11.76 13.92 1986 3.56 560 3.93 6.98 8.04 1998 4.67 1575 6.37 13.09 16.55 1987 3.64 624 3.78 6.59 8.39 1999 5.06 1759 6.16 12.98 20.33 1988 3.67 666 3.84 6.45 8.55 2000 5.01 1994 5.89 12.80 21.96 1989 3.84 717 4.01 7.00 9.37 2001 5.17 2258 6.64 14.10 22.16 1990 4.04 768 3.86 7.32 10.61 2002 5.29 2478 7.04 16.82 23.26 1991 4.03 843 3.98 6.78 10.48 （1）求出该地区关于家庭鸡肉消费需求的如下模型： 01213243ln ln ln ln ln Y X P P P u βββββ=+++++ （2）请分析，鸡肉的家庭消费需求是否受猪肉及牛肉价格的影响。先做回归分析，过程如下：输出结果如下：

GARTEUR 有限元模型修正与确认研究

收稿日期:2003207207;修订日期:2004203225 基金项目:教育部博士学科点专项基金(20010227012)资助项目文章编号:100026893(2004)0420372204 GARTEUR 有限元模型修正与确认研究费庆国,张令弥,郭勤涛 (南京航空航天大学振动工程研究所,江苏南京 210016) Case Study of FE Model Updating and Validation via an Air craft Model Structur e FEI Qing 2guo,Z HANG Ling 2mi,G UO Qin 2tao (Institu te o f Vi brati on Engi neering,Nanjing University of Aeronau tics and Astro nautics,Nanjing 210016,China)摘要:待修正参数的选择以及修正后模型的质量评估是有限元模型修正的两个重要问题。以欧洲学术界广泛采用的GA RTEUR 飞机模型为例,利用基于灵敏度分析的模型修正方法,通过仿真算例研究参数选择对模型修正质量的影响,并以试验数据为目标值对有限元模型进行修正与确认。为全面评估模型的修正质量,引入三级标准对修正后有限元模型进行确认。关键词:固体力学;模型确认;有限元法;模型修正;参数选择中图分类号:O 248121 文献标识码:A Abstr act:Parameter selection and quali ty validation are of g reat i mpo rtance in fini te element model updating.This paper presents so me results which demonstrate the relationship betw een parameter selection and updated model .s quality throu gh si mulation cases.Three q uali ty levels w ith corresponding validation criteria are emplo yed with an emphasis o n updated mod 2el .s predictio n ability.Results of updating based on exper i mental modal test data are sho w n as an application example.A n aircraft test structure,GA RTEUR,which is g enerally utilized in Europe,is employed in bo th the si mulation case and the exper i mental case.Sensi tivity 2based model updating appro ach is applied. Key wor ds:solid mechanics;model validation;finite element method;model updating;parameter selection 在航空工程中,准确的有限元模型对于动态响应预测以及动态设计至关重要。建模过程中的不确定因素,如离散化误差、材料物理参数的不确定性、边界条件的近似等,导致有限元模型必然存在误差。设计规范规定,有限元模型必须通过振动模态试验或者地面共振试验来检验[1]。近30年来,有限元模型修正技术得到了长足的发展[2~6]。根据修正对象的不同可将修正方法分为矩阵型方法和设计参数型方法。后者物理意义明确,更具工程应用价值。本文采用基于灵敏度分析的设计参数型修正方法。基于灵敏度分析的设计参数型修正方法主要包括待修正设计参数选择,灵敏度分析,参数修正以及模型确认等环节。待修正设计参数的选择是模型修正的起始环节。通常,候选参数是有限元模型存在不确定性因素的参数。近20年虽然发展了很多种参数选择或者误差定位的策略与算法,工程应用中仍然难以准确无遗漏地确定误差参数。因此,有必要讨论参数选择对模型修正质量的影响。模型确认是模型修正的检验环节。在当前的研究及工程应用中,通常只要求修正后模型的计算结果能够复现修正过程中利用的试验数据。事实上,为全面评估模型的质量,模型的复现能力与预测能力应予以同等重视[7]。本研究引入了三级质量标准对修正后的有限元模型进行确认。本文采用G ARTE UR 飞机模型为研究对象,通过仿真算例来研究参数选择对模型修正质量的影响,并给出了利用振动模态测试结果对G AR 2TEUR 飞机模型的有限元模型进行修正与确认的结果。 1 模型修正方法与模型确认准则 (1)模型修正方法模型修正可归结为以下的优化问题[8] Min p +R(p )+2 2,R(p )=f E -f A (p )s.t V L [p [V U (1) 其中:p 代表设计参数;f E ,f A 是结构动态特性试验与分析结果;R 代表残差;V L ,V U 是设计参数的下、上限。令设计参数的初始值为p 0,动态特性f 是设计参数的隐函数,其泰勒展开式为第25卷第4期航空学报 Vol 125N o 14 2004年 7月ACT A AERO NA U TICA E T AS TRO NA U TICA SINICA July 2004

板结构有限元分析实例详解

板结构有限元分析实例详解1：带孔平板结构静力分析本节介绍带孔平板结构静力分析问题，同时介绍布尔操作的基本用法。 8.3.1 问题描述与分析有孔的矩形平板，左侧边缘固定，长400mm，宽200 mm，厚度为10 mm，圆孔在板的正中心，半径为40 mm，左侧全约束，右侧边缘均布应力1MPa，如图8.7所示。求板的变形、位移及应力变化情况。（材料的材料属性为：弹性模量为300000 MPa，剪切模量为0.31。）图8.7 带孔的矩形平板由于小孔处边缘不规则，本文采用PLANE82高阶平面单元进行分析。 8.3.2 求解过程 8.3.2.1 定义工作目录及文件名启动ANSYS Mechanical APDL Product Launcher窗口，如图8.8所示。在License下拉选框中选择ANSYS Multiphysics产品，在Working Directory输入栏中输入工作目录：C:\ANSYS12.0 Structural Finite Elements Analysis and Practice\Chapter 8\8-1，在Job Name一栏中输入工作文件名：Chapter8-1。以上参数设置完毕后，单击Run按钮运行ANSYS。

图8.8 ANSYS设置窗口菜单可以先在目标文件位置建立工作目录，然后单击Browse按钮选择工作目录；也可以通过单击Browse按钮选择工作文件名。 8.3.2.2 定义单元类型和材料属性选择Main Menu>Preferences命令，出现Preferences for GUI Filtering对话框，如图8.9所示，在Individual discipline(s) to show in the GUI中勾选Structural，过滤掉ANSYS GUI菜单中与结构分析无关的选项，单击OK按钮关闭该对话框。图8.9 Preferences for GUI Filtering对话框

汽车车架简化模型有限元分析

汽车典型零部件简化模型有限元分析任务1：连杆简化模型的有限元分析 1. 分析任务：对图一所示的连杆的二维简化模型进行有限元分析，确定该设计是否满足结构的强度要求；若强度不够，修改设计直至最大应力减小至材料允许的范围内。在修改结构时，注意不可改变连杆小头衬套的内径和连杆大头的内径，也不可改变连杆各处的厚度和材料。 2. 分析所需数据： a．连杆采用两种材料，连杆本体用的是40Cr结构钢，左侧小头中的衬套用的是铜。 b．连杆杆身和大头的厚度为1.5mm，小头的厚度为3.0mm。注意在杆身和小头的过渡处有R2.0的过渡圆角； c．连杆结构的其它尺寸如图二所示； d．施加在大、小头内壁上的边界条件用于模拟连杆与曲轴及活塞销的连接。假定载荷分布在小头夹角为90o的内壁上，且为锥状分布；约束施加在连杆大头夹角为90o的内壁上； e．40Cr材料的弹性模量：210GPa；泊松比：0.3；屈服极限为：850MPa，设计安全系数为6；铜的弹性模量：120GPa，泊松比：0.33；屈服极限为：250MPa；设计安全系数为4。 3. 完成该分析应掌握的ANSYS技术： a．单元类型的选择；单元的尺寸控制； b．不同厚度和材料的二维实体建模； c．工作平面的灵活应用；

d．按载荷和约束的要求分割线和面； e．模型参数（材料，实常数，单元类型号等） f．粘结、合并等布尔运算操作 g．局部坐标系，旋转节点坐标系； h．线性分布载荷的施加； i．单元网格误差估计； j．Ansys 命令日志文件及其在修改设计中的应用； k．多窗口显示的功能 4. 分析报告内容的基本要求： a.对分析任务的描述；列出分析所需数据： b.利用多窗口显示的功能绘出连杆的实体模型和网格模型，在模型上能反映出连杆各部位材料、厚度的不同； c.绘图反映连杆的边界条件； d.绘出对连杆原设计进行有限元分析后得到的变形图和应力等值线图； e.图示SEPC和SERR并说明有限元分析的建模误差； f.详细说明对不符合设计要求的结构所作的设计修改；及最终符合设计要求的计算结果； g.在分析中遇到的关键问题（在实体建模、网格剖分、边界条件施加等各个步骤中出现的）及解决的办法； h.整理命令日志文件，并在每个语句后添加说明（说明该语句的功能，说明前要加！号）。注意：添加的说明（可以用中文说明）应该反映在建模中的操作步骤而不是简单的ANSYS命令定义。

回归分析方法

第八章回归分析方法当人们对研究对象的内在特性和各因素间的关系有比较充分的认识时，一般用机理分析方法建立数学模型。如果由于客观事物内部规律的复杂性及人们认识程度的限制，无法分析实际对象内在的因果关系，建立合乎机理规律的数学模型，那么通常的办法是搜集大量数据，基于对数据的统计分析去建立模型。本章讨论其中用途非常广泛的一类模型——统计回归模型。回归模型常用来解决预测、控制、生产工艺优化等问题。变量之间的关系可以分为两类：一类叫确定性关系，也叫函数关系，其特征是：一个变量随着其它变量的确定而确定。另一类关系叫相关关系，变量之间的关系很难用一种精确的方法表示出来。例如，通常人的年龄越大血压越高，但人的年龄和血压之间没有确定的数量关系，人的年龄和血压之间的关系就是相关关系。回归分析就是处理变量之间的相关关系的一种数学方法。其解决问题的大致方法、步骤如下：（1）收集一组包含因变量和自变量的数据；（2）选定因变量和自变量之间的模型，即一个数学式子，利用数据按照最小二乘准则计算模型中的系数；（3）利用统计分析方法对不同的模型进行比较，找出与数据拟合得最好的模型；（4）判断得到的模型是否适合于这组数据；（5）利用模型对因变量作出预测或解释。应用统计分析特别是多元统计分析方法一般都要处理大量数据，工作量非常大，所以在计算机普及以前，这些方法大都是停留在理论研究上。运用一般计算语言编程也要占用大量时间，而对于经济管理及社会学等对高级编程语言了解不深的人来说要应用这些统计方法更是不可能。MATLAB 等软件的开发和普及大大减少了对计算机编程的要求，使数据分析方法的广泛应用成为可能。MATLAB 统计工具箱几乎包括了数理统计方面主要的概念、理论、方法和算法。运用MATLAB 统计工具箱，我们可以十分方便地在计算机上进行计算，从而进一步加深理解，同时，其强大的图形功能使得概念、过程和结果可以直观地展现在我们面前。本章内容通常先介绍有关回归分析的数学原理，主要说明建模过程中要做的工作及理由，如模型的假设检验、参数估计等，为了把主要精力集中在应用上，我们略去详细而繁杂的理论。在此基础上再介绍在建模过程中如何有效地使用MATLAB 软件。没有学过这部分数学知识的读者可以不深究其数学原理，只要知道回归分析的目的，按照相应方法通过软件显示的图形或计算所得结果表示什么意思，那么，仍然可以学到用回归模型解决实际问题的基本方法。包括：一元线性回归、多元线性回归、非线性回归、逐步回归等方法以及如何利用MATLAB 软件建立初步的数学模型，如何透过输出结果对模型进行分析和改进，回归模型的应用等。 8.1 一元线性回归分析回归模型可分为线性回归模型和非线性回归模型。非线性回归模型是回归函数关于未知参数具有非线性结构的回归模型。某些非线性回归模型可以化为线性回归模型处理；如果知道函数形式只是要确定其中的参数则是拟合问题，可以使用MATLAB 软件的curvefit 命令或nlinfit 命令拟合得到参数的估计并进行统计分析。本节主要考察线性回归模型。 8.1.1 一元线性回归模型的建立及其MATLAB 实现其中01ββ，是待定系数，对于不同的,x y 是相互独立的随机变量。假设对于x 的n 个值i x ，得到 y 的n 个相应的值i y ，确定01ββ，的方法是根据最小二乘准则，要使取最小值。利用极值必要条件令 01 0,0Q Q ββ??==??，求01ββ，的估计值01??ββ，，从而得到回归直线01 ??y x ββ=+。只不过这个过程可以由软件通过直线拟合完成，而无须进行繁杂的运算。

有限元模型如何查错

有限元模型如何查错作者：PAUL KUROWSKI 在建立有限元模型的过程中很容易出错，如果你知道如何查错，修正这些错误将会变得很简单有限元分析的第一步就是建立被分析对象的数学模型，这要求我们思索建模的理论基础如弹性理论，板的Reissner理论，塑性变形理论等，和考虑问题的其它信息如几何描述、材料特性，约束和荷载等等。分析的目的就是由这些条件，计算得到精确解u_EX并同时得到位移u_EX的应力函数 F(u_EX)如Von Mises应力等。应力函数F (u_EX) 仅仅依赖于数学模型的定义，而与求解该数学问题的数值近似计算方法无关；同时应力函数F(u_EX)也不依赖于网格划分、网格类型和单元尺寸。函数F(u_EX)与模型实体物理性质之间的差异，被称为“模型错误”。下一步就是使用有限元方法去找到精确解u_EX的近似值u_FE。这个过程包括选择网格划分和构件类型，如对二维板用八节点（矩形）单元，依此类推。网格划分＆单元定义被称为有限元的离散化。离散化产生的误差可以被定义为：大部分的分析应该把这个误差控制在10％以内。同时由于建立模型和模型的离散化一定会产生这个误差，正确运用有限元分析就包括对这两类误差进行评估和控制。有限元分析结果中的名义误差＆真实误差是有区别的，最好能够加以区别：名义误差可以比建模误差和离散误差的总和小，二者可能反号而相互抵消。结果的好坏取决于模型是否反应实际（模型误差的大小）和有限元软件在转化过程中的精度控制（离散

化误差的大小）。 WHAT IS MODELING ERROR? 何为模型误差？假设要分析一个支架，我们首先考虑到的问题应该包括：我们想得到什么结果？是最大应力还是最大变形？是固有频率、弯曲刚度、还是温度分布？支架是否处于弹性变形阶段？极限荷载形式有几种情况？如何模拟支撑条件等等。有了一个明确的目标和对我们使用的理论自身局限性的把握，分析者就可以建立模型了。有时这个模型与CAD模型是相似的，但相当多的情况是，为了简化网格的划分，我们有必要修改模型的拓扑描述。部分建模的过程包括以下一些问题：用壳单元模拟薄壁墙体，对对称性、反对称性或两者的运用，是否考虑细部及忽略不重要的特征等。比如，选用壳单元而不用实体单元意味着我们考虑到数学模型和相应的有限元软件的运作方式而作出了一个重要的决定。当(研究对象的)拓扑描述已经比较理想后，我们还需要对材料属性(选择线弹性、弹塑性或其他)、荷载及支撑条件进行理想的简化。我们认为这些简化精确反应了所需模型的重要数据，而建模当中的一些重要决策有时并未过多的考虑这些(方面)。简化了的模型经常是概念错误的，一个检验模型是否不合理的方法是其解析解对应的应变能是否无穷大或趋近于零；另一个方法是对应于数学模型的我们感兴趣的数据在结果没有得到体现。很多分析者认为一个有效的网格生成器可以生成高质量的网格并降低模型误差，其实不尽然，模型是在网格划分前假定的，因此，最合理的网格划分也无法修正一个简化不合理的数学模型。 A SYSTEMATIC APPROACH 一个系统的方法确保模型误差较小的唯一方式是把所需研究的数据放在对模型假设不敏感之处。类似地，通过把所需研究的数据放在对离散不敏感之处（不敏感的表现是：结果对更细的网格划分或更大的p值并不发生明显的改变），以减少离散误差。举个例子：比如说我们对一块简支板沿着边缘方向的剪力感兴趣，那么经典的克西霍夫板模型(Kirchhoff’s plate)是不可用的，可以通过一个Reissner模型或一个全3D的弹性模型轻而易举地检验出来。一个关于板弯曲的Reissner模型假设所有平面内位移沿厚度方向呈线性变化、剪应变沿厚度方向保持不变。若采用更厚的板的话会迫使人们去置疑简支的意义、同时会置疑是否可以给出一

数学建模五步法与灵敏度分析

灵敏度分析简介：研究与分析一个系统（或模型）的状态或输出变化对系统参数或周围条件变化的敏感程度的方法。在最优化方法中经常利用灵敏度分析来研究原始数据不准确或发生变化时最优解的稳定性。通过灵敏度分析还可以决定哪些参数对系统或模型有较大的影响。因此，灵敏度分析几乎在所有的运筹学方法中以及在对各种方案进行评价时都是很重要的。用途：主要用于模型检验和推广。简单来说就是改变模型原有的假设条件之后，所得到的结果会发生多大的变化。举例（建模五步法）：一头猪重200磅，每天增重5磅，饲养每天需花费45美分。猪的市场价格为每磅65美分，但每天下降1美分，求出售猪的最佳时间。建立数学模型的五个步骤： 1.提出问题 2.选择建模方法 3.推到模型的数学表达式 4.求解模型 5.回答问题第一步：提出问题将问题用数学语言表达。例子中包含以下变量：猪的重量w（磅），从现在到出售猪期间经历的时间t（天），t天内饲养猪的花费C（美元），猪的市场价格p（美元/磅），出售生猪所获得的收益R（美元），我们最终要获得的净收益P（美元）。还有一些其他量，如猪的初始重量200磅。（建议先写显而易见的部分）猪从200磅按每天5磅增加（w磅）=（200磅）+（5磅/天）*（t天）饲养每天花费45美分（C美元）=（0.45美元/天）*（t天）价格65美分按每天1美分下降（p美元/磅）=（0.65美元/磅）-（0.01美元/磅）*（t天）生猪收益（R美元）=（p美元/磅）*（w磅）净利润（P美元）=（R美元）-（C美元）用数学语言总结和表达如下：参数设定： t=时间（天）

w=猪的重量（磅） p=猪的价格（美元/磅） C=饲养t天的花费（美元） R=出售猪的收益（美元） P=净收益（美元）假设： w=200+5t C=0.45t p=0.65-0.01t R=p*w P=R-C t>=0 目标：求P的最大值第二步：选择建模方法本例采用单变量最优化问题或极大—极小化问题第三步：推导模型的数学表达式子 P=R-C （1） R=p*w （2） C=0.45t （3）得到R=p*w-0.45t p=0.65-0.01t （4） w=200+5t （5）得到P=（0.65-0.01t)(200+5t)-0.45t 令y=P是需最大化的目标变量，x=t是自变量，现在我们将问题转化为集合S={x:x>=0}上求函数的最大值： y=f（x）=（0.65-0.01x)(200+5x)-0.45x （1-1）第四步：求解模型用第二步中确定的数学方法解出步骤三。例子中，要求（1-1）式中定义的y=f （x）在区间x>=0上求最大值。下图给出了（1-1）的图像和导数（应用几何画板绘制）。在x=8为全局极大值点，此时f（8）=133.20。因此（8，133.20）为f在整个实轴上的全局极大值点，同时也是区间x>=0上的最大值点。第五步：回答问题根据第四步，8天后出售生猪的净收益最大，可以获得净收益133.20美元。只要第一步中的假设成立，这一结果正确。

结构有限元及其应用软件

中国海洋大学本科生课程大纲课程属性：公共基础/通识教育/学科基础/专业知识/工作技能，课程性质：必修、选修一、课程介绍 1.课程描述（中英文）：本课程是一门重要的结构计算分析课程，通过多媒体教学和上机练习，系统学习结构有限元FEM的基本原理和方法，熟悉掌握通用有限元应用软件ANSYS进行结构静力和动力分析的方法和步骤，并初步掌握使用ANSYS进行海工典型结构强度计算的方法。 Structural finite element method and its application software is an important course of structural calculation and analysis. Through multimedia teaching and computer practice, the basic principles and methods of Finite Element Method (FEM) are learned systematically. The general finite element application software ANSYS for the methods and procedures of structural static and dynamic analysis are mastered.At the same time, the strength calculation method of typical ocean engineering structures using ANSYS is preliminarily mastered. 2.设计思路：有限元方法是一种现代设计方法，应用于结构设计中，是一种具有重要经济意义和巨大潜力的先进结构设计技术。因此选择该课程作为结构设计方面的一门必修课程，主要介绍结构有限元的基本原理和方法，还选择了通用的有限元软件ANSYS进行示例分析。包括要求掌握有限元法的基本思想和基本原理、平面刚架结构的有限元法、弹

机械结构有限元分析

机械结构有限元分析有限元分析软件ANSYS在机械设计中的应用摘要：在机械设计中运用ANSYS软件进行有限元分析是今后机械设计发展的必然趋势，将有限元方法引入到机械设计课程教学中，让学生参与如何用有限元法来求解一些典型零件的应力，并将有限元结果与教材上的理论结果进行对照。这种新的教学方法可以大大提高学生的学习兴趣，增强学生对专业知识的理解和掌握，同时还可以培养学生的动手能力。在机械设计课程教学中具有很强的实用价值。关键词：机械设计有限元 Ansys 前言：机械设计课程是一门专业基础课，其中很多教学内容都涉及到如何求取零件的应力问题，比如齿轮、v带、螺栓等零件。在传统的教学过程中，都是根据零件的具体受力情况按材料力学中相应的计算公式来求解。比如，在求解齿轮的接触应力时，是把齿轮啮合转化为两圆柱体的接触，再用公式求解。这些公式本身就比较复杂，还要引入各种修正参数，因此我们在学习这些内容时普遍反映公式难记，学习起来枯燥乏味，而且很吃力。近年来有限元法在结构分析中应用越来越广泛，因此如果能将这种方法运用到机械设计课程中，求解一些典型零件的应力应变，并将分析结果和教材上的理论结果进行对比，那么无论是对于提高学生学习的热情和积极性，增强对重点、难点知识的理解程度，还是加强学生的计算机水平都是一件非常有益的事情。由于直齿圆柱齿轮的接触强度计算是机械设计课程中的一个重要内容，齿轮强度的计算也是课程中工作量最繁琐的部分。下面就以渐开线直齿圆柱齿轮的齿根弯曲疲劳强度的计算为例，探讨在机械设计课程中用ANSYS软件进行计算机辅助教学的步骤和方法，简述如何将有限元方法应用到这门课程的教学中。 1．传统的直齿圆柱齿轮齿根弯曲疲劳强度的计算传统方法把轮齿看作宽度为b的矩形截面的悬臂梁。因此齿根处为危险剖面，它可用30。切线法确定。如图l所示。作与轮齿对称中心线成30。角并与齿根过渡曲线相切的切线，通过两切点作平行与齿轮轴线的剖面，即齿根危险剖面。理论上载荷应由同时啮合的多对齿分担，但为简化计算，通常假设全部载荷作用于齿顶来进行分析，另用重合度系数E对齿根弯曲应力予以修正。由材料力学弯曲应力计算方法求得齿根最大弯曲应力为：

数学建模对偶问题和灵敏度分析资料讲解

数学建模对偶问题和灵敏度分析

对偶问题例题1：某养鸡场所用的混合饲料由n 种天然饲料配合而成。要求在这批配合饲料中必须含有m 种不同的营养成分，且第i 种营养成分的含量不低于bi 。已知第i 种营养成分在每单位第j 种天然饲料中的含量为a ij ，每单位第j 天然饲料的价格为c j 。试问，应如何对这n 种饲料配方，使这批饲料的费用最小? 解设x j 为第j 种天然饲料的用量。显然，a ij x j 即为所用第j 种天然饲料中第i 种营养成分的含量，1n ij j j a x =∑为这批混合饲料中第i 种营养成分的总含量；它不应低于bi 。于是，我们得下列线性规划模型(1—1)： 1 min n j j j f c x ==∑ 1 1,,..01,,n ij j i j j a x b i m s t x j n =?≥=???≥=? ∑ 现设想有一个饲料加工厂欲把这m 种营养成分分别制成m 种营养丸。设第i 种营养丸的价格为ui(i ＝1，…，m)。则养鸡场采购一个单位的第j 种天然饲料，就相当于对这m 种营养丸分别采购数量a 1j ，…a mj ，所化费用为1m ij i i a u =∑养鸡场自然希望在用营养丸代替天然饲料时，在价格上能相对地比较便宜，故而饲料加工厂为了能与天然饲料供应者竞争，在制订价格时必然满足下述条件： 1 1, ,m ij i j i a u c j n =≤=∑ 另一方面，养鸡场如果全部采购营养丸来代替天然饲料进行配料，则第i 种营养丸就需采购bi 个单位，所化费用为b i u i ，总费用为z=∑b i u i

回归分析在数学建模中的应用

摘要回归分析和方差分析是探究和处理相关关系的两个重要的分支,其中回归分析方法是预测方面最常用的数学方法,它是利用统计数据来确定变量之间的关系,并且依据这种关系来预测未来的发展趋势。本文主要介绍了一元线性回归分析方法和多元线性回归分析方法的一般思想方法和一般步骤,并且用它们来研究和分析我们在生活中常遇到的一些难以用函数形式确定的变量之间的关系。在解决的过程中,建立回归方程,再通过该回归方程进行预测。关键词:多元线性回归分析;参数估计;F检验

回归分析在数学建模中的应用 Abstract Regression analysis and analysis of variance is the inquiry and processing of the correlation between two important branches, wherein the regression analysis method is the most commonly used mathematical prediction method, it is the use of statistical data to determine the relationship between the variables, and based on this relationship predict future trends. introduces a linear regression analysis and multiple linear regression analysis method general way of thinking and the general steps, and use them to research and analysis that we encounter in our life, are difficult to determine as a function relationship between the variables in the solving process, the regression equation is established by the regression equation to predict. Keywords:Multiple linear regression analysis; parameter estimation;inspection II

结构有限元分析的形状处理方法_杜平安

结构有限元分析的形状处理方法杜平安摘要　介绍结构形状处理的各种方法,包括类型简化、细节简化、形式变换、局部结构和利用对称性等。关键词　形状处理　有限元分析　建模 Abstract The processing method is intro-duced in the paper ,including ty pe simplifica tion 、details simplifica tio n 、fo rm tra nsfo rmatio n 、local structure a nd symm etry utiliza tion . Key words Shape processing Finite element analysis Modelling 收稿日期:1999-08-18 1　结构类型简化根据结构形状、载荷和约束条件的特点,结构类型可分为空间问题、平面问题、轴对称问题、板壳问题和杆件问题等。其中平面问题和轴对称问题的几何模型是一平面图形,在平面上划分网格比在空间内划分要容易得多,单元数量也少得多。因此将空间问题作适当近似,使其按平面问题来处理,则可使分析过程大为简化。在图1a 中,计算轮毂与轴过盈配合的接触压力时,由于辐孔尺寸较小且远离接触面,因此可以不考虑辐孔而将轮毂简化为轴对称结构。同样,在计算图1b 中螺栓与螺母螺纹面上的接触压力时,由于螺旋升角较小,也可以不考虑升角的影响,而将螺栓与螺母简化为轴对称结构。图1　结构类型简化结构 2　结构细节简化细节是结构中相对尺寸很小的局部,如倒圆、倒角、退刀槽和加工凸台等。根据网格划分特点,一条直线或曲线至少要划分一个单元边;一个平面或曲面至少要划分一个单元面;一个圆至少要用三个单元边离散,因此几何模型中的细节将限制细节处及其附近的网格大小,从而影响整个结构的网格分布和增加网格数量。图2是有无细节时自动划分出的网格,从中可以看出细节对网格划分的影响。图2　细节对网格划分的影响因此,建立几何模型时应尽量忽略一些不必要的细节。在静力分析中,高应力区域中的细节会引起应力集中,细节大小和形状对应力影响很大,这些细节不能忽略。而处于结构低应力区的细节一般可以忽略。在动力计算中,由于结构固有频率和模态振型主要取决于结构的质量分布和刚度,因此细节一般可以忽略。在热分析中,细节不会在结构中引起局部高温,这时也可以考虑较少的细节。 3　结构形式变换有些结构尽管形状不是很复杂,但划分网格却很困难。如果对结构形式作适当变换,则可使网格划分变得容易,划分出的单元更少。例如图3a 所示的带肋板,划分网格时需要用板单元和梁单元组合,且两类单元为偏心连接,自动分网难以满足这种要求。如果将带肋板变换为平板(图3b 所示),则在平板上划分网格要容易得多。由于带肋板用于焊接而成支撑箱式立柱,其特性要求主要是刚度,因此可按等刚度条件作为变换 · 26·《机械与电子》2000(1)

有限元+螺栓简化

1 概述螺栓是机载设备设计中常用的联接件之一。其具有结构简单,拆装方便,调整容易等优点,被广泛应用于航空、航天、汽车以及各种工程结构之中。在航空机载环境下，由于振动冲击的影响，设备往往产生较大的过载，对作为紧固件的螺栓带来强度高要求。螺栓是否满足强度要求，关系到机载设备的稳定性和安全性。传统力学的解析方法对螺栓进行强度校核，主要是运用力的分解和平移原理，解力学平衡方程，借助理论和经验公式，理想化和公式化。没有考虑到连接部件整体性、力的传递途径、部件的局部细节(如应力集中、应力分布)等等。通过有限元法，整体建模，局部细化，可以弥补传统力学解析的缺陷。用有限元分析软件MSC.Patran/MSC.Nastran提供的特殊单元来模拟螺栓连接，过程更方便，计算更精确，结果更可靠。因此，有限元在螺栓强度校核中的应用越来越广泛。 2 有限元模型的建立对于螺栓的模拟，有多种模拟方法，如多点约束单元法和梁元法等。多点约束单元法(MPC)即采用特殊单元RBE2来模拟螺栓连接。在螺栓连接处，设置其中一节点为从节点(Dependent)，另外一个节点为主节点(Independent)。主从节点之间位移约束关系使得从节点跟随主节点位移变化。比例因子选为1,使从节点和主节点位移变化协调一致，从而模拟实际工作状态下，螺栓对法兰的连接紧固作用。梁元法模拟即采用两节点梁单元Beam，其能承受拉伸、剪切、扭转。通过参数设置，使梁元与螺栓几何属性一致。本文分别用算例来说明这两种方法的可行性。 2.1 几何模型如图1所示组合装配体，底部约束。两圆筒连接法兰通过8颗螺栓固定。端面受联合载荷作用。