24小时热门版块排行榜    

查看: 6389  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Ergon

金虫 (著名写手)

[求助] 剪切应变的输出问题—— 已有1人参与

三斜晶体的变形命令用的是
fix        3 all deform 1 x erate 1e-6 remap x units box  
盒子如逾期值变形了,但是剪切应变shear-strain怎么得到?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

fenqin

铜虫 (小有名气)

楼主的问题解决了吗
3楼2015-08-03 15:11:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

Ergon

金虫 (著名写手)

我自己回复一下吧。在手册里有个剪切应变的简单介绍,貌似是x上的变化除以y的长度就是剪切度。
在前边的介绍章节内我发现xy貌似就是x的变化。
因而觉得可以用变量应用。
不知准确与否
认真生活
2楼2015-07-29 11:14:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ergon

金虫 (著名写手)

引用回帖:
3楼: Originally posted by fenqin at 2015-08-03 15:11:27
楼主的问题解决了吗

应该算是解决了。我在二楼已经回复。
我还是黏贴一下脚本吧。主要是变量的应用。
# Uniaxial Tensile Deformation
run 0

variable tmp equal "lx"
variable L0 equal ${tmp}
variable L1 equal "xy"
variable strain equal "v_L1/v_L0"
variable p1 equal "v_strain"
variable p2 equal "-pxx/10000*1.01325"
variable p3 equal "-pyy/10000*1.01325"
variable p4 equal "-pzz/10000*1.01325"
variable p5 equal "lx"
variable p6 equal "ly"
variable p7 equal "lz"
variable p8 equal "temp"
variable t2 equal "epair"
variable t3 equal "ebond"
variable t4 equal "eangle"
variable t5 equal "edihed"
variable t6 equal "-press/10000*1.01325"


fix                1 all nvt/sllod temp 300 300 50 tchain 1
#fix                1 all npt temp 300 300 50 x 0 0 1000 z 0 0 1000 drag 2
fix                2 all deform 1 xy erate 5e-6 units box remap v

compute csym all centro/atom fcc
compute peratom all pe/atom

dump     2 all custom 250 dump.shear.* id type xs ys zs c_csym c_peratom fx fy fz
fix def1 all print 100 "${p1} ${p2} ${p3} ${t6} ${p4} ${p5} ${p6} ${p7} ${p8}" file shear1.txt screen no
fix def2 all print 100 "${p1} ${t2} ${t3} ${t4} ${t5}" file shear2.txt screen no
thermo_style        custom step temp pxx pyy pzz lx ly lz epair ebond eangle edihed
thermo          100
timestep        4
reset_timestep        0
run                17180
unfix 1
unfix 2
unfix def1
unfix def2
认真生活
4楼2015-08-03 16:47:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fenqin

铜虫 (小有名气)

引用回帖:
4楼: Originally posted by Ergon at 2015-08-03 16:47:25
应该算是解决了。我在二楼已经回复。
我还是黏贴一下脚本吧。主要是变量的应用。
# Uniaxial Tensile Deformation
run 0

variable tmp equal "lx"
variable L0 equal ${tmp}
variable L1 equal  ...

谢谢回复,我想问下,在剪切拉伸之前,有做过优化弛豫吗,我看脚本中没有体现
5楼2015-08-05 08:55:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂院校信息 +4 CX 330 2026-03-21 4/200 2026-03-21 23:48 by ms629
[考研] 资源与环境 调剂申请(333分) +5 holy J 2026-03-21 5/250 2026-03-21 22:42 by Catalysis25
[考研] 280求调剂 +11 咕噜晓晓 2026-03-18 12/600 2026-03-21 22:40 by ACS Nano——
[考研] 333求调剂 +5 87639 2026-03-21 7/350 2026-03-21 19:31 by ColorlessPI
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 085601调剂 358分 +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 一志愿武理材料工程348求调剂 +3  ̄^ ̄゜汗 2026-03-19 4/200 2026-03-20 21:01 by zhukairuo
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
信息提示
请填处理意见