Sorry for the self publicity, but for the especific example of $gl_n(k)$, you can view it as $gl_n(k)\cong sl_n(k)\times k$ and we did some work for trivial central extensions that allow you to produce, in particular, all Lie bialgebra extensions on $gl_n$ from data in $sl_n$, independently of the fact that they are comming from Belavin-Drinfel'd classification (they did the complex factorizable case) or not.
If you denote $t$ is a generator of the abelian factor in $gl_n$, then all Lie bialgebra structures on $gl_n$ are of the form
$$\delta (t)=0$$
$$\delta (x)=\delta_{sl_n}(x)+D(x)\wedge t$$
($x\in sl_n$) where $\delta_{sl_n}$ is a Lie bialgebra structure on $sl_n$ and $D$ is a derivation of $sl_n$ that is also a coderivation w.r.t. $\delta_{sl_n}$. Equivalently,
$$D(x)=[H,x]$$
for some $H\in sl_n$ such that $\delta_{sl_n}(H)=0$. Also equivalently, if $\delta_{sl_n}(x)=ad_x(r)$, then $H$ must satisfy $ad_H(r)=0$.
For example, you can take $\delta_{sl_n}\equiv 0$ and pick $H_0\in sl_n$, then
$$\delta(t)=0$$
$$\delta(x)=[H_0,x]\wedge t$$
is a Lie bialgebra structure on $gl_n$.
If $\delta_{sl_n}$ is one of the structures comming from Belavin-Drinfel'd list, then $D=[H,-]$ with $H$ in the same Cartan subalgebra from Belavin-Drinfel'd statement, with some additional conditions. You can take a look at https://arxiv.org/pdf/1110.1072.pdf
We also did some work when the underlying Lie algebra is 2-step nilpotent, generalizing the results for the Heisenberg. https://arxiv.org/pdf/1607.00300.pdf (by the way, the results for the Heisenberg were published, you can also find references there)