 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1851 | 回復(fù): 3 | |||
| 【懸賞金幣】回答本帖問題，作者nulixiaoda將贈(zèng)送您 15 個(gè)金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
nulixiaoda新蟲 (初入文壇)
|
[求助]
lammps拉伸時(shí)變形特別大了還沒有拉斷 已有1人參與
|
||
|
我在用lammps進(jìn)行水合物的拉伸，但拉伸時(shí)三個(gè)尺寸都發(fā)成變形，而且水合物變形特別大了還是不會(huì)斷裂。導(dǎo)出來的應(yīng)力應(yīng)變曲線也是不規(guī)律。請(qǐng)問各位大佬知不知道這是什么原因呢，小弟在這跪求了，謝謝各位！下面是我的in文件。 #水合物拉伸 #初始化 echo screen units real dimension 3 boundary p p p atom_style full #勢(shì)函數(shù) pair_style lj/cut/coul/long 14 kspace_style ewald/disp 1.0e-4 bond_style harmonic angle_style harmonic dihedral_style harmonic pair_modify mix arithmetic #定義文件 read_data "shw.data" #公共近鄰 neighbor 0.3 bin neigh_modify every 2 delay 10 check yes page 100000 velocity all create 273 12345 mom yes rot no #定義變量 計(jì)算值 #應(yīng)力 variable stress equal "-pxx/10000" #應(yīng)變 variable l_x equal lx variable lx0 equal ${l_x} variable strain equal (lx-v_lx0)/v_lx0 #馳豫 timestep 0.001 fix 1 all npt temp 243 243 0.1 z 50 50 0.1 y 50 50 0.1 drag 0.1 thermo 1000 thermo_style custom step lx ly lz press pxx pyy pzz pe temp etotal thermo_modify lost ignore run 20000 #拉伸 reset_timestep 0 unfix 1 fix 1 all npt temp 243 243 0.1 y 50 50 0.1 z 50 50 0.1 drag 0.1 fix 2 all deform 10 x erate 0.001 units box fix 3 all print 100 "${strain} ${stress}" file strain_stress.txt screen no thermo 1000 thermo_style custom step lx ly lz press pxx pyy pzz pe temp etotal thermo_modify lost ignore dump traj all custom 500 shw.xyz id type x y z dump_modify traj sort id timestep 0.001 run 200000 發(fā)自小木蟲Android客戶端 |
回復(fù)此樓» 猜你喜歡
265求調(diào)劑
已經(jīng)有11人回復(fù)
-
材料與化工（0856）304求B區(qū)調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿天津大學(xué)材料與化工275求調(diào)劑
已經(jīng)有19人回復(fù)
-
材料工程307，求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿安徽大學(xué)材料工程專碩313分，求調(diào)劑的學(xué)校
已經(jīng)有7人回復(fù)
-
一志愿山東大學(xué)105500藥學(xué)專碩，總分302求調(diào)劑
已經(jīng)有13人回復(fù)
-
313分生物學(xué)求調(diào)劑
已經(jīng)有7人回復(fù)
-
315食品工程求調(diào)劑
已經(jīng)有7人回復(fù)
-
化工0817調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿211，0860總分286食品工程專業(yè)求調(diào)劑
已經(jīng)有8人回復(fù)
白馬踏平川
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 121.4
- 帖子: 37
- 在線: 3.1小時(shí)
- 蟲號(hào): 29139620
- 注冊(cè): 2022-03-14
- 專業(yè): 原子和分子物理
|
您好我想請(qǐng)教您幾個(gè)問題，以下是我的in文件 #npt弛豫 units real timestep 1 #npt弛豫 fix NPT all npt temp 270 270 100 x 98 98 1000 y 98 98 1000 z 0 0 1000 thermo_style custom time atoms step temp press vol lx ly lz enthalpy pe ke etotal density pxx pyy pzz pxy pxz pyz thermo ${thermo} dump 1 all xyz ${thermo} npt_${name}.xyz dump_modify 1 element O H C fix pt all print 200 "${istep} ${etol} ${han} ${T0} ${V} ${PE} ${P0}" screen no file npt_${name}.dat run ${runmd} undump 1 unfix pt unfix NPT reset_timestep 0 #stress-strain variable tmpz equal "lz" variable Lz0 equal ${tmpz} variable strainz equal "(lz - v_Lz0)/v_Lz0" variable stressz equal "-pzz/10000" #z向拉伸 fix 2 all npt temp 270 270 100 x 98 98 1000 y 98 98 1000 fix 3 all deform 100 z erate 0.00000001 remap x units box thermo_style custom time atoms step temp press vol lx ly lz density pxx pyy pzz pxy pxz pyz thermo ${thermo} dump 1 all xyz ${thermo} tension_${name}.xyz dump_modify 1 element O H C #保存應(yīng)力應(yīng)變數(shù)據(jù) fix def3 all print 100 "${strainx} ${stressx} ${strainy} ${stressy} ${strainz} ${stressz}" screen no file stress-strain_${name}.dat run ${runten} 1.您所說的x方向拉伸，在NPT時(shí)，使x方向的壓強(qiáng)調(diào)節(jié)為零，是在弛豫階段設(shè)定為0嗎？我是z方向拉伸，這么設(shè)置對(duì)嗎？但是我這么設(shè)置應(yīng)力依舊不是從0開始的，我檢查了log文件發(fā)現(xiàn)我的模擬過程控不住壓，壓力波動(dòng)巨大，是否是這個(gè)原因那？應(yīng)該怎么解決？ 2.我的應(yīng)力應(yīng)變文件的屈服強(qiáng)度很大，有8GPa，這明顯是有問題的，這是什么原因那？也是因?yàn)榭刂撇蛔毫?dǎo)致的嗎？ 非常希望能得到您的指點(diǎn)，感謝。 |
3樓2024-01-11 17:06:33
2樓2022-07-04 17:19:20
白馬踏平川
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 121.4
- 帖子: 37
- 在線: 3.1小時(shí)
- 蟲號(hào): 29139620
- 注冊(cè): 2022-03-14
- 專業(yè): 原子和分子物理
4樓2024-01-11 17:13:36
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化工0817調(diào)劑 +7 | 燦若星晨 2026-03-10 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 282分材料專業(yè)求調(diào)劑院校 +15 | 楓橋ZL 2026-03-09 | 21/1050 |
|
|
|
[考研] 一志愿天津大學(xué)，英一數(shù)二305分求調(diào)劑，四六級(jí)已過 +7 | 小小番的茄 2026-03-09 | 7/350 |
|
|
|
[考研] 一志愿天大化工（085600）調(diào)劑總分338 +5 | 蔡大美女 2026-03-09 | 5/250 |
|
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +6 | mingmingry 2026-03-09 | 7/350 |
|
|
|
[考博] 求材料讀博院校 +9 | yanglei131 2026-03-08 | 9/450 |
|
|
|
[碩博家園] 木蟲好像不熱鬧了，是不是？ +4 | 偏振片 2026-03-10 | 4/200 |
|
|
|
[考研] 化工284求調(diào)劑 +12 | 小黑想變白 2026-03-07 | 12/600 |
|
|
|
[考研] 293一志愿華東理工 0817化學(xué)工程與技術(shù) 調(diào)劑 +5 | fjj0912 2026-03-07 | 5/250 |
|
|
|
[考研] 一志愿211 085600 280數(shù)二英二求調(diào)劑 +3 | 月山斜 2026-03-06 | 3/150 |
|
|
|
[考研] 267化工調(diào)劑求助 +7 | 聰少OZ 2026-03-04 | 7/350 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)0856材料與化工，前三科206，總分283，求調(diào)劑 +7 | 26考研求調(diào)劑 2026-03-06 | 7/350 |
|
|
|
[考研] 材料277分求調(diào)劑 +13 | 飯飯星球 2026-03-04 | 14/700 |
|
|
|
[考研] 復(fù)試調(diào)劑 +7 | 呼呼？~+123456 2026-03-05 | 10/500 |
|
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-05 | 5/250 |
|
|
|
[考研] 材料考研339求調(diào)劑 +3 | Karry*^_^* 2026-03-04 | 3/150 |
|
|
|
[考研] 274環(huán)境工程求調(diào)劑 +6 | 扶柳盈江 2026-03-05 | 6/300 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，總分320分，一志愿中南民族大學(xué) +3 | 19286712440 2026-03-04 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +4 | 小賽不吃香菜 2026-03-04 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +6 | 等春來， 2026-03-04 | 6/300 |
|