版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

zhangjiongl

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 26.8
帖子: 8
在线: 12.5小时
虫号: 1321391
注册: 2011-06-12
专业: 固体力学

[求助] 拉伸应力应变曲线怎么里面还有负的

照着手册上有个改的， in 文件
boundary p p p
units       metal
atom_style atomic
read_data data.pi
neighbor 2.0 bin
neigh_modify delay 1 check yes

# LJ potentials
pair_style eam
pair_coeff * * Cu_u3.eam
# ------------------------- SETTINGS ---------------------------------
compute csym all centro/atom fcc
compute peratom all pe/atom
######################################
# EQUILIBRATION
reset_timestep 0
timestep 0.001
velocity all create 300 12345 mom yes rot no
fix 1 all npt temp 300 300 1 iso 0 0 1 drag 1
# Set thermo output
thermo 1000
thermo_style custom step lx ly lz press pxx pyy pzz pe temp
# Run for at least 10 picosecond (assuming 1 fs timestep)
run 20000
unfix 1
# Store final cell length for strain calculations
variable tmp equal "lx"
variable L0 equal ${tmp}
print "Initial Length, L0: ${L0}"
######################################
# DEFORMATION
reset_timestep 0
fix  1 all npt temp 300 300 1 y 0 0 1 z 0 0 1 drag 1
variable srate equal 1.0e10
variable srate1 equal "v_srate / 1.0e12"
fix  2 all deform 1 x erate ${srate1} units box remap x
# Output strain and stress info to file
# for units metal, pressure is in [bars] = 100 [kPa] = 1/10000 [GPa]
# p2, p3, p4 are in GPa
variable strain equal "(lx - v_L0)/v_L0"
variable p1 equal "v_strain"
variable p2 equal "-pxx/10000"
variable p3 equal "-pyy/10000"
variable p4 equal "-pzz/10000"
fix def1 all print 100 "${p1} ${p2} ${p3} ${p4}" file Al_SC_100.def1.txt screen no
# Use cfg for AtomEye
dump 1 all cfg 250 dump.tensile_*.cfg id type xs ys zs c_csym c_peratom fx fy fz
dump_modify 1 element Al
# Display thermo
thermo  1000
#thermo_style custom step v_strain temp v_pxx v_pyy v_pzz ke pe press
thermo_style custom step v_strain temp v_p2 v_p3 v_p4 ke pe press
run  20000
data 文件在附件
stress-strain.jpg

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : data.pi

2013-03-24 12:30:47, 5.3 K

1楼 2013-03-24 12:30:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

应助: 101 (高中生)
金币: 2939.2
红花: 6
帖子: 670
在线: 255.3小时
虫号: 2121784
注册: 2012-11-12
专业: 动力学与控制

【答案】应助回帖

感谢参与，应助指数 +1

你的应力是怎么计算的？另外，你的模型视图是怎样的？

赞一下

回复此楼

2楼2013-03-25 14:40:55

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 zhangjiongl 的主题更新