Clumped beam with cubic thermal load

A beam fixed at both its ends is subject to one-step cubic thermal load:

T=50(x23+x22+x2+1)T = 50(x_2^3+x_2^2+x_2+1)

The material properties are taken to be E = 10,000 MPa,
v = 0.3 and alpha = 0.00001 oC-1.

The numerical results are compared with the reference solution.

Reference

de Medeiros et al. (2004) “The method of fundamental solutions with dual reciprocity for some problems in elasticity”