Here is a solution using complex numbers geometry.
Let us "install" the given figure into a coordinate system with origin in the center of the pentagon and main vertex $A_0$ in $(1,0)$ as shown on Fig. 1, where :
$$b:=e^{2i \pi/40}\tag{1}$$
The polar angle of $b^3$ is $3 \pi/20 \ ( = 27°)$, the angular value we should obtain for $\angle BA_0C$.

Fig. 1.
Let us convert point $B$ into its complex form. Let us use for that relationship $\vec{A_2B}=R(\vec{A_2A_1})$ where $R$ is the $+\pi/2$ rotation :
$$B \leftrightarrow \beta:=\underbrace{b^{16}}_{A_2} + \underbrace{i}_R\underbrace{(b^8-b^{16})}_{A_1-A_2}$$
Or, under a different form :
$$B \longleftrightarrow \beta= b^{16} + b^{10}(b^8-b^{16}) = b^{16}+b^{18}+b^6$$
Remark : we have used various relationships such as
$$b^{10}=i, \ \ b^{20}=-1, \ \ b^{40}=1\tag{2}$$
For similar reasons, relationship $\vec{A_2C}=R^{-1}(\vec{A_2A_3})$ gives :
$$C \longleftrightarrow \gamma:=b^{16} - i(b^{24}-b^{16})$$
With simplifications :
$$C \longleftrightarrow \gamma= b^{16}+b^{14}-b^6$$
Therefore, vectors $\vec{A_0B}$ and $\vec{A_0C}$ have the complex representation :
$$\begin{cases}\vec{A_0B} \longleftrightarrow (\beta-1)\\\vec{A_0C} \longleftrightarrow (\gamma-1)\end{cases}$$
We still have to express with complex numbers that angle $(\vec{A_0B},\vec{A_0C})$ has value $3 \pi/20$ :
$$(\beta-1)b^3 \stackrel{?}{=} k(\gamma - 1)\tag{3}$$
where $k$ is a real number.
Otherwise said, we have to show that complex number
$$z:=\frac{(\beta-1)b^3}{(\gamma - 1)} \stackrel{?}{\in} \mathbb{R} \tag{4}$$
First solution : condition above can be expressed in the following way :
$$ z \stackrel{?}{=} \bar{z}\tag{5}$$
Let us give explicit expressions for $z$ and $\bar{z}$:
$$\begin{cases}z&=&\frac{(b^{16}+b^{18}+b^6-1)b^3}{(b^{16}+b^{14}-b^{6}-1)}&=&\frac{N}{D}\\
\bar{z}&=&\frac{(-b+b^{19}-b^{11}+b^{17})}{(-b^4-b^6+b^{14}-1)}&=&\frac{N'}{D'}\end{cases}$$
Remark : the expression of $\bar{z}$ is obtained by taking the conjugation of all the terms in $z$ knowing that the conjugate of $b^k$ is $b^{40-k}=-b^{20-k}$.
Relationship (4) amounts to the equivalent question for the following polynomial $$\underbrace{ND'-N'D}_{P(b)} \stackrel{?}{=} 0.\tag{6}$$
Expanding the LHS of (6):
$$P(b)=-b^{47} - b^{45}- b^{35} - b^{33} - b^{25} + b^{23} - b^{15} - b^{13} + b^7 + b^3 \stackrel{?}{=}0$$
which is indeed true due to a remarkable cancellation of all its terms ($-b^{47}$ with $b^7$, $-b^{45}$ with $- b^{25}$, etc.) (refer to formulas (2) above) $\square$
Second solution
Based on a simplification of the expression of $z$ using SAGE (a public domain Computer Algebra System) with the following code :
x = polygen(ZZ, 'x')
mp = x^16 - x^12 + x^8 -x^4 + 1
K.<b> = NumberField(mp)
z = (b^18+b^16+b^6-1)*b^3/(b^16+b^14-b^6-1)
z = z.coordinates_in_terms_of_powers()
print(q)
where the field extension $K.<b>$ (alias $\mathbb{Q}[b]$) is classically described using the minimal polynomial of $b$, i.e., $mp(x) = x^{16} - x^{12} + x^8 - x^4 + 1$ (see remark 1 below).
Running this program gives the result :
$$z= -b^{13} + b^7$$
which is indeed a real number :
$$z=-2i b^{10}\left( \frac{b^{3}-b^{-3}}{2i}\right)=2 \sin(3 \pi/40) \color{red}{\in \mathbb{R}}$$
Remarks :
As $b$ is a primitive $40$th root of unity, its minimal polynomial is the so-called cyclotomic polynomial of order $40$, denoted $\Phi_{40}(x)$. Its coefficients can be obtained from the SAGE instructions below :
from sage.rings.polynomial.cyclotomic import cyclotomic_coeffs
print(cyclotomic_coeffs(40))
This second solution has been inspired by Steven Stadnicki's answer.
PM 2Ring has given a similar SAGE program that I find excellent (more concise than mine) ; it can be uploaded here. I reproduce it here:
K.<b> = CyclotomicField(40)
print(K)
print(K.defining_polynomial())
z = (b^18+b^16+b^6-1)*b^3/(b^16+b^14-b^6-1)
print(z)
q = SR(z):#type conversion to complex
print(q)
print(q.real(), q.imag())