1、2014 高教社杯全国大学生数学建模竞赛高教社杯全国大学生数学建模竞赛 承承 诺诺 书书 我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参 赛规则（以下简称为 “竞赛章程和参赛规则” ， 可从全国大学生数学建模竞赛网站下载） 。 我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网 上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。 我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或 其他公开的资料（包括网上查到的资料） ，必须按照规定的参考文献的表述方式在正文 引用处和参考文献中明确列出。 我们郑重承诺，严格遵守竞
2、赛章程和参赛规则，以保证竞赛的公正、公平性。如有 违反竞赛章程和参赛规则的行为，我们将受到严肃处理。 我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展 示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等） 。 我们参赛选择的题号是（从 A/B/C/D 中选择一项填写） ： B 我们的报名参赛队号为（8 位数字组成的编号） ： 所属学校（请填写完整的全名） ： 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： （论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容 请仔细核对，提交后将不再允许做任
3、何修改。如填写错误，论文可能被取消评奖资格。） 日期： 2014 年 9 月 15 日 赛区评阅编号（由赛区组委会评阅前进行编号）： 2014 高教社杯全国大学生数学建模竞赛高教社杯全国大学生数学建模竞赛 编编 号号 专专 用用 页页 赛区评阅编号（由赛区组委会评阅前进行编号）： 赛区评阅记录（可供赛区评阅时使用）： 评 阅 人 评 分 备 注 全国统一编号（由赛区组委会送交全国前编号）： 全国评阅编号（由全国组委会评阅前进行编号）： 1 创意平板折叠桌 摘要 目前住宅空间的紧张导致越来越多的折叠家具的出现。某公司设计制作了一款折叠 桌以满足市场需要。以此折叠桌为背景提出了三个问题，本文运用几
4、何知识、非线性约 束优化模型等方法成功解决了这三个问题， 得到了折叠桌动态过程的描述方程以及在给 定条件下怎样选择最优设计加工参数， 并针对任意形状的桌面边缘线等给出了我们的设 计。 针对问题一，根据木板尺寸、木条宽度，首先确定木条根数为 19 根，接着，根据桌 子是前后左右对称的结构，我们只以桌子的四分之一为研究对象，运用空间几何的相关 知识关系，推导并建立了几何模型。接着用 MATLAB 软件编程，绘制出折叠桌动态变化 过程图。然后求出折叠桌各木条相对桌面的角度、各木条长度、各木条的开槽长度等数 据，相关结果见表 1。然后建立相应的三维坐标系，求出桌角各端点坐标，绘出桌角边 缘线曲线图，并
5、用 MATLAB 工具箱作拟合，求出桌角边缘线的函数关系式，并对拟合效 果做分析（见表 3） 。 针对问题二，在折叠桌高度、桌面直径已知情况下，综合考虑桌子稳固性、加工方 便、用材最少三个方面因素，我们运用材料力学等相关知识，对折叠桌作受力分析，确 定稳固性、 加工方便、 用材最少三个方面因素间的相互制约关系， 建立非线性优化模型。 用 lingo 软件编程，求出对于高 70 cm，桌面直径 80 cm 的折叠桌，平板尺寸 、钢筋位置在桌腿上距离铰链 46.13cm 处、各木条的开槽长 度（见表 3） 、最长木条（桌脚）与水平面夹角。 针对问题三，对任意给出的桌面边缘线（f(x)，不妨假定曲线
6、是对称的（否则，桌 子的稳定性难以保证） ，将对称轴上 n 等份，依照等份点沿着木板较长方向平行的方向 下料，则这些点即是铰接处到木板中垂线（相对于木板长方向）的距离。然后修改问题 二建立的优化模型，用 lingo 软件编程，得到最优设计加工参数（平板尺寸、钢筋位置、 开槽长度等） 。最后，我们根据所建立的模型，设计了一个桌面边缘线为椭圆的折叠桌， 并且给出了 8 个动态变化过程图（见图 10）和其具体设计加工参数（见表 5） 。 最后，对所建立的模型和求解方法的优缺点给出了客观的评价，并指出了改进的方 法。 关键字：折叠桌 曲线拟合 非线性优化模型 受力分析 2 一、 问题重述 1.1 引言
7、 创意平板折叠桌注重于表达木制品的优雅和设计师所想要强调的自动化与 功能性。为了增大有效使用面积。设计师以长方形木板的宽为直径截取了一个圆 形作为桌面，又将木板剩余的面积切割成了若干个长短不一的木条，每根木条的 长度为平板宽到圆上一点的距离，分别用两根钢筋贯穿两侧的木条，使用者只需 提起木板的两侧，便可以在重力的作用下达到自动升起的效果，相互对称的木条 宛如下垂的桌布，精密的制作工艺配以质朴的木材，让这件工艺品看起来就像是 工业革命时期的机器。 1.2 问题的提出 围绕创意平板折叠桌的动态变化过程、设计加工参数，本文依次提出如下问 题： （1） 给定长方形平板尺寸 （120 cm 50 cm
8、3 cm） ， 每根木条宽度 （2.5 cm） ， 连接桌腿木条的钢筋的位置，折叠后桌子的高度（53 cm） 。要求建立模型描述此 折叠桌的动态变化过程， 并在此基础上给出此折叠桌的设计加工参数和桌脚边缘 线的数学描述。 （2）折叠桌的设计应做到产品稳固性好、加工方便、用材最少。对于任意给定 的折叠桌高度和圆形桌面直径的设计要求， 讨论长方形平板材料和折叠桌的最优 设计加工参数，例如，平板尺寸、钢筋位置、开槽长度等。对于桌高 70 cm，桌 面直径 80 cm 的情形，确定最优设计加工参数。 （3）给出软件设计的数学模型，可以根据客户任意设定的折叠桌高度、桌面边 缘线的形状大小和桌脚边缘线的大
9、致形状， 给出所需平板材料的形状尺寸和切实 可行的最优设计加工参数，使得生产的折叠桌尽可能接近客户所期望的形状，并 根据所建立的模型给出几个设计的创意平板折叠桌。 要求给出相应的设计加工参 数，画出至少 8 张动态变化过程的示意图。 一、 模型假设 (1) 忽略实际加工误差对设计的影响； (2) 木条与圆桌面之间的交接处缝隙较小，可忽略； (3) 钢筋强度足够大，不弯曲； (4) 假设地面平整。 3 三、符号说明 符号符号 意义意义 D 木条宽度（cm） 缝宽 L 木板长度（cm） W 木板宽度（cm） N 第n根木条 T 木条根数 木板从外起第1个木条的长度(cm) 木板从外起第n个木条的长
10、度(cm) H 桌子高度（cm） R 桌子半径（cm） R 桌子直径（cm） 桌子厚度(cm) 第n根木条到木板边沿的距离(cm) 第n根木条顶点位置到圆面轴线径向距离(cm) 第n根木条与水平面的夹角(度) 第n根木条开槽长度（cm） 四、问题分析 4.1 问题一分析 题目要求建立模型描述折叠桌的动态变化图，由于在折叠时用力大小的不 同， 我们不能描述在某一时刻折叠桌的具体形态，但我们可以用每根木条的角度 4 变化来描述折叠桌的动态变化。首先，我们知道折叠桌前后左右对称，我们可以 运用几何知识求出四分之一木条的角度变化。最后，根据初始时刻和最终形态两 种状态求出桌腿木条开槽的长度。 4.2
11、问题二分析 题目要求从折叠桌的稳固性好、加工方便、用材最少三个角度，确定设计加 工参数。我们可以从应力、支撑面积考虑稳固性，从开槽长度考虑加工方便，从 木板长度考虑用材最少。而它们之间又是相互制约，我们需要确定最优设计加工 参数，可以建立非线性规划模型，用 lingo 软件来求解最优设计加工参数（平板 尺寸、钢筋位置、开槽长度等） ，这里以合力的方向（斜向上）与最长木条（桌 腿）的夹角方向最小为目标函数，以木条所承受应力小于木条的许用应力、支撑 面积大于桌面面积、木条的开槽长度小于木条本身长为约束条件。 4.3 问题三分析 题目要求制作软件的意思就是客户给定折叠桌高度、桌面边缘线的形状大小 和
12、桌脚边缘线的大致形状，将这些信息输入程序就得到客户想要的桌子。我们在 求解最优设计加工参数时，自行给定桌面边缘线形状（椭圆、相交圆等），桌脚 边缘线形状，折叠桌高度，应用第二问的非线性规划模型，用 MATLAB 软件绘制 折叠桌截面图，得到自己设计的创意平板折叠桌。 问题三流程图： 五、模型建立和解决 已知 f(x)、g(x)、h、w 5 5.1 问题一的模型建立和解决 5.1.1 模型的准备 (1)符号说明 为求出各木条角度关系，现引入下列符号： ：木板从外起第n个木条的长度(cm) ：第n个木条到木板边沿的距离 ：第n个木条与桌面铰接处到桌面轴线距离 ：第n个木条与第n-1个木条桌面铰接处
13、到桌面轴线距离差 ：第n个木条与桌面的夹角 （2）木条数的确定 根据题目意思，长方形平板尺寸，宽50 cm，每根木条宽2.5 cm，知道木条 数越多，桌子越不易松动，即稳固性更好，最大根数为根，考虑木条间 的间隙和刀片的厚度，定为19根，此时，缝宽为： （3）模型近似 从折叠桌实物可以看出，桌面并非为标准的圆面，圆面边上是锯齿形状，考 虑到锯齿长度和圆半径的差异，我们假定圆为过木条中点的圆，在作示意简图和 实际计算时，都以木条端点中点为木条与桌面接触点。 另外，折叠桌以材料最省为设计原则，在木板尺寸一定情况下，应该做到桌 面尽可能大，这里我们取木板宽度为桌面直径。 5.1.2 模型的建立 为帮
14、助理解，我们做折叠桌子两个最长脚（即在未折叠时的木板的同一侧最 长木条）示意图，如图1所示： h-3 A 点 D 点 6 图图 1 1 折叠桌子两个最长脚截面图折叠桌子两个最长脚截面图 (其中A点为最长木条一端到水平面的距离，由于桌实际高度包括桌面厚度3cm， 则A点到水平面距离要减去3cm) 其中 为57cm，因为木板厚度为3cm,有AD为两倍厚度，因为 则知为57cm。记 下面，我们作出平板俯视示意图，如下图2所示 图图2 2 平板俯视示意图平板俯视示意图 对于第n个木条到木板边沿的距离，应该包括（n-1）条缝宽，（n-1）根木条 长度以及它自身一半的长度，则有： 从几何关系上，应用勾股定
15、理可以得出： B 点C 点 E 点 第n根木条 第n-1根木条 0.5（1） 0.5（2） 7 则第n个木条与第n-1个木条顶点位置到圆面轴线径向距离差： 第n根木条长度： 为了求解木条旋转角度，我们沿着钢筋的角度，作出折叠凳示意简图，如图3 所示： 图图3 3 折叠桌示意简图折叠桌示意简图 由上图知 同理可得递推公式，即每根木条旋转角度： 0.50.5h 0.5（3） 0.5（4） 8 （由图3知，可能为负值，说明） 开槽长度 综合以上所分析，可建立如下几何模型： 5.1.3 模型的解决 （1）动态变化过程 动态变化过程：由于用力大小未知，折叠桌与时间的关系不能确定，我们只 能确定桌子从平板
16、到折叠完成后这一过程中，任一角度的桌角位置， （程序见附 录 problem1_3.m）例如当最长木条转过、，通过程序可以得到 各木条相对桌面旋转角度，如表 1 所示： 表表1 1最长木条转过最长木条转过、时各木条转动角度时各木条转动角度 夹角为60夹角为65夹角为70 第1根606570 第2根71.510676.821982.0272 第3根79.72884.982890.063 第4根85.97791.041495.8979 第5根90.765395.6054100.2279 第6根94.383599.0138103.1289 第7根97.0267101.484105.7333 第8根9
17、2.8285103.1591107.2893 第9根99.8766104.1306108.1893 0.5（6） 0.5（7） （5） 9 （2）长槽长度、木条长度、旋转角度 根据以上建立的模型，运用 MATLAB 软件，编程计算每根木条长度、旋转角 度、长槽长度结果如下表 2 所示： 表表2 2 木条长度、旋转角度、长槽长度木条长度、旋转角度、长槽长度 第1根第2根第3根第4根第5根第6根第7根第8根第9根第10根 111.1111.38旋转角度 73.719 85.833 93.737 99.39 103.54 38.765 37.338 36.287 35.563 106.59 108.
18、78 110.25 35.1435 14.793 16.164 17.128 17.70217.892 木条长度 52.089 46.609 43.154 40.65 12.994卡槽长度04.5018 7.9434 10.73 从表 1 可以看出，第一根木条卡槽长度为 0cm,符合实际。 下面我们绘制木条长度（如图4所示），开槽长度（如图5所示）： 02468101214161820 0 10 20 30 40 50 60 序 号 木条长度（单位：cm） 02468101214161820 0 10 20 30 40 50 60 序 号 开槽长度（单位：cm） 图图4 4 木条长度图木条长度
19、图 图图5 5 开槽长度柱形图开槽长度柱形图 （3）桌脚边缘线的描述 为形象描述桌脚边缘线，可以用MATLAB绘图，因此，首先建立三维坐标系， 我们以一个桌角为坐标原点，两桌角（平板状态时为异侧木条）连线为x轴，另 两桌角（平板状态时为同侧木条）连线为y轴，竖直方向为z轴，如图6所示 10 图图6 6 坐标示意图坐标示意图 a. x坐标的确定 考虑到编程的需要，这里直接以数组的形式表示木条桌脚x坐标，记为： 因为每根木条的长度都垂直于x轴（如坐标中红线所示） ，可以得到： b. y坐标的确定 将桌子投影到xoy平面，根据几何关系可以得到木条桌脚y坐标 z y x O 0.5（8） 11 其中：
20、 c. z坐标的确定 将桌子投影到zoy平面，根据几何关系可以得到木条桌脚z坐标 其中： 综合以上分析，运用MATLAB编程（详见problem1_3.m） ，绘制桌角边缘线如图7： 0 10 20 30 40 50 0 5 10 15 0 5 10 15 20 x y z 图图7 7 桌角边缘线桌角边缘线 0.5（9） 0.5（10） 12 为了更精确的描述桌角边缘线，我们可以调用MATLAB拟合工具箱，用多项式 拟合得到桌角边缘线函数和拟合图形，如图8所示： 图图8 8 桌角边缘线函数和拟合图形桌角边缘线函数和拟合图形 拟合函数： Linear model Poly33: f(x,y) =
21、 p00 + p10*x + p01*y + p20*x2 + p11*x*y + p02*y2 + p30*x3 + p21*x2*y + p12*x*y2 + p03*y3 Coefficients (with 95% confidence bounds): p00 = -5.983e-007 (-0.0002297, 0.0002285) p10 = 1.153 (1.15, 1.155) p01 = 0.0173 (0.01688, 0.01771) p20 = -0.02427 (-0.02431, -0.02422) p11 = 0.02388 (0.02366, 0.02411)
22、 p02 = 0.006276 (0.006155, 0.006397) p30 = -4.035e-018 (-1.075e-007, 1.075e-007) p21 = -0.0005028 (-0.0005075, -0.0004981) p12 = -4.395e-017 (-6.067e-007, 6.067e-007) p03 = -3.146e-005 (-4.434e-005, -1.859e-005) 我们还可以得到拟合效果的分析，如表3所示： 表表3 3 拟合效果的分析拟合效果的分析 拟合类 型 误差平方和 复相关系 数 自由度 协方差 均方根误 差 三 维 拟 合 1.0
23、00 9 0 当误差平方和和均方根误差越小，复相关系数越接近于1时标明拟合的越好。 由上表可知，误差平方和为，均方根误差为0，都很小，复相关系数为1， 说明拟合效果很好。 5.2 问题二的模型建立和解决 5.2.1模型准备 13 (1)符号说明 dd:木条厚度 :木条宽度 :木条根数 :木条间的缝隙 s：钢筋位置到桌面圆心的径向距离 H:钢筋位置到桌面的径向距离 :木板长度 :支撑面积 （2）参数确定 木条根数(取整)： 缝隙: 按照问题一同样的处理方法，我们可以得到： 第n个木条与桌面铰接处到桌面轴线距离： 则第n个木条与第n-1个第n个木条与桌面铰接处到桌面轴线距离： 再由几何关系可以得到
24、第一根木条与水平方向夹角： 0.5（11） 0.5（12） 14 钢筋位置到桌面的径向距离H: 每根木条旋转角度: 接着，和第一问相同的处理方法，我们可以得出开槽长度： 5.2.2模型的建立 （1）目标函数： 我们知道，钢条对每根木条都有作用力，当桌子上有物品时，该作用力表现为支 持力，方向朝上，为帮助理解，我们作出桌子受力示意图如图9所示： 图图9 9 桌子受力示意图桌子受力示意图 0.5（13） 0.5（14） 0.5（15） 0.5（16） 15 （粗体表示矢量） 而可由每根木条受到的钢筋对它的作用力的分解再加和得到： (其中为钢筋对木条的作用力，我们知道该作用力大小相等，这里用表示)
25、于是，我们可以得到合力的方向与竖直方向的夹角： 为了使桌子稳固，合力的方向与桌腿方向（斜向上）应该尽量靠近，也就是它们 之间的夹角，我们可以以此为目标函数： （2）约束条件： 应力约束 剪应力是指物体由于外因（受力、湿度变化等）而变形时，在物体内各部分 之间产生相互作用的内力，以抵抗这种外因的作用，并力图使物体从形变后的位 置回复到形变前的位置。 抗拉强度即表征材料最大均匀塑性变形的抗力，拉伸试样在承受最大拉应力 之前，变形是均匀一致的，但超出之后，材料开始出现缩颈现象，即产生集中变 形，对于没有（或很小）均匀塑性变形的脆性材料，它反映了材料的断裂抗力。 在这里我们知道受力最大的是四个桌角，如
26、果四个桌角能够承受剪应力和抗 拉强度，则可以说明桌子是稳定的。 根据定义，我们可以得到剪应力和抗拉强度的数学表达式。 剪应力： （17） （18） 16 抗拉强度： 其中为木条宽度，dd 为木条厚度。在要判定零件或构件受载后的工作应力过 高或过低，需要预先确定一个衡量的标准，这个标准就是许用应力。所以应该有 剪应力小于许用剪应力，抗拉强度小于许用抗拉强度。即： 其中为许用剪应力，为许用抗拉强度。 支撑面积 桌面支撑面积指桌子四条腿所在点按直线连接的形成几何的面积（这里的支 撑面积为矩形） ，支撑面积越大，桌子稳固性越好。在这里我们以桌面支撑面积 大于桌面面积为约束条件。 支撑面积: 应该满足：
27、 即： 长度限制 结合实际情况，第 n 根木条的开槽长度不可能比木条本身长，则有： 即： （19） （20） （21） （22） （23） 17 综合以上分析，可以得到问题二的动态规划模型如下： 5.2.3模型的解决 运用 lingo 软件，编写程序（详见 problem2） ，我们可以得到各木条与桌面 夹角及各木条开槽长度如下表 4 所示： 表表4 4 各木条与桌面夹角及各木条开槽长度各木条与桌面夹角及各木条开槽长度 开槽长度（cm）相对桌面角度（rad） 第1根01.255484 第2根4.0185631.370394 第3根9.8730591.520353 第4根14.506131.62
28、5663 第5根18.272271.703454 第6根21.304161.761349 第7根23.67041.804035 第8根25.412211.833991 第9根26.555571.853039 第10根27.116051.862204 由表 4 知，第 1 根木条即桌脚与水平面夹角为 1.255484rad,对应为 5.3 问题三的模型建立和解决 考虑实际情况，桌子堆放物品时，桌面各点承受力相同，所以桌子应该是前 后、左右对称,这里和问题一、问题二一样，我们仅研究四分之一桌子即可。 对于客户给定的桌面边缘线函数 y=f(x),应该满足 f(x)是关于 y 轴对称的。 (1)参数的
29、确定 （24） 18 木条根数 t： （取整） （其中 d 为木条宽度） 第n个木条到木板边沿的距离： 第 n 根木条到桌面轴线的距离： 从上图知，第 n 根木条到桌面轴线的距离即为木条（上图黑线）与椭圆交点。 在这里，设计加工参数我们依然按照第二问求最优的，即产品稳固性好、加工方 便、用材最少，因此，借用问题二模型，我们可以求出此条件下的折叠桌平板尺 寸、钢筋位置、开槽长度、桌角角度。 现根据我们所建立的模型给出几个自己设计的创意平板折叠桌。 （1）椭圆桌 我们只需要研究四分之一桌面，这里取椭圆第一象限部分： 我们把函数代入 MATLAB 程序，可以得到动态变化过程的示意图如下图 10： x
30、 y F(x) 19 图图1010 椭圆形折叠桌椭圆形折叠桌8 8张动态图张动态图 20 我们利用和问题一类似的解决方法，可以得出最长木条（桌角）相对桌面不 同角度时，其余木条的旋转角度，如表 1。 接着，我们运行 lingo（problem3）程序（见附录） ，求最优设计参数。从运 行结果我们可以得到椭圆桌各木条与桌面夹角及各木条开槽长度，如表 5： 表表5 5 椭圆桌各木条与桌面夹角及各木条开槽长度椭圆桌各木条与桌面夹角及各木条开槽长度 开槽长度（cm） 相对桌面角度（rad） 第1根01.214031 第2根3.191291.424758 第3根5.4305221.553027 第4根7
31、.0406861.636302 第5根8.1022121.687448 第6根8.6483811.712669 第7根8.6930481.7147 第8根8.2372261.693751 第9根7.2695771.647554 以及木板 六、模型评价及改进 6.1 模型评价 6.1.1 问题一模型评价 问题一建立的是几何模型，运用了勾股定理的数学知识，用简单的方法去解 决复杂的问题，简明易懂。本模型主要解决了木条长度、开槽长度、木条旋转角 度，及折叠桌的动态变化过程。该模型基于严密的数学推导，求解过程严瑾，结 果可信度高，说服力强。 6.1.1 问题二模型评价 问题二模型是基于问题一的，进一步
32、的理论推导，从材料力学角度，围绕稳 固性、用材最少、加工方便，对木条作受力分析，解决了在桌子高度、半径确定 情况下，桌脚与水平面的夹角，以及开槽长度。模型理论严谨，假设大胆合理。 6.1.3 问题三模型评价 21 问题三模型综合了问题一的几何模型和问题二的非线性规划模型，巧妙使用 用 MATLAB 软件、lingo 软件编写程序，可以根据客户提供的信息，设计他们自 己的创意平板折叠桌。 6.2 模型改进 由于题目信息量不足，对于不同材料，木条截面承受轴力 N、弯矩 M 共同作 用的强度是不同的，如果是在实际生活中，这些量是已知的，所以对于问题二的 应力约束条件，现在做已下改进，使之更完善。 根
33、据材料力学的假定8，单一均质材料矩形截面杆构件截面承受轴力 N、弯 矩 M 共同作用的强度条件如下: 式中为木条征截面最大正应力;N 为木条正截面所受轴力;A 为木条正截面 面积;M 为木条正截面所受弯矩;为木条抗弯截面系数;为材料许 用应力(假设许用拉应力与许用压应力相等)。 考虑右边的不等式部分，将不等式两边都同时除以，可得: 由于 式中:为构件正截面轴力承载力;为构件正截面抗弯承载 力。 （25） （26） （27） （28） 22 故式(26)可写为: 而构件截面同时承受轴力 N、弯矩 M,剪力 v 时，截面上的正应力与剪应力计算如 下: 式中: 为构件正截面剪应力;v 为构件正截面所
34、受剪力。 求得最大主应力与最小主应力为: 根据第三强度理论： 可得： 将式(30)、式(31）入式(34)，可得: 不等式左边第一个括号内两项分子分母同时乘以，第二个括号内一项分子分 母同时乘以，可得: （29） （30） （31） （32） （33） （34） （35） （36） 23 式中，为构件斜截面抗剪承载力，。不等式两边同时除以， 得: 由式(29)、式(37)可见，单一均质材料矩形截面杆件截面承受轴力、弯矩、 剪力等复合作用时，截面的强度条件可用下式夫示: 对于拉弯剪构件,易损性系数为 ,其值越大则构件越容易受 损。 七、参考文献 （37） （38） 24 1 汪晓银， 周保平，
35、数学建模与数学实验 （第二版） ， 北京： 科学出版式， 2012.8 2 汪晓银，邹庭荣，周保平，数学软件与数学实验（第二版） ，北京：科学出 版式，2012.8 3 姜启源，叶其孝，数学建模，北京：机械工业出版社，2009.8。 4 同济大学数学系， 高等数学 （第二版） 上册， 上海： 同济大学出版社， 2009.10 5 薛定宇， 陈阳泉， 高等数学问题的 MATLAB 求解， 北京： 清华大学出版社,2008 6 刘鸿文，材料力学（第 5 版） ，北京：高等教育出版社，2010.6 7 胡运康，景荣春，理论力学，北京：高等教育出版社，2006.5 8 黄靓，王鉴，陈永亮，李登, 一种简化的结构鲁棒性量化方法,工程力学， 第 30 卷第 10 期，文章编号:1000-4750(2013)10-0046-08，2013.10 源程序引索源程序引索 问题一源程序问题一源程序 MATLAB 程序： problem1_1.m problem1_2.m problem1_3.m 问题二源程序问题二源程序 Lingo 程序： problem2.lg4 问题三源程序问题三源程序 MATLAB 程序： Problem3_1 Problem3_2 Lingo 程序： Problem3.lg4 注：注：m 文件是文件是 Matlab 程序，程序，lg4 文件是文件是 Lingo 程序程序
