动力总成优化代码 （.m文件）

动力总成优化代码 ，有matlab优化工具箱对动力总成六自由度进行优化

%以下是解耦优化的目标函数：
function y=fun_2（x）
%%%%%%%%%%%%---动力总成质量---%%%%%%%%%%%
m=2000;%kg
%%%%%%%%%%%%%%---惯性积---%%%%%%%%%%%%
Jx=705.6;Jy=258.7;Jz=842.8;Jxy=0;Jyz=0;Jzx=0;%kg*m*m
%%%%%%%%%%%%---各个悬置点坐标---%%%%%%%%%%
X1=-0.62; Y1=0.7;  Z1=0.33;
X2=-0.62; Y2=0;    Z2=0.33;
X3=-0.62; Y3=-0.7; Z3=0.33;
X4=-0.62; Y4=-0.7; Z4=-0.33;
X5=-0.62; Y5=0;    Z5=-0.33;
X6=-0.62; Y6=0.7;  Z6=-0.33;%mm在下面将位移坐标变成m
%%%%%%%%%%%---各个悬置的三向刚度---%%%%%%%%%
% K1x=3920000; K1y=2764000; K1z=970000;
% K2x=3920000; K2y=2764000; K2z=970000;
% K3x=3920000; K3y=2764000; K3z=970000;
% K4x=3920000; K4y=2764000; K4z=970000;
% K5x=3920000; K5y=2764000; K5z=970000;
% K6x=3920000; K6y=2764000; K6z=970000;%N/m
%%%%%%%%%%%%---各个悬置的安装角度---%%%%%%%%%%
%%%%%%%%%%%%---悬置1安装角度---%%%%%%%%%%%%%%
ALFA1u=0; BETA1u=90*pi/180;GAMA1u=90*pi/180;
ALFA1v=90*pi/180;BETA1v=0;GAMA1v=90*pi/180;
ALFA1w=90*pi/180;BETA1w=90*pi/180;GAMA1w=0;
%%%%%%%%%%%%%---悬置2安装角度---%%%%%%%%%%%%
ALFA2u=0;BETA2u=90*pi/180;GAMA2u=90*pi/180;
ALFA2v=90*pi/180;BETA2v=0;GAMA2v=90*pi/180;
ALFA2w=90*pi/180;BETA2w=90*pi/180;GAMA2w=0;
%%%%%%%%%%%%%---悬置3安装角度---%%%%%%%%%%%%
ALFA3u=0;BETA3u=90*pi/180;GAMA3u=90*pi/180;
ALFA3v=90*pi/180;BETA3v=0;GAMA3v=90*pi/180;
ALFA3w=90*pi/180;BETA3w=90*pi/180;GAMA3w=0;
%%%%%%%%%%%%%---悬置4安装角度---%%%%%%%%%%%%%
ALFA4u=0;BETA4u=90*pi/180;GAMA4u=90*pi/180;
ALFA4v=90*pi/180;BETA4v=0;GAMA4v=90*pi/180;
ALFA4w=90*pi/180;BETA4w=90*pi/180;GAMA4w=0;
%%%%%%%%%%%%%---悬置5安装角度---%%%%%%%%%%%%%
ALFA5u=0;BETA5u=90*pi/180;GAMA5u=90*pi/180;
ALFA5v=90*pi/180;BETA5v=0;GAMA5v=90*pi/180;
ALFA5w=90*pi/180;BETA5w=90*pi/180;GAMA5w=0;
%%%%%%%%%%%%%---悬置6安装角度---%%%%%%%%%%%%%
ALFA6u=0;BETA6u=90*pi/180;GAMA6u=90*pi/180;
ALFA6v=90*pi/180;BETA6v=0;GAMA6v=90*pi/180;
ALFA6w=90*pi/180;BETA6w=90*pi/180;GAMA6w=0;
%%%%%%%%%%%%%---安装角矩阵Tj---%%%%%%%%%%%%%
T1=[cos（ALFA1u）cos（BETA1u）cos（GAMA1u）;cos（ALFA1v）cos（BETA1v）cos（GAMA1v）;cos（ALFA1w）cos（BETA1w）cos（GAMA1w）];
T2=[cos（ALFA2u）cos（BETA2u）cos（GAMA2u）;cos（ALFA2v）cos（BETA2v）cos（GAMA2v）;cos（ALFA2w）cos（BETA2w）cos（GAMA2w）];
T3=[cos（ALFA3u）cos（BETA3u）cos（GAMA3u）;cos（ALFA3v）cos（BETA3v）cos（GAMA3v）;cos（ALFA3w）cos（BETA3w）cos（GAMA3w）];
T4=[cos（ALFA4u）cos（BETA4u）cos（GAMA4u）;cos（ALFA4v）cos（BETA4v）cos（GAMA4v）;cos（ALFA4w）cos（BETA4w）cos（GAMA4w）];
T5=[cos（ALFA5u）cos（BETA5u）cos（GAMA5u）;cos（ALFA5v）cos（BETA5v）cos（GAMA5v）;cos（ALFA5w）cos（BETA5w）cos（GAMA5w）];
T6=[cos（ALFA6u）cos（BETA6u）cos（GAMA6u）;cos（ALFA6v）cos（BETA6v）cos（GAMA6v）;cos（ALFA6w）cos（BETA6w）cos（GAMA6w）];
%%%%%%%%%%%%%%---参数输入完毕---%%%%%%%%%%%%
%%%%%%%%%%%====================%%%%%%%%%%
%%%%%%%%%%%---能量表达式中的Fj和kj---%%%%%%%%%%
B1=[1000Z1-Y1;010-Z10X1;001Y1-X10];
B2=[1000Z2-Y2;010-Z20X2;001Y2-X20];
B3=[1000Z3-Y3;010-Z30X3;001Y3-X30];
B4=[1000Z4-Y4;010-Z40X4;001Y4-X40];
B5=[1000Z5-Y5;010-Z50X5;001Y5-X50];
B6=[1000Z6-Y6;010-Z60X6;001Y6-X60];
%%%此处将位移坐标转变成m
k1=[x（1）00;0x（2）0;00x（3

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       8426  2019-04-13 12:47  youhua.m

----------- ---------  ---------- -----  ----

                 8426                    1