function [g,ceg]=Constr(x)
%% 变量统一化
x=ones(1,6);
x1=x(1);x2=x(2);x3=x(3);x4=x(4);x5=x(5);x6=x(6);
%% 常数变量
% x1=100;x2=24;x3=5;x4=320;x5=100;x6=164.21;
u=3; % 传动比
P=20; % 输入功率,kw
n=960; % 小齿轮转速,r/min
K=1.3; % 载荷系数,
Z1=188; % 弹性影响系数，查看机械设计
Z2=2.5; % 区域系数，查看机械设计
sigmaF1=261; % 许用弯曲应力,mpa，参考机械设计
sigmaF2=213; % 许用弯曲应力,mpa，参考机械设计
sigmaH=550; % 许用接触疲劳强度，Mpa
T=9550000*P/n; % 小齿轮传递扭矩,N.m
alpha=20*pi/180; % 压力角
Fn=2*T/(x2*x3*cos(alpha)); % 小齿轮的法向载荷
E=200000; % 弹性模量
J=pi*x5^4/64; % 轴的惯性矩
ymax=Fn*x4^3/(48*E*J); % 轴的最大刚度
%% 计算齿形系数、应力修正系数分别于齿数的关系
%Yfa1=0.169+0.006666*x2-0.0000854*x2^2;
%Yfa2=0.2824+0.0003539*(x2*u)-0.000001576*(u*x2)^2;
Yfa1=26.57*x2^(-1.201)+2.073; % 小齿轮齿形系数与齿数的关系
Yfa2=26.57*(u*x2)^(-1.201)+2.073; % 大齿轮齿形系数与齿数的关系
Ysa1=-1.582*x2^(-0.3366)+2.126; % 小齿轮应力修正系数与齿数的关系
Ysa2=-1.582*(u*x2)^(-0.3366)+2.126; % 大齿轮应力修正系数与齿数的关系
%% 计算啮合角及重合度系数,参考机械设计第九版
alpha1=acos(x2*cos(alpha)/(x2+2*1));
alpha2=acos(u*x2*cos(alpha)/(u*x2+2*1));
epsilong=(x2*(tan(alpha1)-tan(alpha))+u*x2*(tan(alpha2)-tan(alpha)))/(2*pi);
Z=sqrt((4-epsilong)/3); % 接触疲劳强度重合度系数
Y=0.25+0.75/epsilong; % 弯曲疲劳强度的重合度系数
%% 目标函数的所有约束条件
g1=17-x2;
g2=0.7-x1/(x2*x3); % 齿宽系数条件
g3=x1/(x2*x3)-1.4; % 齿宽系数条件
g4=2-x3;
g5=x2*x3-200; % 工艺装备条件
g6=60-x5;
g7=x5-150;
g8=100-x6;
g9=x6-200;
g10=x1+0.2*x6+40-x4; % 轴的支撑跨距条件(线性条件)
g11=Z*Z1*Z2*sqrt(2*K*T*(u+1)/(x1*x2^2*x3^2*u))-sigmaH;
g12=2*K*T*Yfa1*Ysa1*Y/(x1*x3^2*x2)-sigmaF1; % 齿根弯曲疲劳强度条件
g13=2*K*T*Yfa2*Ysa2*Y/(x1*x3^2*x2)-sigmaF2; % 齿根弯曲疲劳强度条件
g14=ymax-0.003*x4;
g=[g1;g2;g3;g4;g5;g6;g7;g8;g9;g10;g11;g12;g13;g14];
ceg=[];
