创意平板折叠桌问题的模型设计与优化技术.doc

创意平板折叠桌问题的模型设计与优化一、摘要本文在充分考虑实际设计需求的基础上，讨论了某公司生产的创意平板折叠桌的动态变化过程和一定条件下最优加工参数的设计问题。通过建立空间直角坐标系进行几何分析，构造非线性规划模型，并利用Matlab和Lingo软件编程求解，得出各种条件下的设计参数结果。在问题一中，本文从桌子的稳固性出发，从物理学的角度，根据受力分析，寻找稳固性条件下的约束条件，构建非线性规划模型，并利用Lingo求得单侧20根桌腿情况下的开槽长度、桌腿边缘线等参数，在此基础上描述了折叠桌折叠运动的动态过程。由于问题一采用了构建非线性规划的方法，因此在求解第二题的过程中，本文依然通过寻找约束条件和修改目标函数来优化模型，根据给出的桌面直径和桌子的高度，可以用Matlab求出各个需求的加工参数，所以此模型能够很好地满足设计者及生产者的需求。根据题目中给定的桌面高度70cm和桌面直径80cm的条件，用Matlab编程求解得出从外侧第2根木条到第10根槽长分别为11.86 cm，18.94 cm，24.54 cm，28.98 cm，32.42 cm，34.93 cm，36.58 cm，37.40cm。针对问题三，根据客户要求的折叠桌高度，桌面边缘线的形状大小和桌脚边缘线的大致形状，结合实际情况，发现现实生活中主要的桌面形状是分偶数边正多边形和椭圆形两种情况考虑，通过几何图形分析，分别建立非线性规划模型，根据题意寻找约束条件，优化模型，用Matlab编程求解。 关键词：非线性规划 几何分析 受力分析 空间直角坐标系二、问题重述某公司生产一种可折叠的桌子，桌面呈圆形，桌腿随着铰链的活动可以平摊成一张平板。桌腿由若干根木条组成，分成两组，每组各用一根钢筋将木条连接，钢筋两端分别固定在桌腿各组最外侧的两根木条上，并且沿木条有空槽以保证滑动的自由度。本文根据给定的各种数据研究以下几个问题：1. 给定长方形平板尺寸为120 cm 50 cm 3 cm，每根木条宽2.5 cm，连接桌腿木条的钢筋固定在桌腿最外侧木条的中心位置，折叠后桌子的高度为53 cm。试建立模型描述此折叠桌的动态变化过程，在此基础上给出此折叠桌的设计加工参数（例如，桌腿木条开槽的长度等）和桌脚边缘线（图4中红色曲线）的数学描述。2. 折叠桌的设计应做到产品稳固性好、加工方便、用材最少。对于任意给定的折叠桌高度和圆形桌面直径的设计要求，讨论长方形平板材料和折叠桌的最优设计加工参数，例如，平板尺寸、钢筋位置、开槽长度等。对于桌高70 cm，桌面直径80 cm的情形，确定最优设计加工参数。3. 公司计划开发一种折叠桌设计软件，根据客户任意设定的折叠桌高度、桌面边缘线的形状大小和桌脚边缘线的大致形状，给出所需平板材料的形状尺寸和切实可行的最优设计加工参数，使得生产的折叠桌尽可能接近客户所期望的形状。你们团队的任务是帮助给出这一软件设计的数学模型，并根据所建立的模型给出几个你们自己设计的创意平板折叠桌。要求给出相应的设计加工参数，画出至少8张动态变化过程的示意图。三、问题分析针对问题一，题目给定长方形平板尺寸为120 cm 50 cm 3 cm，每根木条宽2.5 cm，折叠后桌子的高度为53cm，在此问中为简便计算，我们暂不考虑桌面平板的厚度和桌腿之间的缝隙，按照最大20根桌腿的情况进行计算。同时对桌腿进行受力分析得知当时，桌子的稳固性最强，此时可以以的补角的余弦值最大为目标，建立非线性规划模型，在求解最外侧桌腿的参数的基础上，利用递推公式编程逐步得出20条桌腿的设计参数和折叠桌整体变化趋势。针对问题二，我们从问题中的用材最省、稳固性高、加工方便三方面来考虑，建立非线性规划模型。首先，从用材最省的角度出发，我们以平板体积最小为目标函数，其次，要满足桌子的稳定性，由问题一的分析，我们可以知道，在折叠情况下，当最外侧桌腿桌脚到桌面直径的平面距离要和桌面半径的距离保持一致，此时桌子的稳固性最高。然后，我们从加工方便的角度考虑，开槽长度越短，加工就相对方便一点，开槽长度最短也是我们的目标。最后，综合上面三个角度的考虑，建立非线性规划的数学模型，任意给出，桌子的高度和桌面的直径，都能够求得理想的加工参数。针对问题三，结合实际生活，我们只考虑桌面边缘形状为正多边形和椭圆形的两种情况，然后在问题二的基础上，对模型进行延伸和扩展，根据设计者的目的及要求，得到理想的桌子的加工参数。四、模型假设1.假设木条之间没有缝隙且木条表面足够光滑。2.假设桌子折叠展开的时候受力均衡。3.桌子一直处于水平状态。五、符号说明：分别为平板木块的长宽高：每条桌腿钢筋到相应桌腿顶点的距离：每条桌腿与圆形桌面下底面夹角（左侧）：桌腿编号：圆形桌面半径：各桌腿顶点处到过桌面中心垂直于桌腿的直径的距离：折叠桌的高度：桌腿长度 ：桌子平铺时各桌腿中钢筋到铰链的距离：折叠后各桌腿中钢筋到铰链距离：桌腿宽度：桌腿条数（单侧）：钢筋到相应桌腿的长度占整条桌腿的长度比例：偶数正多边形条数：偶数正多边形边长：桌面圆心与第条桌腿所对应的铰链的连线与桌面半径的夹角:第条桌腿的开槽长度：分别为椭圆形桌面长半轴和短半轴长六、模型的建立与求解1.问题一在第一题中，给定长方形平板尺寸为120 cm 50 cm 3 cm，每根木条宽2.5 cm，连接桌腿木条的钢筋固定在桌腿最外侧木条的中心位置，折叠后桌子的高度为53 cm。根据题意得一侧木条的最大数量为根，我们姑且不考虑缝隙，按最大根数计算。由于整个桌子是对称的，可以以桌子侧视图的二分之一来考虑，如图1设桌腿与桌面底部的夹角（左侧 ）为，从外侧向内第根桌腿的夹角为，桌腿钢筋到与桌面结合处的距离ED为，外侧向内第根桌腿的该距离为。图 1设桌子折叠开的高度为（不包含木板的3cm厚度），各桌腿与桌面结合处到圆形桌面直径的距离为，圆形桌面的半径为。实际上，只有四个最外侧桌腿接触到地面，为了保证桌子的稳固性，必须有约束条件：同时，对折叠桌桌腿进行受力分析，根据受力分析，可知桌腿张力越大，则夹角越小，即的补角的余弦值越大，也可表示为目标函数的增大。又根据题意和勾股定理，得约束条件。综上所述建立非线性规划模型如下：为了方便运算我们对桌腿进行编号，如图2，从最外侧的桌腿开始编号依次为1到10。图 2设桌腿长度为，根据题意有由此我们可以将，用多项式表示出来，第1条桌腿与桌面结合处到原面直径的距离：第2条桌腿与桌面结合处到原面直径的距离：第3条桌腿与桌面结合处到原面直径的距离：以此类推，第n条桌腿与桌面结合处到原面直径的距离为从而我们也可以将夹角间接表示出来：以此类推根据以上分析建立方程组如下： 利用Lingo编程求解最外侧木条的参数，结果如表1：表 1变量数值5053.656.3525.79然后我们根据所得结果，以与圆形桌面圆心相对应的地面点为原点，坐标系如图3所示。图 3另外，槽长等于。 根据递推公式编程（Matlab）得出的结果如表2所示：表 226.82506.35001.941853.6500025.238712.61241,708447.38764.676125.000916.28871.562343.71138.114625.164618.94921.456341.050810.938925.477320.98391.376939.016113.286225.828022.56261.316937.437415.215626.159323.77651.272036.223516.760826.439224.67941.239335,320617.943626.649225.30461.217134.695418.778826.778825.67241.204324.327619.2762最后根据我们建立的空间直角坐标系，给出各个桌脚的坐标后，得出边缘线图像如图4和图5所示： 图 4 图 5运用Matlab编程，用坐标表示折叠桌的动态变化趋势见图6。图 62.问题二设桌子的半径为，桌腿的宽度为，桌腿的条数为（单侧桌腿数）,桌腿的编号为。现在将桌腿按照前面所述的规则，从外到内依次编号，由于桌子是对称的，所以我们只考虑单侧桌腿的二分之一的桌子情况。根据问题一的分析，我们分析出桌腿数量的确定如下：根据问题的要求，首先考虑桌子的用材较少，以平板体积最小为目标函数，即，题目中对体积的影响较小，且根据实际生活经验，桌子的厚度差别不大，为了计算方便，我们直接将设定为。其次，当时，桌子的稳固性最高。综上，建立模型如下： 根据此模型，给出设计者想要的桌子半径和桌子高度，我们可以得到平板尺寸，然后考虑桌子加工要方便的角度，即开槽长度要最短。假设为平铺时钢筋到铰链的长度，为折叠时钢筋到铰链的长度，为开槽长度，为钢筋到桌脚的长度占整条桌腿长度的比例，出于安全和稳固性的考虑，开槽长度应该小于槽底到桌脚的长度，综上，建立模型如下： 根据此模型，结合模型，我们可以用求解出开槽的位置以及开槽的长度，根据题意给出的数据，利用求解的具体参数，求解 ，开槽总长度为，其它参数参见表3和表4，具体程序详见附录2表3158.002cm79.9cm3cm4.7cm17cm39.6272cm表4桌腿编号123456789(cm)0 11.8618.9424.5428.9832.4234.9336.5837.403.问题三问题三在问题二模型的基础上，进行相应的改进和扩展，根据现实需要和实际情况，我们只考虑桌面边缘线形状为正多边形和椭圆形的情况，正多边形情况下只考虑桌面边缘形状中两条边与平板边界重合的正多边形。(1)桌子边缘线形状数为多边形的情况当桌面边缘线为偶数正多边形时，给定任意折叠高度、桌面边缘线的形状大小（即给定边长长度和边数），便可求出平板材料的形状尺寸大小和切实可行的最优设计加工参数同问题二的求解类似假设正多边形的变数为，边长为，为桌面圆心与第条桌腿所对应的铰链的连线与桌面半径的夹角，则桌腿数量的确定，由桌腿宽度和木板宽度决定，桌腿数量的确定如下：根据用材最少，桌子稳固性较高，建立模型如下： 根据加工方便，建立模型如下：综上，结合模型和模型即为桌面边缘形状偶数正多边形的模型，给出相应的多边形的边长和桌子高度，即可得到优化模型后的桌子参数。(2)桌子边缘线形状为椭圆形的情况当桌面边缘线椭圆时，给定任意折叠高度、桌面边缘线的形状大小（即给定短半轴、长半轴大小），便可求出平板材料的形状尺寸大小和切实可行的最优设计加工参数同问题二的求解类似。设椭圆形桌面长半轴和短半轴长分别为,桌腿宽度为，桌腿条数为，从最外侧桌腿到中间桌腿依次编号依次为那么三者有以下关系：当为整数时，必定使折叠后的桌面上最外侧桌腿靠外的边与与椭圆相切，若是去掉桌腿，将桌腿在铰链处桌面下底面各突出的端点相连，便形成一个椭圆。当不为整数时，对向下取整，即，将桌腿在铰链处桌面下底面各突出的端点相连，也可以形成一个椭圆。设圆形桌面最外侧桌腿铰链到椭圆桌面长轴的距离为，用勾股定理可得：，设为从外侧看过去对圆形桌面上桌腿铰链到圆形桌面直径的距离,每根桌腿从最外侧开始依次编号为，可以求得：下面先从折叠桌的稳固性分析桌角的位置。由证明的模型定理可知，要使折叠桌稳定性好，则桌子的四条支撑腿所构成的结构面成正方形，即方程，可满足桌子的稳固性。在稳固性最好的基础上再从用材最少，列出下列非线性规划模型求出平板尺寸和圆形桌面最外侧桌腿铰链到圆形桌面直径的距离。 根据上述非线性规划模型，任意给定折叠桌高度和桌面长半轴和短半轴长长度，用Matlab编程便可计算出平板尺寸和椭圆桌面最外侧桌腿铰链到椭圆直长半轴的距离，桌腿宽度为及桌腿条数为。当给定折叠桌高度为和桌面长半轴和短半轴长长度和，用Matlab算出此时的平板尺寸、椭圆桌面最外侧桌腿铰链到的椭圆直长半轴距离，桌腿宽度为及桌腿条数为（程序见附录3，计算出相应参数结果如表5：表5 平板尺寸参数表161.52cm58.00cm3.00cm2.90cm2024.46cm加工方便程度主要受到开槽长度大小的影响，在稳固性最好、材料最省的基础上再以开槽长度最小为为目标，满足一定约束下列出非线性规划模型，求出开槽长度和钢筋位置。为计算方便，我们令椭圆的外切圆的半径近似等于桌面圆心与相对应铰链的长度，非线性规划模型如下： 根据上述非线性规划模型，在稳固性最好、材料最省的基础上，任意给定折叠桌高度和圆形桌面直径，便可先求出相应的最优的平板尺寸、圆形桌面最外侧桌腿铰链到圆形桌面直径的距离，桌腿宽度为及桌腿条数为，对于当给定折叠桌高度为70cm和桌面长半轴和短半轴长长度45cm和30cm,用MATLAB编程（程序见附录3）求解，求得，即最外侧钢筋位置距该桌腿端点距离占该条桌腿总长的比例为0.48。此时各条桌腿中开槽长度下表6：表6 各条桌腿开槽长度桌腿编号12345678910(cm)0 10.2918.5725.4131.0035.5239.0841.7843.6744.78在折叠桌高度由0cm变化到70cm时，取高度分别为3cm、12cm、22cm、31cm、41cm、50cm、60cm、70cm，用MATLAB编程即可分别画出给定高度的桌子状态图如图7和图8所示：图7图8七、模型分析与评价1.优点(1)本文模型建立具有坚实可靠的数学理论基础，易于证明其合理性。(2)模型方程通俗易懂，简单易行，易于推广到实际生活中。2.缺点(1)模型一建立的非线性规划模型没有考虑缝隙对折叠动态变化的影响。(2)模型三只考虑了偶数正多边形和椭圆形状两种情形，没有完全考虑实际生活中所有设计者的需求。(3)模型的实用性有待加强。八、参考文献1杨文茂，李全英，空间解析几何，武汉：武汉大学出版社，2006.2胡良剑，孙晓君，MATLAB数学实验，北京：高等教育出版社，2006. 九、附录附录1：Lingo程序min=cos(H/L);sqrt(L2 - H2) + a =r;a + L = 60; !平板长度为60;r = sqrt(252 + a2);H = 50; !桌面高度;求解加工参数function zhuozi = jgcs( a1,n,h )%n为木块的总数，d1为桌腿钢筋到桌面木块顶点的距离%angel1为第一根桌腿与桌面形成的夹角%求桌子各技术参数的值 l1=60-a1; angle1=asin(h/l1); d1=30-a1/2; zhuozi(1).d=d1; zhuozi(1).a=a1; zhuozi(1).angle=acos(h/l1)+pi/2; zhuozi(1).l=l1; zhuozi(1).t=l1*cos(angle1); zhuozi(1).caochang=0; zhuozi(1).h=0; r=(sqrt(625+a12); x=0; for i=2:n zhuozi(i).a=real(sqrt(r2-(2.5*(11-i)2);%求桌面木块最外侧顶点到直径的距离 zhuozi(i).l=60-zhuozi(i).a;%每条桌腿长度 zhuozi(i).d = real(sqrt(zhuozi(i).a - a1)2 + (l1/2)2 - 2*(zhuozi(i).a - a1)*(l1/2)*cos(angle1); %完全收缩时桌腿钢筋处到桌面木块最外侧顶点的距离 zhuozi(i).caochang = zhuozi(i).d + zhuozi(i).a - (l1/2) - a1; %槽长 zhuozi(i).angle=acos(zhuozi(i).a - a1)2 + zhuozi(i).d2 -(l1/2)2)/(2*(zhuozi(i).a - a1)*zhuozi(i).d); %求桌腿和桌面形成的角度 zhuozi(i).angle=acos(zhuozi(i).a - a1)2 + zhuozi(i).d2 -d12)/(2*(zhuozi(i).a - a1)*zhuozi(i).d); if (zhuozi(i).anglepi/2) zhuozi(i).h=h-cos(zhuozi(i).angle-pi/2)*zhuozi(i).l; %桌腿的高度 zhuozi(i).t=zhuozi(i).a+real(sqrt(zhuozi(i).l2-(h-zhuozi(i).h)2); %求桌腿到直径的距离 else zhuozi(i).h=h-sin(zhuozi(i).angle)*zhuozi(i).l; %桌腿的高度 zhuozi(i).t=zhuozi(i).a-real(sqrt(zhuozi(i).l2-(h-zhuozi(i).h)2); end endend求解边缘线matlab程序x1 = w,fliplr(w); y1 = -25:2.5:-2.5,2.5:2.5:25; z1 = -fix(g,fliplr(g);%fix消除误差带来的影响plot3(x1,y1,z1);grid on;变化趋势%变化形态S=1for h=0:12.5:50 a1=6.35; l1=60-a1; angle1=asin(h/l1); d1=30-a1/2; zhuozi(1).d=d1; zhuozi(1).a=a1; zhuozi(1).angle=acos(h/l1)+pi/2; zhuozi(1).l=l1; zhuozi(1).t=l1*cos(angle1); zhuozi(1).caochang=0; zhuozi(1).h=0; r=(sqrt(625+a12); x=0; for i=2:10 zhuozi(i).a=real(sqrt(r2-(2.5*(11-i)2);%求桌面木块最外侧顶点到直径的距离 zhuozi(i).l=60-zhuozi(i).a;%每条桌腿长度 zhuozi(i).d = real(sqrt(zhuozi(i).a - a1)2 + (l1/2)2 - 2*(zhuozi(i).a - a1)*(l1/2)*cos(angle1); %完全收缩时桌腿钢筋处到桌面木块最外侧顶点的距离 zhuozi(i).caochang = zhuozi(i).d + zhuozi(i).a - (l1/2) - a1; %槽长 zhuozi(i).angle=acos(zhuozi(i).a - a1)2 + zhuozi(i).d2 -(l1/2)2)/(2*(zhuozi(i).a - a1)*zhuozi(i).d); %求桌腿和桌面形成的角度 zhuozi(i).angle=acos(zhuozi(i).a - a1)2 + zhuozi(i).d2 -d12)/(2*(zhuozi(i).a - a1)*zhuozi(i).d); if (zhuozi(i).anglepi/2) zhuozi(i).h=h-cos(zhuozi(i).angle-pi/2)*zhuozi(i).l; %桌腿的高度 zhuozi(i).t=zhuozi(i).a+real(sqrt(zhuozi(i).l2-(h-zhuozi(i).h)2); %求桌腿到直径的距离 else zhuozi(i).h=h-sin(zhuozi(i).angle)*zhuozi(i).l; %桌腿的高度 zhuozi(i).t=zhuozi(i).a-real(sqrt(zhuozi(i).l2-(h-zhuozi(i).h)2); end end x1= -zhuozi.t,fliplr(-zhuozi.t); y1=25:-2.5:2.5,-2.5:-2.5:-25; z1=fix(zhuozi.h,fliplr(zhuozi.h); x2=-zhuozi.a,fliplr(-zhuozi.a); y2=y1; z2=h*ones(1,20); if S = 1 subplot(3,2,1:2); view(-10,75); else subplot(3,2,S+1); view(-25,25); end title(当桌子高度为,num2str(h), 时); hold on for i = 1:20 plot3(x1(i) x2(i),y1(i) y2(i),z1(i) z2(i) end plot3(x1,y1,z1,r); plot3(x2,y2,z2); for i = 1:20 plot3(-x1(i) -x2(i),y1(i) y2(i),z1(i) z2(i) end plot3(-x1,y1,z1,r); plot3(-x2,y2,z2); %桌面与桌腿的连线 plot3(x2(1) -x2(1),y2(1) y2(1),z2(1) z2(1); plot3(x2(20) -x1(20),y2(20) y1(20),z2(20) z2(20); %钢筋 x3 = -(sqrt(zhuozi(1).l/2)2 - (h/2)2) + zhuozi(1).a); y3 = 22.5; z3 = h/2; plot3(x3 x3,y3 -y3,z3 z3,k); plot3(-x3 -x3,y3 -y3,z3 z3,k); grid on hold off S=S+1 End附录2：tich = input(桌高:);R = input(桌面直径：);c = 3; %平板厚度i = 1; %记录坐标for d=1:0.1:10 n1 = R/d; if(n1 = floor(R/d) n = n1-1; else n = floor(R/d); end syms a u1=solve(sqrt(a/2 - sqrt(R/2)2 - (n/2*d)2)2 - (h - c)2) + sqrt(R/2)2 - (n/2*d)2) - n/2*d,a); u1 = vpa(u1); u(i) = eval(u1(1); i = i + 1;endu2 = find(u(1,:) = min(u); while(1) a = u(max(u2); d = 1 + (max(u2) - 1)*0.1; n1 = R/d; if(n1 = floor(R/d) n = n1-1; else n = floor(R/d); end b1 = n*d; %平板宽度 k = sqrt(R/2)2 - (n/2*d)2); %最外侧桌面弦长 o = asin(h - c)/(a/2 - k); %最外侧桌腿与桌面夹角 j = 1; for v = 0:0.00001:1 s8 = 0; for i = 1:ceil(n/2) b5(i) = sqrt(R/2)2 - (n/2 - i + 1)*d)2); %桌面弦长 f5(i) = a/2 - b5(i); %每条桌腿长度 s5(i) = f5(i) - v*f5(1); %开槽顶部到桌腿顶部长度 if i = 1 s6(1) = (1 - v)*f5(1); %开槽底部到桌腿顶部长度 else s6(i) = sqrt(1 - v)*f5(1)2 + (b5(i) - b5(1)2 - 2*(b5(i) - b5(1)*(1 - v)*f5(1)*cos(o); end s7(i) = s6(i) - s5(i); %开槽长度 if s6(i) - s5(i) f5(1)*v s8 = inf; else s8 = s8 + s7(i); %开槽总长 end end if(n/2 = floor(n/2) s8 = 4*s8; else s8 = 4*s8 - 2*s7(i); end s9(j) = s8; j = j+1; end v = (find(min(s9) = s9(1,:)-1)*0.00001; %比率 v = v(1); if inf = s8 break; else u2 = u2 + 1; endendh1 = h;N = 4; %子图个数I = 1; %第i个子图for h = 3 :1:h1 h2 = 3:(h1-3)/(N-1):h1; %取出图像 figure(1); if floor(h2(I) = h figure(2); subplot(2,2,I); view(25,25); if I N %防止越界 I = I + 1; end title(高度为,num2str(h); end o = asin(h - c)/(a/2 - k); %最外侧桌腿与桌面夹角 for i = 1:ceil(n/2) b(i) = sqrt(R/2)2 - (n/2 - i + 1)*d)2); %桌面弦长 f(i) = a/2 - b(i); %每条桌腿长度 s(i) = f(i) - v*f(1); %开槽顶部到桌腿顶部长度 if i = 1 s1(1) = (1 - v)*f(1); %开槽底部到桌腿顶部长度 else s1(i) = sqrt(1 - v)*f(1)2 + (b(i) - b(1)2 - 2*(b(i) - b(1)*(1 - v)*f(1)*cos(o); end s2(i) = s1(i) - s(i); %开槽长度 if i = 1 q(i) = pi - o; else q(i)=acos(b(i) - k)2 + s1(i)2 -s1(1)2)/(2*(b(i) - k)*s1(i); %角度 end if h = 3 %处理0时角度 if(q(i) pi/2) g(i) = cos(q(i) - pi/2)*f(i); %每条桌腿的高度 w(i) = b(i) + sqrt(f(i)2 - g(i)2); %横坐标 else g(i) = sin(q(i)*f(i); %每条桌腿的高度 w(i) = b(i) - sqrt(f(i)2 - g(i)2); end end x = w,fliplr(w); y = -b1/2:d:-d/2,d/2:d:b1/2; z = -fix(g,fliplr(g); x1 = b,fliplr(b); y1 = y; z1 = zeros(1,2*ceil(n/2); n1 = 2*ceil(n/2); cla axis(-a/2 a/2 -b1/2 b1/2 -h1 h1/2); axis equal set(gcf, position, 200 200 800 500); hold on %画出桌面 for i = 1:n1 plot3(x(i) x1(i),y(i) y1(i),z(i) z1(i) end plot3(x,y,z,r); plot3(x1,y1,z1); for i = 1:n1 plot3(-x(i) -x1(i),y(i) y1(i),z(i) z1(i) end plot3(-x,y,z,r); plot3(-x1,y1,z1); %桌面连线 plot3(x1(1) -x1(1),y1(1) y1(1),z1(1) z1(1); plot3(x1(n1) -x1(n1),y1(n1) y1(n1),z1(n1) z1(n1); %钢筋 x3 = (1 - v)*(w(1) - k) + k; y3 = b1/2; z3 = -fix(1 - v)*g(1); plot3(x3 x3,y3 -y3,z3 z3,k); plot3(-x3 -x3,y3 -y3,z3 z3,k); view(20,25); pause(0.05);enddisp(平板长度：); adisp(平板宽度：); b1disp(平板厚度); cdisp(桌腿宽度：); ddisp(单侧桌腿条数：); ndisp(最外侧桌腿钢筋位置比率); vdisp(最外侧桌腿钢筋位置距铰链距离); f(1)*(1 - v)disp(每条桌腿长); fdisp(每条桌腿开槽顶端距铰链距离); sdisp(每条桌腿铰链距平板中线距离); bdisp(开槽长度：); s2disp(开槽总长度：); min(s9)toc附录3：tich = input(桌高:);R = input(桌面短轴：);A = input(桌面长轴：);c = 3; %平板厚度i = 1; %记录坐标for d=1:0.1:10 n1 = R/d; if(n1 = floor(R/d) n = n1-1; else n = floor(R/d); end syms a u1=solve(sqrt(a/2 - sqrt(A/2)2 - (A/R*(n/2)*d)2)2 - (h - c)2) + sqrt(A/2)2 - (A/R*(n/2)*d)2) - n/2*d,a); u1 = vpa(u1); u(i) = eval(u1(1); i = i + 1;endu2 = find(u(1,:) = min(u);while(1) a = u(max(u2); d = 1 + (max(u2) - 1)*0.1; n1 = R/d; if(n1 = floor(R/d) n = n1-1; else n = floor(R/d); end b1 = n*d; %平板宽度 k = sqrt(A/2)2 - (A/R*(n/2)*d)2); %最外侧桌面弦长 o = asin(h - c)/(a/2 - k); %最外侧桌腿与桌面夹角 j = 1; for v = 0:0.00001:1 s8 = 0; for i = 1:ceil(n/2) b5(i) = sqrt(A/2)2 - (A/R*(n/2 - i + 1)*d)2); %桌面弦长 f5(i) = a/2 - b5(i); %每条桌腿长度 s5(i) = f5(i) - v*f5(1); %开槽顶部到桌腿顶部长度 if i = 1 s6(1) = (1 - v)*f5(1); %开槽底部到桌腿顶部长度 else s6(i) = sqrt(1 - v)*f5(1)2 + (b5(i) - b5(1)2 - 2*(b5(i) - b5(1)*(1 - v)*f5(1)*cos(o); end s7(i) = s6(i) - s5(i); %开槽长度 if s6(i) - s5(i) f5(1)*v s8 = inf; else s8 = s8 + s7(i); %开槽总长 end end if(n/2 = floor(n/2) s8 = 4*s8; else s8 = 4*s8 - 2*s7(i); end s9(j) = s8; j = j+1; end v = (find(min(s9) = s9(1,:)-1)*0.00001; %比率 v = v(1); if inf = s8 break; else u2 = u2 + 1; endendh1 = h;N = 8; %子图个数I = 1; %第i个子图for h = 3 :1:h1 h2 = 3:(h1-3)/(N-1):h1; %取出图像 figure(1); if floor(h2(I) = h if I = N/2 %分开画图 figure(2); subplot(2,2,I); else figure(3); subplot(2,2,I-N/2); end view(25,25); if I N %防止越界 I = I + 1; end title(高度为,num2str(h); end o = asin(h - c)/(a/2 - k); %最外侧桌腿与桌面夹角 for i = 1:ceil(n/2) b(i)