The easiest way of solving such problems is to use complex phasors. The total complex impedance is
$$Z_T=R_1+R_2||j\omega L=R_1+\frac{j\omega R_2L}{R_2+j\omega L}=877.18\cdot e^{30.67^{\circ}}\Omega$$
with \$\omega=2\pi\cdot 10,000 Hz\$, \$R_1=470\Omega\$, \$R_2=988\Omega\$, and \$L=10mH\$ (as shown in your circuit diagram).
This gives for the total current
$$I_s=\frac{V_0}{Z_T}=9.12\cdot e^{-30.67^{\circ}} mA$$
with \$V_0=8V\$.
The voltage across \$R_1\$ is
$$V_1=I_sR_1=4.28\cdot e^{-30.67^{\circ}} V$$
The voltage across \$R_2\$ and \$L\$ is
$$V_2=I_s(R_2||j\omega L)=I_s\frac{j\omega R_2L}{R_2+j\omega L}=4.83\cdot e^{26.88^{\circ}} V\quad(=V_0-V_1)$$
The current through \$R_2\$ is
$$I_2=\frac{V_2}{R_2}=4.89\cdot e^{26.88^{\circ}} mA$$
The current through the inductor is
$$I_L=\frac{V_2}{j\omega L}=7.70\cdot e^{-63.12^{\circ}} mA\quad(=I_s-I_2)$$
I realize that there are quite a few differences between these results and your calculations and measurements. I just used the values you gave in the circuit diagram, and I did my best to avoid any errors.
Any help would be appreciated
You can measure phase-neutral voltage and if the supply and load are balanced then you can infer line voltage by: -
line volts = phase volts\$\times\sqrt3\$
This should be one of the easier measurements yet you say you can calculate "real power, reactive power and apparent with ease". This does make me think that you have used one of these measurements, and the RMS measurement of current and back-calculated phase/line voltage and that's when you are seeing a discrepency.
If this is so then I suspect your current measurement may be at fault either through incorrect use of a current transformer or some ratio being incorrect. Or, it could be the \$\sqrt3\$ thing mentioned above.
If you need any further help with this you should consider detailing how you make the other measurements and why you believe voltage to be incorrectly calculated/measured.
Best Answer
Start by minimizing the problem to firstly this: -
Then, because the negative half cycle is a mirror image of the positive half cycle, it will have the same RMS value, thus we end up with only needing to analyse this section to get RMS: -
Then it breaks down into 4 sections that you need to find the square of the RMS for each: -
So, find those individual squares then multiply each by their time duration and the final RMS value of your waveform is this: -
$$\sqrt{\dfrac{d_1V_1^2 + d_2V_2^2 + d_3V_3^2 + d_4V_4^2}{d_1+d_2+d_3+d_4}}$$