YouTube精彩教程分享~
里面涉及到的in文件等信息可以点击下方阅读原文进入官方网站获取（打不开的同学可能需要翻墙）~
in文件内容：
###################################################################################################################setupMDcalculationtocalibratedislocationvelocityasafunctionofappliedstressatfinitetemperature##SG-August14th,##SG-testedonAugust15th,##reference:Grohetal,IJP##TD-RevisedforNickel.March12,##BH-Revisedforageneralmaterial.January13,##BH-Revisedfordifferentpost-processingmethod.February15,#########################################################################################################################################################################################################################################Edgedislocation#1.statictoobtainthecorestructure#2.equilibratetemperatureatK(nve+rescaling)#3.applyaconstantforcetocalculatethedislocationvelocity.Forceisappliedonarigidbody#####################################################################################################################################################################################################################################Usersetvariables#initTemp:desiredtemperature#sigma:shearstressinbarunits#material:materialname#atom_file:nameofthefilecontainingtheatompositions##Miscellaneouscontrolvariables#equilTime:Numberofincrementtoequilibratethetemperature#runTime:Numberofincrementtocalibratethetemperature#energyConv:EnergyconversiontermfromeVtoJtoconvertatomicforcetostress##output:#dump.minimize:atomsbeforeandafterminimization#dump.equilibration:atomsduringtemperatureequilibration#dump.shear:wrappedatomsinthesimulationbox(perfecttolookatthevelocity)#dump.shear.unwrap:unwrappedatoms(tocheckthateverythingisgoingwell)#outputfile(log.lammps):monitorthetemperatureduringtheMDrun####################################################################################################################VariabledefinitionsvariableinitTempequal.#desiredtemperaturevariablesigmaequal.#appliedstressinbarvariablematerialstringTa#materialsymbolvariableatom_filestringatoms.bcc.edge.pad#theconfigurationwasgeneratedbySGwiththepreprocessordislocation.f90variableequilTimeequal#numberofincrementtoequilibratethetemperaturevariablerunTimeequal0#numberofincrementtocalibratethevelocityvariableenergyConvequal.7#conversionfactordimension3boundarypspunitsmetalatom_styleatomicread_data${atom_file}#storeinitialpositionofbottomandtopplanesalongyvariabletmp0equal"ylo+14."variableylo0equal${tmp0}variabletmp1equal"yhi-14.5"variableyhi0equal${tmp1}#variablefordumpingvariableymid1equal"0.5*ylo+0.5*yhi-30.0"variableymiddlenegativeequal${ymid1}variableymid2equal"0.5*(ylo+yhi)+30.0"variableymiddlepositiveequal${ymid2}#definepotentialpair_stylemeampair_coeff**${material}.library.meam${material}${material}.parameter.meam${material}neighbor2.0binneigh_modifydelay5#definitionoftheupperandlowerblocksregionupperblockINFINF${yhi0}INFINFINFunitsboxregionlowerblockINFINFINF${ylo0}INFINFunitsbox#definitionofthegroupgroupupperregionuppergrouplowerregionlowergroupmobilesubtractallupperlower#fixtopandbottongroupfix1lowersetforce0.00.0NULLfix2uppersetforce0.00.0NULL#