24小时热门版块排行榜    

查看: 2382  |  回复: 13

yangwenbin99

铜虫 (小有名气)

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

如图,想在上层加上载荷,结果成这种样子,我用的是fix upper aveforxe NULL NULL -10.0,结果感觉是加速度过大,把-10.0改为-0.1时,经过一段时间是可以达到平衡的,但是载荷大了就不行。希望各位朋友给我出个主意,让我能加上载荷
加载荷过程中产生加速度,模型冲烂,怎么才能正常的加上较大载荷
W17%UH~5`FP8M9U7RPZ2KVM.jpg
回复此楼
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jRl3mE6ddZGq 2026-08-03 4/200 2026-08-04 04:54 by k7OM8YghWbkC
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 04:50 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:41 by k7OM8YghWbkC
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 9/450 2026-08-04 03:56 by kf2781974
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见