现代设计方法计算题.docx
- 文档编号:9573504
- 上传时间:2023-02-05
- 格式:DOCX
- 页数:47
- 大小:45.02KB
现代设计方法计算题.docx
《现代设计方法计算题.docx》由会员分享,可在线阅读,更多相关《现代设计方法计算题.docx(47页珍藏版)》请在冰豆网上搜索。
现代设计方法计算题
1.试用进退法确定函数 f(x)=𝑥2 ‒ 3𝑥 + 5的一维优化初始区间[a,b],给定初始
点𝑥0=-1,初始步长 h=1。
解:
𝑥1=𝑥0=-1,𝑓1=f(𝑥1)=9
𝑥2=𝑥0+h=-1+1=0,𝑓2=f(𝑥2)=5
比较𝑓1,𝑓2,由于𝑓1>𝑓2,作前进计算:
𝑥3=𝑥0+2h=-1+2=1,𝑓3=f(𝑥3)=3
比较𝑓2,𝑓3,由于𝑓2>𝑓3,再作前进计算:
𝑥1=𝑥2 = 0,𝑓1=𝑓2=5
𝑥2=𝑥3=1,𝑓2=𝑓3=3
𝑥3=𝑥0+4h=-1+4=3,𝑓3=f(𝑥3)=5
由于𝑓2 < 𝑓3,可知初始区间已经找到,即[a,b]=[0,3]。
2.设某种单元的可靠度𝑅0(t)=𝑒 ‒ 𝜆𝑡,其中𝜆=0.001/h,试求出:
(1)由这种单元组成的二单元串联系统,二单元并联系统及 2/3(G)表决系
统的平均寿命;
(2)当 t=100h、500h、1000h 时,一单元、二单元串联、二单元并联及 2/3(G)
表决系统的可靠度,并加以比较。
解:
(1)一个单元与系统的平均寿命分别为:
𝜃单=1/𝜆 =1000h𝜃2串=1/2𝜆=500h
𝜃2并=3/2𝜆=1500h𝜃2/3(𝐺)=5/6𝜆=833.3h
(2)当 t=100h 时,一个单元与系统的可靠度分别为:
2
𝑅单=𝑒 ‒ 0.001 × 100=0.905𝑅2串=𝑅单=𝑒 ‒ 0.2=0.819
𝑅2并=1 - (1 - 𝑅单)2=1-(1 ‒ 𝑒 ‒ 0.1)2=0.991
𝑅2/3(𝐺)=3𝑅单2-2𝑅单3=0.975
当 t=500h 时,一个单元与系统的可靠度分别为:
2
𝑅单=𝑒 ‒ 0.001 × 500=0.6065𝑅2串=𝑅单=𝑒 ‒ 0.5 × 2=0.3678
𝑅2并=1-(1 ‒ 𝑅单)2=1-(1 ‒ 0.6065)2=0.8452
𝑅2/3(𝐺)=3𝑅单2-2𝑅单3=0.6575
当 t=1000h 时,一个单元与系统的可靠度分别为:
3.已知约束优化问题:
minf(x)=𝑥1+3𝑥2
2
𝑅单=𝑒 ‒ 0.001 × 1000=0.368𝑅2串=𝑅单=𝑒 ‒ 2=0.135
𝑅2并=1-(1 ‒ 𝑅单)2=1-(1 ‒ 𝑒 ‒ 1)2=0.600𝑅2/3(𝐺)=3𝑅单2 ‒ 2𝑅单3=0.306
从计算结果可以看出:
(1)一个单元的可靠度高于二单元串联系统的可靠度,但低于二单元并联
系统的可靠度;
(2)2/3(G)系统的平均寿命为一个单元的平均寿命的 5/6倍,明显低
于一个单元的平均寿命。
22
S.t.𝑥1 + 𝑥2-1≥0
试写出内点罚函数与外点罚函数的表示式。
解:
内点罚函数:
22
1
当-𝑥1 ‒ 𝑥2+1≤0,Φ(x,𝑟𝑘)=𝑥1+3𝑥2
外点罚函数:
22
当-𝑥1-𝑥2+1>0,Φ(x,𝑟𝑘)=𝑥1+3𝑥2+𝑟𝑘( - 𝑥1 - 𝑥2 + 1)
22
2
maxF(X)=0.5𝑥1 𝑥2
s.t.𝑥 1+3𝑥1𝑥2≤400
4.现在要用钢板制作一个有盖的长方本储水箱,要求各边长均不超过 20 厘米,
且长度为宽度的 2 倍,试确定三边长度值,使该储水箱的容积最大,要求其表
面积不超过 400 平方厘米。
解:
(1)建立数学模型 用复合形法迭代 3 次。
取储水箱长和高为设计变量𝑥1,𝑥2,则其宽 0.5𝑥1,数学模型为
2
2
0≤𝑥1≤200≤𝑥2≤20
(2)用复合形法求解
求得的近似结果为𝑋 ={𝑥1,𝑥2}
∗
𝑇 =
{11.5.7.7}𝑇F(X*)=509
1 已知右上图所示等腰直角三角形的单元刚度矩阵为:
[𝐾](𝑒)=
𝐸𝑡
4
[
3
1
3
对
‒ 2
0
2
称
‒ 1
‒ 1
0
1
‒ 1
‒ 1
0
1
1
0
‒ 2
0
0
0
2
]
右图所示薄板结构中节点 2 处所受载荷以
及材料的弹性模量和板厚分别为:
𝐹2 = 100KN,E=2×107N/cm,t=0.1cm
求节点 2 处的各位移分量。
33
𝐸𝑡
𝑡
[
2 0 𝑢2 ‒ 70.7 × 103
5.用梯度法求下列无约束优化问题:
MinF(X)=𝑥1+4𝑥2,设初始点取为𝑋(0)=
{2,2} ,以梯度模为终止迭代准则,其收敛精度为 5.
𝑢2=-7.07×10 ‒ 2cm,𝑣2 =‒ 0.1414𝑐𝑚
22
𝑇
(1) 求初始点梯度▽F(X)
▽F(X)={2𝑋1,8𝑋2} ▽F(𝑋
𝑇
(0)
)={4.16}𝑇
(2)第一次搜索
(0)
(0)
(0)
𝑇
𝛼(0) = 2.157
𝑋
(1)=𝑋(0)+𝛼(0)𝑆(0)={1.476, ‒ 0.923}𝑇
▽F(𝑥
(1)
)= {2.952, ‒ 0.738}T
丨▽F(𝑋
(1)
)丨=3.043<5.0
最优值 F(𝑋 )=2.21
故满足要求,停止迭代。
最优点 X*={1.476, ‒ 0.0923}
∗
𝑇
6.节点和单元划分如图示的由两根杆组成的平面刚架结构,在节点 3 处作用大
小为 F 的集中载荷,两单元在局部坐标系下的单元刚度矩阵相同,即
[
2
0
0
0
‒ 2
[𝑘]
(1) = [𝑘]
(2) = a 0
0
1
1
0
‒ 1
‒ 1
0
1
1
0
‒ 1
‒ 1
0
0
0
2
0
‒ 2
‒ 2
‒ 1
‒ 1
0
3
1
0
‒ 1
‒ 1
‒ 2
1
3
]
yF
2
i
j
3
l
i
x
其中,a 为常数。
试引人支承条件写出总体平衡方程。
解:
先求单元在总体坐标系下的单元刚度矩阵。
单元
(2)在总体坐标系下的单元刚度矩阵与在局部坐标系下的单元刚度矩阵相
同,即
2
0
[𝐾]
(2) = [𝑘]
(2) = 𝑎[ 00
‒ 2
0
0
1
1
0
‒ 1
‒ 1
0
1
1
0
‒ 1
‒ 1
0
0
0
2
0
‒ 2
‒ 2
‒ 1
‒ 1
0
3
1
0
‒ 1
‒ 1
‒ 2
1
3 ]
单元
(1)的坐标转换矩阵中的 β=-90°,单元
(1)的坐标转换矩阵为
[T]=
]
[
cos ( ‒ 90°) sin ( ‒ 90°) 0
‒ sin ( ‒ 90°) cos ( ‒ 90°) 0
0 0 1
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
cos ( ‒ 90°) sin ( ‒ 90°) 0
‒ sin ( ‒ 90°) cos ( ‒ 90°) 0
0 0 1 =
[
0
1
0
0
0
0
‒ 1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
0
0
0
0
‒ 1
0
0
0
0
0
0
0
1
]
所以单元
(1)在总体坐标系下的单元刚度矩阵为:
[𝐾]
(1)=[𝑇]𝑇[𝐾]
(1)[𝑇]
0
‒ 1
0
0
0
= 0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
‒ 1
0
0
0
0
1
0
0
0
0
0
0
0
1 a
[
2
0
0
0
‒ 2
0
0
1
1
0
‒ 1
‒ 1
0
1
1
0
‒ 1
‒ 1
0
0
0
2
0
‒ 2
‒ 2
‒ 1
‒ 1
0
3
1
0
‒ 1
‒ 1
‒ 2
1
3
][
0
1
0
0
0
0
‒ 1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
1
0
0
0
0
‒ 1
0
0
0
0
0
0
0
1
]
=
[
1
0
1
‒ 1
0
‒ 1
0
2
0
2
0
0
1
0
1
‒ 1
0
‒ 1
‒ 1
2
‒ 1
3
0
1
0
0
0
0
2
2
‒ 1
0
‒ 1
1
2
3
]
单元局部编码和总体编码的对应关系为:
单元
(1)ij→12
单元
(2)ij→23
[
单元刚度矩阵中子块对应关系为;[𝐾]
(1)
𝑘22𝑘23
[
]
(2)
𝑘32𝑘33
𝑘11 𝑘12
= 𝑘21 𝑘22
]
(1)
,[𝑘]
(2)
=
所以总体刚度矩阵为:
[K]=
[
1 1
(1)
(1)
(2)
(2)
(2) 2
]
=a
[
1 0 1
0 2 0
1 0 1
‒ 1 1 ‒ 1
0 0 0
‒ 1 0 ‒ 1
0 0 0
0 0 0
0 0 0
‒ 1 0 ‒ 1
2 0 0
‒ 1 0 ‒ 1
5 0 1
0 3 3
1 3 4
0 0 0
‒ 2 ‒ 1 ‒ 1
0 ‒ 1 ‒ 1
0 0 0
0 0 0
0 0 0
0 ‒ 2 0
0 ‒ 1 ‒ 1
0 ‒ 1 ‒ 1
2 0 ‒ 2
0 3 1
‒ 2 1 3
]
节点的位移矢量为:
{𝑢1𝑣1𝜃1𝑢2𝑣2𝜃2𝑢3𝑣3𝜃3}
约束条件为:
𝑢1=0,𝑣1=0,𝜃1=0
作用到结构上的外力为:
𝐹3𝑦=_F
所以引入支承条件的平衡方程为:
𝑇
a
[
1 0 0
0 1 0
0 0 1
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
5 0 1
0 3 3
1 3 4
0 0 0
‒ 2 ‒ 1 ‒ 1
0 ‒ 1 ‒ 1
0 0 0
0 0 0
0 0 0
0 ‒ 2 0
0 ‒ 1 ‒ 1
0 ‒ 1 ‒ 1
2 0 ‒ 2
0 3 1
‒ 2 1 3
𝑢1
𝑣1
𝑣2 0
𝜃3 = 0
0
0
0
0
0
0
_𝐹
7 一组实验数据如下,试用抛物线插值方法计算 X=92 和 X=198 处的 Y 值。
𝑋𝑖
𝑌𝑖
90 100 110 120 130 140 150
0.68 0.74 0.79 0.83 0.86 0.89 0.92
【参考答案】
抛物线插值公式为:
(𝑥 ‒ 𝑥2)(𝑥 ‒ 𝑥3)
Y(x)=(𝑥1 ‒ 𝑥2)(𝑥1 ‒ 𝑥3)𝑦1
+
(𝑥 ‒ 𝑥1)(𝑥 ‒ 𝑥3)
(𝑥2 ‒ 𝑥1)(𝑥2 ‒ 𝑥3) 𝑦2
(𝑥 ‒ 𝑥1)(𝑥 ‒ 𝑥2)
𝑦
当 x=136 时
∵ x∈(130,140),|136-130|>|136-140|
∴选择插值节点:
(𝑥1,𝑦1)=(130,0.86),(𝑥2,𝑦2)=(140,0.89),(𝑥3,𝑦3)=(150,0.92)
将以上数据和 x=136 代入抛物线插值公式,得 x-136 时 y 值为:
(136 ‒ 140)(136 ‒ 150)
(136 ‒ 130)(136 ‒ 150)
Y(136)=
(130 ‒ 140)
(130 ‒ 150) ×0.86+(140 ‒
130)(140 ‒
150)×0.89+
(136 ‒ 130) (136 ‒ 140)
(150 ‒ 130)
(150 ‒ 140)×0.92=0.878
8 如图所示的平面刚架,由两个单元
(1)和
(2)组成,两单元的长度和载面
尺寸及材料特性相同,单元
(1)的局部坐标正方向为沿轴线方向节点 1 指向节
点 2,单元
(2)的局部坐标正方向为沿轴线方向由节点 3 指向节点 1,在局部
坐标系下每个单元的刚度矩阵为
100
0126
64
‒ 100
0‒ 12‒ 6
062
(1) 求刚架总体刚度矩阵[K]。
(2) 引入支撑条件,写出平衡方程。
平面刚架的坐标转换矩阵为。
‒ 1 0 0
0 ‒ 12 6
0 ‒ 6 2
1 0 0
0 12 6
0 6 4
]
𝑇[𝑒]=
[
cos 𝑎 𝑠𝑖𝑛𝑎 0
‒ 𝑠𝑖𝑛𝑎 𝑐𝑜𝑠𝑎 0
0 0 1
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
𝑐𝑜𝑠𝑎 𝑠𝑖𝑛𝑎 0
‒ 𝑠𝑖𝑛𝑎 𝑐𝑜𝑠𝑎 0
0 0
由局部坐标系与总体坐标系的关系知:
单元
(1)a=0,单元
(2)a=
[
所以𝑇
(1)
=[1]
0 1 0
‒ 1 0 0
0 0 1
0 0 0
0 0 0
𝑇
(2)= 0 0 0
0 0 0
0 0 0
0 0 0
0 1 0
‒ 1 0 0
0 0 1
]
在整体坐标系下,单元的刚度矩阵为:
[𝑘]
(𝑒)
= [𝑇(𝑒)][K](𝑒)[𝑇(𝑒)]
]
∴[𝐾]
(1)
=[𝑇
(1) 𝑒
[𝐾]
(1)[𝑇
(1)]-[𝐾]
(1)
[𝐾]
(2)-[𝑇
(2)]𝑇[𝐾]
(2)[𝑇
(2)]-[𝐾]
(2)-
𝐴
2
[
12 0 6
0 1 0
‒ 6 0 4
‒ 12 0 6
0 ‒ 1 0
‒ 6 0 2
12 0 6
0 ‒ 1 0
6 0 2
12 0 ‒ 6
0 1 0
‒ 6 0 4
]
单元
(1)局部码对应的总码为,2,单元
(2)局部码对应的总码为 3.1
9.已知某零件的工作应力和材料强度均服从指数分布,且强度和应力的均值分
别为µ𝑟=210Mpa 和µ𝑠=160Mpa,试确定零件的可靠度。
零件的工作应力和材料强度均服从指数分布,
且µ𝑟=210MPA;µ𝑠 =160Mpa
1 1
∴𝜆𝑠=µ𝑠,𝜆𝑟=µ𝑟
𝜆𝑠
µ𝑟
R=𝜆𝑠 + 𝜆𝑟=µ𝑟 + µ𝑠
210
=
210 + 160=0.5675676
将下列实验测试数据拟合成 y=ax 形式的经验公式。
(计算过程中保留小数点后
该零件的可靠度为:
R=0.5675676
𝑏
两位)
𝑥𝑖
𝑦𝑖
1.18 1.58 2.40 3.00 3.80
2.76 3.29 4.23 4.83 5.57
将 Y=axb 两边取对数,得:
lny=lna+blnx
令 U=lnyA=lnaB=bV=lnx,则原式变为:
U:
A+BV
将表数据取对数:
𝑉𝑖
𝑈𝑖
0.17 0.46 0.88 1.10 1.34
1.02 1.19 1.44 1.57 1.72
按以上的𝑉𝑖,𝑈𝑖进行最小二乘拟合得
{5A+(
5 5
)B=
(
5 5 5
)A+( )B=
5𝐴 + 3. 95𝐵 = 6. 94
代入数据得{3. 95𝐴 + 4. 02𝐵 = 6
a=𝑒 =2.5 b=B=0.6
求解得 A=0.92B=0.60
𝐴
拟合的经验公式为 y=2.5𝑥
0.6
10.如图所示的平面刚架,由两个单元
(1)和
(2)组成,两单元的长度和载
面尺寸及材料特性相同,单元
(1)的局部坐标正方向为沿轴线方向节点 1 指向
节点 2,单元
(2)的局部坐标正方向为沿轴线方向由节点 3 指向节点 1,在局
部坐标系下每个单元的刚度矩阵为.
[𝑘]
(1) = [𝑘]
1 0 0
0 12 6
6 4
0 0
0 ‒ 12 ‒ 6
0 6 2
‒ 1 0 0
0 ‒ 12 6
0 ‒ 6 2
1 0 0
0 12 6
0 6 4
]
(1) 求刚架总体刚度矩阵[K]。
(2) 引入支撑条件,写出平衡方程。
解. 平面刚架的坐标转换矩阵为
𝑇[𝑒]=
[
cos 𝑎 𝑠𝑖𝑛𝑎 0
‒ 𝑠𝑖𝑛𝑎 𝑐𝑜𝑠𝑎 0
0 0 1
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0
𝑐𝑜𝑠𝑎 𝑠𝑖𝑛𝑎 0
‒ 𝑠𝑖𝑛𝑎 𝑐𝑜𝑠𝑎 0
0 0 1
]
由局部坐标系与总体坐标系的关系知:
单元
(1)a=0,单元
(2)a=90°
[
所以
𝑇
(1)=[1]
0 1 0
‒ 1 0 0
0 0 1
0 0 0
0 0 0
𝑇
(2)= 0 0 0
0 0 0
0 0 0
0 0 0
0 1 0
‒ 1 0 0
0 0 1
]
在整体坐标系下,单元的刚度矩阵为:
[𝑘]
(𝑒)
= [𝑇(𝑒)]𝑇[K](𝑒)[𝑇(𝑒)]
]
∴[𝐾]
(1)
=[𝑇
(1) 𝑇
[𝐾]
(1)[𝑇
(1)]=[𝐾]
(1)
[𝐾]
(2)-[𝑇
(2)]𝑇[𝐾]
(2)[𝑇
(2)]=[𝐾]
(2)=
𝐴
2
[
12 0 ‒ 6
0 1 0
‒ 6 0 4
‒ 12 0 6
0 ‒ 1 0
‒ 6 0 2
‒ 12 0 ‒ 6
0 ‒ 1 0
6 0 2
12 0 ‒ 6
0 1 0
‒ 6 0 4
]
单元
(1)局部码对应的总码为,2,单元
(2)局部码对应的总码为 3,1
∴[𝐾]
(1)
𝑘11 𝑘12
(1)
[𝐾]
(2)
𝑘33 𝑘31
(2)
所以按照刚度集成法,可得出总体风度矩阵为:
[𝐾]22
(1)
[K]=
[
[𝐾]11
(1) + [𝐾]11
(2)
[𝐾]21
(1)
[𝐾]31
(2)
[𝐾]12
(1) [𝐾]13
(2)
0
0 [𝐾]33
(2)
]
𝐴
=2
[
13 0 ‒ 6
0 13 6
‒ 6 6 8
‒ 1 0 0
0 ‒ 12 ‒ 6
0 6 2
‒ 12 0 ‒ 6
0 ‒ 1 0
6 0 2
‒ 1 0 0
0 ‒ 12 6
0 ‒ 6 2
1 0 0
0 12 6
0 6 4
0 0 0
0 0 0
0 0 0
‒ 12 0 6
0 ‒ 1 0
‒ 6 0 2
0 0 0
0 0 0
0 0 0
12 0 ‒ 6
0 ‒ 1 0
‒ 6 0 4
]
由于只有在节点 2 处作用有沿 Y 轴负方向的外载荷 F=100N,所以节点载荷矢量
为{F}={0,0,0,0 - 100,0,0,0,0}
𝑇
支撑条件为𝑢3=𝑣3=𝜃3=0,所以等式右端的力矢量无须修改,矩阵[K]中 7 至 9 的
各行各列修改成除主对角线元素为 1 外,其余各元素均为零 1 分
总体平衡方程为
𝐴
2
[
13 0 ‒ 6
0 13 6
‒ 6 6 8
‒ 1 0 0
0 ‒ 12 ‒ 6
0 6 2
0 0 0
0 0 0
0 0 0
𝑢1
‒ 1 0 0 0 0 0 𝑣1
0 ‒ 12 6 0 0 0 𝜃1 0
0 ‒ 6 2 0 0 0 0
1 0 0 0 0 0 𝑢2 0
0 12 6 0 0 0 𝑣2 0
0 6 4 0 0 ‒ 100
0 0 0 1 0 0
3
0 0 0 0 0 1 𝑣3 0
𝜃3 = 0
11.某机电系统由 10 台相同设备组成,各设备可靠度为 0.9,若该系统至少有
7 台设备正常运行就可以保证整个系统正常工作,试求该系统的可靠度。
解:
该系统的每台设备或是正常工作或是发生故障,其失效数为正整数。
因此
是离散型随机变量,且服从二项分布。
∴系统的可靠度由下式计算
∑𝑖 = 0𝐶𝑁𝑅𝑁 ‒ 1
R(r)=
𝑟
F
由题意知:
r=3,N=10,R=0.9,F=1-0.9=0.1,则
10!
R(3)=(10 ‒ 0)!
× 0!
×0.9
(10 ‒ 0)
10!
×0.10+(10 ‒ 1)!
× 1!
×0.9(10 ‒ 1)×0.11
+
10!
10!
(10 ‒ 2)!
× 2!
×0.9 (10 ‒ 2) 0.12+(10 ‒ 3)!
× 3!
×0.9(10 ‒ 3)×0.13≈0.9872
该系统的可靠度为 0.9872
𝑋𝑖
1.20
2.40
3.20
4.54
5.82
𝑌𝑖
7.63
16.41
24.19
40.65
60.41
001
201
12. 已知△ABC= 111 ,将该三角形沿 X 方向移动 1 个单元,沿 Y 方向移动
2 个单位后,再放大一倍,求变换后△ABC 各顶点的坐标。
这是一组合变换,先求出组合变换的变换矩阵。
沿 X 方向移动 1 个单位,沿 Y 方向移动 2 个单位,变换矩阵为
[
100
010
𝑇1= 121
]
放大一倍,变换矩阵为:
2 0 0
0 2 0
𝑇2= 0 0 1
]
所以组合变换矩阵为:
T=
1 0 0 2 0 0 2 0 0
0 1 0 0 2 0 0 2 0
𝑇1𝑇2= 1 2 1 0 0 1 = 2 4 1
]
所以变换后的三角形顶点的坐标短阵为:
ˊBˊCˊ=ABC*T=
[
0 0 1 2 0 0 2 4 1
2 0 1 0 2 0 6 4 1
1 1 1 2 4 1 = 4 6 1
]
变换后的三角形各顶点的坐标为:
A(2,4),B(6,4),C(4,6)
13.用最小二乘法将下列数据拟合成 Y=𝑎0 + 𝑎1𝑥 + 𝑎2𝑥
2
形式的经验公式。
(计算
结果中保留两位小数)
乘法拟合思想得:
{
𝑎0
𝑎0
5𝑎0 + 𝑎1
5
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 现代 设计 方法 算题