current.c 6.6 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228
  1. #include <math.h>
  2. #include "bsp/adc.h"
  3. #include "bsp/pwm.h"
  4. #include "foc/motor/current.h"
  5. #include "foc/core/PMSM_FOC_Core.h"
  6. #include "libs/utils.h"
  7. #include "libs/logger.h"
  8. #include "math/fix_math.h"
  9. static current_samp_t g_cs;
  10. #define NB_OFFSET_SAMPLES 32
  11. #define Rvbus 0.0005f
  12. #define Gvbus (13.1f) //母线电流的运放
  13. #define Rds_Defualt 0.005f//欧
  14. #define Gmos (1.7f)//mos 电流的运放
  15. #define Sample_R Rds_Defualt
  16. #define Lower_Pass_p 0.2f
  17. #define VBUS_VOL(adc) (((float)(adc)) * 3.3f / 4096.0f / Gvbus)
  18. #define MOSds_VOL(adc) (((float)(adc)) * 3.3f / 4096.0f / Gmos)
  19. #define current_i(v, r) ((v)/(r))
  20. /*
  21. va = iv*kv - ia*k
  22. vb = iv*kv - ib*k
  23. vc = iv*kv - ic*k -> vb+vc = 2iv*kv - (ib+ic)*k ->vb+vc = 2iv*kv + ia*k -->va+vb+vc =
  24. ia = -(ib+ic)
  25. */
  26. void phase_current_init(void) {
  27. current_samp_t *cs = &g_cs;
  28. cs->offset_sample_count = NB_OFFSET_SAMPLES + 1;
  29. cs->adc_ia = 0;
  30. cs->adc_ib = 0;
  31. cs->adc_ic = 0;
  32. }
  33. void phase_current_start_cali(void){
  34. g_cs.adc_offset_a = 0;
  35. g_cs.adc_offset_b = 0;
  36. g_cs.adc_offset_c = 0;
  37. phase_current_init();
  38. g_cs.is_calibrating_offset = true;
  39. g_cs.c_phases = PHASE_AB;
  40. adc_current_sample_config(g_cs.c_phases);
  41. }
  42. void phase_current_wait_cali(void) {
  43. while(g_cs.is_calibrating_offset) {
  44. wdog_reload();
  45. }
  46. }
  47. bool phase_current_offset(void) {
  48. current_samp_t *cs = &g_cs;
  49. if (!cs->is_calibrating_offset) {
  50. return false;
  51. }
  52. s32 phase_current1 = 0 , phase_current2 = 0;
  53. adc_phase_current_read(cs->c_phases, &phase_current1, &phase_current2);
  54. if (cs->offset_sample_count == (NB_OFFSET_SAMPLES + 1)) {
  55. cs->offset_sample_count --;
  56. return true;
  57. }
  58. if (cs->offset_sample_count > 0) {
  59. cs->offset_sample_count--;
  60. if (cs->c_phases == PHASE_AB && cs->offset_sample_count >= 0) {
  61. cs->adc_offset_a += phase_current1;
  62. cs->adc_offset_b += phase_current2;
  63. if (cs->offset_sample_count == 0) {
  64. cs->adc_offset_a = cs->adc_offset_a / NB_OFFSET_SAMPLES;
  65. cs->adc_offset_b = cs->adc_offset_b / NB_OFFSET_SAMPLES;
  66. }
  67. }
  68. if (cs->c_phases == PHASE_BC && cs->offset_sample_count >= 0) {
  69. cs->adc_offset_c += phase_current2;
  70. if (cs->offset_sample_count == 0) {
  71. cs->adc_offset_c = cs->adc_offset_c / NB_OFFSET_SAMPLES;
  72. }
  73. }
  74. }else {
  75. if (cs->c_phases == PHASE_AB) {
  76. cs->c_phases = PHASE_BC;
  77. phase_current_init();
  78. adc_current_sample_config(cs->c_phases);
  79. }else {
  80. cs->is_calibrating_offset = false;
  81. sys_debug("offset %d, %d, %d\n", g_cs.adc_offset_a, g_cs.adc_offset_b, g_cs.adc_offset_c);
  82. }
  83. }
  84. return true;
  85. }
  86. void phase_current_get(float *iABC){
  87. current_samp_t *cs = &g_cs;
  88. s32 phase_current1, phase_current2;
  89. adc_phase_current_read(cs->c_phases, &phase_current1, &phase_current2);
  90. if (cs->c_phases == PHASE_AB) {
  91. /* Current on Phase C is not accessible */
  92. /* Ia = PhaseAOffset - ADC converted value) */
  93. cs->adc_ia = (phase_current1 - cs->adc_offset_a);
  94. cs->adc_ib = (phase_current2 - cs->adc_offset_b);
  95. cs->adc_ic = -(cs->adc_ia + cs->adc_ib);
  96. }else if (cs->c_phases == PHASE_BC) {
  97. /* Current on Phase A is not accessible */
  98. /* Ib = PhaseBOffset - ADC converted value) */
  99. cs->adc_ib = (phase_current1 - cs->adc_offset_b);
  100. cs->adc_ic = (phase_current2 - cs->adc_offset_c);
  101. cs->adc_ia = -(cs->adc_ib + cs->adc_ic);
  102. }else if (cs->c_phases == PHASE_AC) {
  103. /* Current on Phase B is not accessible */
  104. /* Ia = PhaseAOffset - ADC converted value) */
  105. cs->adc_ia = (phase_current1 - cs->adc_offset_a);
  106. cs->adc_ic = (phase_current2 - cs->adc_offset_c);
  107. cs->adc_ib = -(cs->adc_ia + cs->adc_ic);
  108. }
  109. iABC[0] = -cs->adc_ia * ADC_TO_CURR_ceof;
  110. iABC[1] = -cs->adc_ib * ADC_TO_CURR_ceof;
  111. iABC[2] = -cs->adc_ic * ADC_TO_CURR_ceof;
  112. }
  113. #if 0
  114. void phase_current_point(void *p){
  115. FOC_OutP *out = p;
  116. current_samp_t *cs = &g_cs;
  117. u32 low_side_low_duty = FOC_PWM_Half_Period - out->n_lowDuty;
  118. u32 low_side_mid_duty = FOC_PWM_Half_Period - out->n_midDuty;
  119. cs->sector = out->n_Sector;
  120. out->n_Sample1 = FOC_PWM_Half_Period + 1;
  121. out->n_Sample2 = FOC_PWM_Half_Period + 1;
  122. s16 potest = 7;
  123. /*底边开mos的时间是2倍的 low_side_low_duty(一个周期)*/
  124. if (low_side_low_duty * 2 >= TSampleMIN) { //可以采样
  125. if (low_side_low_duty >= (TADC + TDead)) {//可以在pwm的中心点采样
  126. out->n_Sample1 = FOC_PWM_Half_Period - 1;
  127. cs->sector = SECTOR_1;
  128. potest = 1;
  129. }else {
  130. u32 Samp_p = out->n_lowDuty + TSampleBefore;
  131. if (Samp_p >= FOC_PWM_Half_Period) { //需要在pwm中心点过后采样,需要配置PWM0模式
  132. out->n_Sample2 = ( 2u * FOC_PWM_Half_Period ) - Samp_p - (uint16_t) 1;
  133. potest = 2;
  134. }else {
  135. out->n_Sample1 = Samp_p;
  136. potest = 3;
  137. }
  138. }
  139. }else if (low_side_mid_duty * 2 >= TSampleMIN){
  140. if (low_side_mid_duty >= (TADC + TDead)) {//可以在pwm的中心点采样
  141. out->n_Sample1 = FOC_PWM_Half_Period - 1;
  142. potest = 4;
  143. }else {
  144. u32 Samp_p = out->n_midDuty + TSampleBefore;
  145. if (Samp_p >= FOC_PWM_Half_Period) { //需要在pwm中心点过后采样,需要配置PWM0模式
  146. out->n_Sample2 = ( 2u * FOC_PWM_Half_Period ) - Samp_p - (uint16_t) 1;
  147. potest = 5;
  148. }else {
  149. out->n_Sample1 = Samp_p;
  150. potest = 6;
  151. }
  152. }
  153. }
  154. plot_1data16(potest);
  155. out->n_Sector = cs->sector;
  156. }
  157. #else
  158. static __inline__ s16 get_mid_duty(u32 max_hi_side, FOC_OutP *out) {
  159. u32 min_low_side = FOC_PWM_Half_Period - max_hi_side;
  160. if (min_low_side >= MAX(TSampleBefore, TADC)) { //可以采样
  161. if ((min_low_side > TSampleBefore) && (min_low_side > TADC)) {
  162. out->n_Sample1 = FOC_PWM_Half_Period - 1;
  163. return 0;
  164. }else {
  165. u32 Samp_p = max_hi_side + TSampleBefore;
  166. if (Samp_p >= FOC_PWM_Half_Period) { //需要在pwm中心点过后采样,需要配置PWM0模式
  167. out->n_Sample2 = ( 2u * FOC_PWM_Half_Period ) - Samp_p - (uint16_t) 1;
  168. return 1;
  169. }else {
  170. out->n_Sample1 = Samp_p;
  171. return 2;
  172. }
  173. }
  174. }
  175. return -1;
  176. }
  177. void phase_current_point(void *p){
  178. FOC_OutP *out = p;
  179. current_samp_t *cs = &g_cs;
  180. u32 max_hi_side = MAX(out->n_Duty[0], MAX(out->n_Duty[1], out->n_Duty[2]));
  181. out->n_Sample1 = FOC_PWM_Half_Period + 1;
  182. out->n_Sample2 = FOC_PWM_Half_Period + 1;
  183. if (get_mid_duty(max_hi_side, out) == -1) {
  184. u32 mid_hi_side;
  185. out->n_Sample1 = FOC_PWM_Half_Period + 1;
  186. out->n_Sample2 = FOC_PWM_Half_Period + 1;
  187. if (out->n_Duty[0] == max_hi_side) {
  188. //sample B, C
  189. mid_hi_side = MAX(out->n_Duty[1], out->n_Duty[2]);
  190. cs->c_phases = PHASE_BC;
  191. }else if (out->n_Duty[1] == max_hi_side) {
  192. //sample A, C
  193. mid_hi_side = MAX(out->n_Duty[0], out->n_Duty[2]);
  194. cs->c_phases = PHASE_AC;
  195. }else {
  196. //sample A, B
  197. mid_hi_side = MAX(out->n_Duty[0], out->n_Duty[1]);
  198. cs->c_phases = PHASE_AB;
  199. }
  200. get_mid_duty(mid_hi_side, out);
  201. }else {
  202. cs->c_phases = PHASE_AB;
  203. }
  204. out->n_CPhases = cs->c_phases;
  205. }
  206. #endif
  207. void phase_current_adc_triger(void){
  208. adc_enable_ext_trigger();
  209. }