24小时热门版块排行榜    

查看: 5474  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 束卷卷 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于豆爷回答的JTJC与%2F数量 +3 yang182083 2026-08-06 4/200 2026-08-07 11:07 by 虫友是什么虫
[基金申请] 基金中了 +11 laoda193707 2026-08-06 11/550 2026-08-07 10:54 by TT和NN
[基金申请] filecode +14 等待解的谜 2026-08-06 18/900 2026-08-07 10:13 by naalan7001
[基金申请] 听说今天filecode变了 +21 布布和一二 2026-08-06 44/2200 2026-08-07 10:11 by 布布和一二
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 6/300 2026-08-06 20:19 by 苏知砚
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 12/600 2026-08-06 15:11 by eddyin
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[教师之家] 咨询面上基金 +4 李长云 2026-07-31 7/350 2026-08-06 11:20 by 李长云
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +22 低垂的野花 2026-07-31 30/1500 2026-08-05 18:03 by 低垂的野花
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了,论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见