1 00:00:00,900 --> 00:00:07,950 Welcome to the South Asian Development of Vehicle Model in Matlab simulation, in this situation, we 2 00:00:07,980 --> 00:00:13,320 will develop the mathematical model of vehicle by solving the differential equation with the help of 3 00:00:13,320 --> 00:00:14,300 Matlab simulation. 4 00:00:14,730 --> 00:00:21,030 If you are already about the using of Matlab simulation for solving the differential equation, then 5 00:00:21,030 --> 00:00:21,380 it's fine. 6 00:00:21,420 --> 00:00:28,260 Otherwise we will also go step by step procedure for this particular example of vehicle model we have 7 00:00:28,260 --> 00:00:34,830 already discussed about the dynamic equation of the vehicle and in decision value the simulation for 8 00:00:34,830 --> 00:00:37,200 developing the Mathematica model. 9 00:00:39,700 --> 00:00:50,110 We will use this equation they applied for is equal to A that is M into my DETI with the velocity plus 10 00:00:50,110 --> 00:00:53,150 opposing force of aerodynamic drag. 11 00:00:53,830 --> 00:01:01,780 So this opposing force and plus if it is a vehicle, is at incline angle B to then the component of 12 00:01:01,780 --> 00:01:08,080 gravitational force is Modisane Beita to these three force and opposing the supply force. 13 00:01:08,080 --> 00:01:11,740 And it will result in equilibrium speed for applying for. 14 00:01:11,770 --> 00:01:16,930 So it will determine the stability, steady state speed and also give us the idea about the response, 15 00:01:17,110 --> 00:01:22,610 how fast or slow we will let you the steady state speed for the Applied Force. 16 00:01:23,110 --> 00:01:25,630 So this is the system level simulation here. 17 00:01:25,750 --> 00:01:34,160 Our main input parameters that weight and inertia of the vehicle and supply force. 18 00:01:34,570 --> 00:01:40,340 So this is called the system simulation because we use the general assistance simulation to represent 19 00:01:40,340 --> 00:01:44,640 the mechanical, electrical or even magnetic or chemical system. 20 00:01:45,370 --> 00:01:51,970 Next, higher in the simulation level at the unit level simulation that we used up in that simulation 21 00:01:51,970 --> 00:02:00,000 of different units like the motor converter than engine or attacks and different segments, level simulation. 22 00:02:00,010 --> 00:02:02,440 And then we join this simulation together. 23 00:02:02,740 --> 00:02:03,520 This is next step. 24 00:02:03,520 --> 00:02:09,880 We are not going to discuss that a unit level and furthermore, for the higher prediction and to predict 25 00:02:09,880 --> 00:02:15,790 the more information we use, the component level simulation simulation of Las Vegas, simulation of 26 00:02:15,790 --> 00:02:18,320 joint, and then we can work together. 27 00:02:18,640 --> 00:02:27,310 So because this is the multidisciplinary course, we use the first system level simulation here we are 28 00:02:27,310 --> 00:02:31,440 designing the simulation of vehicle simulation models vehicle. 29 00:02:31,960 --> 00:02:39,120 So this differential equation will be simulated in Matlab if we use this equation. 30 00:02:39,150 --> 00:02:43,030 So but it is equal to supply force. 31 00:02:43,180 --> 00:02:51,370 Applied Force minus 40 is that dynamic drag force and minus Modisane Beta. 32 00:02:51,820 --> 00:02:57,220 This can be easily simulated with the help of Simulant Library of. 33 00:02:58,360 --> 00:02:58,930 Let's see. 34 00:03:00,370 --> 00:03:07,200 Here is the model I developed for the vehicle so we can see the, the. 35 00:03:08,580 --> 00:03:11,740 And platforms and component design. 36 00:03:12,480 --> 00:03:20,160 So if you recall this equation at this point, so we are having the question and to be reminded is equal 37 00:03:20,160 --> 00:03:25,580 to supply force F that is force. 38 00:03:25,590 --> 00:03:34,230 We apply opposing force represented by minus sign F. Lee and M.G. S. Bita. 39 00:03:34,440 --> 00:03:37,070 So opposing force that F3 and M.G. to. 40 00:03:37,110 --> 00:03:43,940 So we had m divided is equal to F minus F and minus Modisane Bita. 41 00:03:44,220 --> 00:03:46,780 So we'll see what this FBN, what it means and B. 42 00:03:46,980 --> 00:03:54,260 But at this point what we need to realize, the equation that we discussed is realized at this point 43 00:03:54,450 --> 00:04:00,240 m divided is equal to F minus 70 minus Modisane. 44 00:04:01,080 --> 00:04:07,500 Now in itself a question of the sometimes we use the negative side to represent it is opposite to the 45 00:04:07,500 --> 00:04:09,290 direction of velocity. 46 00:04:09,720 --> 00:04:17,230 But to understand by assuming the direction of velocity in the positive direction, it will oppose SOYER 47 00:04:17,280 --> 00:04:18,180 negative SANWAL. 48 00:04:20,610 --> 00:04:30,840 So after this point, if we move forward, then if we want to integrate it and before if we want to 49 00:04:30,840 --> 00:04:37,360 divide by the mass with the gain block, you can find the gain block and commonly used block yet again 50 00:04:37,380 --> 00:04:38,230 you can it. 51 00:04:38,580 --> 00:04:42,270 And this gain block, you can set whatever you want. 52 00:04:42,810 --> 00:04:46,020 So for more clarity, I read the same quote. 53 00:04:46,020 --> 00:04:48,310 One by 200, 200 is the math. 54 00:04:48,330 --> 00:04:51,330 So here we have assumed the vehicle math is 200 cagy. 55 00:04:51,930 --> 00:04:56,310 And if we divide this M divided by 200. 56 00:04:56,340 --> 00:04:58,240 So what will remain at the site? 57 00:04:58,710 --> 00:05:00,780 So we will be by DETI. 58 00:05:01,050 --> 00:05:10,260 So here we are having DBI did not feel the integration of divisibility then integration block will be 59 00:05:10,260 --> 00:05:11,750 used that also you can find. 60 00:05:11,760 --> 00:05:12,000 Yes. 61 00:05:12,330 --> 00:05:16,980 And disintegrator block we can use. 62 00:05:18,300 --> 00:05:25,980 So integrator is used so that at the end of integrated output will be regex. 63 00:05:26,580 --> 00:05:29,010 So integration up divided this way. 64 00:05:29,580 --> 00:05:36,390 So what we want, we want to predict if we apply a certain amount of force, what will be the speed. 65 00:05:36,660 --> 00:05:42,200 Yet force is also constant, steady state, but in actual spirit it will be in actual performance. 66 00:05:42,210 --> 00:05:42,990 It will be ready. 67 00:05:42,990 --> 00:05:45,270 We'll see later how we can do this. 68 00:05:45,270 --> 00:05:46,920 But this is the vehicle simulation. 69 00:05:46,920 --> 00:05:53,600 We can do this for the applied for what we've done speed yet spellbind meters per second. 70 00:05:54,120 --> 00:06:02,460 And if you remember the equation of the aerodynamic force, then we are having this out in component 71 00:06:02,730 --> 00:06:10,470 Pantelleria, concentric Koshien, much density of air and this is for one by two multiplied by zero 72 00:06:10,470 --> 00:06:18,830 point one by two, the frontal area, a direct quotient and mass density of at rule. 73 00:06:19,920 --> 00:06:24,800 Yet we use the typical values one point two zero point four and one point two. 74 00:06:25,170 --> 00:06:32,430 This is representing one by two entities representing the Met function that is squared squared off the 75 00:06:32,580 --> 00:06:33,300 velocity. 76 00:06:33,630 --> 00:06:37,690 So whatever input it will be had and it will output. 77 00:06:37,950 --> 00:06:42,450 So yeah, we need to the squared output because we also need to. 78 00:06:44,130 --> 00:06:54,030 Use the sign off velocity so that if velocity is reversed, the direction of force, FDR would also 79 00:06:54,090 --> 00:07:00,450 reverse, but because here we are not going to reverse the direction we're applying the negative force. 80 00:07:00,460 --> 00:07:04,790 So for the moment, we can let it only to function as a square. 81 00:07:04,800 --> 00:07:10,920 You can find the math function from the SIMULANT library and met operator. 82 00:07:10,950 --> 00:07:12,780 You can find this math function. 83 00:07:13,200 --> 00:07:19,950 And if you open this parameter, you can select the function you want and magnitude square that can 84 00:07:20,160 --> 00:07:26,270 use the input magnitude will be scared and the square will be the output at this terminal. 85 00:07:26,610 --> 00:07:34,980 So you can see this equation of Apne V squared multiplied by one by two, multiplied by one point two 86 00:07:34,980 --> 00:07:42,000 multiplied by zero point for the and multiply by frontal area one point two that this is representing 87 00:07:42,000 --> 00:07:48,570 the FBI aerodynamic drag force that is opposing to the Applied Force. 88 00:07:50,930 --> 00:07:55,880 Now, come to this spot, this is and this would be to say a few. 89 00:07:57,410 --> 00:08:06,430 Go and read the William M Jerrell Agenda Masses 200 Kaiji into nine point eight one is the constant 90 00:08:07,100 --> 00:08:11,050 so and Rita is currently we can see it from here. 91 00:08:11,450 --> 00:08:18,990 Currently it is set to zero and sine function we can use from the Met operator. 92 00:08:29,340 --> 00:08:37,740 This same function can be used from Meritocrat, so this will become wita into to sign to sign Bitta 93 00:08:37,950 --> 00:08:40,150 and sign it and multiply by M.G.. 94 00:08:40,170 --> 00:08:43,860 So this will be this input becomes M.G. sign Bita. 95 00:08:44,610 --> 00:08:50,250 So this is our complete simulation model of the vehicle having the weight mass. 96 00:08:50,610 --> 00:08:55,950 And if we are having the radius, then we can also supplied the talk instead of force. 97 00:08:55,950 --> 00:09:02,880 That Burnsy not to simulate this model, we need to just run it for the particular model. 98 00:09:02,880 --> 00:09:06,670 Second, by default, it will be 10 seconds in your model. 99 00:09:06,900 --> 00:09:10,310 You need to include the second because of the high inertia. 100 00:09:10,320 --> 00:09:19,830 It will take too much time in physical seconds so our vehicle cannot reach to speed at 13 kilometers 101 00:09:19,830 --> 00:09:22,960 per hour in the fraction of second in one second or two seconds. 102 00:09:23,340 --> 00:09:31,200 So we need to simulate for certain amount of time support and setting the two under second and fourth 103 00:09:31,200 --> 00:09:34,230 applied force in Newton and Newton. 104 00:09:35,640 --> 00:09:43,380 If I simulate, we will see this, this represents the speed and meter per second and final speed, 105 00:09:44,010 --> 00:09:45,010 let's say, is it? 106 00:09:45,030 --> 00:09:48,060 The statistics are not to check this. 107 00:09:48,070 --> 00:09:56,510 We need to check the performance so you can see it is not showing the exponential plot. 108 00:09:56,730 --> 00:09:58,280 And this is still increasing. 109 00:09:58,290 --> 00:10:00,190 It is not steady state speed. 110 00:10:00,210 --> 00:10:02,000 We need to increase the second. 111 00:10:02,010 --> 00:10:07,320 So we will do the two things for refining this data point. 112 00:10:07,710 --> 00:10:11,850 We need to we know that this is exponential response. 113 00:10:12,180 --> 00:10:19,410 But for the defense, you can go to the simulation model configuration parameter and under the data 114 00:10:19,410 --> 00:10:20,430 input expert. 115 00:10:23,060 --> 00:10:29,660 You can see the zero option under the safe option output option if, in fact, you need to increase 116 00:10:29,660 --> 00:10:32,360 the defense factor if output is not. 117 00:10:32,870 --> 00:10:39,230 So when you are using the same Willink model component, not the discrete and not using any electrical 118 00:10:39,860 --> 00:10:43,750 sympathectomy library component of the simulation will be very fast. 119 00:10:43,760 --> 00:10:46,030 But you need to use this different factor. 120 00:10:46,040 --> 00:10:51,080 You can increase to 10, 20 or I used 50 for refinement. 121 00:10:51,080 --> 00:10:56,440 Even if you design functional sine wave, sometime you will find the distortion. 122 00:10:56,630 --> 00:11:00,760 We often think that this refined factor and let's run it again. 123 00:11:02,720 --> 00:11:13,000 So we see this is now our actual response, exponential response of differential equation. 124 00:11:13,010 --> 00:11:20,930 And we can see the steady state speed is eighteen point six three meters per second. 125 00:11:21,950 --> 00:11:26,750 So this is what our we model will give us the exact idea. 126 00:11:26,750 --> 00:11:33,590 If we apply the constant photon at this rate, speed will increase and what the final steady state of 127 00:11:33,590 --> 00:11:37,010 the vehicle for the particular amount of force and. 128 00:11:38,550 --> 00:11:45,690 We will have also at how much time, what will be just be retained by the weekend, we can convert it 129 00:11:45,690 --> 00:11:48,460 in cell phone from the meters per second to kilometer. 130 00:11:48,580 --> 00:11:59,220 But with this game clock, so if you have the speed in kilometer per hour, we will have around sixty 131 00:11:59,790 --> 00:12:00,180 five. 132 00:12:03,950 --> 00:12:04,790 You can see. 133 00:12:07,120 --> 00:12:15,760 Sixty seven kilometer per hour, final speed, maximum speed, steady state speed, and this development, 134 00:12:16,690 --> 00:12:26,530 we can say, will also make considering the Guarascio one, what we can see the girls drive six three 135 00:12:26,560 --> 00:12:27,870 eight eight PM. 136 00:12:28,450 --> 00:12:30,820 You can see the evolution, but Mynott. 137 00:12:31,120 --> 00:12:37,470 So to convert the minute into second, it is multiplied by 60 because of the evolution minute. 138 00:12:37,480 --> 00:12:38,650 So it will become the evolution. 139 00:12:38,650 --> 00:12:45,310 But second and one, the evolution representa to buy at meters to represent one. 140 00:12:45,310 --> 00:12:51,850 The evolution is divided by two by two into three point one, four into Arietta areas that we assume. 141 00:12:52,450 --> 00:12:57,020 So it is around zero five to seven, nine is around 11 each. 142 00:12:57,550 --> 00:13:07,960 We can see it is we are doing the mopper or such kind of vehicle because we are a motorcycle as the 143 00:13:07,960 --> 00:13:09,760 radius around soundin each. 144 00:13:11,710 --> 00:13:12,670 So we are not. 145 00:13:13,410 --> 00:13:21,250 So if you are using the model bike, then accordingly you can set the radius for real speed and you 146 00:13:21,250 --> 00:13:22,210 can use it. 147 00:13:24,120 --> 00:13:31,200 So we will also have the distance covered because this is the velocity of the military, so if you want 148 00:13:31,200 --> 00:13:39,090 to know how much distance travel at the end of this time, we can also see distance travel is around 149 00:13:39,090 --> 00:13:48,330 three thousand in between 3000 to 30000 meter administrator, three point five km. 150 00:13:50,070 --> 00:13:51,750 In 200 second. 151 00:14:00,670 --> 00:14:06,550 So this is all output, we will help for this input, and this is called a system level model of the 152 00:14:07,000 --> 00:14:07,640 vehicle. 153 00:14:08,620 --> 00:14:12,430 Now we can also talk component instead of force if we want. 154 00:14:12,880 --> 00:14:16,440 So what we need to do, we need to use the gain. 155 00:14:17,020 --> 00:14:19,950 So what is the talk I read again? 156 00:14:20,680 --> 00:14:24,600 So it will remain force if it is dark, a new transmitter. 157 00:14:25,360 --> 00:14:31,180 So I am changing this to talk about and become. 158 00:14:32,530 --> 00:14:41,420 Back and force is equal to talk radio so we can talk is equal to force, multiplied by radius. 159 00:14:41,740 --> 00:14:48,160 So if I use this is one by a radius is zero point two seven nine. 160 00:14:50,450 --> 00:14:58,040 But this is talk, this is for now, if we apply the ten hundred Newton Newton of talk, then we can 161 00:14:58,040 --> 00:15:05,360 see the performance speed is 120 greater than 120. 162 00:15:07,100 --> 00:15:11,230 Steady state speed is 120 seven km. 163 00:15:11,260 --> 00:15:13,400 But this is inevitable. 164 00:15:13,640 --> 00:15:21,890 Second will be so you can use the display block for all this also so you can measure the speed and relative 165 00:15:22,220 --> 00:15:25,490 and distance travel for the airplane. 166 00:15:26,690 --> 00:15:32,690 So this is how we can also include the angle to observe the performance. 167 00:15:33,080 --> 00:15:39,310 So this is how we can simulate the vehicle model with the help of Matlab blocks. 168 00:15:40,250 --> 00:15:48,560 Now, what we need to do is you need to store the characteristics by doing the experiment in this simulation 169 00:15:48,560 --> 00:16:00,800 model so that you need to do is starting from the minimum task that is, say, five Newton Mitka dissimulate 170 00:16:00,800 --> 00:16:10,250 eight and record the steady state speed intermap kilometer, but not ensure that it is steady state 171 00:16:10,250 --> 00:16:16,390 speed, because when you reduce to talk, it will take more time to read the steady state speed they 172 00:16:16,400 --> 00:16:16,760 produce. 173 00:16:16,760 --> 00:16:22,360 Taking longer than your simulation time, then it may not be steady state speed. 174 00:16:22,760 --> 00:16:24,620 So yet I need to return. 175 00:16:24,620 --> 00:16:25,990 It is kilometer per. 176 00:16:27,490 --> 00:16:27,960 Ah. 177 00:16:33,460 --> 00:16:34,420 And this is. 178 00:16:37,340 --> 00:16:37,730 And we. 179 00:16:40,100 --> 00:16:45,250 So let's observe it is state or not, you can say this is not a steady state speed. 180 00:16:45,650 --> 00:16:52,730 When we are playing the less amount of risk, it will take more time to fly under so we can have steady 181 00:16:52,730 --> 00:16:53,390 state speed. 182 00:16:53,810 --> 00:16:57,500 So you can see now the steady state speed. 183 00:16:57,510 --> 00:16:58,730 So it is fine. 184 00:16:59,150 --> 00:17:00,170 So what we want. 185 00:17:00,590 --> 00:17:01,430 We want to. 186 00:17:03,440 --> 00:17:08,300 The talk and measure the speed in KM. 187 00:17:09,390 --> 00:17:12,750 But I will use this table. 188 00:17:18,740 --> 00:17:19,780 You can meet the. 189 00:17:23,180 --> 00:17:28,440 So for the five, the spirit is twenty eight point four. 190 00:17:29,030 --> 00:17:33,290 Let's do the experiment one more time, increased by five to four. 191 00:17:33,290 --> 00:17:37,610 The ten Newton meters statistics is forty 192 00:17:40,400 --> 00:17:44,480 four point forty point one six. 193 00:17:45,660 --> 00:17:52,810 Likewise, you need to plug up to the hundred with the increment of five. 194 00:17:53,550 --> 00:17:53,970 So. 195 00:17:56,820 --> 00:18:02,970 For the 15 days, forty nine point one, so when you plot these characteristics, toxified character, 196 00:18:03,180 --> 00:18:06,650 it called the toxified characteristic of a vehicle, we're measuring the statistics. 197 00:18:06,860 --> 00:18:13,800 We will use this characteristic for merging the different system level simulation, electrical with 198 00:18:14,430 --> 00:18:17,130 this mechanical model of the vehicle. 199 00:18:18,240 --> 00:18:20,190 So this talk is generated. 200 00:18:20,460 --> 00:18:28,200 So, you know, this experiment and Shari'ah result of the task with characteristic of the vehicle in 201 00:18:28,200 --> 00:18:30,700 terms of local people as well as in terms of plot. 202 00:18:30,720 --> 00:18:33,240 So it will be exponentially increasing. 203 00:18:33,240 --> 00:18:36,450 You will realize that it will increase exponentially. 204 00:18:36,780 --> 00:18:44,970 But this is called the speed characteristic of our speed makeup that can be predicted by this model. 205 00:18:47,600 --> 00:18:51,050 So this is the output in terms of speed. 206 00:18:52,030 --> 00:18:58,420 And acceleration performance also we can predict and distance travel, but main important thing, we 207 00:18:58,420 --> 00:19:04,660 are interested in the actual speed and the acceleration of the vehicle and this talk is generated from 208 00:19:04,660 --> 00:19:06,250 the electric motor. 209 00:19:06,610 --> 00:19:14,560 So in that section of the electric drive, we will discuss about the aspect of electrical, motor and 210 00:19:14,560 --> 00:19:21,010 system level simulation of electric motor that can generate this talk and this talk. 211 00:19:21,310 --> 00:19:22,980 When applied to the vehicle. 212 00:19:23,050 --> 00:19:31,320 We can predict the performance of the system of combining the electrical as well as the mechanical systems. 213 00:19:31,360 --> 00:19:37,760 When we supply the voltage to the electrical system, it will generate the thought that we will see 214 00:19:37,800 --> 00:19:43,510 the motor and this is applied to the vehicle model and it will result in speed. 215 00:19:43,870 --> 00:19:48,670 So when we supply the voltage, we will get the speed and the performance that. 216 00:19:49,910 --> 00:19:57,410 Hope you got that basic idea about the system level simulation of vehicle with the help of Matlab simulation 217 00:19:57,740 --> 00:19:58,250 model.