Znn3bq.jpeg
²é¿´: 2650  |  »Ø¸´: 12
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

kite175

Òø³æ (³õÈëÎÄ̳)


[½»Á÷] ÓÐÈË×öÁÑÎÆÄ£ÄâµÄÂ𣿿´¿´ÎÒÕâ¸öÁÑÎÆµÄinÎļþÓÐʲôÎÊÌâ

# 3d Morse crack simulation


units          metal
boundary        s s p
atom_style        atomic


# create geometry

lattice                fcc 4.05  orient x 1 0 0  orient y 0 1 0 orient z 0 0 1
region                box block 0 60 0 40 -0.25 0.25
create_box        5 box
create_atoms        1 box

mass                1 1.0
mass                2 1.0
mass                3 1.0
mass                4 1.0
mass                5 1.0


# Morse potentials

pair_style        morse 2.5
pair_coeff        * * 1.0 1.0 2.5
neighbor        0.3 bin
neigh_modify        every 5 delay 0 check yes

# define groups

region                1 block INF INF INF 2 INF INF
group                lower region 1
region                2 block INF INF 38 INF INF INF
group                upper region 2
group                boundary union lower upper
group                mobile subtract all boundary


region                leftupper block INF 20 20 INF INF INF
region                leftlower block INF 20 INF 20 INF INF

region          bonds   block  INF 20  19.75 20.25 -0.25 0.25
delete_atoms    region  bonds

group                leftupper region leftupper
group                leftlower region leftlower


group                halfupper1 region halfupper1
group                halfupper2 region halfupper2

set                group leftupper type 3
set                group leftlower type 4


set                group lower type 1
set                group upper type 2

neigh_modify        exclude type 3 4
neigh_modify        exclude type 1 2


# ÄÜÁ¿×îС»¯
min_style sd
minimize       1.0e-12 1.0e-12  10000 10000

# initial velocities ³ÚÔ¥

compute                  new mobile temp
velocity        mobile create 0.001 887723 temp new
timestep        0.001
thermo                200
run 10000

velocity        upper  set 0.0 0.01 0.0
velocity        lower  set 0.0 -0.01 0.0
velocity        mobile ramp vy 0.01 0.01 y 20 38 sum yes

# fixes
fix            1 all nve  
fix            2 lower setforce  0.0 NULL 0.0
fix            3 upper setforce  0.0 NULL 0.0
# run
thermo                200
thermo_modify        temp new

dump                1 all atom 1000 dump.crack.lammpstrj
run                300000
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

» ÇÀ½ð±ÒÀ²£¡»ØÌû¾Í¿ÉÒԵõ½:

²é¿´È«²¿É¢½ðÌù

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

tfl03

ľ³æ (Ö°Òµ×÷¼Ò)


¡ï ¡ï ¡ï ¡ï
Сľ³æ: ½ð±Ò+0.5, ¸ø¸öºì°ü£¬Ð»Ð»»ØÌû
zh1987hs: ½ð±Ò+3, лл 2012-10-01 11:02:11
# 3d Morse crack simulation


units          metal
boundary        s s p
atom_style        atomic


# create geometry

lattice                fcc 4.05  orient x 1 0 0  orient y 0 1 0 orient z 0 0 1
region                box block 0 60 0 40 -0.25 0.25
create_box        5 box
create_atoms        1 box

mass                1 1.0
mass                2 1.0
mass                3 1.0
mass                4 1.0
mass                5 1.0


# Morse potentials

pair_style        morse 2.5
pair_coeff        * * 1.0 1.0 2.5
neighbor        0.3 bin
neigh_modify        every 5 delay 0 check yes

# define groups

region                1 block INF INF INF 2 INF INF
group                lower region 1
region                2 block INF INF 38 INF INF INF
group                upper region 2
group                boundary union lower upper
group                mobile subtract all boundary


region                leftupper block INF 20 20 INF INF INF
region                leftlower block INF 20 INF 20 INF INF

region          bonds   block  INF 20  19.75 20.25 -0.25 0.25
delete_atoms    region  bonds

group                leftupper region leftupper
group                leftlower region leftlower


group                halfupper1 region halfupper1
group                halfupper2 region halfupper2

set                group leftupper type 3
set                group leftlower type 4


set                group lower type 1
set                group upper type 2

neigh_modify        exclude type 3 4
neigh_modify        exclude type 1 2


# ÄÜÁ¿×îС»¯
min_style sd
minimize       1.0e-12 1.0e-12  10000 10000

# initial velocities ³ÚÔ¥

compute                  new mobile temp
velocity        mobile create 0.001 887723 temp new
timestep        0.001
thermo                200
run 10000

velocity        upper  set 0.0 0.01 0.0
velocity        lower  set 0.0 -0.01 0.0
velocity        mobile ramp vy 0.01 0.01 y 20 38 sum yes

# fixes
fix            1 all nve  
fix            2 lower setforce  0.0 NULL 0.0
fix            3 upper setforce  0.0 NULL 0.0
# run
thermo                200
thermo_modify        temp new

dump                1 all atom 1000 dump.crack.lammpstrj
run                300000
6Â¥2012-09-29 11:51:41
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 13 ¸ö»Ø´ð

kite175

