24小时热门版块排行榜    

查看: 5251  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 束卷卷 的主题更新
信息提示
请填处理意见