应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 加载荷过程中产生加速度，模型冲烂，怎么才能正常的加上较大载荷
52/1返回列表
查看: 2017  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yangwenbin99

铜虫 (小有名气)

[求助] 加载荷过程中产生加速度，模型冲烂，怎么才能正常的加上较大载荷

如图，想在上层加上载荷，结果成这种样子，我用的是fix upper aveforxe NULL NULL -10.0，结果感觉是加速度过大，把-10.0改为-0.1时，经过一段时间是可以达到平衡的，但是载荷大了就不行。希望各位朋友给我出个主意，让我能加上载荷
加载荷过程中产生加速度，模型冲烂，怎么才能正常的加上较大载荷
W17%UH~5`FP8M9U7RPZ2KVM.jpg
回复此楼

» 猜你喜欢
1楼 2013-08-28 22:10:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangwenbin99

铜虫 (小有名气)
引用回帖:
10楼: Originally posted by ddll_slsl at 2013-09-12 09:24:04
有原子轨迹图吗？ 你的顶层绑定了。 飞出去的原子应该是来自下层， 单个速度大， 但没有到太夸张的程度。 建议用多种压力，作一个势态曲线。

我不太懂，用多种压力是什么意思
赞一下 回复此楼
11楼2013-09-12 10:06:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

一般应该对upper的自由度做个限制， 贴出你的in文件看看。
赞一下 回复此楼
3楼2013-09-10 15:11:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangwenbin99

铜虫 (小有名气)
引用回帖:
3楼: Originally posted by ddll_slsl at 2013-09-10 15:11:16
一般应该对upper的自由度做个限制， 贴出你的in文件看看。

# 3-d LJ flow simulation

dimension        3
boundary        p p p

units  metal
lattice bcc 2.8664
atom_style        full
neighbor        2.0 bin
neigh_modify         every 1 delay 5 check yes



bond_style harmonic
angle_style harmonic
read_data layer(13).lammps05
#read_restart restart.equil
mass 1 15.9994
mass 2 1.008
mass 3 55.847

bond_coeff 1 25.6971 0.9572
angle_coeff 1 2.385 104.5


pair_style hybrid eam/fs lj/cut 8.6
pair_coeff 1 1 lj/cut 0.006596 3.1507 8.6
pair_coeff 2 2 lj/cut 0.0 0.0 8.6
pair_coeff 3 3 lj/cut 0.2589 4.3 8.6
pair_coeff 1 2 lj/cut 0.0 1.57535 8.6
pair_coeff 2 3 lj/cut 0.0 2.15 8.6
pair_coeff 1 3 lj/cut 0.04132 3.72535 8.6
pair_coeff * * eam/fs Fe_mm.eam.fs NULL NULL Fe


region   lob  block INF INF INF INF INF 5.7328 side in units box
group    lob  region lob
region   lohwc   block INF INF INF INF 5.7328 14.332 side in units box
group    lohwc   region lohwc
region   londc    block INF INF INF INF 14.332 33.664 side in units box
group    londc    region londc
region   water   block INF INF INF INF 33.664 49.9583 side in units box
group    water   region water
region   upb  block INF INF INF INF 72.8895 84.1223 side in units box
group    upb  region upb
region   uphwc block INF INF INF INF 64.2903 72.8895 side in units box
group    uphwc region uphwc
region   upndc block INF INF INF INF 49.9583 64.2903 side in units box
group    upndc region upndc
group    hwc    union lohwc uphwc
group    ndc   union upndc londc
group    mobile   union londc upndc water lohwc uphwc
group    boundary union lob upb
group    upper union upb uphwc upndc
group    upmov union upb uphwc
timestep 0.0001

minimize 1.0e-4 1.0e-6 100 1000
reset_timestep 0

#initial velocity

timestep 0.0001
velocity all create 333 482748  
fix 1 upb rigid group 1 upb  

fix 2 mobile nve
compute t_h mobile temp
fix 3 mobile temp/rescale 1 333 333.1 0.02 1.0
fix_modify 3 temp t_h

fix 5 upb aveforce NULL  NULL  -0.25675
fix 6 lob setforce 0.0 0.0 0.0

compute temp_mobile mobile temp

thermo 100
thermo_style custom step temp etotal c_temp_mobile
dump 3 all atom 100 output.pfx.lammpstrj
run  700000

unfix 1
fix 8 upb move linear  0.5 0.0 0.0
compute t_b boundary temp
fix 9 boundary temp/rescale 1 333 333.1 0.02 1.0
fix_modify 9 temp t_b
run 800000 upto

compute s upb stress/atom
compute str upb reduce sum c_s[1]

compute minx_upb upb reduce min x
compute maxx_upb upb reduce max x

compute miny_upb upb reduce min y
compute maxy_upb upb reduce max y

compute minz_upb upb reduce min z
compute maxz_upb upb reduce max z

variable lx_upb  equal c_maxx_upb-c_minx_upb
variable ly_upb  equal c_maxy_upb-c_miny_upb
variable lz_upb  equal c_maxz_upb-c_minz_upb

variable vol_upb equal v_lx_upb*v_ly_upb*v_lz_upb
variable s_yz    equal v_ly_upb*v_lz_upb

variable stress  equal c_str/v_vol_upb

variable friction equal v_stress*v_s_yz

fix 10 upb ave/time 1 800000 1000000   v_friction file upb.temp

thermo 100
thermo_style custom step temp etotal  c_str v_stress v_vol_upb v_s_yz v_friction
dump 4 all atom 100 output1.pfx.lammpstrj

run 1000000 upto
write_restart restart.equil
赞一下 回复此楼
4楼2013-09-10 15:43:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

应该是你的作用力太大造成的。 把你的作用力转换成应力看看？ 查一下log文件，其能量，力是否出现过大的情况。
需要先确定你要施加多大的应力。
赞一下 回复此楼
5楼2013-09-11 09:30:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