Òø³æ (³õÈëÎÄ̳)


ÔõôûÓÐÈ˻شðѽ£¿¾À½áѽ
2Â¥2012-09-05 10:11:52
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

1779064387

гæ (³õÈëÎÄ̳)


¡ï ¡ï ¡ï
Сľ³æ: ½ð±Ò+0.5, ¸ø¸öºì°ü£¬Ð»Ð»»ØÌû
jiaoyixiong: ½ð±Ò+2, ¹ÄÀø½»Á÷ 2012-09-27 19:49:16
½¨Òé°ÑÄãµÄdataÎļþÌùÒ»ÏÂ~
3Â¥2012-09-27 15:19:20
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
¡ï ¡ï ¡ï
Сľ³æ: ½ð±Ò+0.5, ¸ø¸öºì°ü£¬Ð»Ð»»ØÌû
chaizhm: ½ð±Ò+2, лл~ 2012-09-28 09:38:30
ÒýÓûØÌû:
2Â¥: Originally posted by kite175 at 2012-09-05 10:11:52
ÔõôûÓÐÈ˻شðѽ£¿¾À½áѽ

²»ÖªÎªºÎ¾À½á£¿

ÏëÖªµÀÄãµÄinÎļþÓкÎÎÊÌ⣬×î¼òµ¥µÄ°ì·¨¾ÍÊÇÄãÔËÐÐһϰ¡£¬
ÓÐʲôÎÊÌ⣬Á¢¿Ì¾ÍÖªµÀÁË¡£
Èç¹û˳ÀûÔËÐУ¬ÄǾÍûʲôÎÊÌ⣬¶àºÃ°¡¡£
4Â¥2012-09-27 19:51:35
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] Çóµ÷¼Á +4 ¿àÃüÈË¡£¡£¡£ 2026-04-18 4/200 2026-04-19 02:01 by ÑÌÓêÁ÷ÑÄ
[¿¼ÑÐ] 300Çóµ÷¼Á +12 ³Èa777 2026-04-15 12/600 2026-04-18 23:51 by ·²¡Çé
[¿¼ÑÐ] ½ÓÊÜÈκε÷¼Á +6 Ò²¾ÍÊÇÀõ×Ó 2026-04-17 7/350 2026-04-18 17:20 by º­ÖñÁõ
[¿¼ÑÐ] 22408 312Çóµ÷¼Á +24 ÃÅ·ÃþÃþ 2026-04-14 26/1300 2026-04-18 13:04 by wunaiy88
[¿¼ÑÐ] ¼±Ðèµ÷¼Á +9 ¾ø²»·ÅÆú22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[¿¼ÑÐ] »¯¹¤Ñ§Ë¶294·Ö£¬Çóµ¼Ê¦ÊÕÁô +33 yzyzx 2026-04-12 37/1850 2026-04-17 23:00 by wunaiy88
[¿¼ÑÐ] Ò»Ö¾Ô¸»ªÖÐũҵ071010£¬320Çóµ÷¼Á +17 À§À§À§À§À¤À¤ 2026-04-14 19/950 2026-04-17 20:08 by ¹ØÒ»ÕµµÆcd
[Óлú½»Á÷] ¶þ±½¼×ͪËáÀàÑÜÉúÎï 50+3 С°×°®Ö÷ÈË 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[¿¼ÑÐ] 322Çóµ÷¼Á +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[¿¼ÑÐ] 291Çóµ÷¼Á +11 ¹ØÒä±±. 2026-04-14 11/550 2026-04-16 15:18 by jiahl2024
[¿¼ÑÐ] ҩѧÇóµ÷¼Á +14 à¶¹þ¼ÓÓÍ 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[¿¼ÑÐ] 289 ·Ö105500ҩѧר˶Çóµ÷¼Á(ÕÒBÇøÑ§Ð£) +4 °×ÔÆ123456789 2026-04-13 4/200 2026-04-16 00:18 by ·ÛÇßÈô³¾
[¿¼ÑÐ] Çóµ÷¼ÁѧУ +14 ²»»á³ÔÈâ 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[¿¼ÑÐ] ͨÐŹ¤³ÌÇóµ÷¼Á£¡£¡£¡ +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by ѧԱJpLReM
[¿¼ÑÐ] 105500ҩѧÇóµ÷¼Á +4 x_skys 2026-04-12 4/200 2026-04-14 13:37 by rndfc
[¿¼ÑÐ] ¿¼ÑÐÓ¢Ò»ÊýÒ»338·Ö +9 ³¤½­´óѧ¶«Ð£Çø 2026-04-13 10/500 2026-04-14 00:41 by Íõ¬Bè±
[¿¼ÑÐ] 2026˶ʿµ÷¼Á_Äܶ¯_ºÓÄÏũҵ´óѧ +4 ºÓÄÏũҵ´óѧ-ÄÜ 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖÐÄÏ´óѧ 0855 »úе 286 Çóµ÷¼Á +11 ²»»á³ÔÈâ 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[¿¼ÑÐ] 302Çóµ÷¼Á +10 Ò×£¡? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[¿¼ÑÐ] 297¹¤¿Æ£¬Çóµ÷¼Á? +13 ºÓÄÏũҵ´óѧ-ÄÜ 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û