Questions about setting quasi-static displacement load on the boundary

Hello,

I’m trying to set the a displacement load u_x = 1 on the right surface of a cubic domain, the problem is that I need to apply the load step by step (u_x=steps*0.01, steps=100) until u_x= 1. then I plan to use the quasi-static scheme to implement.

I looked through the tutorial https://jsdokken.com/dolfinx-tutorial/chapter2/hyperelasticity.html, but this part code refers to the traction applying on it

tval0 = -1.5
for n in range(1, 10):
    T.value[2] = n * tval0
    num_its, converged = solver.solve(u)
    assert(converged)
    u.x.scatter_forward()

Could some one help me on this? thanks!