发明名称 一种三轴向编织复合材料格栅力学模量的预测方法
摘要 一种三轴向编织复合材料格栅力学模量的预测方法,它有四大步骤:一、根据纤维束的三轴编织方式,选择最小的重复性单元作为代表性体积元,由此确定其胞体单元;二、根据外载荷施加方式以及步骤一中的胞体单元,对胞体单元内的纤维束进行受力分析,从而,建立三轴编织复合材料格栅胞体单元纤维织布的细观力学模型,确定胞体单元内纤维织布的总应变余能U<sup>*</sup>,并利用最小势能原理,求解胞体单元纤维织布的内力;三、根据卡式定理或单位载荷法求解胞体单元纤维织布的变形,再根据应力‑应变本构方程,得到三轴编织复合材料格栅织布的力学模量;四、根据混合定理,得到三轴编织复合材料格栅的弹性模量。
申请公布号 CN103455712B 申请公布日期 2016.08.10
申请号 CN201310363477.1 申请日期 2013.08.20
申请人 北京航空航天大学 发明人 熊峻江;白江波;满孜郁;刘猛
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种三轴向编织复合材料格栅力学模量的预测方法,其特征在于:该方法具体步骤如下:步骤一、根据纤维束的三轴编织方式,选择最小的重复性单元作为代表性体积元,由此确定其胞体单元;将三轴向纤维束理想化成正弦的曲梁,并根据纤维束中心线位置,分别建立0°、‑60°和60°纤维束坐标系,于是,得到三轴向纤维束的中心线Z坐标表达式:<maths num="0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>Z</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mi>&pi;</mi><mi>h</mi></mrow><mrow><mn>2</mn><mi>L</mi></mrow></mfrac><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mfrac><mi>&pi;</mi><mi>L</mi></mfrac><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>Z</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mi>&pi;</mi><mi>h</mi></mrow><mrow><mn>2</mn><mi>L</mi></mrow></mfrac><mi>c</mi><mi>o</mi><mi>s</mi><mrow><mo>(</mo><mfrac><mi>&pi;</mi><mi>L</mi></mfrac><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>Z</mi><mn>3</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mi>&pi;</mi><mi>h</mi></mrow><mrow><mn>2</mn><mi>L</mi></mrow></mfrac><mi>s</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mfrac><mi>&pi;</mi><mi>L</mi></mfrac><msub><mi>x</mi><mn>3</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000985639380000011.GIF" wi="1206" he="407" /></maths>式中h和L分别为纤维束截面高度和纤维束长度;步骤二、根据外载荷施加方式以及步骤一中的胞体单元,对胞体单元内的纤维束进行受力分析,从而,建立三轴编织复合材料格栅胞体单元纤维织布的细观力学模型,确定胞体单元内纤维织布的总应变余能U<sup>*</sup>,并利用最小势能原理,求解胞体单元纤维织布的内力;胞体单元的总应变余能U<sup>*</sup>的表达式为<maths num="0002"><math><![CDATA[<mrow><msup><mi>U</mi><mo>*</mo></msup><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msubsup><mi>U</mi><mi>i</mi><mo>*</mo></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000985639380000012.GIF" wi="1107" he="127" /></maths>式中,<maths num="0003"><math><![CDATA[<mrow><msubsup><mi>U</mi><mi>i</mi><mo>*</mo></msubsup><mo>=</mo><msup><mfrac><mn>1</mn><mrow><mi>E</mi><mi>I</mi></mrow></mfrac><mfrac><mi>L</mi><munderover><mo>&Integral;</mo><mn>0</mn><mn>2</mn></munderover></mfrac></msup><msup><mi>M</mi><mn>2</mn></msup><mi>d</mi><mi>x</mi><mo>+</mo><msup><mfrac><mn>1</mn><mrow><mi>E</mi><mi>A</mi></mrow></mfrac><mfrac><mi>L</mi><munderover><mo>&Integral;</mo><mn>0</mn><mn>2</mn></munderover></mfrac></msup><msup><mi>N</mi><mn>2</mn></msup><mi>d</mi><mi>x</mi><mo>+</mo><msup><mfrac><mn>1</mn><mrow><msub><mi>GI</mi><mi>p</mi></msub></mrow></mfrac><mfrac><mi>L</mi><munderover><mo>&Integral;</mo><mn>0</mn><mn>2</mn></munderover></mfrac></msup><msup><mi>T</mi><mn>2</mn></msup><mi>d</mi><mi>x</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000985639380000013.GIF" wi="1410" he="181" /></maths>其中,M、N和T分别为弯矩、轴力和扭矩;I、A和I<sub>p</sub>分别为纤维束截面惯性矩、面积和极惯性矩;E为纤维束的弹性模量;根据最小势能原理,确定胞体单元各个内力;步骤三、根据卡式定理或单位载荷法求解胞体单元纤维织布的变形,再根据应力‑应变本构方程,得到三轴编织复合材料格栅织布的力学模量;由卡式定理确定的系统在外载荷包括拉伸、压缩和剪切载荷作用下三轴编织复合材料织布的变形Δ为<maths num="0004"><math><![CDATA[<mrow><mi>&Delta;</mi><mo>=</mo><mfrac><mrow><mo>&part;</mo><msup><mi>U</mi><mo>*</mo></msup></mrow><mrow><mo>&part;</mo><mi>P</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000985639380000021.GIF" wi="908" he="126" /></maths>其中P代表所受的外载荷;根据应力和应变关系表达式,得到纤维织布的拉伸、压缩和剪切模量公式:<maths num="0005"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>E</mi><mi>t</mi></msub><mo>=</mo><mfrac><msub><mi>&sigma;</mi><mrow><mi>t</mi><mi>x</mi></mrow></msub><msub><mi>&epsiv;</mi><mrow><mi>t</mi><mi>x</mi></mrow></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>E</mi><mi>c</mi></msub><mo>=</mo><mfrac><msub><mi>&sigma;</mi><mrow><mi>c</mi><mi>x</mi></mrow></msub><msub><mi>&epsiv;</mi><mrow><mi>c</mi><mi>x</mi></mrow></msub></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mi>G</mi><mo>=</mo><mfrac><mi>&tau;</mi><mi>&gamma;</mi></mfrac></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000985639380000022.GIF" wi="1109" he="400" /></maths>式中E<sub>t</sub>、E<sub>c</sub>和G分别表示纤维织布的拉伸、压缩和剪切模量,σ<sub>tx</sub>、σ<sub>cx</sub>和τ分别表示拉伸、压缩和剪切应力,ε<sub>tx</sub>、ε<sub>cx</sub>和γ分别表示拉伸、压缩和剪切应变;步骤四、根据混合定理,得到三轴编织复合材料格栅的弹性模量三轴编织复合材料格栅的弹性模量为<maths num="0006"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>E</mi><mrow><mi>t</mi><mi>l</mi><mi>a</mi></mrow></msub><mo>=</mo><msub><mi>E</mi><mi>t</mi></msub><msub><mi>V</mi><mi>f</mi></msub><mo>+</mo><msub><mi>E</mi><mi>m</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>V</mi><mi>f</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>E</mi><mrow><mi>c</mi><mi>l</mi><mi>a</mi></mrow></msub><mo>=</mo><msub><mi>E</mi><mi>c</mi></msub><msub><mi>V</mi><mi>f</mi></msub><mo>+</mo><msub><mi>E</mi><mi>m</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>V</mi><mi>f</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>G</mi><mrow><mi>l</mi><mi>a</mi></mrow></msub><mo>=</mo><msub><mi>GV</mi><mi>f</mi></msub><mo>+</mo><msub><mi>G</mi><mi>m</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>V</mi><mi>f</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>.</mo></mrow>]]></math><img file="FDA0000985639380000023.GIF" wi="547" he="236" /></maths>
地址 100191 北京市海淀区学院路37号