CN108038318B - 变截面金属点阵结构初始刚度及塑性破坏强度计算算法 - Google Patents
变截面金属点阵结构初始刚度及塑性破坏强度计算算法 Download PDFInfo
- Publication number
- CN108038318B CN108038318B CN201711374124.6A CN201711374124A CN108038318B CN 108038318 B CN108038318 B CN 108038318B CN 201711374124 A CN201711374124 A CN 201711374124A CN 108038318 B CN108038318 B CN 108038318B
- Authority
- CN
- China
- Prior art keywords
- unit cell
- bending moment
- tangential force
- lattice structure
- variable cross
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种变截面金属点阵结构初始刚度及塑性破坏强度计算算法,包括以下步骤:s1:通过建立平面坐标系,并结合BCC变截面金属点阵单胞杆的剖面尺寸得到BCC变截面金属点阵单胞杆截面的半径表达式;s2:通过建立空间坐标系并结合胡克定律及弯压组合变形公式得到所述单胞杆的结点所受的空间切向力、弯矩与空间位移的关系式;s3:通过能量守恒定律,结合所述空间切向力、弯矩与空间位移的关系式,得到一个关于所述单胞杆的结点处的切向力的一元二次方程,解该方程得到不含所述结点的空间位移的切向力表达式,并结合所述结点处的切向力与轴向力、弯矩的关系得到所述轴向力和弯矩的表达式;s4:依据胡克定律,结合步骤s3得出的所述切向力、弯矩、轴向力的表达式,得到BCC变截面金属点阵结构的初始刚度及塑性破坏强度。
Description
技术领域
本发明涉及新型变截面金属点阵结构初始刚度及塑性破坏强度计算算法,具体地,涉及一种通过变截面金属点阵结构的剖面尺寸、结点空间位移及受力,结合材料力学胡克定律、功能原理、平面几何等最终计算出变截面金属点阵结构的初始刚度及塑性破坏强度。
背景技术
目前,金属点阵结构的研究对象以体心立方(Body-centered Cubic;BCC)点阵结构为主,因为这类结构是由结点和结点间连接杆件单元组成的按一定规则重复排列构成的空间桁架结构,所以其构型简单且具有各向同性的特点,能较好地适应SLM成型工艺,制备方便可靠、破坏形式单一方便观察与分析。
但是BCC点阵结构的力学性能远远差于其他具有基础拓扑结构的点阵结构诸如BCCZ、FCCZ等,因此探讨从根本上解决BCC点阵现存的轻质化与承载能力间此消彼长的矛盾具有重要的意义。
申请号为CN201510229244.1的专利所公开的一种BCC变截面金属点阵结构,其截面直径是由变截面杆件在载荷作用下需满足的内力条件来精确确定的,实现了构成变截面杆件的材料根据变截面杆件内的应力的梯度增减,提高了材料的利用率,减少了应力集中,在保持BCC金属点阵结构现存优势的情况下同时提高了结构的比强度和比刚度。但是,这种新型的BCC变截面金属点阵结构的相关数学模型并未建立,不能预测其初始刚度及塑性破坏强度,阻碍了对该结构力学性能进一步的优化。而现存的BCC点阵结构的初始刚度及塑性破坏强度的数学模型主要针对的是具有匀质杆的BCC点阵结构,不适用于具有变截面杆的BCC金属点阵结构。
因此,本发明针对BCC变截面金属点阵结构初始刚度及塑性破坏强度所建立的数学模型能填补目前缺乏BCC变截面金属点阵结构理论预测模型的空白,能为该结构力学性能的进一步优化提供理论支持。
发明内容
有鉴于此,本发明提供一种通过变截面金属点阵结构的剖面尺寸、结点空间位移及受力,结合材料力学胡克定律、功能原理、平面几何等最终计算出变截面金属点阵结构的初始刚度及塑性破坏强度。
本发明的变截面金属点阵结构初始刚度及塑性破坏强度计算算法,包括以下步骤:s1:通过建立平面坐标系,并结合BCC变截面金属点阵单胞杆的剖面尺寸得到BCC变截面金属点阵单胞杆截面的半径表达式;s2:通过建立空间坐标系并结合胡克定律及弯压组合变形公式得到所述单胞杆的结点所受的空间切向力、弯矩与空间位移的关系式;s3:通过能量守恒定律,结合所述空间切向力、弯矩与空间位移的关系式,得到一个关于所述单胞杆的结点处的切向力的一元二次方程,解该方程得到不含所述结点的空间位移的切向力表达式,并结合所述结点处的切向力与轴向力、弯矩的关系得到所述轴向力和弯矩的表达式;s4:依据胡克定律,结合步骤s3得出的所述切向力、弯矩、轴向力的表达式,得到BCC变截面金属点阵结构的初始刚度及塑性破坏强度;
进一步的,在步骤s1中,所述半径表达式为:
进一步的,在步骤s2中,所述空间切向力、弯矩与空间位移的关系式为:
进一步的,在步骤s3中,所述切向力的一元二次方程为:
C9F1 2+C10F1+C11=0;
C9=-2R(2R2C8+C1);
且所述切向力F1的表达式为:
其中,
C12=2R2C8+C1;
所述轴向力N1及弯矩M1的表达式分别为:
σz为点阵单胞结构受到的压应力。
进一步的,在所述步骤s4中,所述初始刚度EplBCC为:
进一步的,所述塑性破坏强度σplBCC为:
本发明的有益效果:
1、本发明建立了BCC变截面金属点阵单胞的杆半径表达式,可通过控制表达式中的参数准确的控制杆的截面变化。
2.本发明结合功能原理及胡克定律等方法能准确的理论计算出BCC变截面金属点阵结构的初始刚度及塑性破坏强度。
附图说明
下面结合附图和实施例对本发明作进一步描述。
图1为BCC变截面金属点阵单胞及截面示意图;
图2为BCC变截面金属点阵单胞杆的示意图;
图3为BCC变截面金属点阵单胞的受力示意图;
图4为BCC变截面金属点阵单胞杆的空间受力及变形示意图;
图5为S7点在平面o'x'y'内的位移示意图;
图6为BCC变截面金属点阵单胞杆在平面o”x”y”内的受力及变形示意图;
图7为变密度与定密度样件理论计算及实验的应力应变对比图。
具体实施方式
1.变截面杆结构设计
选取单胞为边长是L的立方体进行BCC变截面点阵结构的初始刚度及塑性破坏强度的计算,如图1所示。取单胞中所有杆的尺寸是一样的,本文取杆lS7S9来分析。图2为变密度杆截面尺寸示意图,变密度杆通过半径为R、跨度为杆长l的圆弧线来控制半径均匀变化,通过几何关系可以得到圆弧半径R的表达式为:
以变密度杆其中一端面中点为原点,轴向及横向分别为x、y轴建立图4所示平面坐标系oxy,则圆弧圆心M点的坐标值为(l/2,-R-R1),并且圆弧上任意一点的坐标值符合下面关系:
则变密度杆在平面坐标系oxy下任意横截面的半径r(x)为:
圆心M到圆弧上任意一点的矢量与x轴的夹角α取值范围为α0~α1,且α1=π-α0。根据直角三角形的定理与几何关系可求得:
2变密度点阵结构力学性能分析
2.1初始刚度EplBCC
如图3所示,当BCC单胞受到压应力σz作用时,结点可以在空间中任意移动。由于整个单胞关于S9点中心对称,各杆力在S9点相互抵消,故S9点的节点位移为零。以S7S9杆为例,设S7S9杆与平面S5S6S7S8的夹角为θ,单胞底面对角线S5S7与边线S7S8的夹角为θ',S7点受到轴向力N1、切向力F1及弯矩M1的作用后在空间坐标系o'x'y'z'中三个方向的位移分别为(u,ν,w),如图4、5所示。
变形示意图见图6,以等效悬臂梁的固定点S9点为原点,杆的轴向及横向分别为x轴、y轴建立平面坐标系o”x”y”,杆受到轴向力N1、切向力F1及弯矩M1的作用后产生的平面轴向位移为λ1、挠度为ω1、截面转角为β1,图中实线部分为杆的初始状态,虚线部分为杆变形以后的状态。
基于以上的分析,通过S7点在空间坐标系o'x'y'z'中三个方向的位移(u,ν,w)可得杆的轴向位移为:
根据材料力学中的胡克定律可知在轴向力N1的作用下,变截面杆的轴向位移为:
结合式(5)和式(6)可得杆的轴向力N1为:
通过S7点在空间坐标系o'x'y'z'中三个方向的位移(u,ν,w)可得杆切向位移为:
在平面坐标系o”x”y”下,由于杆不仅受到切向力F1的作用,还受到弯矩M1的作用,因此根据结构弯压组合变形的方程可得杆在切向力及弯矩的作用下切向位移为:
另外,因为对称的关系,各结点处的弯曲力矩相同,所以结点S7、S9处的弯矩MS7、MS9相等且为M1,且
结合式(8)、(9)、(10)可得杆件所受的切向力F1及弯矩M1分别为:
其中,
通过S7点在空间坐标系o'x'y'z'中所受的切向力F1及弯矩M1可得S7在z'向所受的力Fz'为:
Fz'=F1 cosθ-N1 sinθ (12)
另外,由于单胞底面四个杆件的端点S5、S6、S7、S8同等的受到压应力σz的作用,因此各端点z'向受力Fz'均为σzL2/4。
由前述分析可知,杆的变形主要为轴向变形与横向弯曲变形,其弹性应变能包括轴向伸缩能与弯曲应变能,由功能原理与胡克定律可知,S7S9杆受到轴向力N1、切向力F1及弯矩M1作用时的应变能U为:
一个单胞里面有八根杆件,所以单胞的总应变能UBCT为单杆应变能的8倍既是UBCT=8U。另,BCC单胞z'方向上的压应力σz对单胞做的功Uw为:
Uw=2σzL2w (14)
根据能量守恒,单胞的应变能与外力给单胞做的功相等,即UBCT=Uw,结合式(13)、(14)得:
结合式(7)、式(11)可将F1的表达式可变为:
由式(12)及Fz'=σzL2/4得N1的表达式为:
结合式(10)、(15)、(16)及(17)可消掉u、ω得到关于F1的一元二次方程为:
C9F1 2+C10F1+C11=0 (18)
其中,
C9=-2R(2R2C8+C1);
通过式(18)求得切向力F1得:
其中,
C12=2R2C8+C1;
所以结合式(10)、(17)、(19)可得杆件所受的轴向力N1及弯矩M1分别为:
由于单胞具有对称性,因此单胞在z'方向上的总位移为杆端点S7在z'方向上的位移w的两倍,结合式(16)、(27)、(19)以及(20)可得BCC单胞在z'方向上的应变为:
本文着眼于分析受到z'向加载的BCC点阵单胞塑性破坏前的力学性能,即其弹性阶段内的力学性能,故其应力应变满足胡克定律。结合式(21),可得BCC单胞在受z'方向上应变时,单胞的初始刚度EplBCC为:
2.2塑性破坏强度σplBCC
BCC变密度单胞杆的塑形破坏主要考虑弯矩的作用,忽略轴向力及屈服面剪切的影响,当变密度杆的弯矩M1为杆的极限弯矩Mu时,此时点阵块所受到的压应力σz为结构的塑性破坏强度σplBCC。在求杆的极限弯矩Mu时,近似的认为变密度杆的弯矩M1的作用点在结点处,极限弯矩Mu的表达式中的半径取为R2。则杆的极限弯矩Mu为:
式中,σs为母体材料的屈服极限,为944MPa;
当M1=Mu时,σz为结构的塑性破坏强度σplBCC,即结构的塑性破坏强度σplBCC为:
3实验验证
图7中σGBCC、σ*GBCC和EGBCC分别为BCC变密度结构通过实验求出的屈服强度、抗压强度以及初始刚度。σplGBCC及EplGBCC为BCC变密度结构理论计算的塑形破坏强度和初始刚度,BCC变密度结构通过实验得出的初始刚度为406.2±0.5Mpa,比理论计算出的数值468.47Mpa要低13.3±0.1%;通过实验得出的塑形破坏强度为16.6±0.3Mpa,比理论计算出的数值19.14Mpa低13.3±0.6%,体现了本计算算法的可行性。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.变截面金属点阵结构初始刚度及塑性破坏强度计算方法,其特征在于,包括以下步骤:
s1:通过建立平面坐标系,并结合体心立方点阵结构变截面金属点阵单胞杆的剖面尺寸得到体心立方点阵结构变截面金属点阵单胞杆截面的半径表达式;
s2:通过建立空间坐标系并结合胡克定律及弯压组合变形公式得到所述单胞杆的结点所受的空间切向力、弯矩与空间位移的关系式;
s3:通过能量守恒定律,结合所述空间切向力、弯矩与空间位移的关系式,得到一个关于所述单胞杆的结点处的切向力的一元二次方程,解该方程得到不含所述结点的空间位移的切向力表达式,并结合所述结点处的切向力与轴向力、弯矩的关系得到所述轴向力和弯矩的表达式;
s4:依据胡克定律,结合步骤s3得出的所述切向力、弯矩、轴向力的表达式,得到体心立方点阵结构变截面金属点阵结构的初始刚度及塑性破坏强度;步骤s1中,所述半径表达式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711374124.6A CN108038318B (zh) | 2017-12-19 | 2017-12-19 | 变截面金属点阵结构初始刚度及塑性破坏强度计算算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711374124.6A CN108038318B (zh) | 2017-12-19 | 2017-12-19 | 变截面金属点阵结构初始刚度及塑性破坏强度计算算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108038318A CN108038318A (zh) | 2018-05-15 |
CN108038318B true CN108038318B (zh) | 2020-10-30 |
Family
ID=62099894
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711374124.6A Active CN108038318B (zh) | 2017-12-19 | 2017-12-19 | 变截面金属点阵结构初始刚度及塑性破坏强度计算算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108038318B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108824634B (zh) * | 2018-06-14 | 2020-10-27 | 西北工业大学 | 一种轻质空间点阵结构 |
CN109084169A (zh) * | 2018-09-17 | 2018-12-25 | 中南大学 | 一种点阵结构及点阵零件 |
CN109829261A (zh) * | 2019-04-03 | 2019-05-31 | 大连理工大学 | 一种由变截面圆管组成的轻质点阵结构及其优化设计方法 |
CN110064755B (zh) * | 2019-04-22 | 2021-07-20 | 北京空间飞行器总体设计部 | 3d打印轻量化后埋件以及后埋方法 |
CN111310296B (zh) * | 2019-10-12 | 2024-03-29 | 南京航空航天大学 | 一种应用于复杂曲面的变密度金属板材点阵芯体设计方法 |
CN111895015B (zh) * | 2020-07-03 | 2022-07-08 | 重庆大学 | 基于增材制造的变体梯度点阵结构 |
CN112268092B (zh) * | 2020-09-30 | 2022-05-27 | 重庆大学 | 一种点阵流变智能结构及其制作方法 |
CN112800553B (zh) * | 2021-01-04 | 2024-01-30 | 中国科学院力学研究所 | 多级可控渐进吸能点阵结构 |
CN115249862A (zh) * | 2021-04-26 | 2022-10-28 | 北京航空航天大学 | 用于车辆的动力电池防护系统、动力电池总成和车辆 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104715091A (zh) * | 2013-12-16 | 2015-06-17 | 华中科技大学 | 一种铝合金周期性点阵多孔结构的快速成形制造方法 |
CN105020566A (zh) * | 2015-05-07 | 2015-11-04 | 重庆大学 | 变截面金属点阵结构及其加工方法 |
-
2017
- 2017-12-19 CN CN201711374124.6A patent/CN108038318B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104715091A (zh) * | 2013-12-16 | 2015-06-17 | 华中科技大学 | 一种铝合金周期性点阵多孔结构的快速成形制造方法 |
CN105020566A (zh) * | 2015-05-07 | 2015-11-04 | 重庆大学 | 变截面金属点阵结构及其加工方法 |
Non-Patent Citations (2)
Title |
---|
Mechanical Properties of Ti-6Al-4V Selectively Laser Melted Parts with Body-Centred-Cubic Lattices of Varying cell size;I. Maskery等;《Experimental Mechanics》;20150930;第55卷;第1261-1272页 * |
金属三维点阵结构拓扑构型研究及应用现状综述;易长炎;《功能材料》;20171031;第10055-10065页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108038318A (zh) | 2018-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108038318B (zh) | 变截面金属点阵结构初始刚度及塑性破坏强度计算算法 | |
Ushijima et al. | Prediction of the mechanical properties of micro-lattice structures subjected to multi-axial loading | |
Kubiak | Static and dynamic buckling of thin-walled plate structures | |
Chen et al. | Elastic buckling and static bending of shear deformable functionally graded porous beam | |
Hyun et al. | Simulated properties of Kagomé and tetragonal truss core panels | |
Carrera et al. | Computations and evaluations of higher-order theories for free vibration analysis of beams | |
Zuo et al. | Evolutionary topology optimization of structures with multiple displacement and frequency constraints | |
CN110516317B (zh) | 一种嵌套式类蜂窝夹层结构 | |
Niu et al. | Directional mechanical properties and wave propagation directionality of Kagome honeycomb structures | |
CN110287637A (zh) | 一种弹塑性屈曲承载力计算方法 | |
Hualin et al. | An equivalent continuum method of lattice structures | |
Fazli et al. | Design of tensegrity structures for supporting deployable mesh antennas | |
Veenendaal et al. | Design and optimization of fabric‐formed beams and trusses: evolutionary algorithms and form‐finding | |
Hwang et al. | Investigation for the bending modes of a semi-circular pyramidal kagome sandwich structure and the bending load calculation | |
Zuo et al. | Variable cross-section rectangular beam and sensitivity analysis for lightweight design of bus frame | |
Kang et al. | Dynamical modeling and non-planar coupled behavior of inclined CFRP cables under simultaneous internal and external resonances | |
Galliot et al. | Structural behavior of symmetric spindle-shaped Tensairity girders with reinforced chord coupling | |
Baradaran et al. | Application of an improved genetic algorithm for optimal design of planar steel frames | |
Dahmani et al. | Lateral torsional buckling response of steel beam with different boundary conditions and loading | |
Savio et al. | Design method and taxonomy of optimized regular cellular structures for additive manufacturing technologies | |
Lee et al. | Performance of quadruplex module tensegrities using new pin-jointed connections | |
Pekbey et al. | Buckling optimization of composite columns with variable thickness | |
Šakalys et al. | Numerical investigation on web crippling behaviour of cold-formed C-section beam with vertical stiffeners | |
Shrestha et al. | Crane hook analysis for different cross-section using ANSYS | |
Schmidt et al. | Architectural columns as a result of shape optimization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |