fix 5 all deform 1000x erate 0.002 units box #設置沿x方向拉伸，應變率為0.002（1/ps）
fix ID group-ID deform N parameter args ... keyword value ...
N是每N步實施一次變形，如果這個數越大，那么每N步變形就越大，如果這個數較小，比如說是1，那么每跑一步就變形一次，這樣模擬的更細，精度更高。總變形是由后面的參數確定的，N影響不了總變形。
fix 1 all deform 100 x erate 0.0001 units box remap x
run 40000
deform后面的N為1或100的區別在於，N為1時，是每運行1步，就變形1次；而N為100時，就是第一步變形了，第2-100步沒有變形，第101步變形，第102-200步沒有變形。
0.0001：應變速率為0.0001.
總應變是由erate=0.0001/ps、timestep=0.001ps、run 300000步三者決定，總應變為0.03，與deform后的N值無關！
upto關鍵字表示從當前時間步長開始直到指定的時間步長執行運行。例如。如果當前時間步長為10,000並且使用“運行100000 upto”，則將運行額外的90,000個時間步長。
（1）N=1時，意味着每運行一步，就變形一次，總共運行300000步要達到變形量0.03，也就是說每次運行的應變為10^-7（0.03/300000），即每運行一步，應變增加10^-7，運行到10000步時，總應變為10^-3;則在step=1~10000的過程中，總應變依次為10^-7,2*10^-7,3*10^-7.........,10^-3
(2)N =10000時，意味着每運行10000步，就變形一次，總共運行300000步要達到變形量0.03，亦即在每個10000th變形一次，且應變為10^-3;;則在step=1~10000的過程中，總應變依次為0,0 0,...，0,10^-3，step達到10000時，應變一下子就由0增加大10^-3。
上述分析說明，實際上N=1和N=10000時，模型的實際變形速率確實是不一致的，N=1是每一步變形，變形速率為10^-4/ps; N=10000時在10000步之前都保持0應變，而在每個10000 th時，直接變形達到規定值,變形速率為1/ps。N=10000時的實際變形速率是N=1的10000倍。
erate value = R
R = engineering strain rate (1/time units)
trate value = R
R = true strain rate (1/time units)
- change in length = deltaL = r * time
- engineering strain = deltaL/L = r * time/L
- engineering strain rate = strain per time = r/L
- true strain = ln(1+ engineering strain) = ln(1+ r*time/L)
- L = length of specimen in loaded direction; r = rate of loading
- true (effective) strain rate varies with time = d(true strain)/dt = [ln(1+r*time2/L) - ln(1+r*time1/L)]/(time2-time1)
- 剪切應變是無單位的，定義為偏移量/長度，其中長度是垂直於剪切方向的盒子長度（例如，用於xy變形的y盒子長度），偏移量是剪切方向上的位移距離（例如，用於xy變形的x方向） 從不受約束的方向。