24小时热门版块排行榜

返回列表

honghaier90

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 208.6
散金: 10
红花: 1
帖子: 52
在线: 46小时
虫号: 2581277
注册: 2013-08-05
性别: MM
专业: 高分子合成化学

[求助] 杆件材料的单轴拉伸问题

我在做复合材料的单轴拉伸时，遇到了以下问题：
1、当加载速度很快时，变形区域的温度急剧升高，我需要做温度控制吗？
2、我应用了  velocity right set 0.1 0.0 0.0
                  velocity mobile ramp vx 0.0 0.1 x 2 198 sum yes
加载方式，我不理解的是：这种加载方式是首先给mobile区域一个0—0.1线性增加的速度之后，这个速度指令就立即消失,然后速度随着被加载部位的速度的变化而变化？还是mobile区域的速度一直是0——0.3线性变化？
3、在计算应变时，我需要首先得到材料的原长，这个原长是取材料能量最小化之前的长度，还是能量最小化后的速度？

以下是我的in文件，希望求得大神指点。

dimension 2
boundary s p p

atom_style atomic
neighbor 0.3 bin
neigh_modify every 1 delay 0 check yes
lattice hex 0.91648

read_data    data.new3

pair_style lj/cut 1.5
pair_coeff 1 1 1.0 1.0  1.5
pair_coeff    2 2 100.0  1.0  1.5
pair_coeff    1 2 1.0  1.0 1.5

region left block INF 2    INF  INF  INF  INF
region right  block 198 INF INF  INF  INF  INF
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 100 1000

variable       tmp equal  lx
variable       L  equal ${tmp}
print          "initial  length , L:${L}"

#deformation
reset_timestep  0
velocity       left set 0.0 0.0 0.0
compute          new mobile temp
compute       p all pressure thermo_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       strain equal "(lx - v_L)/v_L"
thermo_style custom step temp vol lx ly lz pxx v_x c_strex v_strain
thermo_modify temp new

velocity mobile create 0.01 887723 temp new
velocity right  set 0.1 0.0 0.0
velocity mobile ramp vx 0.0 0.1 x 2 198 sum yes



fix 1 all nve
fix 2 boundary setforce    0.0 0.0 0.0
fix          3 mobile temp/rescale 100  0  0.05  0.002  0.1

timestep 0.001
thermo 200

dump 1 all atom 200 dump.right.lammpstrj
run 20000
请求各位大神的帮助。以下是我的应力——应变曲线。

JTO9O_]3]$[ZW~OCY@MF~QR.jpg

回复此楼