Nghiên cứu thiết kế động cơ không đồng bộ 3 pha Rôto lồng sóc rãnh sâu có tính đến hiệu ứng mặt ngoài

BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THỦY LỢI LÊ QUỐC DUY NGHIÊN CỨU THIẾT KẾ ĐỘNG CƠ KHÔNG ĐỒNG BỘ PHA RƠTO LỒNG SĨC RÃNH SÂU CĨ TÍNH ĐẾN HIỆU ỨNG MẶT NGỒI TRONG THANH DẪN RƠTO LUẬN VĂN THẠC SĨ NINH THUẬN, NĂM 2017 BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THỦY LỢI LÊ QUỐC DUY NGHIÊN CỨU THIẾT KẾ ĐỘNG CƠ KHÔNG ĐỒNG BỘ PHA RƠTO LỒNG SĨC RÃNH SÂU CĨ TÍNH ĐẾN HIỆU ỨNG MẶT NGỒI TRONG THANH DẪN RƠTO Chuyên ngành: KỸ THUẬT ĐIỆN Mã số: 60520202 NGƯỜI HƯỚNG DẪN KHOA HỌC TS: Lê Quang Cường NINH THUẬN, NĂM 2017 LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu thân Các kết nghiên cứu kết luận luận văn trung thực, không chép từ nguồn hình thức Việc tham khảo nguồn tài liệu thực trích dẫn ghi nguồn tài liệu tham khảo quy định Tác giả luận văn Lê Quốc Duy i LỜI CÁM ƠN Để có kết hồn thành tốt thời gian học tập Trường Đại Học Thủy lợi em xin chân thành cảm ơn đến tất quý thầy cô khoa Năng Lượng Trường Đại Học Thủy lợi Trên thực tế khơng có thành cơng mà không gắn liền với hỗ trợ, giúp đỡ dù hay nhiều, dù trực tiếp hay gián tiếp Trong suốt khoảng thời gian học tập tai trường em nhận nhiều quan tâm, giúp đỡ q thầy, gia đình bạn bè Với lòng biết ơn sâu sắc nhất, em xin gửi đến q thầy, khoa Năng Lượng, Trường Đại học Thủy Lợi Trong suốt thời gian học tập truyền đạt kiến thức thiết thực cần thiết ngành nghề điện mà em theo học, với hướng dẫn dạy bảo tận tình thầy suốt tháng làm luận văn Đặc biệt em xin gửi lời cảm ơn đến GVHD thầy TS Lê Quang Cường giúp đỡ em nhiều trình học tập luận văn em Do không tránh khỏi thiếu sót điều chắn em mong nhận ngững ý kiến đóng góp quý báo quý thầy cô bạn học lớp để kiến thức em lĩnh vực hoàn thiện Em gởi lời cảm ơn chân thành lời tri ân sâu sắc đến quý thầy cô khoa Năng Lượng, Trường Đại học Thủy Lợi tạo điều kiện cho em hồn thành khóa học, trình độ lý luận kinh nghiệm thực tiễn hạn chế nên khơng tránh khỏi thiếu sót em mong nhận gớp ý quý thầy cô để em học thêm nhiều kinh nghiệm hoàn thành tốt Em xin chân thành cám ơn ii MỤC LỤC DANH MỤC CÁC HÌNH ẢNH v MỞ ĐẦU vii CHƯƠNG 1.1 ĐẠI CƯƠNG VỀ MÁY ĐIỆN KHÔNG ĐỒNG BỘ Kết cấu, đại lượng định mức công dụng .1 1.1.1 Lõi sắt: phần dẫn từ 1.1.2 Dây quấn: 1.1.3 Vỏ máy: gồm thân máy, nắp máy chân đế 1.1.4 Lõi thép 1.1.5 Dây quấn rôto: 1.2 Quan hệ điện từ máy điện không đồng 1.2.1 Máy điện không đồng làm việc rôto đứng yên 1.2.2 Máy điện không đồng làm việc rôto quay .10 1.3 Các chế độ làm việc, giản đồ lượng, đồ thị vectơ máy điện không đồng 15 1.4 Các đặc tính máy điện khơng đồng 18 1.4.1 Các đặc tính làm việc động điện không đồng 18 1.4.2 Đặc tính moment 19 1.4.3 Tổn hao đặc tính hiệu suất động h = f (P2) 19 1.4.4 Đặc tính hệ số cơng suất cosφ = f (P2) 19 1.4.5 Năng lực tải 20 1.5 Hiệu ứng mặt ngồi dẫn rơto lịng sóc 20 CHƯƠNG THIẾT KẾ MÁY ĐIỆN KHÔNG ĐỒNG BỘ 21 2.1 Những vấn đề chung thiết kế động không đồng 21 2.2 Tính tốn thiết kế máy điện khơng đồng 27 2.2.1 Xác định kích thước chủ yếu 27 2.2.2 Thiết kế stato 29 2.2.3 Thiết kế rôto 31 2.2.4 Tính tốn điện từ 34 2.2.5 Tham số động định mức 44 2.2.6 Tính tốn tổn hao 49 iii 2.2.7 Tính tốn đặc tính làm việc 51 2.2.8 Tính tốn đặc tính khởi động 54 CHƯƠNG KHẢO SÁT ĐẶC TÍNH ĐỘNG CỦA ĐỘNG CƠ KHƠNG ĐỒNG BỘ PHA RƠTO LỒNG SĨC RÃNH SÂU KHI TÍNH ĐẾN HIỆU ỨNG MẶT NGỒI TRONG THANH DẪN RÔTO 59 3.1 Xây dựng module tính tốn ảnh hưởng hiệu ứng tần số đến điện trở điện kháng rơto lồng sóc rãnh sâu 59 3.1.1 Sự thay đổi tham số tượng hiệu ứng mặt ngồi dịng điện .59 3.1.2 Phương pháp tính toán ảnh hưởng hiệu ứng tần số dẫn rôto 63 3.2 Mô khảo sát đặc tính động động khơng đồng tổng trở rơto thay đổi q trình làm việc .69 3.2.1 Công nghê xây dựng S-Funtion SIMULINK 69 3.2.2 Mơ hình tốn học 72 3.2.3 Xây dựng mơ hình DKB Simulink sử dụng hàm S-function lập trình cho DKB .78 KẾT LUẬN VÀ KIẾN NGHỊ VỀ NHỮNG NGHIÊN CỨU TIẾP THEO 83 TÀI LIỆU THAM KHẢO 84 PHỤ LỤC 85 iv DANH MỤC CÁC HÌNH ẢNH Hình 1.1 Lá thép kĩ thuật điện Hình 1.2 Dây quấn stator Hình 1.3 Stator máy điện không đồng Hình 1.4 Lá thép rơto kĩ thuật điện không đồng Hình 1.5 Rơto (a) sơ đồ mạch điện (b) rôto dây quấn .3 Hình 1.6 Dây quấn rơto lồng sóc (a) rơto lồng sóc rãnh chéo (b) Hình 1.7 Đồ thị véc tơ máy điện điện không đồng ngắn mạch 10 Hình 1.8 Mạch điện thay máy không đồng rôto đứng yên 10 Hình 1.9 Giản đồ thay hình T 14 Hình 1.10 Giản đồ lượng động điện 15 Hình 1.11 Giản đồ lượng máy phát điện không đồng 17 Hình 1.12 Đồ thị vec tơ máy điện không đồng 17 Hình 1.13 Đồ thị vectơ giản đồ lượng máy điện không đồng chế độ hãm điện từ .18 Hình 2.1 Dạng rãnh rơto loại thường 33 Hình 2.2 Kích thước rãnh cách điện 37 Hình 2.3 Kích thước rãnh rơto vịng ngắn mạch 40 Hình 2.4 Đồ thị biểu thị đặc tính làm việc động khơng đồng rơto lồng sóc 53 Hình 3.1 kích thước rãnh rơto dung để xác định chiều sâu qui đổi hx, hr 61 Hình 3.2 Đường cong  , '  f ( ) 61 Hình 3.3 Sơ đồ thay cuộn dây rôto chia nhiều lớp .63 Hình 3.4 Bố trí khơng gian tổng quát 72 Hình 3.5 Mơ hình hệ thống xây dựng 78 v DANH MỤC BẢNG BIỂU Bảng 2.1 Trị số KD 27 Bảng 2.2 Hiệu suất cosϕ dãy động điện không đồng 3K 28 Bảng 2.3 Xác định bước stato .31 Bảng 2.4 Số liệu đặc tính làm việc .52 vi MỞ ĐẦU Tính cấp thiết Đề tài: Trong loại động động khơng đồng có kết cấu đơn giản nhất, đặc biệt động rơto lồng sóc Do đơn giản nên giá rẻ, làm việc tin cậy, thao tác vận hành thuận tiện Ngày nay, với phát triển biến tần nên việc điều khiển Động không đồng thực dễ ràng Do nguyên nhân động không đồng sử dụng rộng rãi Động không đồng bao gồm động rôto dây quấn động rơto lồng sóc Các động khơng đồng rơto dây quấn rơto lồng sóc đơn giản khơng sử dụng rộng rãi mơ men mở máy nhỏ dòng điện mở máy lại lớn Ngồi loại rơto dây quấn phải có vành trượt, cấu tạo phức tạp, địi hỏi phải có điện trở mở máy nên làm việc đảm bảo, điều khiển mở máy phức tạp Trên thực tế ngưởi ta sử dụng rộng rãi động rôto lồng sóc rãnh kép rãnh sâu để cải thiện đặc tính mở máy Khi thiết kế tính tốn động khơng đồng bộ, việc khảo sát tính tốn đặc tính động thường nghiên cứu phương pháp giải tích, tính tốn đặc tính động dự sơ đồ thay thế, khảo sát đặt tính theo độ trượt Trong q trình tính tốn ta coi điện kháng tản stato rơto, từ hóa khơng phụ thuộc vào bão hịa lõi thép, cịn tổn thất từ thơng tản sóng điều hịa bậc cao khơng, bỏ qua phụ thuộc tổng trở rôto vào hệ số trượt Thực chất thông số động thay đổi tượng bão hòa mạch từ, hiệu ứng tần số rôto Với phát triển công nghệ thông tin, vận dụng để tính tốn khảo sát đặc tính động động có tính đến thay đổi tham số động theo thời gian, nhờ tính tốn tiêu động động bước thiết kế Chính em xin chọn đề tài “Nghiên cứu thiết kế động khơng đồng pha rơto lồng sóc rãnh sâu có đến tính đến hiệu ứng mặt ngồi dẫn rơto” vii Mục đích Đề tài: Thiết kế động khơng đồng rơto lồng sóc rãnh sâu nhằm cải thiện đặc tính khởi động, xây dựng module tính tốn mơ hình khảo sát đặc tính động động tham số động thay đổi q trình làm việc, từ tìm tiêu động động trình thiết kế Đối tượng phạm vi nghiên cứu Động không đồng pha rơto lồng sóc rãnh sâu Cách tiếp cận phương pháp nghiên cứu: Nghiên cứu lý thuyết, tính tốn thiết kế, lập trình giải tốn máy tính, mơ thực nghiệm mơ hình mơ để khảo sát đặc tính động động không đồng Cấu trúc luận văn Luận văn gồm chương, bảng 22 hình trình bày 82 trang với nội dung đại cương máy điện, thiết kế máy điện không đồng Xây dựng module tính tốn ảnh hưởng hiệu ứng tần số đến điện trở điện kháng rơto lồng sóc rãnh sâu Mơ khảo sát đặc tính động ĐCKĐB tổng trở rơto thay đổi q trình làm việc kết luận luận văn viii #if defined(MDL_SET_INPUT_PORT_WIDTH) && defined(MATLAB_MEX_FILE) static void mdlSetInputPortWidth(SimStruct *S, int portIndex, int width) { } /* end mdlSetInputPortWidth */ #endif /* MDL_SET_INPUT_PORT_WIDTH */s #define MDL_SET_OUTPUT_PORT_WIDTH /* Change to #undef to remove function */ #if defined(MDL_SET_OUTPUT_PORT_WIDTH) && defined(MATLAB_MEX_FILE) static void mdlSetOutputPortWidth(SimStruct *S, int portIndex, int width) { } /* end mdlSetOutputPortWidth */ #endif /* MDL_SET_OUTPUT_PORT_WIDTH */ #undef MDL_SET_INPUT_PORT_DIMENSION_INFO /* Change to #define to add function */ #if defined(MDL_SET_INPUT_PORT_DIMENSION_INFO) && defined(MATLAB_MEX_FILE) static void mdlSetInputPortDimensionInfo(SimStruct *S, int_T portIndex, const DimsInfo_T *dimsInfo) 90 { } /* mdlSetInputPortDimensionInfo */ #endif /* MDL_SET_INPUT_PORT_DIMENSION_INFO */ #undef MDL_SET_OUTPUT_PORT_DIMENSION_INFO /* Change to #define to add function*/ #if defined(MDL_SET_OUTPUT_PORT_DIMENSION_INFO) && defined(MATLAB_MEX_FILE) static void mdlSetOutputPortDimensionInfo(SimStruct *S, int_T portIndex, const DimsInfo_T *dimsInfo) { } /* mdlSetOutputPortDimensionInfo */ #endif /* MDL_SET_OUTPUT_PORT_DIMENSION_INFO */ #undef MDL_SET_DEFAULT_PORT_DIMENSION_INFO /* Change to #define to add fcn */ #if defined(MDL_SET_DEFAULT_PORT_DIMENSION_INFO) defined(MATLAB_MEX_FILE) static void mdlSetDefaultPortDimensionInfo(SimStruct *S) { } /* mdlSetDefaultPortDimensionInfo */ #endif /* MDL_SET_DEFAULT_PORT_DIMENSION_INFO */ #define MDL_SET_INPUT_PORT_SAMPLE_TIME 91 && #if defined(MDL_SET_INPUT_PORT_SAMPLE_TIME) && defined(MATLAB_MEX_FILE) static void mdlSetInputPortSampleTime(SimStruct *S, int_T portIdx, real_T sampleTime, real_T offsetTime) { } /* end mdlSetInputPortSampleTime */ #endif /* MDL_SET_INPUT_PORT_SAMPLE_TIME */ #define MDL_SET_OUTPUT_PORT_SAMPLE_TIME #if defined(MDL_SET_OUTPUT_PORT_SAMPLE_TIME) defined(MATLAB_MEX_FILE) static void mdlSetOutputPortSampleTime(SimStruct *S, int_T portIdx, real_T sampleTime, real_T offsetTime) { } /* end mdlSetOutputPortSampleTime */ #endif /* MDL_SET_OUTPUT_PORT_SAMPLE_TIME */ static void mdlInitializeSampleTimes(SimStruct *S) { /* Register one pair for each sample time */ 92 && ssSetSampleTime(S, 0, CONTINUOUS_SAMPLE_TIME); ssSetOffsetTime(S, 0, 0.0); } /* end mdlInitializeSampleTimes */ #define MDL_SET_INPUT_PORT_DATA_TYPE /* Change to #undef to remove function */ #if defined(MDL_SET_INPUT_PORT_DATA_TYPE) && defined(MATLAB_MEX_FILE) static void mdlSetInputPortDataType(SimStruct *S, int portIndex,DTypeId dType) { } /* mdlSetInputPortDataType */ #endif /* MDL_SET_INPUT_PORT_DATA_TYPE */ #define MDL_SET_OUTPUT_PORT_DATA_TYPE /* Change to #undef to remove function */ #if defined(MDL_SET_OUTPUT_PORT_DATA_TYPE) && defined(MATLAB_MEX_FILE) static void mdlSetOutputPortDataType(SimStruct *S,int portIndex,DTypeId dType) { } /* mdlSetOutputPortDataType */ #endif /* MDL_SET_OUTPUT_PORT_DATA_TYPE */ #define MDL_SET_DEFAULT_PORT_DATA_TYPES #undef to remove function*/ 93 /* Change to #if defined(MDL_SET_DEFAULT_PORT_DATA_TYPES) && defined(MATLAB_MEX_FILE) static void mdlSetDefaultPortDataTypes(SimStruct *S) { } /* mdlSetDefaultPortDataTypes */ #endif /* MDL_SET_DEFAULT_PORT_DATA_TYPES */ #define MDL_SET_INPUT_PORT_COMPLEX_SIGNAL /* Change to #undef to remove */ #if defined(MDL_SET_INPUT_PORT_COMPLEX_SIGNAL) && defined(MATLAB_MEX_FILE) static void mdlSetInputPortComplexSignal(SimStruct *S, int portIndex, CSignal_T cSignalSetting) { } /* mdlSetInputPortComplexSignal */ #endif /* MDL_SET_INPUT_PORT_COMPLEX_SIGNAL */ #define MDL_SET_OUTPUT_PORT_COMPLEX_SIGNAL /* Change to #undef to remove */ #if defined(MDL_SET_OUTPUT_PORT_COMPLEX_SIGNAL) defined(MATLAB_MEX_FILE) static void mdlSetOutputPortComplexSignal(SimStruct *S, int portIndex, CSignal_T cSignalSetting) 94 && { } /* mdlSetOutputPortComplexSignal */ #endif /* MDL_SET_OUTPUT_PORT_COMPLEX_SIGNAL */ #define MDL_SET_DEFAULT_PORT_COMPLEX_SIGNALS /* Change to #undef to remove */ #if defined(MDL_SET_DEFAULT_PORT_COMPLEX_SIGNALS) && defined(MATLAB_MEX_FILE) static void mdlSetDefaultPortComplexSignals(SimStruct *S) { } /* mdlSetDefaultPortComplexSignals */ #endif /* MDL_SET_DEFAULT_PORT_COMPLEX_SIGNALS */ #define MDL_SET_WORK_WIDTHS /* Change to #undef to remove function */ #if defined(MDL_SET_WORK_WIDTHS) && defined(MATLAB_MEX_FILE) static void mdlSetWorkWidths(SimStruct *S) { } #endif /* MDL_SET_WORK_WIDTHS */ #define MDL_INITIALIZE_CONDITIONS /* Change to #undef to remove function */ #if defined(MDL_INITIALIZE_CONDITIONS) static void mdlInitializeConditions(SimStruct *S) 95 { int_T i; real_T *xcont = ssGetContStates(S); int_T nCStates = ssGetNumContStates(S); for (i = 0; i < nCStates; i++) { *xcont++ = 0.0; } } #endif /* MDL_INITIALIZE_CONDITIONS */ #define MDL_START /* Change to #undef to remove function */ #if defined(MDL_START) static void mdlStart(SimStruct *S) { } #endif /* /* Define MDL_START */ to indicate that this S-Function has the mdlG[S]etSimState mothods */ #define MDL_SIM_STATE /* Change to #undef to remove this function */ #if defined(MDL_SIM_STATE) static mxArray* mdlGetSimState(SimStruct* S) { 96 } static void mdlSetSimState(SimStruct* S, const mxArray* inSimState) { } #endif /* MDL_SIM_STATE */ #define MDL_GET_TIME_OF_NEXT_VAR_HIT /* Change to #undef to remove function */ #if defined(MDL_GET_TIME_OF_NEXT_VAR_HIT) && (defined(MATLAB_MEX_FILE) || \ defined(NRT)) static void mdlGetTimeOfNextVarHit(SimStruct *S) { time_T timeOfNextHit = ssGetT(S) /* + offset */ ; ssSetTNext(S, timeOfNextHit); } #endif /* MDL_GET_TIME_OF_NEXT_VAR_HIT */ #define MDL_ZERO_CROSSINGS /* Change to #undef to remove function */ #if defined(MDL_ZERO_CROSSINGS) (defined(MATLAB_MEX_FILE) || defined(NRT)) static void mdlZeroCrossings(SimStruct *S) { 97 && } #endif /* MDL_ZERO_CROSSINGS */ static void mdlOutputs(SimStruct *S, int_T tid) { real_T *y0 = ssGetOutputPortRealSignal(S,0); real_T *y1 = ssGetOutputPortRealSignal(S,1); real_T *y2 = ssGetOutputPortRealSignal(S,2); real_T *y3 = ssGetOutputPortRealSignal(S,3); real_T *y4 = ssGetOutputPortRealSignal(S,4); real_T *y5 = ssGetOutputPortRealSignal(S,5); real_T *m = ssGetOutputPortRealSignal(S,6); real_T *x = ssGetContStates(S); InputRealPtrsType u3Ptrs = ssGetInputPortRealSignalPtrs(S,3); double fas, fbs, fcs, far1, fbr, fcr, lm, ls, lr, zp; double k, kk, ias, ibs, ics, iar, ibr,icr; fas=x[0]; fbs=x[1]; fcs=x[2]; far1=x[3]; fbr=x[4]; fcr=x[5]; 98 lm=mxGetPr(plm)[0]; ls=mxGetPr(plos)[0]+lm; lr=lor[0]+lm; zp=mxGetPr(pzp)[0]; k=2.0*lr*ls-4.0*lm*lm+lr*lm+lm*ls; kk=-lm*lm*ls-lr*ls*lm+2.0*lr*lr*ls+lr*lr*lm+4.0*lm*lm*lm5.0*lr*lm*lm; y0[0]=(2.0*lr*ls+lm*ls-2.0*lm*lm-lr*lm)*fas/(k*(lslm))+lm*(lr-lm)*fbs/(k*(ls-lm))+lm*(lr-lm)*fcs/(k*(lslm))-2.0*lm*far1/k+lm*fbr/k+lm*fcr/k; ias=y0[0]; y1[0]=lm*(lr-lm)*fas/(k*(ls-lm))+(2.0*lr*ls+lm*ls2.0*lm*lm-lr*lm)*fbs/(k*(ls-lm))+lm*(lr-lm)*fcs/(k*(lslm))+lm*far1/k-2.0*lm*fbr/k+lm*fcr/k; ibs=y1[0]; y2[0]=lm*(-lm+lr)*fas/(k*(ls-lm))+lm*(-lm+lr)*fbs/(k*(lslm))+(2.0*lr*ls+lm*ls-2.0*lm*lm-lr*lm)*fcs/(k*(lslm))+lm*far1/k+lm*fbr/k-2.0*lm*fcr/k; ics=y2[0]; y3[0]=-2.0*lm*fas/k+lm*fbs/k+lm*fcs/k+(2.0*lr*ls+lr*lmlm*ls-2.0*lm*lm)*far1/kk+lm*(ls-lm)*fbr/kk+lm*(lslm)*fcr/kk; iar=y3[0]; 99 y4[0]=lm*fas/k-2.0*lm*fbs/k+lm*fcs/k+lm*(lslm)*far1/(k*(lr-lm))+(2.0*lr*ls+lr*lm-lm*ls2.0*lm*lm)*fbr/(k*(lr-lm))+lm*(ls-lm)*fcr/kk; ibr=y4[0]; y5[0]=lm*fas/k+lm*fbs/k-2.0*lm*fcs/k+lm*(lslm)*far1/kk+lm*(ls-lm)*fbr/kk+(2.0*lr*ls+lr*lm-lm*ls2.0*lm*lm)*fcr/kk; icr=y5[0]; m[0]=sqrt(3.0)/2.0 * zp * lm * ((ias*icr+ibs*iar+ics*ibr)(ias*ibr+ibs*icr+ics*iar)); } /* end mdlOutputs */ #define MDL_UPDATE /* Change to #undef to remove function */ #if defined(MDL_UPDATE) static void mdlUpdate(SimStruct *S, int_T tid) { } #endif /* MDL_UPDATE */ #define MDL_DERIVATIVES /* Change to #undef to remove function */ #if defined(MDL_DERIVATIVES) static void mdlDerivatives(SimStruct *S) { real_T *dx = ssGetdX(S); 100 real_T *x = ssGetContStates(S); InputRealPtrsType u0Ptrs = u1Ptrs = u2Ptrs = u3Ptrs = u4Ptrs = u5Ptrs = ssGetInputPortRealSignalPtrs(S,0); InputRealPtrsType ssGetInputPortRealSignalPtrs(S,1); InputRealPtrsType ssGetInputPortRealSignalPtrs(S,2); InputRealPtrsType ssGetInputPortRealSignalPtrs(S,3); InputRealPtrsType ssGetInputPortRealSignalPtrs(S,4); InputRealPtrsType ssGetInputPortRealSignalPtrs(S,5); double fas, fbs, fcs, far1, fbr, fcr, uas, ubs, ucs, wr; double rs, rr; double k, kk; double lm,ls,lr; fas=x[0]; fbs=x[1]; fcs=x[2]; far1=x[3]; fbr=x[4]; fcr=x[5]; uas=usa[0]; 101 ubs=usb[0]; ucs=usc[0]; wr=w[0]; rs=mxGetPr(prs)[0]; rr=rrr[0]; lm=mxGetPr(plm)[0]; ls=mxGetPr(plos)[0]+lm; lr=lor[0]+lm; k=2.0*lr*ls-4.0*lm*lm+lr*lm+lm*ls; kk=-lm*lm*ls-lr*ls*lm+2.0*lr*lr*ls+lr*lr*lm+4.0*lm*lm*lm5.0*lr*lm*lm; dx[0]=-rs*(2.0*lr*ls+lm*ls-2.0*lm*lm-lr*lm)*fas/(k*(lslm))-rs*lm*(lr-lm)*fbs/(k*(ls-lm))-rs*lm*(lrlm)*fcs/(k*(ls-lm))+2.0*rs*lm*far1/k-rs*lm*fbr/krs*lm*fcr/k+uas; dx[1]=-rs*lm*(lr-lm)*fas/(k*(ls-lm))-rs*(2.0*lr*ls+lm*ls2.0*lm*lm-lr*lm)*fbs/(k*(ls-lm))-rs*lm*(lr-lm)*fcs/(k*(lslm))-rs*lm*far1/k+2.0*rs*lm*fbr/k-rs*lm*fcr/k+ubs; dx[2]=-rs*lm*(-lm+lr)*fas/(k*(ls-lm))-rs*lm*(lm+lr)*fbs/(k*(ls-lm))-rs*(2.0*lr*ls+lm*ls-2.0*lm*lmlr*lm)*fcs/(k*(ls-lm))-rs*lm*far1/krs*lm*fbr/k+2.0*rs*lm*fcr/k+ucs; dx[3]=2.0*rr*lm*fas/k-rr*lm*fbs/k-rr*lm*fcs/krr*(2.0*lr*ls+lr*lm-lm*ls-2.0*lm*lm)*far1/kk+(-rr*lm*(lslm)/kk-(1.0/3.0)*wr*sqrt(3.0))*fbr+(-rr*lm*(lslm)/kk+(1.0/3.0)*wr*sqrt(3.0))*fcr; 102 dx[4]=-rr*lm*fas/k+2.0*rr*lm*fbs/k-rr*lm*fcs/k+(rr*lm*(ls-lm)/(k*(-lm+lr))+(1.0/3.0)*wr*sqrt(3.0))*far1rr*(2.0*lr*ls+lr*lm-lm*ls-2.0*lm*lm)*fbr/(k*(-lm+lr))+(rr*lm*(ls-lm)/kk-(1.0/3.0)*wr*sqrt(3.0))*fcr; dx[5]=-rr*lm*fas/k-rr*lm*fbs/k+2.0*rr*lm*fcs/k+(rr*lm*(ls-lm)/kk-(1.0/3.0)*wr*sqrt(3.0))*far1+(-rr*lm*(lslm)/kk+(1.0/3.0)*wr*sqrt(3.0))*fbr-rr*(2.0*lr*ls+lr*lmlm*ls-2.0*lm*lm)*fcr/kk; } #endif /* MDL_DERIVATIVES */ static void mdlTerminate(SimStruct *S) { } #define MDL_RTW /* Change to #undef to remove function */ #if defined(MDL_RTW) && defined(MATLAB_MEX_FILE) static void mdlRTW(SimStruct *S) { } #endif /* MDL_RTW */ /*=============================* * Required S-function trailer * *=============================*/ #ifdef MATLAB_MEX_FILE /* Is this file being compiled as a MEX-file? */ 103 #include "simulink.c" /* MEX-file interface mechanism */ #else #include "cg_sfun.h" /* Code generation registration function */ #endif 104 ... theo thời gian, nhờ tính tốn tiêu động động bước thiết kế Chính em xin chọn đề tài ? ?Nghiên cứu thiết kế động không đồng pha rơto lồng sóc rãnh sâu có đến tính đến hiệu ứng mặt ngồi dẫn rơto” vii...BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ NÔNG NGHIỆP VÀ PTNT TRƯỜNG ĐẠI HỌC THỦY LỢI LÊ QUỐC DUY NGHIÊN CỨU THIẾT KẾ ĐỘNG CƠ KHÔNG ĐỒNG BỘ PHA RƠTO LỒNG SĨC RÃNH SÂU CĨ TÍNH ĐẾN HIỆU ỨNG MẶT NGỒI TRONG... tìm tiêu động động trình thiết kế Đối tượng phạm vi nghiên cứu Động khơng đồng pha rơto lồng sóc rãnh sâu Cách tiếp cận phương pháp nghiên cứu: Nghiên cứu lý thuyết, tính tốn thiết kế, lập trình