【LAMMPS学习】八、基础知识(5.9)LAMMPS 近场动力学

8. 基础知识

此部分描述了如何使用 LAMMPS 为用户和开发人员执行各种任务。术语表页面还列出了 MD 术语,以及相应 LAMMPS 手册页的链接。 LAMMPS 源代码分发的 examples 目录中包含的示例输入脚本以及示例脚本页面上突出显示的示例输入脚本还展示了如何设置和运行各种模拟。

8.1.通用基础知识

8.2. 设置入门

8.3. 分析入门

8.4. 力场入门

8.5. 软件包入门

8.5.1.有限尺寸球形和非球形粒子

8.5.2. 粒度模型

8.5.3.体粒子

8.5.4.黏合粒子(BMP)模型

8.5.5.极化模型

8.5.6. 绝热核/壳模型

8.5.7.Drude感应偶极子 

8.5.8. LAMMPS 中热化 Drude 振荡器教程

8.5.9. LAMMPS 近场动力学

本指南基于 Michael L. Parks、Pablo Seleson、Steven J. Plimpton、Richard B. Lehoucq 和 Stewart A. Silling 撰写的桑迪亚报告 2010-5549。


概述

近场动力学是经典连续介质力学的非局域扩展。离散近场动力学模型具有与分子动力学模型相同的计算结构。本 Howto 简要概述了连续体的近场动力学模型,然后讨论了如何在 LAMMPS 中离散化近场动力学模型,如原始文章 (Parks) 中所述。还包括一个带有注释的示例问题。


快速入门

近场动力学样式包含在可选的 PERI package包中。如果您的 LAMMPS 可执行文件尚未包含 PERI 包,您可以查看 PERI packagebuild instructions for packagesPERI package,了解如何在从源代码编译自定义版本的 LAMMPS 时启用该包。

这是设置近场动力学模拟的最小示例。

units         si
boundary      s s s
lattice       sc 0.0005
atom_style    peri
atom_modify   map array
neighbor      0.0010 bin
region        target cylinder y 0.0 0.0 0.0050 -0.0050 0.0 units box
create_box    1 target
create_atoms  1 region targetpair_style    peri/pmb
pair_coeff    * * 1.6863e22 0.0015001 0.0005 0.25
set           group all density 2200
set           group all volume 1.25e-10
velocity      all set 0.0 0.0 0.0 sum no units box
fix           1 all nve
compute       1 all damage/atom
timestep      1.0e-7

关于此输入示例的一些注释:

  • 近场动力学模拟通常使用 SI 单位

  • 粒子必须在简单的立方晶格上创建

  • 需要使用 atom style peri

  • 索引粒子需要原子图

  • 计算邻居列表时使用的 skin distance应根据您选择的模拟参数进行适当定义。 skin distance应设置为一个值,使得近场动力学地平线加上  skin distance大于两个键合粒子之间的最大可能距离(在它们的键断裂之前)。这里设置为0.001米。

  • 需要peridynamics pair style 。目前可用的选项有:peri/eps、peri/lps、peri/pmb 和 peri/ves。模型参数使用pair_coeff 命令设置。

  • 必须定义每个颗粒的质量密度和体积分数。这是通过密度和体积的两个设置命令来完成的。对于简单的立方晶格,粒子的体积应等于晶格常数的立方,此处为 𝑉𝑖=Δ𝑥3 。

  • 使用速度命令时,所有粒子最初都处于静止状态

  • 使用普通速度-Verlet时间积分器,其在代数上等价于中心时间差,但在数值上更稳定

  • 您可以使用计算损伤/原子来计算每个粒子位置的损伤

  • 最后,使用 timestep 命令将时间步设置为 0.1 微秒。


连续体的近场动力学模型

以下并不是近场动力学的完整概述,而是仅讨论我们在 LAMMPS 中实现的模型特有的细节。有关近场动力学理论的更多信息,读者可以参考(Silling 2007)。首先,我们定义将使用的符号。

基本符号

在近场动力学文献中,通常使用以下符号约定。参考配置中给定点的位置是 x 。让 𝑢(𝑥,𝑡) 和 𝑦(𝑥,𝑡) 分别表示点 𝑥 在时间 𝑡 的位移和位置。将两个粘合点 x 和 x′ 的相对位置和位移向量分别定义为 𝜉=x′−x 和 𝜂=u(x′,𝑡)−u(x,𝑡) 。我们在这里注意到 𝜂 与时间相关,而 𝜉 则不然。由此可见,当前配置中两个绑定点的相对位置可以写为 𝜉+𝜂=𝑦(𝑥′,𝑡)−𝑦(𝑥,𝑡) 

