发明名称 一种求解由地基变形引起的机床加工误差的方法
摘要 一种求解由地基变形引起的机床加工误差的方法,属于机床研究技术领域。其针对地基不均匀沉降对机床加工精度的影响进行研究,以大型龙门式机床为研究对象,假设机床底座由于地基不均匀沉降导致两底座产生沉降差,由底座沉降差引起机床整体的结构变形使刀尖点产生位置和角度的偏移,从而影响加工精度,对机床和地基整体模型进行简化,提供了底座的沉降差对机床加工精度影响的建模方法,可为重型龙门式机床的设计、结构优化及误差补偿提供理论与技术指导。
申请公布号 CN104122839A 申请公布日期 2014.10.29
申请号 CN201410332520.2 申请日期 2014.07.11
申请人 北京工业大学 发明人 刘志峰;张柯;王建华;杨恩芝;宁越;田杨
分类号 G05B19/404(2006.01)I 主分类号 G05B19/404(2006.01)I
代理机构 北京思海天达知识产权代理有限公司 11203 代理人 纪佳
主权项 一种求解由地基变形引起的机床加工误差的方法,,步骤包括:简化机床模型、简化地基变形模型、应用材料力学进行受力分析、刀尖点横向偏移、刀尖点纵向偏移、刀尖点转角偏差;其特征在于:具体步骤如下:S1简化机床模型把重型龙门机床结构简化为框架结构,提取重型龙门机床尺寸、材料参数和截面特性,获得机床材料的弹性模量E,横梁长度l,立柱高度h和刀杆的长度l<sub>0</sub>,并由截面特性分别算得横梁和立柱的惯性矩I<sub>l</sub>和I<sub>h</sub>。S2地基变形模型把地基的变形简化为不同高度的两个平台,机床的两根导轨分别坐落在地基不同高度平台上,两个平台之间的沉降差为Δh。S3应用材料力学进行受力分析S3.1假设立柱为刚体求横梁变形假设立柱为刚性不变形,当地基变形时产生沉降差,使机床横梁弯曲,根据材料力学可求得横梁一端B处的挠度ω<sub>l</sub>和转角θ<sub>l</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>l</mi></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mfrac><mi>F</mi><mn>6</mn></mfrac><msup><mi>x</mi><mn>3</mn></msup><mo>-</mo><mfrac><mi>Fl</mi><mn>2</mn></mfrac><msup><mi>x</mi><mn>2</mn></msup></mrow><msub><mi>EI</mi><mi>l</mi></msub></mfrac></mrow>]]></math><img file="FDA0000536673520000011.GIF" wi="467" he="225" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mi>l</mi></msub><mo>=</mo><msubsup><mi>&omega;</mi><mi>l</mi><mo>&prime;</mo></msubsup><mo>=</mo><mo>-</mo><mfrac><mrow><mfrac><mi>F</mi><mn>2</mn></mfrac><msup><mi>x</mi><mn>2</mn></msup><mo>-</mo><mi>Flx</mi></mrow><msub><mi>EI</mi><mi>l</mi></msub></mfrac></mrow>]]></math><img file="FDA0000536673520000012.GIF" wi="521" he="222" /></maths>式中:ω<sub>l</sub>为横梁一端处的挠度,F为使横梁挠曲的力,l为机床横梁长度,x为刀杆所在横梁位置,E为机床材料的弹性模量,I<sub>l</sub>为横梁惯性矩,θ<sub>l</sub>为横梁弯曲的转角,ω′<sub>l</sub>为横梁挠度的导数;以及地基沉降差为y<sub>1</sub>时使横梁变形所需要的力F的值:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mo>-</mo><mfrac><mrow><mn>3</mn><mi>E</mi><msub><mi>I</mi><mi>l</mi></msub><msub><mi>y</mi><mn>1</mn></msub></mrow><msup><mi>l</mi><mn>3</mn></msup></mfrac></mrow>]]></math><img file="FDA0000536673520000013.GIF" wi="330" he="148" /></maths>式中:E为机床材料的弹性模量,I<sub>l</sub>为横梁惯性矩,y<sub>1</sub>为地基沉降差,为刀杆长度;S3.2假设横梁为刚性求立柱的变形假设横梁为刚性不变形,当地基变形时产生沉降差,使机床的立柱发生弯曲,根据材料力学可求得横梁所在立柱位置A点处即机床左立柱与横梁的交点的挠度ω<sub>h</sub>和转角θ<sub>h</sub>:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>h</mi></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mi>Fl</mi><msup><mi>h</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><mi>E</mi><msub><mi>I</mi><mi>h</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000536673520000021.GIF" wi="319" he="177" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mi>h</mi></msub><mo>=</mo><msubsup><mi>&omega;</mi><mi>h</mi><mo>&prime;</mo></msubsup><mo>=</mo><mo>-</mo><mfrac><mi>Flh</mi><msub><mi>EI</mi><mi>h</mi></msub></mfrac></mrow>]]></math><img file="FDA0000536673520000022.GIF" wi="380" he="159" /></maths>式中:ω<sub>h</sub>为立柱一端处的挠度,F为使横梁挠曲的力,l为机床横梁长度,h为横梁所在立柱位置,E为机床材料的弹性模量,I<sub>h</sub>为立柱惯性矩,θ<sub>h</sub>为立柱弯曲的转角,ω′<sub>h</sub>为立柱挠度的导数。和立柱A点出现的转角引起B点处即机床右立柱与横梁的交点的下沉量y<sub>2</sub>:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><msub><mi>&theta;</mi><mi>h</mi></msub><mi>l</mi><mo>=</mo><mo>-</mo><mfrac><mrow><mi>F</mi><msup><mi>l</mi><mn>2</mn></msup><mi>h</mi></mrow><msub><mi>EI</mi><mi>h</mi></msub></mfrac></mrow>]]></math><img file="FDA0000536673520000023.GIF" wi="412" he="167" /></maths>式中:y<sub>2</sub>为立柱的下沉量,θ<sub>h</sub>为立柱弯曲的转角,F为使横梁挠曲的力,l为机床横梁长度,h为横梁所在立柱位置,E为机床材料的弹性模量,I<sub>h</sub>为立柱惯性矩。S3.3横梁和立柱的叠加变形由于立柱的弯曲,A点出现转角,导致横梁跟随转动而使横梁B点出现y方向的位移,当假设立柱为刚体时,横梁一端收到力F的左右也会出现弯曲变形,因此横梁B点的下沉量Δh为立柱和横梁变形的叠加:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mo>-</mo><mi>&Delta;h</mi><mo>=</mo><msub><mi>y</mi><mn>1</mn></msub><mo>+</mo><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><mo>-</mo><mfrac><msup><mi>Fl</mi><mn>3</mn></msup><mrow><mn>3</mn><msub><mi>EI</mi><mi>l</mi></msub></mrow></mfrac><mo>-</mo><mfrac><mrow><msup><mi>Fl</mi><mn>2</mn></msup><mi>h</mi></mrow><msub><mi>EI</mi><mi>h</mi></msub></mfrac></mrow>]]></math><img file="FDA0000536673520000024.GIF" wi="770" he="179" /></maths>式中:y<sub>1</sub>、y<sub>2</sub>分别为横梁和立柱弯曲时地基的下沉量,F为使横梁挠曲的力,l为机床横梁长度,h为横梁所在立柱位置,E为机床材料的弹性模量,I<sub>l</sub>为横梁惯性矩,I<sub>h</sub>为立柱惯性矩。然后可用地基沉降差表示出横梁端点B处所受力F的大小,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mfrac><mi>&Delta;hE</mi><mrow><mfrac><msup><mi>l</mi><mn>3</mn></msup><mrow><mn>3</mn><msub><mi>I</mi><mi>l</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mi>l</mi><mn>2</mn></msup><mi>h</mi></mrow><msub><mi>I</mi><mi>h</mi></msub></mfrac></mrow></mfrac></mrow>]]></math><img file="FDA0000536673520000031.GIF" wi="358" he="252" /></maths>式中:Δh为横梁和立柱同时弯曲时地基的下沉量,E为机床材料的弹性模量,l为机床横梁长度,h为横梁所在立柱位置,I<sub>l</sub>为横梁惯性矩,I<sub>h</sub>为立柱惯性矩。并求出横梁不同位置处x、y方向的偏移量以及转角公式:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><msub><mi>&omega;</mi><mi>h</mi></msub><mo>=</mo><mfrac><mrow><mi>Fl</mi><msup><mi>h</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><msub><mi>EI</mi><mi>h</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>&Delta;hl</mi><msup><mi>h</mi><mn>2</mn></msup></mrow><mrow><msub><mrow><mn>2</mn><mi>I</mi></mrow><mi>h</mi></msub><mrow><mo>(</mo><mfrac><msup><mi>l</mi><mn>3</mn></msup><mrow><mn>3</mn><msub><mi>I</mi><mi>l</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mi>l</mi><mn>2</mn></msup><mi>h</mi></mrow><msub><mi>I</mi><mi>h</mi></msub></mfrac><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000536673520000032.GIF" wi="773" he="244" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><msub><mi>&theta;</mi><mi>h</mi></msub><mi>x</mi><mo>+</mo><msub><mi>&omega;</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mfrac><mi>F</mi><mn>6</mn></mfrac><msup><mi>x</mi><mn>3</mn></msup><mo>-</mo><mfrac><mi>Fl</mi><mn>2</mn></mfrac><msup><mi>x</mi><mn>2</mn></msup></mrow><msub><mi>EI</mi><mi>l</mi></msub></mfrac></mrow>]]></math><img file="FDA0000536673520000033.GIF" wi="700" he="222" /></maths><maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mn>1</mn></msub><mo>=</mo><msub><mi>&theta;</mi><mi>l</mi></msub><mo>+</mo><msub><mi>&theta;</mi><mi>h</mi></msub><mo>=</mo><mo>-</mo><mfrac><mrow><mfrac><mi>F</mi><mn>2</mn></mfrac><msup><mi>x</mi><mn>2</mn></msup><mo>-</mo><mi>Flx</mi></mrow><msub><mi>EI</mi><mi>l</mi></msub></mfrac><mo>-</mo><mfrac><mi>Flh</mi><msub><mi>EI</mi><mi>h</mi></msub></mfrac></mrow>]]></math><img file="FDA0000536673520000034.GIF" wi="775" he="226" /></maths>式中:ω<sub>h</sub>为立柱一端处的挠度,ω<sub>l</sub>为横梁一端处的挠度,F为使横梁挠曲的力,θ<sub>h</sub>为立柱弯曲的转角,θ<sub>l</sub>为横梁弯曲的转角,x为刀杆所在横梁位置,E为机床材料的弹性模量,l为机床横梁长度,h为横梁所在立柱位置,I<sub>l</sub>为横梁惯性矩,I<sub>h</sub>为立柱惯性矩,Δh为横梁和立柱同时弯曲时地基的下沉量。S3.4考虑刀杆尺寸的影响刀杆在横梁上滑动,求出刀尖点相对于刀杆与横梁焦点处的偏移,叠加上一步求得的横梁不同位置处x方向的偏移量算出刀尖点的偏移量:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msub><mi>x</mi><mn>2</mn></msub><mo>=</mo><mi>&theta;</mi><msub><mi>l</mi><mn>0</mn></msub><mo>=</mo><mo>-</mo><mrow><mo>(</mo><mfrac><mrow><mfrac><mi>F</mi><mn>2</mn></mfrac><msup><mi>x</mi><mn>2</mn></msup><mo>-</mo><mi>Flx</mi></mrow><msub><mi>EI</mi><mi>l</mi></msub></mfrac><mo>+</mo><mfrac><mi>Flh</mi><msub><mi>EI</mi><mi>h</mi></msub></mfrac><mo>)</mo></mrow><msub><mi>l</mi><mn>0</mn></msub></mrow>]]></math><img file="FDA0000536673520000035.GIF" wi="684" he="258" /></maths>式中:θ<sub>1</sub>为横梁弯曲的转角,l<sub>0</sub>为刀杆的长度,F为使横梁挠曲的力,x为刀杆所在横梁位置,E为机床材料的弹性模量,l为机床横梁长度,h为横梁所在立柱位置,I<sub>h</sub>为立柱惯性矩。S4刀尖点横向偏移根据横梁和立柱的叠加变形以及刀杆尺寸的影响,最后得出刀尖点横向偏移量公式:<maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><msub><mi>x</mi><mn>0</mn></msub><mo>=</mo><mfrac><mrow><mi>&Delta;hl</mi><msup><mi>h</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><msub><mi>I</mi><mi>h</mi></msub><mrow><mo>(</mo><mfrac><msup><mi>l</mi><mn>3</mn></msup><mrow><mn>3</mn><msub><mi>I</mi><mi>l</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mi>l</mi><mn>2</mn></msup><mi>h</mi></mrow><msub><mi>I</mi><mi>h</mi></msub></mfrac><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><mi>&Delta;h</mi><mrow><mfrac><msup><mi>l</mi><mn>3</mn></msup><mrow><mn>3</mn><msub><mi>I</mi><mi>l</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mi>l</mi><mn>2</mn></msup><mi>h</mi></mrow><msub><mi>I</mi><mi>h</mi></msub></mfrac></mrow></mfrac><mo>&CenterDot;</mo><mrow><mo>(</mo><mfrac><mrow><mfrac><msup><mi>x</mi><mn>2</mn></msup><mn>2</mn></mfrac><mo>-</mo><mi>lx</mi></mrow><msub><mi>I</mi><mi>l</mi></msub></mfrac><mo>+</mo><mfrac><mi>lh</mi><msub><mi>I</mi><mi>h</mi></msub></mfrac><mo>)</mo></mrow><msub><mi>l</mi><mn>0</mn></msub></mrow>]]></math><img file="FDA0000536673520000041.GIF" wi="1199" he="344" /></maths>式中:Δh为地基沉降差,l为机床横梁长度,h为机床横梁所在立柱高度,I<sub>h</sub>为立柱的惯性矩,I<sub>l</sub>为横梁惯性矩,x为刀杆所在横梁位置,l<sub>0</sub>为刀杆长度。S5刀尖点纵向偏移根据横梁和立柱的叠加变形以及刀杆尺寸的影响,最后得出刀尖点纵向偏移量公式:<maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><mi>y</mi><mo>=</mo><mo>-</mo><mfrac><mi>&Delta;h</mi><mrow><mfrac><msup><mi>l</mi><mn>3</mn></msup><mrow><mn>3</mn><msub><mi>I</mi><mi>l</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mi>l</mi><mn>2</mn></msup><mi>h</mi></mrow><msub><mi>I</mi><mi>h</mi></msub></mfrac></mrow></mfrac><mo>&CenterDot;</mo><mfrac><mrow><mfrac><mn>1</mn><mn>6</mn></mfrac><msup><mi>x</mi><mn>3</mn></msup><mo>-</mo><mfrac><mi>l</mi><mn>2</mn></mfrac><msup><mi>x</mi><mn>2</mn></msup></mrow><msub><mi>I</mi><mi>l</mi></msub></mfrac></mrow>]]></math><img file="FDA0000536673520000042.GIF" wi="731" he="340" /></maths>式中:Δh为地基沉降差,l为机床横梁长度,h为机床横梁所在立柱高度,I<sub>h</sub>为立柱的惯性矩,I<sub>l</sub>为横梁惯性矩,x为刀杆所在横梁位置。S6刀尖点转角偏差根据横梁和立柱的叠加变形以及刀杆尺寸的影响,最后得出刀尖点转角偏差公式:<maths num="0015" id="cmaths0015"><math><![CDATA[<mrow><mi>&theta;</mi><mo>=</mo><mo>-</mo><mfrac><mi>&Delta;h</mi><mrow><mfrac><msup><mi>l</mi><mn>3</mn></msup><mrow><mn>3</mn><msub><mi>I</mi><mi>l</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mi>l</mi><mn>2</mn></msup><mi>h</mi></mrow><msub><mi>I</mi><mi>h</mi></msub></mfrac></mrow></mfrac><mo>&CenterDot;</mo><mrow><mo>(</mo><mfrac><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>x</mi><mn>2</mn></msup><mo>-</mo><mi>lx</mi></mrow><msub><mi>I</mi><mi>l</mi></msub></mfrac><mo>+</mo><mfrac><mi>lh</mi><msub><mi>I</mi><mi>h</mi></msub></mfrac><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000536673520000043.GIF" wi="814" he="338" /></maths>式中:Δh为地基沉降差,l为机床横梁长度,h为机床横梁所在立柱高度,I<sub>h</sub>为立柱的惯性矩,I<sub>l</sub>为横梁惯性矩,x为刀杆所在横梁位置。
地址 100124 北京市朝阳区平乐园100号