Nghiên cứu, thiết kế động cơ cảm ứng startor dây quấn kép

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC QUY NHƠN BÙI ANH TÀI NGHIÊN CỨU, THIẾT KẾ ĐỘNG CƠ CẢM ỨNG STARTOR DÂY QUẤN KÉP LUẬN VĂN THẠC SĨ KỸ THUẬT ĐIỆN Bình Định - Năm 2022 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƢỜNG ĐẠI HỌC QUY NHƠN BÙI ANH TÀI NGHIÊN CỨU, THIẾT KẾ ĐỘNG CƠ CẢM ỨNG STARTOR DÂY QUẤN KÉP Ngành: Kỹ thuật điện Mã số: 8520201 Ngƣời hƣớng dẫn: PGS.TS Đoàn Đức Tùng LỜI CAM ĐOAN Tôi xin cam đoan kết khoa học đƣợc trình bày luận văn thành nghiên cứu thân suốt thời gian thực đề tài chƣa xuất công bố tác giả khác Các kết đạt đƣợc xác trung thực Tác giả luận văn Bùi Anh Tài LỜI CẢM ƠN Để hoàn thành cơng trình này, tơi xin bày tỏ lịng biết ơn sâu sắc kính trọng đến PGS.TS Đồn Đức Tùng trực tiếp hƣớng dẫn, định hƣớng khoa học trình nghiên cứu Thầy dành nhiều thời gian, bảo hỗ trợ nhiều cho suốt trình thực đề tài Tác giả xin trân trọng cảm ơn Lãnh đạo Trƣờng Đại học Quy Nhơn, Phòng Đào tạo Sau Đại học, Khoa Kỹ thuật & Công nghệ tạo điều kiện thuận lợi cho học viên suốt trình học tập nghiên cứu Tôi xin chân thành cảm ơn lãnh đạo Trƣờng hỗ trợ, tạo điều tạo điều kiện cơng tác để tơi có thời gia học tập Tác giả luận văn Bùi Anh Tài MỤC LỤC LỜI CAM ĐOAN LỜI CẢM ƠN MỤC LỤC DANH MỤC TỪ VIẾT TẮT DANH MỤC BẢNG BIỂU DANH MỤC HÌNH ẢNH MỞ ĐẦU CHƢƠNG TỔNG QUAN VỀ MÁY ĐIỆN CẢM ỨNG STATOR DÂY QUẤN KÉP 1.1 Giới thiệu 1.2 Cấu tạo, nguyên lý máy điện cảm ứng stator dây quấn kép 1.2.1 Kết cấu stator máy điện 1.2.2 Kết cấu rotor máy điện 1.2.3 Nguyên lý hoạt động 1.3 Phân loại máy điện cảm ứng stator dây quấn kép 12 1.4 Các chế độ vận hành 13 1.4.1 Chế độ vận hành với nguồn đồng 13 1.4.2 Chế độ vận hành với nguồn không đồng 13 1.5 Kết luận chƣơng 13 CHƢƠNG PHƢƠNG PHÁP THIẾT KẾ MÁY ĐIỆN CẢM ỨNG STATOR DÂY QUẤN KÉP 15 2.1 Tổng quan 15 2.2 Cơ sở lý tuyết thiết kế máy điện cảm ứng stator dây quấn kép 15 2.3 Thiết kế tự động phần khung máy điện cảm ứng stator dây quấn kép 18 2.3.1 Cơ sở tiền đề 18 2.3.2 Nội dung tốn thiết kế tự động động không đồng 18 2.4 Cơ sở liệu bảo đảm công việc tự động thiết kế máy điện không đồng bộ…… 21 2.4.1 Xác định khích thƣớc chủ yếu 21 2.4.2 Kết thiết kế 25 2.5 Thiết kế tính tốn lại thông số máy cảm ứng stator dây quấn kép 27 2.5.1 Cuộn dây Stator 27 2.5.2 Mật độ từ thơng khe hở khơng khí 28 2.5.3 Số vòng dây pha 32 2.5.4 Ƣớc tính cơngsuất 33 2.6 Ƣớc lƣợng thông số máy điện 35 2.6.1 Điện cảm từ hóa L m2 35 2.6.2 Điện cảm rò stator 40 2.6.3 Điện cảm rò rotor 43 2.6.4 Điện trở rotor rb 44 2.6.5 Điện trở vòng ngắn mạch re 44 2.6.6 Điện trở Stator rs 44 2.7 Mơ hình tốn học DSIM 45 2.7.1 Mạch tƣơng đƣơng DSIM d-q 45 2.7.2 Mơ hình Simulink 46 2.8 Kết luận chƣơng 49 CHƢƠNG PHÂN TÍCH VÀ ĐÁNH GIÁ PHƢƠNG PHÁP THIẾT KẾ MÁY ĐIỆN CẢM ỨNG STATOR DÂY QUẤN KÉP 51 3.1 Giới thiệu 51 3.2 Phần mềm ANSYS 52 3.3 Thiết kế động với ANSYS RMxprt Maxwell 54 3.4 Xây dựng mơ hình 2D MAXWELL 55 3.5 Mơ hình FEM cho máy cực 12 cực tiêu chuẩn 56 3.5.1 Mô máy điện dây quấn kép dựa phần mềm Ansys Maxwell 56 3.5.2 Thiết lập dây quấn 59 3.5.3 Đặc tính hiệu suất dƣới dạng đồ thị 60 3.5.4 Mơ hình FEM cho DSIM 64 3.5.5 Kết mô cho chế độ hoạt động khác DSIM .67 3.6 Kết luận chƣơng 71 KẾT LUẬN VÀ KIẾN NGHỊ 73 DANH MỤC TÀI LIỆU THAM KHẢO 74 PHỤ LỤC 76 DANH MỤC TỪ VIẾT TẮT DSP: Bộ xử lý tín hiệu kỹ thuật số DSIG: Máy phát điện cảm ứng cuộn dây Stator kép DSIM: Máy điện cảm ứng Stator dây quấn kép DANH MỤC BẢNG BIỂU Bảng 2.1.Cơ sở liệu hệ số KD theo số đôi cực 21 Bảng 2.2.Cơ sở liệu xác định đƣờng kính Dn theo chiều cao tâm trục 21 Bảng 2.3 Cơ sở liệu xác định hệ số công suất hiệu suất theo công suất động .21 Bảng 2.4.Cơ sở liệu xác định hệ số KE theo đƣờng kính ngồi Stator 22 Bảng 2.5.Cơ sở liệu xác định A B 23 Bảng 2.6.Cơ sở liệu xác định hệ số từ thông tản 24 Bảng 2.7.Cơ sở liệu xác định tích số AJ 24 Bảng 2.8.Sức từ động trung bình rãnh stator 25 Bảng 2.9.Các thơng số stator 26 Bảng 2.10 Các thơng số rotor 27 DANH MỤC HÌNH ẢNH Hình 1.1 Vỏ máy động Hình 1.2 Lõi thép stato động Hình 1.3 Dây quấn stato động Hình 1.4 Rotor trục động Hình 1.5 Nguyên lý DSIM 10 Hình 1.6 Mặt cắt DSIM 10 Hình 2.1 Thuật toán thiết kế máy điện cảm ứng stator dây quấn kép 17 Hình 2.2 Quy trình thiết kế phần khung máy điện KĐB rotor lồng sóc 19 Hình 2.3 Cơ sở liệu cho thiết kế tự động 25 Hình 2.4 Đặc tính làm việc động 26 Hình 2.5 Sơ đồ dây quấn độc lập với số cực khác 28 Hình 2.6 Mật độ từ thơng theo hàm K1 để xác định điều kiện hoạt động máy điện 32 Hình 2.7 Cấu hình đầu cuối cuộn dây 41 Hình 2.8 Hệ tọa độ khơng gian dq 45 Hình 2.9 Mạch d-q tƣơng đƣơng DSIM 46 Hình 3.1 Tích hợp tam giác hình học vật lý 55 Hình 3.2 Các loại máy điện 56 Hình 3.3 Giao diện làm việc 57 Hình 3.4 Thơng số stator 57 Hình 3.5 Thơng số rãnh stator 58 x=['3.3/26: Buoc rang roto: t2=', num2str(t2),'(cm)'];disp(x); %So bo dinh chieu rong rang roto Bz2=Bz1; bpz2=Bdelta*t2/(Bz2*kc); x=['3.4/27: So bo chieu rong rang roto: bpz2=', num2str(bpz2),'(cm)'];disp(x); %Duong kinh truc roto Dt=0.3*D; x=['3.5/28: Duong kinh truc roto: Dt=', num2str(Dt),'(cm)'];disp(x); %Dong dien dan roto load('TrakI.mat'); %trang244 for i=1:7 if(TrakI(8-i,1)cosphi) kImax=TrakI(i,2); cosphimax=TrakI(i,1); break; end end kI=kImin+(kImax-kImin)*(cosphi-cosphimin)/(cosphimax-cosphimin); Itd=kI*I1*6*w1*kd/Z2; I2=Itd; x=['3.6/29: Dong dien dan roto: Itd=I2=', num2str(Itd),'(A)'];disp(x); x=['Voi kI=', num2str(kI),];disp(x); %Dong dien vanh ngan mach Iv=Itd/(2*sin((pi*pp/2)/Z2)); x=['3.7/30: Dong dien vanh ngan mach: Iv=', num2str(Iv),'(A)'];disp(x); J2=3; Std=Itd/J2; x=['Tiet dien dan bang nhom: Std=', num2str(Std),'(mm2)'];disp(x); %Tiet dien vanh ngan mach Jv=2.5; Sv=Iv/Jv; x=['3.8/32: Tiet dien vanh ngan mach: Sv=', num2str(Sv),'(mm2)'];disp(x); %Kich thuoc roto va vanh ngan mach %xem Trang 248 if h>160 b42=1.5; h42=0.5; Isroto=false; for d1r=5.5:0.1:7 if Isroto==true break; end for hr2=20:0.1:40 d2r=d1r; dr=d1r; h12r=hr2-h42-dr; Sr2=pi*dr*dr/4+h12r*dr; %CT/34 bz2=pi*(Dro-2*h42/10-4*(h12r/10+dr/10)/3)/Z2-dr/10; %CT/36 hg2=(Dro-Dt)/2-hr2/10+dr/60; %CT/37 Srss=Sr2-Std; Srss=abs(Srss); %lay tri tuyet doi if Srss=(Z2*2/pp) Rho1=TraRho1(a+1,i); break; end end k41=1-0.033*(b41/10)*(b41/10)/t1/(delta/10); lamdat1=0.9*t1*(q1*kd)*(q1*kd)*Rho1*k41*Sigma1/(kdelta*(delta/10)); x=['5.11/70: He so tu dan tap stato: lamdat1=', num2str(lamdat1),];disp(x); %He so tu dan tan phan dau noi lamdad1=0.34*q1*(ld1-0.64*beta*To)/l1; x=['5.12/71: He so tu dan tan phan dau noi: lamdad1=', num2str(lamdad1),];disp(x); %He so tu dan tan stato Tonglamda1=lamdar1+lamdat1+lamdad1; x=['5.13/72: He so tu dan tan stato: Tonglamda1=', num2str(Tonglamda1),];disp(x); %Dien khang day quan stato x1=0.158*0.5*(w1/100)*(w1/100)*l1*Tonglamda1/(pp/2*q1); x=['5.14/73: Dien khang day quan stato: x1=', num2str(x1),'(Ohm)'];disp(x); xsao1=x1*I1/220; x=[' Tinh theo don vi tuong doi: xsao1=', num2str(xsao1),];disp(x); %He so tu dan tan ranh roto lamdar2=(h12r+dr/2-dr*0.1)*(1-pi*dr*dr/8/Sr2)*(1-pi*dr*dr/8/Sr2)/(3*dr)+0.66-b42/2/dr+h42/b42; x=['5.15/74: He so tu dan tan ranh roto: lamdar2=', num2str(lamdar2),];disp(x); %He so tu dan tap roto if (Z2/(3*pp))=(Z2/(3*pp)) Sigma2=TraSigma2(2,j); break; end end elseif (Z2/(3*pp))>4 Sigma2=9.15/(Z2/(3*pp))/(Z2/(3*pp)); end Sigma2=Sigma2/100; lamdat2=0.9*t2*(Z2/(3*pp))*(Z2/(3*pp))*Sigma2/(kdelta*delta*0.1); x=['5.16/75: He so tu dan tap roto: lamdat2=', num2str(lamdat2),];disp(x); %He so tu tan phan dau noi lamdad2=(2.3*Dv/10)/(Z2*l2*lamda*lamda)*log10((4.7*Dv/10)/(av/10+2*bv/10)); x=['5.17/76: He so tu dan tan phan dau noi: lamdad2=', num2str(lamdad2),];disp(x); %He so tu dan tan ranh nghieng lamdarn=0.5*lamdat2*(t1/t2)*(t1/t2); x=['5.18/77: He so tu tan ranh nghieng: lamdarn=', num2str(lamdarn),];disp(x); %He so tu tan roto Tonglamda2=lamdar2+lamdat2+lamdad2+lamdarn; x=['5.19/78: He so tu dan tan roto: Tonglamda2=', num2str(Tonglamda2),];disp(x); %Dien khang tan day quan roto x2=7.9*50*l2*Tonglamda2*0.00000001; x=['5.20/79: Dien khang day quan roto: x2=', num2str(x2),'(Ohm)'];disp(x); %Dien khang roto da quy doi xp2=x2*gamma; xpsao2=xp2*I1/220; x=['5.21/80: Dien khang roto da quy doi: xp2=', num2str(xp2),'(Ohm)'];disp(x); x=[' Tinh theo don vi tuong doi: xpsao2=', num2str(xpsao2),];disp(x); %Dien khang ho cam x12=(220-Imu*x1)/Imu; xsao12=x12*I1/220; x=['5.22/81: Dien khang ho cam: x12=', num2str(x12),'(Ohm)'];disp(x); x=[' Tinh theo don vi tuong doi: xsao12=', num2str(xsao12),];disp(x); %Tinh lai kE kEtt=(220-Imu*x1)/220; x=['5.23/82: Tinh lai kE: kEtt=', num2str(kEtt),];disp(x); disp(' Tri so khong sai khac nhieu so voi gia thiet ban dau'); disp(' -'); disp(' CHUONG 6: TON HAO THEP VA TON HAO CO'); disp(' ******************'); %Trong luong rang stato hpz1=hr1/10-d2s/30; Gz1=7.8*Z1*bz1*hpz1*l1*kc*0.001; x=['6.1/83: Trong luong rang stato: Gz1=', num2str(Gz1),'(kg)'];disp(x); %Trong luong gong tu stato Gg1=7.8*l1*Lg1*hg1*pp*kc*0.001; x=['6.2/84: Trong luong gong tu stato: Gg1=', num2str(Gg1),'(kg)'];disp(x); %Ton hao sat loi sat stato %He so gia cong kgc tr.140 PFez1=1.8*2.5*Bz1tt*Bz1tt*Gz1*0.001; PFeg1=1.6*2.5*Bg1tt*Bg1tt*Gg1*0.001; PpFe=PFez1+PFeg1; disp('6.3/85: Ton hao sat loi sat stato:'); x=['Trong rang: PFez1=', num2str(PFez1),'(kW)'];disp(x); x=['Trong gong: PFeg1=', num2str(PFeg1),'(kW)'];disp(x); x=['Trong ca loi sat stato: PpFe=', num2str(PpFe),'(kW)'];disp(x); %Ton hao be mat tren rang roto load('Trabeta0'); for i=1:9 if Trabeta0(i,1)=b41/delta bdmin=Trabeta0(i,1); bdmax=Trabeta0(i+1,1); beta0min=Trabeta0(i,2); beta0max=Trabeta0(i+1,2); end end beta0=(b41/delta-bdmin)*(beta0max-beta0min)/(bdmax-bdmin)+beta0min; B0=beta0*kdelta*Bdelta; pbm=0.5*2*((Z1*n/10000)^1.5)*((10*B0*t1)^2); Pbm=pp*To*((t2-b42/10)/t2)*l2*pbm*(10^-7); x=['6.4/86: Ton hao be mat tren rang roto: Pbm=', num2str(Pbm),'(kW)'];disp(x); %Ton hao dap mach tren rang roto Bdap=(v1*delta/10)*Bz2tt/(2*t2); Gz2=7.8*Z2*(hr2/10-dr/30)*bz2*l2*kc*0.001; Pdm=0.11*((Z1*n/1000*Bdap)^2)*Gz2*0.001; x=['6.5/87: Ton hao dap mach tren rang roto: Pdm=', num2str(Pdm),'(kW)'];disp(x); %Tong ton hao thep PFe=PpFe+Pbm+Pdm; x=['6.6/88: Tong ton hao thep: PFe=', num2str(PFe),'(kW)'];disp(x); %Ton hao co Pco=((n/1000)^2)*((Dn/10)^4)*0.001; x=['6.7/89: Ton hao co: Pco=', num2str(Pco),'(kW)'];disp(x); %Ton hao khong tai P0=PFe+Pco; x=['6.8/90: Ton hao khong tai: P0=', num2str(P0),'(kW)'];disp(x); disp(' '); disp(' CHUONG 7: TINH TOAN DAC TINH LAM VIEC'); disp(' ******************'); disp('Cac thong so tinh toan:'); x=[' r1=', num2str(r1),'(Ohm)'];disp(x); x=[' x1=', num2str(x1),'(Ohm)'];disp(x); x=[' rp2=', num2str(rp2),'(Ohm)'];disp(x); x=[' xp2=', num2str(xp2),'(Ohm)'];disp(x); x=[' x12=', num2str(x12),'(Ohm)'];disp(x); C1=1+x1/x12; Idbx=Imu; Idbr=(PFe*1000+3*Imu*Imu*r1)/(3*220); E1=220-Imu*x1; kItt=6*w1*kd/Z2; Ip2=Itd/kItt; sdm=Ip2*rp2/E1; sm=rp2/(x1/C1+xp2); x=[' C1=', num2str(C1),];disp(x); x=[' Idbx=', num2str(Idbx),'(A)'];disp(x); x=[' Idbr=', num2str(Idbr),'(A)'];disp(x); x=[' E1=', num2str(E1),'(V)'];disp(x); x=[' Ip2=', num2str(Ip2),'(A)'];disp(x); x=['He so truot dinh muc: sdm=', num2str(sdm),];disp(x); x=['He so truot max: sm=', num2str(sm),];disp(x); %Hien thi bang dac tinh lam viec sh=[0.005 0.01 sdm 0.06 sm]; dactinhlamviec=zeros(18,5); for i=1:5 s=sh(1,i); dactinhlamviec(1,i)=C1*C1*((r1/C1) + (rp2/s)); %rns dactinhlamviec(2,i)=C1*C1*(x1/C1 + xp2); %xns dactinhlamviec(3,i)=sqrt(dactinhlamviec(1,i)*dactinhlamviec(1,i) + dactinhlamviec(2,i)*dactinhlamviec(2,i));%Zns dactinhlamviec(4,i)=C1*220/dactinhlamviec(3,i);%I'2 dactinhlamviec(5,i)=dactinhlamviec(1,i)/dactinhlamviec(3,i);%cosfi'2 dactinhlamviec(6,i)=dactinhlamviec(2,i)/dactinhlamviec(3,i);%sinfi'2 dactinhlamviec(7,i)=Idbr + dactinhlamviec(4,i)*dactinhlamviec(5,i)/C1;%I1r dactinhlamviec(8,i)=Idbx + dactinhlamviec(4,i)*dactinhlamviec(6,i)/C1;%I1x dactinhlamviec(9,i)=sqrt(dactinhlamviec(7,i)*dactinhlamviec(7,i) + dactinhlamviec(8,i)*dactinhlamviec(8,i));%I1 dactinhlamviec(10,i)=dactinhlamviec(7,i)/dactinhlamviec(9,i);%cosfi dactinhlamviec(11,i)=3*220*dactinhlamviec(7,i)*10.^-3;%P1 dactinhlamviec(12,i)=3*dactinhlamviec(9,i)*dactinhlamviec(9,i)*r1*10.^-3;%PCu1 dactinhlamviec(13,i)=3*dactinhlamviec(4,i)*dactinhlamviec(4,i)*rp2*10.^-3;%Pcu2 dactinhlamviec(14,i)=0.005*dactinhlamviec(11,i);%Pf dactinhlamviec(15,i)=P0; %P0 dactinhlamviec(16,i)=dactinhlamviec(12,i) + dactinhlamviec(13,i) + dactinhlamviec(14,i) + dactinhlamviec(15,i);%Tong P dactinhlamviec(17,i)=dactinhlamviec(11,i)-dactinhlamviec(16,i);%P2 dactinhlamviec(18,i)=dactinhlamviec(17,i)*100/dactinhlamviec(11,i);%hieu suat end cnames = ['0.005 0.01 ',num2str(sdm),' 0.06 ' num2str(sm),]; rnames = ['rns(ohm) xns(ohm) Zns(ohm) Ip2(A) cosphip2 sinphip2 I1r(A) I1x(A) I11(A) cosphitt P1(kW) Pcu1(kW) Pcu2(kW) Pf(kW) P0(kW) TongP(kW) P2(kW) hstt']; printmat(dactinhlamviec,'Bang dac tinh lam viec cua DC dien KDB roto long soc:',rnames,cnames); %Tinh toan sdm rnsdm=(C1^2)*(r1/C1+rp2/sdm); xnsdm=C1*C1*(x1/C1+xp2); Znsdm=sqrt(rnsdm^2+xnsdm^2); Ip2dm=C1*220/Znsdm; cosphip2dm=rnsdm/Znsdm; sinphip2dm=xnsdm/Znsdm; I1rdm=Idbr+Ip2dm*cosphip2dm/C1; I1xdm=Idbx+Ip2dm*sinphip2dm/C1; I11dm=sqrt(I1rdm^2+I1xdm^2); cosphittdm=I1rdm/I11dm; P1dm=3*220*I1rdm*0.001; Pcu1dm=3*I11dm*I11dm*r1*0.001; Pcu2dm=3*Ip2dm*Ip2dm*rp2*0.001; Pfdm=0.005*P1dm; TongPdm=Pcu1dm+Pcu2dm+Pfdm+P0; P2dm=P1dm-TongPdm; hsttdm=P2dm*100/P1dm; %Tinh toan sm rnsmax=(C1^2)*(r1/C1+rp2/sm); xnsmax=C1*C1*(x1/C1+xp2); Znsmax=sqrt(rnsmax^2+xnsmax^2); Ip2max=C1*220/Znsmax; disp('Cac thong so lam viec:'); x=[' I11dm=', num2str(I11dm),'(A)'];disp(x); x=[' Ip2dm=', num2str(Ip2dm),'(A)'];disp(x); x=[' P2dm=', num2str(P2dm),'(kW)'];disp(x); cosphittdm=cosphittdm*1000;cosphittdm=round(cosphittdm);cosphittdm=cosphittdm/1000; x=[' cosphittdm=', num2str(cosphittdm),];disp(x); hsttdm=hsttdm*100;hsttdm=round(hsttdm);hsttdm=hsttdm/100; x=[' hsttdm=', num2str(hsttdm),'(%)'];disp(x); x=[' Ip2max=', num2str(Ip2max),'(A)'];disp(x); %Boi so momen cuc dai mmax=((Ip2max/Ip2dm)^2)*(sdm/sm); x=['91/ Boi so momen cuc dai: mmax=', num2str(mmax),];disp(x); load('Trammax3K.mat'); a=(pp/2)+1; for i=1:18 if Trammax3K(i,1)==P mmax3K=Trammax3K(i,a); break; end end x=['- Boi so momen cuc dai cua day may 3K: mmax3K=', num2str(mmax3K),];disp(x); if mmax>=mmax3K disp('- Boi so momen cuc dai dat tieu chuan thiet ke'); elseif mmax0 tong=tong+1; end; i=tong; rnsao=(C1^2)*(r1/C1+rp2/sao); xnsao=C1*C1*(x1/C1+xp2); Znsao=sqrt(rnsao^2+xnsao^2); Ip2ao=C1*220/Znsao; cosphip2ao=rnsao/Znsao; sinphip2ao=xnsao/Znsao; I1rao=Idbr+Ip2ao*cosphip2ao/C1; I1xao=Idbx+Ip2ao*sinphip2ao/C1; I11ao=sqrt(I1rao^2+I1xao^2); cosphittao=I1rao/I11ao; P1ao=3*220*I1rao*0.001; Pcu1ao=3*I11ao*I11ao*r1*0.001; Pcu2ao=3*Ip2ao*Ip2ao*rp2*0.001; Pfao=0.005*P1ao; TongPao=Pcu1ao+Pcu2ao+Pfao+P0; P2ao=P1ao-TongPao; hsttao=P2ao*100/P1ao; if hsttao