| 1 | 1/1 | 返回列表 |
| 查看: 5039 | 回復(fù): 0 | ||
honghaier90新蟲 (小有名氣)
|
[求助]
lammps里,關(guān)于應(yīng)力應(yīng)變的問題
|
|
我在用lammps做材料的單軸拉伸和壓時(shí),遇到以下幾個(gè)問題: 1,第0步輸出的應(yīng)力不為零,第0步材料未發(fā)生變形,為什么會(huì)出現(xiàn)應(yīng)力,按理說應(yīng)力也應(yīng)該為零啊。拉壓之前,均做過能量最小化,但是還是會(huì)出現(xiàn)這種情況,這個(gè)該如何處理? 2,對(duì)于同種材料分別做拉伸壓縮,在得出應(yīng)力應(yīng)變曲線后,通過直線擬合,發(fā)現(xiàn)兩條直線斜率不相等,也就是彈性模量不相等,這是為什么? 以下是我做拉伸時(shí)的in文件,懇請(qǐng)大家給于指點(diǎn)。 dimension 2 boundary s s p atom_style angle neighbor 0.3 bin neigh_modify every 10 delay 0 check yes lattice hex 0.91648 read_data data.new1 pair_style lj/cut 1.5 pair_coeff 1 1 1.0 1.0 1.5 bond_style harmonic bond_coeff 1 100 1.12246 angle_style harmonic angle_coeff 1 1 60 region left block INF 10 INF INF INF INF units lattice region right block 190 INF INF INF INF INF units lattice group left region left group right region right group boundary union left right group mobile subtract all boundary thermo_style custom lx #min_style sd #minimize 1.0e-12 1.0e-12 10000 10000 variable tmp equal lx variable L equal ${tmp} print "initial length , L:${L}" #deformation reset_timestep 0 velocity left set 0.0 NULL 0.0 compute new mobile temp variable pressx equal c_p[1] compute s all stress/atom compute strex all reduce sum c_s[1] variable x equal c_strex/vol variable strain1 equal "(lx - 224.6726)/224.6726" #variable strain2 equal "(-ly + 69.107548)/69.107548" thermo_style custom step temp vol lx ly lz pxx v_x c_strex v_strain1 thermo_modify temp new velocity mobile create 0.01 887723 temp new velocity right set 0.3 0.0 0.0 velocity mobile ramp vx 0.0 0.3 x 10 190 sum yes fix 1 all nve fix 2 boundary setforce 0.0 NULL 0.0 fix 3 mobile temp/rescale 100 0 0.05 0.002 0.1 timestep 0.001 thermo 100 dump 1 all atom 100 dump.pull1.lammpstrj run 80000 |

找到一些相關(guān)的精華帖子,希望有用哦~
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +4 |  ̄^ ̄゜汗 2026-03-19 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +10 | 申子申申 2026-03-19 | 16/800 |
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +6 | 皮皮皮123456 2026-03-22 | 6/300 |
|
|
[考研] 298求調(diào)劑 +5 | 上岸6666@ 2026-03-20 | 5/250 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 318求調(diào)劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
[考研] 083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑 +3 | innocenceF 2026-03-17 | 3/150 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計(jì)算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|