24小时热门版块排行榜    

CyRhmU.jpeg
查看: 6218  |  回复: 9

Ergon

金虫 (著名写手)

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

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

» 猜你喜欢

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

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

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的回帖

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的回帖

Ergon

金虫 (著名写手)

引用回帖:
5楼: Originally posted by fenqin at 2015-08-05 08:55:05
谢谢回复,我想问下,在剪切拉伸之前,有做过优化弛豫吗,我看脚本中没有体现...

我只是给你剪切的脚本
说白了就是
fix                2 all deform 1 xy erate 5e-6 units box remap v
其余的和分析输出有关,数据什么的都没给呢,要自己建模根据模拟需求。
认真生活
6楼2015-08-05 11:34:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuihuina

银虫 (小有名气)

你好,我按你的deformation脚本进行剪切拉伸,但是剪切应变率变化范围是(-0.5,0.5),如此这样的循环,这是怎么回事啊?谢谢解答
天道酬勤
7楼2017-04-26 21:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

青丘白浅

新虫 (初入文坛)

【答案】应助回帖

引用回帖:
7楼: Originally posted by cuihuina at 2017-04-26 21:08:56
你好,我按你的deformation脚本进行剪切拉伸,但是剪切应变率变化范围是(-0.5,0.5),如此这样的循环,这是怎么回事啊?谢谢解答

想咨询一下楼楼解决这个问题了吗,能解答一下吗,谢谢
8楼2018-03-31 14:48:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

3477180472

铁虫 (初入文坛)

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

发自小木虫IOS客户端
9楼2020-05-25 00:22:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kkhkkh

新虫 (初入文坛)

想问一下楼主和大家,如果是正交体系,该怎么实现剪切变形呢?(是剪切变形 不是拉伸变形)

发自小木虫Android客户端
10楼2021-03-19 09:00:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Ergon 的主题更新
信息提示
请填处理意见