24小时热门版块排行榜    

查看: 670  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见