二维圆弧阵列结构共面冲击下的平台应力研究_孙德强.pdf

第4 1卷 第3期 陕西科技大学学报 V o l.4 1N o.3 2 0 2 3年6月 J o u r n a l o fS h a a n x iU n i v e r s i t yo fS c i e n c e&T e c h n o l o g y J u n.2 0 2 3*文章编号:2 0 9 6-3 9 8 X(2 0 2 3)0 3-0 1 4 7-0 9二维圆弧阵列结构共面冲击下的平台应力研究孙德强1,本金翠1,张 超2,许亚利3,葛 凤1,常 露1,李芸娴1(1.陕西科技大学 中国轻工业功能印刷与运输包装重点实验室3 S包装新科技研究所 轻化工程国家级实验教学示范中心,陕西 西安 7 1 0 0 2 1;2.西安西电变压器有限责任公司,陕西 西安 7 1 0 0 7 7;3.西安国际港务区新农小学,陕西 西安 7 1 0 0 2 6)摘 要:借助数值模拟和理论分析的方法研究了二维圆弧阵列结构在不同冲击速度下的共面平台应力.利用有限元分析软件AN S Y S/L S-D YNA建立了二维圆弧阵列结构共面冲击下的有限元模型,并用压缩试验验证了有限元模型的可靠性.数值模拟结果表明该结构在低速冲击下呈现出“波浪”形坍塌,中速冲击下呈现出“一”字形加“波浪”形坍塌,高速冲击下表现为稳定的“一”字形坍塌.基于线性回归理论推导出不同变形模式间的临界速度表达式,发现其与相对密度呈线性关系.根据特征单元变形模式,推导出了预测二维圆弧阵列结构在不同冲击速度下平台应力的理论模型,理论值与有限元计算值吻合较好,最大相对误差为9.7 9%.此外,大量参数化计算后发现:同一结构参数下,二维圆弧阵列结构共面平台应力与冲击速度的平方呈线性关系,并随冲击速度的增大而增大;同一冲击速度下,其共面平台应力随着壁厚半径比的增大而增大;同一相对密度下,其共面平台应力随着圆弧角度的增大而增大.关键词:二维圆弧阵列结构;共面;平台应力;有限元模型;理论模型中图分类号:T B 4 8 5.1 文献标志码:AS t u d yo np l a t e a us t r e s s e so f t w o-d i m e n s i o n a l c i r c u l a ra r cl a t t i c e s s t r u c t u r eu n d e r i n-p l a n e i m p a c t l o a d i n g sS UND e-q i a n g1,B E NJ i n-c u i1,Z HANGC h a o2,XUY a-l i3,G EF e n g1,CHANGL u1,L IY u n-x i a n1(1.K e yL a b o r a t o r yo fF u n c t i o n a lP r i n t i n ga n dT r a n s p o r tP a c k a g i n go fC h i n aN a t i o n a lL i g h t I n d u s t r y,3 SR e-s e a r c hI n s t i t u t eo fN o v e lP a c k a g i n gS c i e n c ea n dT e c h n o l o g y,N a t i o n a lD e m o n s t r a t i o nC e n t e r f o rE x p e r i m e n t a lL i g h tC h e m i s t r yE n g i n e e r i n gE d u c a t i o n,S h a a n x iU n i v e r s i t yo fS c i e n c e&T e c h n o l o g y,X i a n7 1 0 0 2 1,C h i n a;2.X i a nX DT r a n s f o r m e rC o.,L t d.,X i a n7 1 0 0 7 7,C h i n a;3.X i n n o n gP r i m a r yS c h o o lo fX i a nI n t e r n a t i o n a lT r a d ea n dL o g i s t i c sP a r k,X i a n7 1 0 0 2 6,C h i n a)A b s t r a c t:T h e i n-p l a n ep l a t e a us t r e s s e s o f t w o-d i m e n s i o n a l c i r c u l a r a r c l a t t i c e s s t r u c t u r eu n d e rd i f f e r e n t i m p a c t v e l o c i t i e sw e r e a n a l y z e d t h r o u g hn u m e r i c a l s i m u l a t i o n s a n d t h e o r e t i c a l a n a l y-s e s.T h e f i n i t ee l e m e n tm o d e l f o r i n-p l a n e i m p a c t a n a l y s i so f t w o-d i m e n s i o n a l c i r c u l a r a r c l a t-t i c e ss t r u c t u r ew a se s t a b l i s h e dw i t ht h es o f t w a r eAN S Y S/L S-D YNA,a n dv e r i f i e dt h er e l i a-b i l i t yt h r o u g ht h ec o m p r e s s i o nt e s tr e s u l t s.N u m e r i c a ls i m u l a t i o n sr e v e a l e dt h e“w a v e”d e-f o r m a t i o nm o d eu n d e r l o wi m p a c tv e l o c i t i y,t h e“”a n d“w a v e”d e f o r m a t i o nm o d e su n d e rm e d i u mi m p a c tv e l o c i t ya n dt h e“”d e f o r m a t i o nm o d eu n d e rh i g h i m p a c tv e l o c i t y.B a s e do nt h e l i n e a r r e g r e s s i o n t h e o r y,t h e e x p r e s s i o no f t h e c r i t i c a l v e l o c i t i e sb e t w e e nd i f f e r e n t d e f o r m-*收稿日期:2 0 2 2-1 0-0 9基金项目:国家自然科学基金项目(5 1 5 7 5 3 2 7);陕西省教育厅重点实验室及基地科研项目(1 6 J S 0 1 4)作者简介:孙德强(1 9 7 6),男,山东沂水人,教授,博士,研究方向:轻质绿色智能共享循环可持续性包装系统的结构、性能与设计DOI:10.19481/ki.issn2096-398x.2023.03.010陕西科技大学学报第4 1卷a t i o nm o d e sw e r e f i t t e d,a n d i tw a s f o u n dt h a t i t i s l i n e a r l yr e l a t e dt ot h er e l a t i v ed e n s i t yo ft w o-d i m e n s i o n a l c i r c u l a ra r c l a t t i c e ss t r u c t u r e.T h e o r e t i c a lm o d e l sw e r ed e r i v e db a s e do nd e-f o r m a t i o nm o d eo f t h e r e p r e s e n t a t i v eu n i t f o rp r e d i c t i n g t h ep l a t e a us t r e s so f t w o-d i m e n s i o n a lc i r c u l a ra r c l a t t i c e s s t r u c t u r eu n d e r d i f f e r e n t i m p a c t v e l o c i t i e s.G o o da g r e e m e n tw a s f o u n db e-t w e e nt h e t h e o r e t i c a l p r e d i c t i o n sa n dt h en u m e r i c a l c a l c u l a t e dr e s u l t sw i t ham a x i m u mr e l a-t i v ee r r o ro f 9.7 9%.I na d d i t i o n,f r o mal o to fp a r a m e t r i cc a l c u l a t i o n s,i tw a sf o u n dt h a t t h ei n-p l a n ep l a t e a us t r e s so f t w o-d i m e n s i o n a l c i r c u l a ra r c l a t t i c e ss t r u c t u r e i s l i n e a r l yr e l a t e dt ot h es q u a r eo f i m p a c tv e l o c i t y,a n di n c r e a s e sw i t ht h ei n c r e a s eo f i m p a c tv e l o c i t y,w h e nt h es t r u c t u r a l p a r a m e t e r sa r ef i x e d.U n d e rt h es a m ei m p a c tv e l o c i t y,t h ei n-p l a n ep l a t e a us t r e s si n c r e a s e sw i t ht h e i n c r e a s eo f t h er a t i oo f c e l lw a l l t h i c k n e s s t or a d i u s.F o r t h es a m er e l a t i v ed e n s i t y,t h e i n-p l a n ep l a t e a us t r e s s i n c r e a s e sw i t ht h e i n c r e a s eo f t h ea r ca n g l e.K e yw o r d s:t w o-d i m e n s i o n a l c i r c u l a r a r c l a t t i c e s s t r u c t u r e;i n-p l a n e;p l a t e a us t r e s s;f i n i t e e l e-m e n tm o d e l;t h e o r e t i c a lm o d e l0 引言蜂窝结构具有轻量化、耐冲击和缓冲吸能性好等特点,在航空航天、汽车防撞、建筑和包装等领域中被制作成各种缓冲吸能装置,来吸收外部冲击载荷产生的能量1-4.相关研究表明蜂窝材料孔穴形状不同,相应的能量吸收性能也会大不相同5,6.W a n g等7研究了七种不同构型蜂窝结构的共面力学 行 为,通 过 比 较 不 同 的 蜂 窝 结 构,得 出K a g o m e型蜂窝结构的力学性能是最优越的.刘颖等8对三角形和四边形蜂窝结构进行了共面冲击有限元计算,发现单位质量的三角形蜂窝比四边形蜂窝具有更高的平台应力和更高的能量吸收能力,并且从压缩变形过程说明相关成因.为了满足不同的设计需求,蜂窝夹芯层的种类变得越来越多.新型蜂窝夹芯层的研究也越来越多.G u o等9在传统的双V形的基础上提出双U形蜂窝,借助理论、数值模拟和试验相结合的方法研究其缓冲性能.孙德强等1 0利用有限元法研究了凹六边形蜂窝的准静态压缩行为,探究了壁厚边长比和扩展角的影响,推导了有关计算公式.Q i等1 1提出一种双圆弧代替凹六边形侧棱结构,极大提高了蜂窝材料的缓冲性能.Z h a n g等1 2以柚子皮的微观结构为蓝本,设计了一种仿柚子皮层级蜂窝结构,这种仿生设计增强了蜂窝结构共面和异面的耐撞性.Y a n g等1 3模仿马蹄的形状设计了一系列马蹄形蜂窝结构,证明了马蹄形状的设计方案可以增大蜂窝的平台应力,提高吸能性能.C h e n等1 4研究了一种新型负泊松比蜂窝结构的准静态力学特性,通过试验和模拟相结合的方式研究结构参数对冲击性能的影响.这些研究表明,改变传统蜂窝孔穴形状可以得到综合力学性能更优的蜂窝结构材料.在对新构型蜂窝结构力学行为,特别是在大变形条件下的力学行为进行研究,已成为研究人员密切关注的一个热点话题.同时,对产生的各种新构型蜂窝的力学参数要从理论上开展相应研究,并通过试验和有限元数值模拟加以验证,比如共面动态平台应力与各结构参数之间的关系.二维圆弧阵列结构是一种新型蜂窝结构材料,本文首先建立其共面冲击分析的有限元模型,并通过压缩试验验证有限元模型的可靠性.通过大量参数化模拟计算,揭示其变形模式并得到不同变形模式间临界速度的计算公式.其次,基于二维圆弧阵列结构的特征单元变形模式,推导该结构在不同冲击速度下的平台应力理论模型,并借助有限元计算结果验证理论模型的正确性.最后,利用理论模型和有限元模拟,分析冲击速度、壁厚半径比和圆弧角度对该结构共面平台应力的影响规律.1 有限元模型及可靠性分析1.1 有限元模型参考S u n等1 5建立有限元模型的思路和方法,利用有限元分析软件AN S Y S/L S-D YNA建立二维圆弧阵列结构有限元模型,如图1所示.其中R为圆弧半径,b为样品厚度(胞元孔深),为圆弧角度,t为圆弧壁厚.二维圆弧阵列结构置于刚性上压板P1和支撑板P2之间,P1只允许在y方向上移动,限制其余自由度,P2的所有自由度均被限制为0,使得P1以恒定速度v向下运动,压缩样品变形.对样品采用5个积分点的B e l y t s c h k o-T s a y壳单元S h e l l 1 6 3进行网格划分;通过收敛性分析,确定壳单元边长为841第3期孙德强等:二维圆弧阵列结构共面冲击下的平台应力研究0.6m m.整个模型定义单面无摩擦自动接触,样品与刚板之间定义面面自动接触,静动摩擦系数分别为0.2和0.1 5.基材选用对应变率不敏感的双线性应变硬化材料,各项力学参数为:杨氏模量Es为6 8.9 7G P a,屈服应力y s为2 9 2M P a,正切模量Et a n为6 8 9.7M P a,泊松比Vs为0.3 5,基材密度为s为27 0 0k g/m3.样品底部所有节点固定,其它节点只能沿共面方向移动.为了消除尺寸效应对其仿真结果的影响,需要对样品施加周期性边界条件,即将二维圆弧阵列结构左右两侧各节点沿x方向的位移、绕y轴与z轴的转动 等自由度约 束为0,如 图1(d)所示.图1 二维圆弧阵列结构共面冲击有限元模型1.2 二维圆弧阵列结构试验样品图2为二维圆弧阵列结构样品.其制作首先是将铝合金圆管采用激光切割,而后利用高强度粘合剂粘合而成.其总体尺寸为:1 0 0mm6 0mm5 0mm.其中R=1 0mm、b=5 0mm、=1 8 0、t=0.5mm.图2 二维圆弧阵列结构试验样品1.3 模型可靠性分析为了验证有限元建模方法的可靠性,需要采用压缩试验来验证.为了获得相应有限元模拟的材料参数,需要对压缩试验中使用试样的基体材料进行拉伸试验.1.3.1 基体材料拉伸试验构成二维圆弧阵列结构的基体材料为6 0 6 3-T 6型铝合金,按照G B/T2 2 8.1-2 0 2 11 6中金属材料室温拉伸试验的规定对基体材料进行拉伸试验.拉伸试样尺寸如图3(a)所示.采用激光切割法制备的实际拉伸试样如图3(b)所示.采用万能材料试验机以2mm/m i n的拉伸速率对三组试样进行拉伸,如图3(d)所示.拉伸后的试样如图3(c)所示.最终得到的应力-应变曲线如图4所示,由图4可以看出6 0 6 3-T 6型铝合金属于典型的理想弹塑性材料,可在有限元模型中采用双线性应变硬化材料模型.对拉伸试验结果进行处理,可得到试验样品基体材料的各项力学参数为:杨氏模量Es为6 9G P a,屈服应力y s为2 1 7M P a,正切模量Et a n为6 7 5M P a,泊松比Vs为0.3 3,基材密度为s为27 0 0k g/m3.图3 基材拉伸试样及试验941陕西科技大学学报第4 1卷图4 试样拉伸应力-应变曲线1.3.2 压缩试验压缩试验采用万能材料试验机进行,试验样品居中放置于上下刚性压板之间,如图5所示.上压板以1 2mm/m i n的速度向下压缩试样,采用相机录制试样变形全过程.图5 二维圆弧阵列结构的压缩试验1.3.3 模型可靠性分析建立与压缩试验样品尺寸和基材参数完全等效的有限元模型,从变形模式和应力-应变曲线两方面对试验结果和仿真结果进行对比分析.(1)变形模式图6表示压缩过程中试样在某一时刻的变形模式以及同一时刻的仿真模拟轮廓.可以看出,试验所得到的变形模式与仿真所得到的变形模式具有很好的一致性.在加载初期,以上部最中间胞元为顶点,出现倒“V”形局部坍塌带(图6(b);然后倒“V”形坍塌带的两个分支向附近单元扩展,引起更多的胞元发生坍塌,坍塌带逐渐变宽(图6(c)(d);最后一直没有完全坍塌的中部胞元发生坍塌后,样品达到密实(图6(d)(e).图6 压缩试验与仿真变形模式对比(2)应力-应变曲线图7表示压缩试验与有限元仿真所得到的应力-应变曲线对比图,两者具有很好的一致性.尽管此试验载荷条件是低速条件下的准静态压缩,但本研究所用基材(6 0 6 3-T 6型铝合金)是应变率不敏感的双线性应变硬化材料,可推出在高速载荷条件下试验和有限元模拟结果也会有很好的吻合性.据此,可证明前述所建有限元模型的可靠性.图7 压缩试验与仿真应力-应变曲线对比2 共面变形模式与临界速度2.1 变形模式与其它二维多孔缓冲材料一样,二维圆弧阵列结构在低速、中速和高速共面冲击载荷作用下,会表现出三种不同的变形模式1 5,1 7.2.1.1 低速冲击载荷下的变形模式图8为二维圆弧阵列结构在低速冲击载荷作用下的变形模式(v=1m/s、R=3mm、=1 2 0 和t=0.1mm).冲击初期,出现了一条“波浪”形局部坍塌带(图8(a);随着上压板继续往下移动,该条051第3期孙德强等:二维圆弧阵列结构共面冲击下的平台应力研究“波浪”形坍塌带向周边扩展而变宽,同时在其上方出现了第二条“波浪”形局部坍塌带(图8(b);随着时间的推移,两条坍塌带中间的单元逐步被压溃,使得两条坍塌带接触并逐渐融合在一起(图8(c);最终样品达到密实(图8(d).图8 二维圆弧阵列结构低速冲击载荷下的变形模式2.1.2 中速冲击载荷下的变形模式图9为二维圆弧阵列结构在中速冲击载荷作用下的变形模式(v=6m/s、R=3mm、=1 2 0 和t=0.1mm).在这一变形模式下,样品顶部第一行单元首先出现“一”字形坍塌带(图9(a);随着上压板的向下移动,中间部分出现一条类似低速变形模式下的“波浪”形局部坍塌带(图9(b);一定变形过后,各坍塌带中间出现更多的局部坍塌带,使得更多的单元局部密实(图9(c);随着上压板继续下移,各局部坍塌带连成一片,直到样品被完全压实(图9(d).可看出这一变形模式下,除了加载初期顶部短时间的“一”字形坍塌带外,与低速冲击载荷下的变形模式相似.图9 二维圆弧阵列结构中速冲击载荷下的变形模式2.1.3 高速冲击载荷下的变形模式图1 0为二维圆弧阵列结构在高速冲击载荷作用下的变形模式(v=7 0m/s、R=3mm、=1 2 0 和t=0.2 5mm).这一变形模式下,惯性效应占据了主导地位,下层单元还未发生明显变形的条件下,顶部首层单元已经完全坍塌,并出现“一”字形坍塌带(图1 0(a);随后“一”字形坍塌带随着上压板的下压,保持同种变形模式持续向下推进(图1 0(b)(c),直至最后完全被压实(图1 0(d).图1 0 二维圆弧阵列结构高速冲击载荷下的变形模式由上述分析可以看出,可以用不同的变形模式来区分冲击速度的高低.中速冲击载荷作用下的变形模式属于低速和高速冲击载荷作用下的混合变形模式,表明在从低速模式向高速模式的过渡,与W e i等1 7就 新 型 星 型 蜂 窝 变 形 模 式 的 结 论 相一致.2.2 临界速度相关研究表明1 7-1 9,相对密度和冲击速度对蜂窝材料的变形模式起着重要作用.为了揭示相对密度和冲击速度对二维圆弧阵列结构变形模式的影响规律,将不同相对密度和冲击速度下的变形模式列于图1 1.其中,二维圆弧阵列结构的相对密度公式为:=t2Rs i n21-c o s2(1)图1 1 二维圆弧阵列结构不同冲击速度下的变形模式由图1 1可以看出,二维圆弧阵列结构不同变形模式间的临界速度与相对密度近似呈线性关系,与W e i等1 7、Q i等1 8和A n等1 9得出的结论一致.根据最小二乘法,通过回归分析得到了不同变形模式间的临界速度经验公式分别为:v1=5 3.0 6+2.7 3(2)v2=1 3 3.8 3+1 3.0 4(3)151陕西科技大学学报第4 1卷式(2)、(3)中:v1和v2分别为低速冲击载荷下的变形模式向中速冲击载荷下的变形模式、中速冲击载荷下的变形模式向高速冲击载荷下的变形模式转变的临界速度.3 平台应力理论计算公式平台应力p的大小决定了蜂窝材料的能量吸收能力,因此对平台应力的研究很有必要.由于在中速冲击载荷下的变形模式除了顶部短时间的“一”字形坍塌带外,更多的表现为低速冲击载荷下的变形模式.参考W e i等1 7的研究方法,可将中速冲击下的平台应力近似等于相应低速冲击下的平台应力.因此本节根据不同的变形模式,在建立该结构准静态平台应力理论模型的基础上,只研究低速和高速冲击下平台应力的理论模型.3.1 准静态平台应力由于v=1m/s时,二维圆弧阵列结构动能占比很小,与Q i等1 1的研究方法一样,可以将v=1m/s视为准静态.在准静态变形模式中提取出一个特征胞元,其变形前后的形状如图1 2所示.其中,特征胞元的原始高度和宽度为H0和L0,密实化之后高度和宽度为HD和LD.图1 2 特征胞元的准静态变形模式根据图1 2(b),可假设只有4段圆弧发生变形(弧N C、弧MC、弧K D和弧L B),其中K、L、M和N分别为弧A D、A B、B C和C D的中点.各段圆弧在发生变形后曲率变化不大,因此可以假设塑性铰只产生于圆弧两端附近1 1,故一个胞元内有8个塑性铰,用空心圆圈表示于图1 2(c)中.各段圆弧发生变形后,塑性铰转动的角度可以用图1 3表示.其中弧L B和弧K D中每个塑性铰转动的角度a=/4,弧N C和弧MC中每个塑性铰转动的角度b=a/2=/8.图1 3 特征胞元准静态变形下塑性铰转动的角度因此,塑性铰转动的总角度和特征胞元的高度变化量H分别为:=4(+)=44+8(4)H=H0-HD=R-Rc o s2-Rc o s4+Rc o s8(5)准静态载荷下可忽略动能的影响1 1,根据塑性耗散理论,外力做功W等于塑性铰产生的塑性耗散EP,可以得出下式:W=0L0bH=EP=48+4MP=32 MP(6)式(6)中:MP为塑性矩,其值为MP=y sb t2/4,原始宽度L0=2Rs i n(/2).根据公式(5)和公式(6)可得到准静态下的平台应力表达式为:0=3 y st21 6R2s i n21-c o s2-c o s4+c o s8(7)3.2 低速冲击时的平台应力假设由样品与冲击端和固定端的接触力标准化后的平台应力分别为1和2,在准静态载荷下1=2.但动态冲击载荷下由于惯性效应,二者不相等,2不受冲击速度影响,其值近似等于准静态平台应力01 8,1 9.按照乔锦秀2 0的方法来研究在低速冲击载荷下(1m/s vv1)的平台应力.选取相邻两层胞元和作为一个特征单元来分析低速冲击载荷下的平台应力.特征单元的变形模式如图1 4所示.特征单元变形前原始高度为H,变形至T1时刻和T2时刻的高度分别为H1和H2.251第3期孙德强等:二维圆弧阵列结构共面冲击下的平台应力研究图1 4 低速冲击载荷下特征单元的变形模式根据动量守恒原理,1和2从T1时刻到T2时刻的冲量差值等于特征单元的动量变化.即:L0bT2T1(1-2)dT=P(8)特征单元从T1到T2时刻动量改变量为:P=P2+P2-P1-P1(9)式(9)中:Pji表示胞元j在i时刻的动量(j=,;i=1,2).由于在T1时刻,胞元发生初始变形,胞元未发生变形(图1 4(b)和图1 4(e);在T2时刻,胞元达到密实,胞元产生与胞元在T1时刻相同的变形(图1 4(c)和图1 4(f).所以,T1时刻胞元动量为0,T2时刻胞元动量等于T1时刻胞元的动量.即:P1=0;P2=P1(1 0)根据公式(9)、(1 0)可以得到:P=P2=m v=2 R t bsv(1 1)式(1 1)中:m为胞元的质量.由于速度是恒定的,因此T1时刻到T2时刻的时间间隔T为:T=T2-T1=H1-H2v=R1-c o s2-c o s4+c o s8v(1 2)由公式(8)、公式(1 1)和公式(1 2)可得:L0b(1-2)R1-c o s2-c o s4+c o s8v=2 R t bsv(1 3)将L0=2Rs i n(/2)、2=0代入公式(1 3),可以得到低速冲击载荷下的平台应力L的理论公式为:L=0+tsRs i n21-c o s2-c o s4+c o s8v2(1 4)3.3 高速冲击时的平台应力与低速冲击下平台应力的推导方法一样2 0,在高速冲击载荷下(vv2),选取相邻两层胞元作为一个特征单元来分析其在共面高速冲击下的平台应力,如图1 5所示.根据公式(8)(1 3)可以得出高速冲击载荷下的平台应力H的理论公式为:H=0+ts2Rs i n21-c o s2v2(1 5)图1 5 高速冲击载荷下特征单元的变形模式3.4 平台应力理论计算公式验证图1 6为二维圆弧阵列结构不同冲击速度下平台应力理论值与仿真值的对比图.351陕西科技大学学报第4 1卷图1 6 平台应力理论值与仿真值的对比由图1 6可以看出,理论值与仿真值吻合较好.并用表1给出了平台应力理论值与仿真值的相对误差,两者的误差1 0%.证明公式(7)、公式(1 4)和公式(1 5)可以较好地预测二维圆弧阵列结构在共面不同冲击速度下的平台应力,同时也证明了前述有限元模型的可靠性.表1 平台应力理论值与仿真值的相对误差相对密度冲击速度/(m/s)仿真值/MP a理论值/MP a相对误差/%0.0 510.1 4 30.1 5 79.7 930.1 4 70.1 5 98.1 61 0 01.4 9 01.5 7 15.4 41 2 52.2 3 32.3 6 65.9 60.0 810.2 2 70.2 4 57.9 330.2 3 70.2 4 95.0 61 0 02.3 8 72.4 2 21.4 71 2 53.5 6 23.6 4 62.3 60.1 310.3 6 30.3 8 66.3 430.3 6 70.3 9 16.5 41 0 03.5 8 13.7 9 96.0 91 2 55.4 3 45.7 1 95.2 44 结果讨论与分析4.1 冲击速度对平台应力的影响图1 7表示二维圆弧阵列结构冲击速度与平台应力的关系曲线.可以明显地看出平台应力理论值与仿真值吻合度较高,都随着冲击速度的增大而增大.且中速冲击下的平台应力与低速冲击和高速冲击下的平台应力可近似为在一条曲线上,可视为与冲击速度的平方呈线性关系.图1 7 冲击速度与平台应力关系曲线4.2 壁厚半径比对平台应力的影响图1 8为圆弧角度相同时,二维圆弧阵列结构壁厚半径比与平台应力关系曲线.可以看出理论值与仿真值具有很好的一致性.同一圆弧角度下,平台应力随着壁厚半径比的增大而增大.这是因为壁厚半径比越大,该结构相对密度越高,越不容易发生屈曲且抗冲击能力也越强.图1 8 壁厚半径比与平台应力关系曲线4.3 圆弧角度对平台应力的影响图1 9为相对密度一致时,圆弧角度与平台应力的关系曲线.图1 9 圆弧角度与平台应力关系曲线451第3期孙德强等:二维圆弧阵列结构共面冲击下的平台应力研究由图1 9可以看出,理论值与仿真值具有很好的一致性,平台应力随着圆弧角度的增大而增大.这是因为在相对密度一致时,随着圆弧角度的增大,二维圆弧阵列结构的壁厚、可变形空间和塑性铰转动的角度进一步增大,因此就需要更大的冲击力来达到与其它较小圆弧角度相同的应变.5 结论本文借助理论、仿真和试验的方法研究了二维圆弧阵列结构共面冲击载荷条件下的平台应力,得到如下结论:(1)二维圆弧阵列结构在低速冲击下的变形模式为“波浪”形坍塌;中速冲击下的变形模式为“一”字形加“波浪”形坍塌;高速冲击下的变形模式为从顶部到底部稳定的“一”字形坍塌;并通过回归分析给出不同变形模式之间的临界速度经验公式,得出临界速度与相对密度呈线性关系.(2)运用塑性耗散理论和动量守恒原理推导出了二维圆弧阵列结构在不同冲击速度下共面平台应力的理论公式,并用有限元计算结果验证了理论公式的正确性.(3)同一结构参数下,二维圆弧阵列结构的共面平台应力与冲击速度的平方呈线性关系,并随冲击速度的增大而增大;冲击速度相同时,其共面平台应力随着壁厚半径比的增大而增大;相对密度一致时,其 共面平台应 力随着圆弧 角度的增 大而增大.参考文献1W a n gH,L uZ,Y a n gZ,e ta l.An o v e lr e-e n t r a n ta u x e t i ch o n e y c o m bw i t he n h a n c e di n-p l a n e i m p a c tr e s i s t a n c eJ.C o m p o s i t eS t r u c t u r e s,2 0 1 9,2 0 8:7 5 8-7 7 0.2F uJ,L i uQ,L i u f uK,e ta l.D e s i g no fb i o n i c-b a m b o ot h i n-w a l l e ds t r u c t u r e sf o re n e r g ya b s o r p t i o nJ.T h i n-W a l l e dS t r u c t u r e s,2 0 1 9,1 3 5:4 0 0-4 1 3.3L i uJY,L i uHT,A nMR.C r u s h i n gb e h a v i o r s o f n o v e l d i-a b o l os h a p e dh o n e y c o m b sw i t he n h a n c e de n e r g ya b s o r p-t i o np e r f o r m a n c eJ.I n t e r n a t i o n a lJ o u r n a lo fM e c h a n i c a lS c i e n c e s,2 0 2 2,2 2 9:1-1 6.4X i eS,C h e nP,W a n gN,e ta l.C r a s h w o r t h i n e s ss t u d yo fc i r c u l a rt u b e ss u b j e c t e dt or a d i a le x t r u s i o nu n d e rq u a s i-s t a t i c l o a d i n gJ.I n t e r n a t i o n a l J o u r n a lo fM e c h a n i c a lS c i-e n c e s,2 0 2 1,1 9 2:1-1 7.5胡玲玲,蒋 玲.胞孔构型对金属蜂窝动态力学性能的影响机理J.爆炸与冲击,2 0 1 4,3 4(1):4 1-4 6.6C o m b e s c u r eC,E l l i o t tRS,T r i a n t a f y l l i d i sN.D e f o r m a t i o np a t t e r n sa n dt h e i r s t a b i l i t y i nf i n i t e l ys t r a i n e dc i r c u l a r c e l lh o n e y c o m b sJ.J o u r n a lo f t h eM e c h a n i c sa n dP h y s i c so fS o l i d s,2 0 2 0,1 4 2:1-1 6.7W a n gAJ,M cD o w e l lD L.I n-p l a n es t i f f n e s sa n dy i e l ds t r e n g t ho fp e r i o d i cm e t a lh o n e y c o m b sJ.J o u r n a lo fE n-g i n e e r i n g M a t e r i a l s a n d T e c h n o l o g y,2 0 0 4,1 2 6(2):1 3 7-1 5 6.8刘 颖,张新春.胞元微拓扑结构对蜂窝材料面内冲击性能的影响J.爆炸与冲击,2 0 0 8,2 8(6):4 9 4-5 0 2.9G u oMF,Y a n gH,M aL.D e s i g na n da n a l y s i so f 2 Dd o u b-l e-U a u x e t i ch o n e y c o m b sJ.T h i n-W a l l e d S t r u c t u r e s,2 0 2 0,1 5 5:1-1 6.1 0孙德强,张小强,王复会,等.凹六边形蜂窝芯材共面准静态力学行为J.陕西科技大学学报(自然科学版),2 0 1 5,3 3(5):1 5 3-1 5 6.1 1Q iC,J i a n gF,R e m e n n i k o vA,e t a l.Q u a s i-s t a t i c c r u s h i n gb e h a v i o ro fn o v e l r e-e n t r a n t c i r c u l a ra u x e t i ch o n e y c o m b sJ.C o m p o s i t e sP a r tB:E n g i n e e r i n g,2 0 2 0,1 9 7:1-1 2.1 2Z h a n gW,Y i nS,Y uTX,e t