本发明涉及新型变截面金属点阵结构初始刚度及塑性破坏强度计算算法，具体地，涉及一种通过变截面金属点阵结构的剖面尺寸、结点空间位移及受力，结合材料力学胡克定律、功能原理、平面几何等最终计算出变截面金属点阵结构的初始刚度及塑性破坏强度。
背景技术：
目前，金属点阵结构的研究对象以体心立方(body-centeredcubic；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中，所述空间切向力、弯矩与空间位移的关系式为：
式中，γ为单胞杆在切向力及弯矩的作用下切向位移；f1为空间切向力，m1为弯矩，b2(α)＝sinα/[rsinα-(r+r1)]4；
进一步的，在步骤s3中，所述切向力的一元二次方程为：
其中，
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的作用下，变截面杆的轴向位移为：
本文中设bi＝bi(α)为关于α的函数，ci为常量；式(11)中e为母体材料的弹性模量；b1(α)＝sinα/[rsinα-(r+r1)]2；
结合式(5)和式(6)可得杆的轴向力n1为：
通过s7点在空间坐标系o'x'y'z'中三个方向的位移(u，ν，w)可得杆切向位移为：
在平面坐标系o”x”y”下，由于杆不仅受到切向力f1的作用，还受到弯矩m1的作用，因此根据结构弯压组合变形的方程可得杆在切向力及弯矩的作用下切向位移为：
式中：b2(α)＝sinα/[rsinα-(r+r1)]4；
另外，因为对称的关系，各结点处的弯曲力矩相同，所以结点s7、s9处的弯矩ms7、ms9相等且为m1，且
结合式(8)、(9)、(10)可得杆件所受的切向力f1及弯矩m1分别为：
其中，
通过s7点在空间坐标系o'x'y'z'中所受的切向力f1及弯矩m1可得s7在z'向所受的力fz'为：
fz'＝f1cosθ-n1sinθ(12)
式中：
另外，由于单胞底面四个杆件的端点s5、s6、s7、s8同等的受到压应力σz的作用，因此各端点z'向受力fz'均为σzl2/4。
由前述分析可知，杆的变形主要为轴向变形与横向弯曲变形，其弹性应变能包括轴向伸缩能与弯曲应变能，由功能原理与胡克定律可知，s7s9杆受到轴向力n1、切向力f1及弯矩m1作用时的应变能u为：
式中：b4(α)＝2rcosα+l；b5(α)＝b42(α)+l2；
一个单胞里面有八根杆件，所以单胞的总应变能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的一元二次方程为：
c9f12+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实验验证
图中σgbcc、σ*gbcc和egbcc分别为bcc变密度结构通过实验求出的屈服强度、抗压强度以及初始刚度。σpigbcc及epigbcc为bcc变密度结构理论计算的塑形破坏强度和初始刚度，bcc变密度结构通过实验得出的初始刚度为406.2±0.5mpa，比理论计算出的数值468.47mpa要低13.3±0.1％；通过实验得出的塑形破坏强度为16.6±0.3mpa，比理论计算出的数值19.14mpa低13.3±0.6％，体现了本计算算法的可行性。
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。