The piezoelectric semiconductor combines the characteristics of both piezoelectric and semiconducting materials, making it highly promising for a variety of engineering applications. The present work adopts a Galerkin approach to derive a general weak form of the coupled governing equations for the analysis of piezoelectric semiconductors. The coupling effect between the strain, the electric field, and the gradients of the concentrations of holes and electrons has been considered. Piezoelectric semiconductors are often deployed in the shape of rods and beams and therefore commonly experience axial and transverse load. Thus, the Euler–Bernoulli beam theory is employed to formulate a numerical approach to analyse the coupling behaviour of such piezoelectric semiconductor beams. In order to satisfy the C1 continuity requirement of Euler–Bernoulli beam theory without introducing additional degrees of freedom, a discretisation with B-Splines is used. Non-mechanical variables, including electrical potential, holes and electrons concentrations, are decomposed into constant, linear and quadratic terms through series expansion along the thickness, which ensures compatibility with beam formulation. To validate the proposed method, an N-type piezoelectric semiconducting beam under a sinusoidal distributed load, for which an analytical solution is available, is initially examined. Subsequently, complex examples including a cantilever beam with generalised loading conditions and a beam with PN-junction are analysed to demonstrate the capabilities of the proposed method in handling more challenging scenarios.