24小时热门版块排行榜    

查看: 671  |  回复: 0

honghaier90

新虫 (小有名气)

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

我在做复合材料的单轴拉伸时,遇到了以下问题:
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
回复此楼
优秀是一种习惯!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 honghaier90 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 我不理解! +18 Edward_pc 2026-08-26 26/1300 2026-09-02 11:01 by shanzhipu
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +4 瞬息宇宙 2026-09-01 8/400 2026-09-02 10:35 by 瞬息宇宙
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
信息提示
请填处理意见