Dear Lammps users,
I used the below command to calculate the local stress of polymer.
fix SS all ave/correlate $s $p $d &
v_pxy v_pyz v_pxz type auto file profile.gk ave running
variable scale equal 1.0/t/{kb}vol$s*dt
variable vxy equal trap(f_SS[3])*${scale}
variable vyz equal trap(f_SS[4])*${scale}
variable vxz equal trap(f_SS[5])*${scale}
thermo_style custom step temp press pxy pyz pxz v_vxy v_vyz v_vxz
run 10000000
variable v equal (v_vxy+v_vyz+v_vxz)/3.0
print “running average viscosity: $v”
How can I calculate the normal stress. Would you give me some suggestion? Thank you for your help.
Best Regards,
Gerry Wong
发送自 Windows 10 版邮件应用