《现代设计方法》实验报告(制造++马龙)精要.doc

重 庆 交 通 大 学 学 生 实 验 报 告 实验课程名称 现代设计方法 开课实验室 计算机中心（南岸校区） 学院 机电与车辆工程学院 年级 2012制造专业一班 学 生 姓 名 马龙 学 号 631224110114 开 课 时 间 2015 至 2016 学年第 一 学期 总 成 绩 教师签名 批改日期 实验项目 机械可靠性设计 实验时间 2015.9.24 实验地点 南岸校区计算机中心 实验性质 (验证性 ?设计性 (综合性 教师评价： 评价教师签名：朱孙科 一：实验目的 利用MATLAB软件，采用中心点法和蒙特卡洛方法求取给定零件的可靠性参数。 二：实验主要内容及过程 （一）实验内容 1、圆截面直杆承受轴向拉力P=100N。设杆的材料的屈服极限fy和直径d为随机变量，其均值和标准差分别为，；，。求杆的可靠指标。为正态分布 ，假设只考虑尺寸系数和表面质量系数对零件疲劳强度的影响，并且假设它们均为正态分布，即，。试采用蒙特卡洛方法确定零件强度的数字特征和分布。 （二）实验过程或计算结果 1 中心点法求杆的可靠指标 1.1 计算结果 （1）轴向拉力P：14000N （2）杆的可靠指标：10.6746 1.2 MATLAB程序 clear; clc; P=14000; muX=[290;30];sigmaX=[25;3]; g=muX(1)-4/pi*P/muX(2)^2; gX=[1;8*P/pi/muX(2)^3]; % bbetaC2=g/norm(gX.*sigmaX) %以应力表示 2 蒙特卡罗方法求零件强度分布 2.1 计算结果 （1）模拟次数N：140000 （2）均值： 328.1687 （3）方差：49.7326 （4）计算时间：0.5148 （5）直方图： 2.2 MATLAB程序 clear; clc; N=input(输入模拟次数N=); t=cputime; %计时起点 P=rand(1,3*N); %产生随机数 x1=norminv(P(1:N),551.4,44.1);%正态逆累积分布函数 x2=norminv(P(N+1:2*N),0.70,0.05);%正态逆累积分布函数 x3=norminv(P(2*N+1:3*N),0.85,0.09);%正态逆累积分布函数 S=x1.*x2.*x3; E=sum(S)/N %求均值 sigma=std(S) %求标准差 sort(S); %排序 hist(S,30) %画直方图 t=cputime-t %计时重点 实验项目 无约束优化问题 实验时间 2015.10.08 实验地点 南岸校区计算机中心 实验性质 (验证性 ?设计性 (综合性 教师评价： 评价教师签名：朱孙科 一：实验目的 利用MATLAB软件求解无约束优化问题 二：实验主要内容及过程 （一）实验内容 1.（一维无约束优化问题）: 采用fminbnd函数，求解一维无约束优化问题f(x)=(x3+cosx+xlogx/ex)，在区间[0,1]中的极小值。 2.（二维无约束优化问题）: 采用fminsearch函数，求解二维无约束优化问题：f(x)=(x14+3x12+x22-2x1-2x2-2x12x2+A)的极小值。 （二）实验过程或计算结果 题目一（一维无约束问题）: 1 函数的极小值：0.3974 2 函数取极小值时的x值：0.5224 3 原函数的曲线图： 4 Matlab程序（需加注释） %求解一维优化问题 clear; clc; fun=inline((x^3+cos(x)+x*log(x))/exp(x),x); %目标函数 x1=0;x2=1; %搜索区间 [xopt,fopt]=fminbnd(fun,x1,x2) %绘制图形 ezplot(fun,[0,10]) title((x^3+cosx+xlogx)/e^x) grid on 题目二（二维无约束问题）: 1 函数的极小值：12.0000 2 函数取极小值时的x1、x2值：1.0000 2.0000 3 Matlab程序（需加注释） %求解二维优化问题 clear; clc; fun=x(1)^4+3*x(1)^2+x(2)^2-2*x(1)-2*x(2)-2*x(1)^2*x(2)+14; x0=[0,0];