近场动力学模型经常使用状态来编写,我们在这里简要描述。出于我们讨论的目的,所有状态都是作用于 𝑅3 中向量的运算符。有关状态的更完整讨论,请参阅(Silling 2007)。矢量状态是其图像是矢量的算子,并且可以被视为二阶张量的推广。类似地,标量状态是一个其图像是标量的算子。特别令人感兴趣的是矢量力状态 \underline{\mathbf{T}}[\mathbf{x}, t]\left\langle\mathbf{x}^{\prime}-\mathbf{x}\right\rangle ,它是矢量 𝑥′−𝑥 到力矢量状态场的映射,单位为每体积平方的力。向量状态运算符 \underline{\mathbf{T}}本身可以是 𝑥 和 𝑡 的函数。本构模型完全包含在 \underline{\mathbf{T}}内。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/833531.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Ubuntu软件中心不显示

装完Ubuntu后没有Software -- 更新apt sudo apt update -- 升级apt sudo apt upgrade -- 重启 sudo systemctl reboot-- 安装snap sudo apt-get install snap -- 安装软件商店 sudo snap install snap-store -- 更新软件商店 sudo snap refresh snap-store安装成功&#xff01…

学QT的第二天~

小黑子鉴别界面 #include "mywidget.h" void MyWidget::bth1() { if(edit3 ->text()"520cxk"&&edit4 ->text()"1314520") { qDebug()<< "你好&#xff0c;真爱粉"; this->close(); } else { speecher->sa…

本地加载hugging face模型:Bert

找了个hf的镜像站&#xff0c;把config.json和pytorch_model.bin两个文件进行下载下来&#xff0c;模型文件uncased_L-12_H-768_A-12.zip下载下来先。 解压模型文件压缩包&#xff0c;把前面下载的两个文件也放进去&#xff0c;总共6个文件。这个文件夹就是代码 tokenizer Be…

【建网护网三十载】 守护不息创新不止,C3安全AI未来!

30年&#xff0c;中国互联网从起步探索到领先全球。1994年4月20日&#xff0c;中国正式开通首条64K的国际专线&#xff0c;标志着我国成功实现与国际互联网的全功能接轨&#xff0c;展开互联网快速发展的三十载。 回望30年&#xff0c;亲历建网&#xff0c;投身建设&#xff0c…

Charles抓包工具

Charles是一个HTTP代理工具&#xff0c;使开发人员能够查看客服端和服务器之间的所有HTTP/ HTTPS/SSL网络请求。 Charles是在PC环境下常用的网络抓包截取工具&#xff0c;在做移动开发时&#xff0c;我们为了调试客户端与服务端的网络通讯协议&#xff0c;常常需要截取网络请求…

【ITK配准】第五期 变换(Transform)(三)

很高兴在雪易的CSDN遇见你 VTK技术爱好者 QQ:870202403 公众号:VTK忠粉 前言 本文分享用于ITK配准中的三维变换(Transform),希望对各位小伙伴有所帮助! 感谢各位小伙伴的点赞+关注,小易会继续努力分享,一起进步! 你的点赞就是我的动力(^U^)ノ~YO 目录 …

大数据集成平台建设方案-word原件资料

基础支撑平台主要承担系统总体架构与各个应用子系统的交互&#xff0c;第三方系统与总体架构的交互。需要满足内部业务在该平台的基础上&#xff0c;实现平台对于子系统的可扩展性。基于以上分析对基础支撑平台&#xff0c;提出了以下要求&#xff1a; (1) 基于平台的基础架构&…

iframe通信postMessage

iframe嵌入页面是实现微前端的方式之一。由于浏览器的跨域限制&#xff0c;iframe与父页面之间的通信变得不是那么容易。postMessage解决了这个问题。从广义上讲&#xff0c;一个窗口可以获得对另一个窗口的引用&#xff08;比如 targetWindow window.parent&#xff09;&…

Mysql:交叉连接、内连接

打开Navicat Premium&#xff0c;新建数据库&#xff0c;创建两个表dept和emp&#xff1a; create table dept(did int(4) primary key,dname varchar(36)); create table emp(id int(4) primary key,name varchar(36),age int(2),did int(4) not null); 插入数据&#xff1a…

spring bean生命周期全部过程

