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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 加载荷过程中产生加速度，模型冲烂，怎么才能正常的加上较大载荷
141/2返回列表12下一页
查看: 2009  |  回复: 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

铜虫 (小有名气)
自己先顶一下
回复此楼
2楼2013-08-29 09:07:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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

yangwenbin99

铜虫 (小有名气)
引用回帖:
5楼: Originally posted by ddll_slsl at 2013-09-11 09:30:33
应该是你的作用力太大造成的。 把你的作用力转换成应力看看？ 查一下log文件，其能量，力是否出现过大的情况。
需要先确定你要施加多大的应力。

怎么转换成应力，我大体上算了一下，这种情况是50Mpa，不算大啊
赞一下 回复此楼
6楼2013-09-11 10:07:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

你的单位， 力是eV/A, 就是1.6nN。 看图片，横截面积大约是1nm边长的正方形。1个单位应力就是1.6GPa。你设10的话就是16GPa了。 查过你的log文件吗？ 其能量，力变化如何？
赞一下 回复此楼
7楼2013-09-11 16:46:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangwenbin99

铜虫 (小有名气)
引用回帖:
7楼: Originally posted by ddll_slsl at 2013-09-11 16:46:49
你的单位， 力是eV/A, 就是1.6nN。 看图片，横截面积大约是1nm边长的正方形。1个单位应力就是1.6GPa。你设10的话就是16GPa了。 查过你的log文件吗？ 其能量，力变化如何？

Step Temp TotEng str stress vol_upb s_yz friction
  800000     333.0119   -15435.106 2.8873201e+008    75718.902    3813.2091    134.83582     10209620
  800100    332.92624   -15433.512 2.8808579e+008    75549.433    3813.2091    134.83582     10186769
  800200    332.92629   -15431.712 2.8683741e+008     75222.05    3813.2091    134.83582     10142626
  800300    332.92633    -15429.31 2.851059e+008     74767.97    3813.2091    134.83582     10081400
  800400    332.93791   -15428.392 2.83436e+008    74330.044    3813.2091    134.83582     10022352
  800500    332.92641   -15429.525 2.8228536e+008    74028.294    3813.2091    134.83582    9981665.4
  800600    332.92645   -15431.839 2.8175973e+008    73890.448    3813.2091    134.83582    9963078.8
  800700    332.92649   -15434.094 2.817182e+008    73879.558    3813.2091    134.83582    9961610.5
  800800    332.92573   -15434.795 2.8215163e+008    73993.224    3813.2091    134.83582    9976936.7
  800900    332.92657   -15433.883 2.8304636e+008    74227.863    3813.2091    134.83582     10008574
  801000    332.92661   -15432.318 2.8395547e+008    74466.274    3813.2091    134.83582     10040721
赞一下 回复此楼
8楼2013-09-12 09:06:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangwenbin99

铜虫 (小有名气)
引用回帖:
7楼: Originally posted by ddll_slsl at 2013-09-11 16:46:49
你的单位， 力是eV/A, 就是1.6nN。 看图片，横截面积大约是1nm边长的正方形。1个单位应力就是1.6GPa。你设10的话就是16GPa了。 查过你的log文件吗？ 其能量，力变化如何？

我看着没问题啊
回复此楼
9楼2013-09-12 09:06:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

有原子轨迹图吗？ 你的顶层绑定了。 飞出去的原子应该是来自下层， 单个速度大， 但没有到太夸张的程度。 建议用多种压力，作一个势态曲线。
赞一下 回复此楼
10楼2013-09-12 09:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangwenbin99 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