| 查看: 4362 | 回复: 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图像 |
回复此楼» 收录本帖的淘帖专辑推荐
 力学性能 |
» 猜你喜欢
291分调剂
已经有9人回复
-
调剂求收留
已经有34人回复
-
291 求调剂
已经有38人回复
-
22408 312求调剂
已经有17人回复
-
一志愿华中农业071010，320求调剂
已经有6人回复
-
290调剂生物0860
已经有41人回复
-
291求调剂
已经有3人回复
-
211本科材料化工求调剂
已经有23人回复
-
山东省基金2026
已经有9人回复
-
药学求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教：赴法签证及语言问题
已经有11人回复
- 公派签证体检求教~~教本人要去西班牙，对此情况知之甚少，愿大神们出手援助
已经有14人回复
- 公派成功，但是接下来的问题也很迷茫，求教大家！
已经有11人回复
- 请教几点PBT机械拉伸疑问
已经有3人回复
- 求教：水电厂某一钢管的腐蚀形貌是怎么形成的？？
已经有7人回复
- 求助各位大虾们，实验堵住了！求教
已经有12人回复
- 求教循环伏安法中的补偿问题
已经有6人回复
- 求教做负极材料电池方面高手
已经有8人回复
- 关于拉伸模量的问题
已经有50人回复
- 蚜虫 GC-EAD交流（求教）
已经有9人回复
- 关于fix nvt里的Tdamp
已经有7人回复
- 使用fix ave/spatial指令分层后，层有编号么
已经有8人回复
- 【求助】强制性精修时用AFIX指令，后面的数字什么意思
已经有9人回复
9楼2014-03-10 09:57:33
★ ★
感谢参与，应助指数 +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
