24小时热门版块排行榜    

查看: 2463  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yangwenbin99

铜虫 (小有名气)

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

如图,想在上层加上载荷,结果成这种样子,我用的是fix upper aveforxe NULL NULL -10.0,结果感觉是加速度过大,把-10.0改为-0.1时,经过一段时间是可以达到平衡的,但是载荷大了就不行。希望各位朋友给我出个主意,让我能加上载荷
加载荷过程中产生加速度,模型冲烂,怎么才能正常的加上较大载荷
W17%UH~5`FP8M9U7RPZ2KVM.jpg
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见