# How can you calculate Peak Pulse Power for a Non-Rectangular Power Pulse

I was reading this Vishay document on how Peak Pulse Power can affect resistors.

There is a flow which is attached in the document describing the process of calculating the Peak Pulse Power for a Non-Rectangular Power Pulse (In my case a discharging capacitor curve). Does anyone know how to transfer to a square wave and calculate the Peak Pulse Power as in the flow attached? • I would use a simulator because they can do this calculation for you. Sep 26 at 16:57

The same document explains this on item 4: So, you need to write the voltage's equation and solve the energy integral:

$$E = \int_{T1}^{T2}p(t) \ dt$$

If you have $$\p(t)\$$ put it directly. If you have voltage instead, use:-

$$E = \int_{T1}^{T2}\frac{v^2(t)}{R} \ dt$$

That'll give you the pulse energy. From here, you can get a rectangular pulse with same peak power/voltage but different duration (TX) to give the same pulse energy.

$$E = P_{pk} \ T_X = \frac{V_{pk}^2}{R} \ T_X$$

Consider the pulse shape below just as an example (yours will be exponential): We can write the following:-

$$v(t) = -\frac{V_{pk}}{T_p} \ t + V_{pk}$$

So the pulse energy is:-

$$E = \int_0^{T_p}\frac{v^2}{R} \ dt = ... = \frac{V_{pk}^2}{R} \ T_p\Big(\frac{1}{3}-1+1\Big)=\frac{V_{pk}^2}{R} \ \frac{T_p}{3}$$

So, as you can see, to have the same pulse energy from a rectangular pulse having the same peak power (or voltage), the duration must be one third i.e. $$\T_X=T_p/3\$$.

Numerical example: For a triangular pulse as shown above, having $$\\mathrm{V_{pk}=200\ V}\$$ and $$\\mathrm{T_{p}=3\ ms}\$$ across $$\\mathrm{R = 100 \ \Omega}\$$ , pulse energy will be $$\\mathrm{E = 0.4 \ J}\$$. A rectangular pulse having the same peak power, $$\\mathrm{P_{pk}=200^2/100=400 \ W}\$$, should have a duration of $$\\mathrm{T_X=E/P_{pk}= 0.4/400=1 \ ms}\$$.

Now using these information you can go through the info on the datasheet.

NOTE: The triangular waveform above is given just as an example for the sake of simplicity i.e. to give an idea. Since the pulse waveform for your experiment will be exponential, the voltage will approach zero asymptotically. As Mattman944 stated in their comment, you can take 3 to 5 times the time constant ($$\\mathrm{\tau=RC}\$$) for the pulse duration.

• OP said it was a discharging capacitor curve, AKA exponential. So, the pulse duration is not obvious. I would use 3 time constants, but it is a judgement call. Sep 26 at 15:17
• @Mattman944 I know. Triangular pulse was just as an example (I should have indicated this, edited my answer accordingly) as it was rather simple for me to show integrals. Sep 26 at 15:24
• Would this link for a pulse-energy-calculator be accurate as an approximation for the discharge curve Pulse Energy? vishay.com/en/resistors/pulse-energy-calculator-noheader Sep 28 at 13:01
• @PDK11 the calculators in the link takes the integral I explained above for you and gives you the pulse energy (haven't checked if the calculators do the job correctly, though). So, once you get the pulse energy (E) from the calculator you can calculate the rectangular pulse duration, Tx, with the same peak voltage/power i.e. E = P.Tx. Sep 28 at 13:52

This is something of a non-answer, alas; but also the best possible answer?

That's a bit contradictory, so let me explain.

When the resistor is thermally limited, the problem is to avoid exceeding maximum peak rated temperature at any point within the device.

Already we see a problem here, because resistors aren't rated by internal temperatures (how would you figure that out, anyway?). If we had the means to measure every point within, that would be one thing; but since we can't, what then?

The next best thing is to simulate the process -- make a thermal model, following reasonable assumptions about materials and geometry. If we can extract such a model, we can simulate any arbitrary pulse ourselves, and determine if it's within ratings or not.

It's possible that a manufacturer might possess such data. Unfortunately, I haven't come across any myself, even for familiar-packaged parts (like TO-220 thin films, which might be particularly easy and worthwhile to characterize, if anything?). So, lacking data to start with... this isn't a very good answer.

It's also possible we can take the pulse rating (energy or power) curves, and turn them into expected, or limiting, transient thermal impedance curves; but I don't have a strong feel on this and would have to do some work to come up with a method. Anyway, it's likely that such curves are themselves guardbanded (i.e. derated from ultimate device capability) by some margin, by the manufacturer already, given that they probably don't have transient thermal resistance data themselves. (But, they might, and just keep it locked away for internal use, or privileged customers; you'd have to ask.)

## Using Transient Thermal Impedance

If you do have a thermal impedance network -- simply set up your power waveform (typically as a current representing power, for voltage as temperature) into the impedance network, and run the simulation.

This is a standard method for semiconductors, which generally do provide transient impedance characteristics (if not necessarily the network parameters, so you may have to do some curve fitting yourself still), and generally do have a single point at which to measure internal temperature (usually by characterizing VGS(th) or body diode VF vs. T).

Here's an exercise where I derived ZθJA for a D2PAK device: ITH is the thermal source, current corresponding to power at the transistor; CTHn are the thermal capacitances or masses in the system, and RTHn are the resistances connecting them. Voltage corresponds to temperature, and VTA is the voltage corresponding to a temperature of ambient (which I don't care about here, so it's set to zero).

As a plotting aid, the E source buffers the junction temperature TJ, rescales it, and M1 is a log function, so that TJNL matches the scale on the datasheet plot. (I can set the time axis to display logarithmic in the viewer; I could set the same vertical functions in the viewer, actually, but this shows them off explicitly.) I overlay a capture from the datasheet using OnTopReplica, and compare simulated to measured until they match: Tweaking the RC values is a bit tedious, but as you can see, the 4th order fit is quite good. Well, it's a bit hard to see what all is going on here, there's transparency everywhere, but the on-top window is set to 50% transparency hence the slight outline of the window, and the datasheet plot extends further leftward but I'm only modeling the rightmost curve (which is ZθJA i.e. the device on a real circuit board in ambient conditions; the other is ZθJC i.e. for the case joined to an infinite heatsink). The background matter, "Zth Sim.sdf", the vertical (V) axis, etc., is the simulation output view in Altium.

This was done to simulate a device under pulsed conditions. Here is the representative schematic: IC and VCE are measured by E and H sources, multiplied by A1 to get voltage representing power, and converted to a current representing power into the simulation network. An integrator also records total energy (here, a bit over 6J after 400ms).

The simulation results: Note that TJ reads in units of kelvin, temperature rise above ambient. During this pulse (an automotive load dump transient), peak power greatly exceeds the continuous rating of the device (a couple watts for a D2PAK on board with no additional heatsinking), but falls well within its headline rating (over 100W); somewhere inbetween, there is a pulse amplitude and duration that will survive, and the challenge was to determine where. As you can see, the immediate heating is strong, but not destructive, 45K -- and the application was approved for production.

In comparison, a simulated square pulse, of peak amplitude and half width (38W for 205ms), overestimates the temperature rise as 59.5 instead of 44K. A triangular approximation (10ms rise, 5ms peak dwell at 38W, 400ms linear ramp down to zero) does however get pretty close at 45.5K.

This pulse is a bit different from an exponential capacitor discharge (it's a truncated exponential), but a triangle width/duration of about 1.5 time constants should get pretty close. (Calculating the exact ratio would be a good exercise for the proverbial student. :) )