| 查看: 4098 | 回复: 10 | |||
[求助]
fix deform 指令拉伸求教 已有1人参与
|
|||
|
面是写的IN文件，但是模拟得出的结果，一直到应变为0.5时才将近拉断，而且感觉其中应力一度还出现了特别大的值，所以有哪位用过fix deform指令拉伸的还望指教一下？ #lammps graphene test units metal boundary p s s atom_style atomic neighbor 2.0 nsq neigh_modify delay 1 check yes read_data data.graphene pair_style tersoff pair_coeff * * SiCGe.tersoff C region 1 block INF 3 INF INF INF INF units box group left region 1 region 2 block 61 INF INF INF INF INF units box group right region 2 group boundary union left right group mobile subtract all boundary velocity mobile create 0.01 4928459 dist gaussian units box velocity boundary set 0.0 0.0 0.0 units box #fix fix 1 all nvt temp 0.01 0.01 0.001 fix 2 boundary setforce 0.0 0.0 0.0 min_style sd minimize 1.0e-12 1.0e-12 1000 1000 #run timestep 0.001 compute 2 all pe/atom compute 3 all reduce sum c_2 variable avepo equal c_3/atoms compute 4 all stress/atom compute strec all reduce sum c_4[1] variable s1 equal c_strec/(vol*10000) #X方向应力 thermo_style custom step temp etotal vol v_avepo v_s1 lx thermo 1000 run 100000 dump 1 all custom 1000 dump.out id type x y z #for loop label forloop variable i loop 400 fix 5 all deform 1000 x erate 0.0015 units box run 1000 variable strainx equal lx/62.480000003-1 #X应变 print "@ @ @ @ : ${strainx} ${s1} " next i jump in.graphene forloop  拉伸后的VMD图像 |
回复此楼» 收录本帖的淘帖专辑推荐
 力学性能 |
» 猜你喜欢
职称评审没过，求安慰
已经有51人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教：赴法签证及语言问题
已经有11人回复
- 公派签证体检求教~~教本人要去西班牙，对此情况知之甚少，愿大神们出手援助
已经有14人回复
- 公派成功，但是接下来的问题也很迷茫，求教大家！
已经有11人回复
- 请教几点PBT机械拉伸疑问
已经有3人回复
- 求教：水电厂某一钢管的腐蚀形貌是怎么形成的？？
已经有7人回复
- 求助各位大虾们，实验堵住了！求教
已经有12人回复
- 求教循环伏安法中的补偿问题
已经有6人回复
- 求教做负极材料电池方面高手
已经有8人回复
- 关于拉伸模量的问题
已经有50人回复
- 蚜虫 GC-EAD交流（求教）
已经有9人回复
- 关于fix nvt里的Tdamp
已经有7人回复
- 使用fix ave/spatial指令分层后，层有编号么
已经有8人回复
- 【求助】强制性精修时用AFIX指令，后面的数字什么意思
已经有9人回复
5楼2013-08-07 11:48:16
8楼2014-03-10 09:32:12
★ ★
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-06-01 10:25:37
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-06-01 10:25:37
|
本帖内容被屏蔽 |
2楼2012-05-31 00:49:08
3楼2012-06-01 09:56:48
yinhanqing
金虫 (著名写手)
- 应助: 14 (小学生)
- 金币: 12742.8
- 散金: 70
- 红花: 4
- 帖子: 1461
- 在线: 323.2小时
- 虫号: 1756621
- 注册: 2012-04-15
- 专业: 固体力学
4楼2012-06-17 20:32:15
6楼2013-08-07 11:49:21
ldx626
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1700.1
- 帖子: 19
- 在线: 26.8小时
- 虫号: 2323076
- 注册: 2013-03-06
- 性别: GG
- 专业: 极端条件下使用的金属材料
7楼2014-03-09 22:35:29
9楼2014-03-10 09:57:33
ldx626
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1700.1
- 帖子: 19
- 在线: 26.8小时
- 虫号: 2323076
- 注册: 2013-03-06
- 性别: GG
- 专业: 极端条件下使用的金属材料
10楼2014-03-10 10:22:39