Dear Shrushti,
it is straightforward to deal with strained graphene. 
In the region you want to make strain, just set the function which modifies the site position (this mainly depends on the amount and type of the strain). Second, once you have shifted the site position you need to shift the hopping as illustrated in the file provided by Mr. Antonio.
here you find a simple script about how to shift the site position in the case of uniaxial strain (adjusted to you case)
def Triaxial_transform(pos, center, I):
    """ I is a parameter that lets us to control the intensity of the strain"""
    x, y = pos  
    cx, cy = center[0], center[1] # the center of the strain
if r<sigma: ux = 2*I * x*y uy = I * (x**2 - y**2) return x + ux, y + uy else: return x, y
Then use this position to set your hopping as illustrated in a previous discussion [Ref-1].

You have to be careful in case you are using a strain in the z axis since you will not be able to plot the current.  It is somehow tricky in this scenario.


I hop this will help
Best wishes

Le jeu. 15 juil. 2021 à 14:19, <> a écrit :
I want to create the graphene strained superlattice-like structure, having uniaxial strain defined at the specified region. Please, someone can suggest to me how to define the strained region and interface between unstrained and strained graphene regions. 
Thanks in Advance