Spring Bean的生命周期包括以下全部过程&#xff1a; 实例化&#xff1a;在Spring容器启动时&#xff0c;根据配置文件或注解等信息创建Bean的实例。属性赋值&#xff1a;如果Bean有属性需要进行初始化&#xff0c;Spring容器会自动为这些属性进行赋值。自定义初始化方法&…

线程理论篇1

本章问题&#xff1a;什么是线程?线程的使用场景&#xff1f;什么是线程池&#xff1f;线程池是如何工作的&#xff1f;线程池共享了哪些资源?线程安全代码怎么写&#xff1f;什么是线程安全? 什么是线程&#xff1f; 线程是为了提高进程的效率。进程的地址空间中保存了cpu…

【mysql】mysql单表查询、多表查询、分组查询、子查询等案例详细解析

✨✨ 欢迎大家来到景天科技苑✨✨ &#x1f388;&#x1f388; 养成好习惯&#xff0c;先赞后看哦~&#x1f388;&#x1f388; &#x1f3c6; 作者简介&#xff1a;景天科技苑 &#x1f3c6;《头衔》&#xff1a;大厂架构师&#xff0c;华为云开发者社区专家博主&#xff0c;…

嵌入式C语言高级教程:实现基于STM32的人工智能语音识别系统

在嵌入式系统中实现语音识别技术可以极大地增强设备的交互性。本教程将指导您如何在STM32微控制器上使用TensorFlow Lite for Microcontrollers实现基本的语音识别功能。 一、开发环境准备 硬件要求 微控制器&#xff1a;STM32F746NG&#xff0c;支持足够的运算能力和内存来…

spring高级篇(七)

1、异常处理 在DispatcherServlet中&#xff0c;doDispatch(HttpServletRequest request, HttpServletResponse response) 方法用于进行任务处理&#xff1a; 在捕获到异常后没有立刻进行处理&#xff0c;而是先用一个局部变量dispatchException进行记录&#xff0c;然后统一由…

VMware虚拟机中ubuntu使用记录(5)—— 如何在ubuntu中安装USB相机ros驱动并获取usb摄像头数据

提示&#xff1a;文章写完后&#xff0c;目录可以自动生成&#xff0c;如何生成可参考右边的帮助文档 文章目录 一、ROS下USB相机驱动1.准备工作(1) 下载驱动(2) 创建ROS工作空间 2. 安装usb_cam驱动(1) 安装usb_cam驱动包(2) 编译代码 3. 修改usb_cam驱动的配置文件(1) 查看US…

掌握TypeScript的非空断言(!)和可选链(?):开发效率翻倍!

引言 标题&#xff1a;掌握TypeScript的非空断言和可选链&#xff1a;开发效率翻倍&#xff01;简短介绍&#xff1a;在TypeScript中&#xff0c;?和!操作符是提高代码安全性和开发效率的强大工具。本文将为你揭示它们的使用方式和最佳实践。 背景知识 易于理解的解释&…

vue2+@tinymce/tinymce-vue富文本编辑器的实现

1.使用的第三方的库 npm install @tinymce/tinymce-vue@3.2.8 tinymce@5.10.9 --save 2.引入应用,直接上代码,支持上传图片,视频 <template><div class="tinymce-container" ref="tinymceContainer"><Eidtor api-key="" v-if…

一路串联电机的绕制原理

这里要说明的一点是 对于一路串联的电机&#xff0c;无论是一把线圈还是两把线圈&#xff0c;出来的都是只有两个线头&#xff0c;可看做一个整体来对待&#xff01; 绕制具体原理 同心式线圈绕制 前面说的都是等距式的 线圈绕制&#xff0c;下面我们讲解一下同心式的绕制办法…

numpy学习笔记8

对不同列标识不同信息进行数据读取时 步骤&#xff1a; 1.根据不同列的数据标识的含义和类型&#xff0c;因此需要自定义数据类型 2.使用自定义的数据类型读取数据 有点类似于数据库读取数据操作 例子: # 1. 以上数据由于不同列数据标识的含义和类型不同,因此需要自定义数据类…

整体意义的构成与构建

输入数据中不同部分之间的复杂关系是指数据中各个元素如何相互作用和相互影响以构成整体意义的能力。 理解上述观点可以从以下几个方面进行&#xff1a; 元素的相互作用&#xff1a;在任何给定的数据集中&#xff0c;单个数据元素往往不是孤立存在的。它们之间存在某种形式的交…