|
|
[求助]
杆件材料的单轴拉伸问题
我在做复合材料的单轴拉伸时,遇到了以下问题:
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 |
|