24小时热门版块排行榜    

Znn3bq.jpeg
查看: 6425  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854求调剂 +12 门路摸摸 2026-04-15 13/650 2026-04-16 00:14 by Art1977
[考研] 急需调剂 +5 绝不放弃22 2026-04-15 5/250 2026-04-15 23:11 by Tide man
[考研] 求调剂 +11 小聂爱学习 2026-04-11 15/750 2026-04-15 21:57 by noqvsozv
[考研] 291求调剂 +7 关忆北. 2026-04-14 7/350 2026-04-15 14:39 by 黑科技花岗岩
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 6/300 2026-04-14 16:40 by zhouxiaoyu
[考研] 一志愿沪9,326求生物学调剂 +10 刘墨墨 2026-04-13 10/500 2026-04-14 15:16 by zs92450
[考研] 本科211,报考085601-310分 +16 ararak 2026-04-13 16/800 2026-04-14 14:55 by Delta2012
[考研] 人工智能320调剂08工类还有机会吗 +18 振—TZ 2026-04-10 19/950 2026-04-14 10:34 by screening
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] 290求调剂 +18 柯淮然 2026-04-12 20/1000 2026-04-13 12:56 by cyh—315
[考研] 0831生医工第一轮调剂失败求助 +12 小熊睿睿_s 2026-04-11 16/800 2026-04-12 16:28 by 钰璞
[考研] 326求调剂 +6 Shansyn 2026-04-10 6/300 2026-04-12 09:46 by hammer3
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 346,工科求调剂 +3 moser233 2026-04-09 3/150 2026-04-11 10:04 by zhq0425
[考研] 一志愿北理工298英一数二已上岸,感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
信息提示
请填处理意见