24小时热门版块排行榜    

Znn3bq.jpeg
查看: 5361  |  回复: 5

束卷卷

新虫 (初入文坛)

[求助] 剪切应力应变曲线怎么输出

根据lammps自带的in.shear剪切模拟,进行模拟后,怎么输出应力应变曲线啊,或者说剪切应力应变曲线在in文件中写的啊,求各位知道,谢谢!
# 3d metal shear simulation

units                metal
boundary        s s p

atom_style        atomic
lattice                fcc 3.52
region                box block 0 16.0 0 10.0 0 2.828427
create_box        3 box

lattice                fcc 3.52 orient        x 1 0 0 orient y 0 1 1 orient z 0 -1 1 origin 0.5 0 0
create_atoms        1 box

pair_style        eam
pair_coeff        * * Ni_u3.eam

neighbor        0.3 bin
neigh_modify        delay 5

region                lower block INF INF INF 0.9 INF INF
region                upper block INF INF 6.1 INF INF INF
group                lower region lower
group                upper region upper
group                boundary union lower upper
group                mobile subtract all boundary

set                group lower type 2
set                group upper type 3

# void

#region                void cylinder z 8 5 2.5 INF INF
#delete_atoms        region void

# temp controllers

compute                new3d mobile temp
compute                new2d mobile temp/partial 0 1 1

# equilibrate

velocity        mobile create 300.0 5812775 temp new3d
fix                1 all nve
fix                2 boundary setforce 0.0 0.0 0.0

fix                3 mobile temp/rescale 10 300.0 300.0 10.0 1.0
fix_modify        3 temp new3d

thermo                25
thermo_modify        temp new3d

timestep        0.001
run                100

# shear

velocity        upper set 1.0 0 0
velocity        mobile ramp vx 0.0 1.0 y 1.4 8.6 sum yes

unfix                3
fix                3 mobile temp/rescale 10 300.0 300.0 10.0 1.0
fix_modify        3 temp new2d

#dump                1 all atom 100 dump.shear

#dump                2 all image 100 image.*.jpg type type &
#                axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0
#dump_modify        2 pad 4

#dump                3 all movie 100 movie.mpg type type &
#                axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0
#dump_modify        3 pad 4

thermo                100
thermo_modify        temp new2d

reset_timestep        0
run                3000
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luyaobao

木虫 (著名写手)

应变我不知道,是不是可以通过上面施加的驱动和盒子尺寸计算,应力如果不需要界面附近的数据可以用compute stress/atom计算后统计平均输出

发自小木虫Android客户端
2楼2018-12-28 12:44:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ruibingqilin

铁虫 (小有名气)

你好,我刚接触lammps,按照例子在熟悉操作。从NIST官网上复制了势函数到text文件里,命名为Al99.eam.alloy,放在计算的文件夹里了,但运行后还是出现Cannot open EAM potential file Al99.eam.alloy 这个错误。请问应该怎么解决? 谢谢!
3楼2019-02-12 11:53:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

束卷卷

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by ruibingqilin at 2019-02-12 04:53:33
你好,我刚接触lammps,按照例子在熟悉操作。从NIST官网上复制了势函数到text文件里,命名为Al99.eam.alloy,放在计算的文件夹里了,但运行后还是出现Cannot open EAM potential file Al99.eam.alloy 这个错误。请问 ...

把势函数文件命名为以.eam结尾的,后面的.alloy删掉,跟lammps自带的shear的例子一样,我是这样的,可以读取势函数
4楼2019-03-05 11:26:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shanks2017

铁虫 (小有名气)

楼主你的问题解决了吗?能不能分享一下方法?谢啦
5楼2019-04-12 09:53:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

嘿咻小尹子

新虫 (小有名气)

楼主请问剪应力输出解决了吗?可否告知是怎么做的?谢谢啦~

发自小木虫Android客户端
6楼2019-09-23 21:43:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 束卷卷 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0831生医工第一轮调剂失败求助 +13 小熊睿睿_s 2026-04-11 17/850 2026-04-16 15:37 by Poly12123
[考研] 271求调剂 +35 2261744733 2026-04-11 37/1850 2026-04-16 13:16 by jiahl2024
[考研] 307中医考研调剂 +5 于以采蘩 2026-04-14 5/250 2026-04-16 10:40 by dick_runner
[考研] 求调剂 +22 MAX怅惘 2026-04-09 24/1200 2026-04-16 10:34 by Espannnnnol
[考研] 求助调剂,跨调 +19 X十甫寸Y 2026-04-11 20/1000 2026-04-15 21:18 by cuisz
[考研] 药学305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-15 13:21 by 西北望—风沙
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[考研] 一志愿沪9,326求生物学调剂 +10 刘墨墨 2026-04-13 10/500 2026-04-14 15:16 by zs92450
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
[考研] 一志愿哈工大 085600 277 12材科基求调剂 5+5 chenny174 2026-04-10 37/1850 2026-04-14 07:39 by Abskk
[考研] 求调剂 +12 璃茉一定上岸 2026-04-10 13/650 2026-04-14 00:08 by Equinoxhua
[考研] 332求调剂 +15 蕉蕉123 2026-04-10 15/750 2026-04-13 23:12 by pies112
[考研] 一志愿双非085400电子信息344 求调剂,对材料和化学方向也感兴趣 +12 无情的小羊 2026-04-09 13/650 2026-04-13 14:17 by 张zhihao
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
[考研] 本人女孩 +7 吼吼, 2026-04-10 9/450 2026-04-11 14:45 by ACS Nano——
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 一志愿京区985,085401,与本科专业一致,电子信息工程, +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
信息提示
请填处理意见