matlab多目标优化计算.ppt

function y=barelefor(k,u) % 计算单元节点力矢量 y=k*u%u 位移 function y=barelestr(k,u,a) % 计算单元应力矢量 y=k*u/a 杆单元 E=210GPa,A=0.003m^2,P=10kN, 右端点位 移 0.002m,L12=1.5m,L23=1m 1 2 P 3 杆单元 e=210e6 ； a=0.003 ； l1=1.5 ； l2=1 ； k1=bareal(e,a,l1) % 单元刚度矩阵 k2=bareal(e,a,l2) 杆单元 k1 = 420000 -420000 -420000 420000 k2 = 630000 -630000 -630000 630000 kk=zeros(3,3); % 整体 kk=barass(kk,k1,1,2) kk=barass(kk,k2,2,3) 杆单元 kk = 420000 -420000 0 -420000 420000 0 0 0 0 kk = 420000 -420000 0 -420000 1050000 -630000 0 -630000 630000 P=10kN, 节点 3 位移 0.002m 0 u2 0.002 f1 -10 f3 k=kk(2,2);k0=kk(2,3); u0=0.0020;f=[-10];% 边界条件 f0=f-k0*u0; u=k\f0% 节点 2 位移 uu=[0;u;u0]; ff=kk*uu% 支反力 杆单元 u = 0.0012 ff = -500.0000 -10.0000 510.0000 u1=[0;uu(2)] f1=barelefor(k1,u1) s1=barelestr(k1,u1,a) u2=[uu(2);uu(3)] f2=barelefor(k2,u2) s2=barelestr(k2,u2,a) f1 = -500.0000 500.0000 s1 = 1.0e+005 * -1.6667 1.6667 f2 = -510.0000 510.0000 s2 = 1.0e+005 * -1.7000 1.7000 杆单元 6 多目标优化问题 多目标优化问题的 MATLAB 函数有 fgoalattain 需确定各分目标的加权系数 需知各分目标的单个的最优值 目标函数的最大值逐次减小 fminimax 6.1 函数 fgoalattain min v s.t. f i (X)- w i v ≤ goal i i =1,2,…,t AX≤b （线性不等式约束） AeqX=beq （线性等式约束） C(X)≤0 （非线性不等式约束条件） Ceq(X)=0 （非线性等式约束） Lb ≤X ≤Ub （边界约束条件） 一、多目标优化问题数学模型 标量变量 各分目标函数 分目标函数的权重 各分目标函数的目标值 二、优化函数使用格式 [x,fval,exitflag,output, grad,hessian]= fgoalattain(@fun,x0,goal,w,A,b,Aeq,beq,Lb,Ub,Nlc,options,P1,P2…) 设置优化选项参数 各分目标权重 各分目标期望值 目标函数在最优解的海色矩阵 返回目标函数在最优解的梯度 优化算法信息的一个数据结构 返回算法的终止标志 返回目标函数的最优值 返回目标函数的最优解 附加参数 非线性约束条件的函数名 设计变量的下界和上界 线性等式约束的常数向量 线性等式约束的系数矩阵 线性不等式约束的常数向量 线性不等式约束的系数矩阵 无定义时以空矩阵 符号“ [ ]” 代替 6.1 函数 fgoalattain 初始点 目标函数文件名 三、例题 设计带式输送机传动装置上的普通 V 带传动。已知电机额定功率 P=4kW ，转 速 n 1 =1440r/min ，传动比 i=3 ，采用 A 型 V 带，每天工作不超过 10 小时。要求 传动结构紧凑（带的根数尽量少，带轮直径和中心距尽量小）。 解 :(1) 建立优化设计的数学模型 ①设计变量： V 带传动的独立设计变量是小带轮直径 d d1 和带的基准长度 L d X=[d d1 ,L d ] T =[x 1 ,x 2 ] T ②目标函数： ? 小带轮直径： ? 中心距 : ? 带的根数 : L 0 0 A 3 K K ) P P ( P K z ) X ( f