24小时热门版块排行榜    

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

Ergon

金虫 (著名写手)

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

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

» 猜你喜欢

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

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

3477180472

铁虫 (初入文坛)

楼主,剪切应力应该怎么输出

发自小木虫IOS客户端
9楼2020-05-25 00:22:42
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297,工科调剂? +10 河南农业大学-能 2026-04-14 10/500 2026-04-15 21:50 by noqvsozv
[考研] 322求调剂 +7 123安康 2026-04-12 14/700 2026-04-15 21:48 by noqvsozv
[考研] 085404 22408 309分求调剂 +9 lzmk 2026-04-14 10/500 2026-04-15 20:02 by 学员JpLReM
[考研] 307中医考研调剂 +4 于以采蘩 2026-04-14 4/200 2026-04-15 19:14 by AN流800
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 农学0904 312求调剂 +4 Say Never 2026-04-11 4/200 2026-04-14 09:10 by zs92450
[考研] 考研英一数一338分 +9 长江大学东校区 2026-04-13 10/500 2026-04-14 00:41 by 王珺璞
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 314求调剂 +18 xhhdjdjsjks 2026-04-09 19/950 2026-04-10 18:53 by HPUCZ
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
信息提示
请填处理意见