| 查看: 5427 | 回复: 0 | |||
honghaier90新虫 (小有名气)
|
[求助]
lammps里，关于应力应变的问题
|
|
我在用lammps做材料的单轴拉伸和压时，遇到以下几个问题： 1，第0步输出的应力不为零，第0步材料未发生变形，为什么会出现应力，按理说应力也应该为零啊。拉压之前，均做过能量最小化，但是还是会出现这种情况，这个该如何处理？ 2，对于同种材料分别做拉伸压缩，在得出应力应变曲线后，通过直线拟合，发现两条直线斜率不相等，也就是弹性模量不相等，这是为什么？ 以下是我做拉伸时的in文件，恳请大家给于指点。 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 |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有9人回复
-
厉娜
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
找到一些相关的精华帖子，希望有用哦~
vasp计算应力--应变曲线时应力值发生突变的问题
已经有8人回复
- 内应力与应变的关系 求助
已经有11人回复
- 求应力-应变曲线分析
已经有8人回复
- 泡沫金属的应力-应变曲线
已经有4人回复
- 怎样输出应力值
已经有6人回复
- 拉伸使用引伸计的问题
已经有16人回复
- 已知材料的弹性模量和真实应力应变曲线，如何求屈服强度？
已经有11人回复
- 用origin画应力应变图中出现的问题
已经有7人回复
- 如何把力-位移曲线转化为应力-应变曲线
已经有11人回复
- 请教如何计算应力-应变曲线下的面积
已经有19人回复
- 铝合金应力—应变曲线
已经有21人回复
- 应力应变拉伸曲线的处理
已经有8人回复
- 拉伸应力应变曲线
已经有26人回复
- 应力应变曲线为什么是这样的？
已经有4人回复
- 关于工程应力应变曲线换算成真实应力应变曲线的问题...
已经有12人回复
- 【求助】低碳钢动态再结晶 应力应变曲线求教
已经有6人回复
- 【交流】为什么此拉伸应力-应变曲线如此怪异？
已经有32人回复
- 【求助】关于钢的应力-应变曲线图
已经有10人回复
- 【求助】MS计算应力应变曲线
已经有31人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
