现代设计方法计算题

1.试用进退法确定函数f（x）=的一维优化初始区间[a,b],给定初始点=-1，初始步长h=1。解：==-1，=f（）=9=+h=-1+1=0，=f（）=5比较，，由于，作前进计算：=+2h=-1+2=1，=f（）=3比较，，由于，再作前进计算：=，==5==1，==3=+4h=-1+4=3，=f（）=5由于，可知初始区间已经找到，即[a,b]=[0,3]。2.设某种单元的可靠度(t)=，其中=0.001/h，试求出：（1）由这种单元组成的二单元串联系统，二单元并联系统及2/3（G）表决系统的平均寿命；（2）当t=100h、500h、1000h时，一单元、二单元串联、二单元并联及2/3(G)表决系统的可靠度，并加以比较。解：（1）一个单元与系统的平均寿命分别为：单=1/=1000h串=1/2=500h并=3/2=1500h（）=5/6=833.3h（2）当t=100h时，一个单元与系统的可靠度分别为：单==0.905串=单==0.819并=（单）=1-（）=0.991=单-单=0.975当t=500h时，一个单元与系统的可靠度分别为：单==0.6065串=单==0.3678并=1-（单）=1-（）=0.8452（）=单-单=0.6575当t=1000h时，一个单元与系统的可靠度分别为：单==0.368串=单==0.135并=1-（单）=1-（）=0.600=单单=0.306从计算结果可以看出：（1）一个单元的可靠度高于二单元串联系统的可靠度，但低于二单元并联系统的可靠度；（2）2/3（G）系统的平均寿命为一个单元的平均寿命的5/6倍，明显低于一个单元的平均寿命。3.已知约束优化问题：minf（x）=+S．t．-1≥0试写出内点罚函数与外点罚函数的表示式。解：内点罚函数：Φ（）=++外点罚函数：当-+1≤0，Φ（x，）=+当--+1＞0，Φ（x，）=++（）4.现在要用钢板制作一个有盖的长方本储水箱，要求各边长均不超过20厘米，且长度为宽度的2倍，试确定三边长度值，使该储水箱的容积最大，要求其表面积不超过400平方厘米。解：（1）建立数学模型用复合形法迭代3次。取储水箱长和高为设计变量，，则其宽0.5，数学模型为maxF（X）=0.5s．t．+3≤4000≤≤200≤≤20（2）用复合形法求解求得的近似结果为∗={}T={77}TF(X*)=5091已知右上图所示等腰直角三角形的单元刚度矩阵为：[K]（e）=E4[对称]右图所示薄板结构中节点2处所受载荷以及材料的弹性模量和板厚分别为：F100KN，E=2×7N/cm，t=0.1cm求节点2处的各位移分量。解得Fx7.7×N,FY=-70.7×NE[]{uv}={7777}u=-7.07×cm，v44cm5.用梯度法求下列无约束优化问题：MinF（X）=+，设初始点取为（）={}，以梯度模为终止迭代准则，其收敛精度为5.（1）求初始点梯度▽F（X）▽F（X）={8}T▽F(={4}T(2)第一次搜索丨▽F（（））丨=16.5，s（）=-▽F（（））/16.5=-{497}Tα7（）=（）+α（）S（）={479}T▽F()={978}T丨▽F()丨=3.043＜5.0故满足要求，停止迭代。最优点X*={47，9}T最优值F（∗）=2.216.节点和单元划分如图示的由两根杆组成的平面刚架结构，在节点3处作用大小为F的集中载荷，两单元在局部坐标系下的单元刚度矩阵相同，即[k][k]a[]其中，a为常数。试引人支承条件写出总体平衡方程。解：先求单元在总体坐标系下的单元刚度矩阵。单元（2）在总体坐标系下的单元刚度矩阵与在局部坐标系下的单元刚度矩阵相同，即iil2yF3jx[K][k]a[]单元（1）的坐标转换矩阵中的β=-90°，单元（1）的坐标转换矩阵为[T]=[cos9°）sin（9°）sin（9°）cos（9°）cos（9°）sin（9°）sin（9°）cos（9°）]=[]所以单元（1）在总体坐标系下的单元刚度矩阵为：[K]=[T]T[K][T]=[]a[][]=[]单元局部编码和总体编码的对应关系为：单元（1）ij→12单元（2）ij→23单元刚度矩阵中子块对应关系为；[K]（）=[kkkk]（），[k]（）=[kkkk]（）所以总体刚度矩阵为：[K]=[k（）k（）k（）k（）k（）k（）k（）k（）]=a[4]节点的位移矢量为：{uvuvuv}T约束条件为：u=0，v=0，=0作用到结构上的外力为：Fy=_F所以引入支承条件的平衡方程为：a[4]{uvuvuv}={_F}7一组实验数据如下，试用抛物线插值方法计算X=92和X=198处的Y值。i90100110120130140150Yi0.680.740.790.830.860.890.92【参考答案】抛物线插值公式为：Y（x）=（xx2）（xx3）（x1x2）（x1x3）y+（xx1）（xx3x2−x1x2x3y+xx1xx2x3−x1x3x2y当x=136时∵x∈（130,140），|136-130|＞|136-140|∴选择插值节点：（，y）=（130,0.86），（，y）=（140,0.89），（，y）=（150,0.92）将以上数据和x=136代入抛物线插值公式，得x-136时y值为：Y（136）=（64）（6）（4）（）×0.86+（6）（6）（4）（4）×0.89+（6）（64）（）（4）×0.92=0.8788如图所示的平面刚架，由两个单元（1）和（2）组成，两单元的长度和载面尺寸及材料特性相同，单元（1）的局部坐标正方向为沿轴线方向节点1指向节点2，单元（2）的局部坐标正方向为沿轴线方向由节点3指向节点1，在局部坐标系下每个单元的刚度矩阵为[k][k]（）A[44]（1）求刚架总体刚度矩阵[K]。（2）引入支撑条件，写出平衡方程。平面刚架的坐标转换矩阵为。T[e]=[cosasinasinacosacosasinasinacosa]-由局部坐标系与总体坐标系的关系知：单元（1）a=0，单元（2）a=所以T=[1]T（）=[]在整体坐标系下，单元的刚度矩阵为：[k]e[Te）][K]e[Te]∴[K]=[T]e[K][T]-[K][K]-[T]T[K][T]-[K]-A[44]单元（1）局部码对应的总码为，2，单元（2）局部码对应的总码为3.19.已知某零件的工作应力和材料强度均服从指数分布，且强度和应力的均值分别为µr=210Mpa和µs=160Mpa，试确定零件的可靠度。零件的工作应力和材料强度均服从指数分布，且µr=210MPA；µs=160Mpa∴s=µs，r=µrR=ssr=µrµrµs=6=0.5675676该零件的可靠度为：R=0.5675676将下列实验测试数据拟合成y=ab形式的经验公式。（计算过程中保留小数点后两位）i1.181.582.403.003.80yi2.763.294.234.835.57将Y=axb两边取对数，得：lny=lna+blnx令U=lnyA=lnaB=bV=lnx，则原式变为：U：A+BV将表数据取对数：Vi0.170.460.881.101.34Ui1.021.191.441.571.72按以上的Vi，Ui进行最小二乘拟合得{5A+（∑Vii=）B=∑Uii=（∑Vii=）A+（∑Vii=）B=∑ViUii=代入数据得{A9B949A4B求解得A=0.92B=0.60a=A=2.5b=B=0.6拟合的经验公式为y=610.如图所示的平面刚架，由两个单元（1）和（2）组成，两单元的长度和载面尺寸及材料特性相同，单元（1）的局部坐标正方向为沿轴线方向节点1指向节点2，单元（2）的局部坐标正方向为沿轴线方向由节点3指向节点1，在局部坐标系下每个单元的刚度矩阵为.[k]（）[k]（）A[44]（1）求刚架总体刚度矩阵[K]。（2）引入支撑条件，写出平衡方程。解.平面刚架的坐标转换矩阵为T[e]=[cosasinasinacosacosasinasinacosa]由局部坐标系与总体坐标系的关系知：单元（1）a=0，单元（2）a=90°所以T=[1]T（）=[]在整体坐标系下，单元的刚度矩阵为：[k]e[Te）]T[K]e[Te]∴[K]=[T]T[K][T]=[K][K]-[T]T[K][T]=[K]=A[44]单元（1）局部码对应的总码为，2，单元（2）局部码对应的总码为3,1∴[K]（）=[kkkk]（）[K]（）=[kkkk]（）所以按照刚度集成法，可得出总体风度矩阵为：[K]=[[K][K][K][K][K][K][K][K]]=A[844]由于只有在节点2处作用有沿Y轴负方向的外载荷F=100N，所以节点载荷矢量为{F}={}T支撑条件为u=v==0，所以等式右端的力矢量无须修改，矩阵[K]中7至9的各行各列修改成除主对角线元素为1外，其余各元素均为零1分总体平衡方程为A[84]{uvuvuv}=[]11.某机电系统由10台相同设备组成，各设备可靠度为0.9，若该系统至少有7台设备正常运行就可以保证整个系统正常工作，试求该系统的可靠度。解：该系统的每台设备或是正常工作或是发生故障，其失效数为正整数。因此是离散型随机变量，且服从二项分布。∴系统的可靠度由下式计算R（r）=∑𝐂𝐍𝐫𝐢=𝐑𝐍F由题意知：r=3，N=10，R=0.9，F=1-0.9=0.1，则R（3）=！！！×𝟗（）×+！！！×𝟗（）×+！！！×𝟗（）×+！𝟑！𝟑！×𝟗（𝟑）×𝟑≈0.9872该系统的可靠度为0.987212.已知△ABC=[]，将该三角形沿X方向移动1个单元，沿Y方向移动2个单位后，再放大一倍，求变换后△ABC各顶点的坐标。这是一组合变换，先求出组合变换的变换矩阵。沿X方向移动1个单位，沿Y方向移动2个单位，变换矩阵为T=[]放大一倍，变换矩阵为：T=[]所以组合变换矩阵为：T=TT=[][]=[4]所以变换后的三角形顶点的坐标短阵为：△AˊBˊCˊ=△ABC*T=[][4]=[444]变换后的三角形各顶点的坐标为：A（2，4），B（6,4），C（4,6）13.用最小二乘法将下列数据拟合成Y=aaa形式的经验公式。（计算结果中保留两位小数）i1.202.403.204.545.82Yi7.6316.4124.1940.6560.41解：拟合公式为f（x）=y=aax+a，表中共有5组数据则m=5由最小二乘法拟合思想得：{aa∑ii=a∑ii=∑yii=a∑ii=a∑i=a∑ii=∑ii=yia∑ii=a∑ii=a∑i4i=∑ii=yi将表中数据代入得：{a7a79a4997a79a9a979a9a79a7解得a=2.31，a=3.00，a=1.20所以，由表中数据拟合成的经验公式为：_y=2.31+3.00x+1.2014.已知目标函数：minF（x）=-4-约束条件：g（x）=--3≤，0g（x）=--2≤0g（x）=+-4≤0g4（x）=≥0试用Kuhn-Tucker条件判断点（）T和（）T是否为该有约束问题的极值点。15.（1）先确定起作用约束将x={}T代入约束方程，只有g（x）=0和g（x）=0所以起作用的约束为g（x）和g（x）(2)根据K-T条件，有▽F（x）=-∑ni=▽g（x）▽F(x)={4}▽g={}▽g={}在x={}T处-{}{4}+{}解得=-2/3=2/3不满足K-T条件，所以x={}T不是约束最优点。将x={}T代入约束方程，只有g和g=0所在起作用的约束方程为g和g（根据K-T条件，有▽F(x)=-∑ni=▽g(x)▽F(x)={4}▽g={}▽g={}在x={}T处-{}={}+{}解得=1/4=1/4满足K-T条件，所以x={}T是约束极值点。16.一厚度为t，边长为l的正方形钢板，其支撑和受力情况如图示，按平面问题进行有限元分析时，划分的单元、结点的局部和整体编码如图所示。两三角形单元在整体坐标系下的单元刚度矩阵为；[k][k]=[对称]试引入支撑条件写出平衡方程。单元节点局部编码与总体编码的对