24小时热门版块排行榜    

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

Ergon

金虫 (著名写手)

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

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

» 猜你喜欢

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

认真生活
已阅   回复此楼   关注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的回帖
查看全部 10 个回答

Ergon

金虫 (著名写手)

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

fenqin

铜虫 (小有名气)

楼主的问题解决了吗
3楼2015-08-03 15:11:27
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂推荐 +8 小聂爱学习 2026-04-14 8/400 2026-04-16 07:22 by 学员JpLReM
[考研] 273求调剂 +6 白居不易. 2026-04-09 8/400 2026-04-15 22:02 by wooluyong
[考研] 291分调剂 +11 上岸小莹加油 2026-04-09 12/600 2026-04-15 21:55 by noqvsozv
[考研] 297,工科调剂? +10 河南农业大学-能 2026-04-14 10/500 2026-04-15 21:50 by noqvsozv
[考研] 药学求调剂 +11 RussHu 2026-04-12 13/650 2026-04-15 19:07 by zhuwenxu
[考研] 药学305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-15 13:21 by 西北望—风沙
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[基金申请] 2026 WR青拔 +3 冬日阳光CAS 2026-04-09 6/300 2026-04-13 18:40 by liuchb715
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 一志愿华中农微生物,288分,三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 材料与化工300求调剂 +39 肖开文 2026-04-09 43/2150 2026-04-12 01:30 by 秋豆菜芽
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 调剂 +6 青灯不负 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 08:57 by zhq0425
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
信息提示
请填处理意见