发明名称 一种基于力反馈技术的毛笔建模方法
摘要 一种基于力反馈技术的毛笔建模方法,属于三维虚拟绘制技术领域,是在三维虚拟绘制中引入力反馈技术,毛笔几何模型将毛笔定义为毛笔中心折线以及过各节点上的不同半径的圆,毛笔笔尖为一个点,当毛笔受力发生变形时,各节点上的圆发生变形,形成椭圆;毛笔力学模型为弹簧振子模型:在毛笔笔头安装一个垂直于纸面的弹簧,弹簧始终沿纸面法向变形。通过力反馈设备向毛笔施压弹簧发生变形产生反馈力,由反馈力计算毛笔变形后的中心折线及变形后各椭圆的长短轴,由此确定毛笔笔道宽度,从而建立毛笔压力与毛笔变形的关系。本发明揭示了三维虚拟绘制过程中力对毛笔变形的影响机理,满足了基于力反馈技术的虚拟绘制要求。
申请公布号 CN103345773A 申请公布日期 2013.10.09
申请号 CN201310248722.4 申请日期 2013.06.20
申请人 大连理工大学 发明人 侯增选;郭超;殷婷婷;尹天慈;闫友森;李鑫;郑栓柱;李俊强;陈广州;杨广卿
分类号 G06T15/10(2006.01)I 主分类号 G06T15/10(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 梅洪玉
主权项 1.一种基于力反馈技术的毛笔建模方法,其特征包括以下步骤:(1)建立基于力反馈技术的毛笔几何模型将毛笔定义为毛笔中心折线以及过各节点上的不同半径Ri的圆,毛笔笔尖为一个点,当毛笔受力发生变形时,各节点上的圆发生变形,形成椭圆;假设该圆在虚拟绘制过程中形状不发生变化;设毛笔模型中心折线的长度为L,在中心折线上从毛笔笔尖到笔头根部依次设立控制节点P<sub>0</sub>,P<sub>1</sub>,P<sub>2</sub>,…,P<sub>9</sub>,每相邻两个控制节点之间由线段相连形成折线,由于毛笔尖部相较于毛笔笔头根部更容易产生变形,因此设定毛笔笔尖到笔头根部每相邻两个控制节点之间的线段长度逐渐增大,呈等差数列分布,公差<img file="FDA00003377333300011.GIF" wi="447" he="131" />其中P<sub>0</sub>P<sub>1</sub>是控制节点P<sub>0</sub>和P<sub>1</sub>之间的线段长度,n为中心折线上线段的个数,n=9;(2)建立毛笔力学模型,毛笔力学模型为弹簧振子模型在毛笔笔头安装一个垂直于纸面的弹簧,弹簧始终沿纸面法向变形,在控制节点P<sub>9</sub>和P<sub>9</sub>在纸面上的投影点<img file="FDA00003377333300012.GIF" wi="70" he="84" />处设置虚拟弹簧,毛笔未受力作用时,<img file="FDA00003377333300013.GIF" wi="83" he="82" />与P<sub>0</sub>重合;当毛笔受力作用时,弹簧受力向下运动,压力F的计算公式为:<maths num="0001"><![CDATA[<math><mrow><mi>F</mi><mo>=</mo><mfrac><mrow><mi>&lambda;</mi><mo>&CenterDot;</mo><mi>H</mi></mrow><mrow><mn>1</mn><mo>+</mo><mi>W</mi></mrow></mfrac><mo>&times;</mo><mi>X</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中λ为压力因子,单位为N/mm;X为毛笔受压接触纸面后向下移动的距离,即弹簧变形量,单位为mm;H为毛笔硬度,与F成正比,变化范围为:0&lt;H≤1;W为毛笔湿度,即毛笔内部剩余水墨量与毛笔所蘸含的总水墨量的比值,变化范围为:0≤W≤1,F随着湿度的减小而增大;(3)确定力与毛笔变形的关系压力使毛笔产生弯曲变形,S<sub>i</sub>为毛笔笔头中心线P<sub>0</sub>到P<sub>i</sub>线段在平面<img file="FDA00003377333300015.GIF" wi="122" he="82" />上的投影,压力越大,S<sub>i</sub>越大,当施加于笔上的压力为0时,笔尖P<sub>0</sub>与纸面刚刚接触,S<sub>i</sub>为0;当压力最大时,毛笔笔头的中心线全部投影在纸面上;设p为压力影响因子,p值为压力与最大输出力的比值,将其转换为角度<img file="FDA00003377333300021.GIF" wi="396" he="105" />力反馈设备能够提供的最大输出力为7.9牛顿,因此p=F/7.9,变化范围为:0≤p≤1;S<sub>i</sub>=(P<sub>0</sub>P<sub>1</sub>+P<sub>1</sub>P<sub>2</sub>+...+P<sub>i-1</sub>P<sub>i</sub>)×cosα         (2)其中P<sub>i-1</sub>P<sub>i</sub>是毛笔中心线控制节点P<sub>i-1</sub>到P<sub>i</sub>之间的线段长度;设毛笔沿运动方向移动微小距离,由x轴正方向逆时针旋转到运动方向的反方向时,所经过的角度为γ角,其取值范围为0≤γ<2π,建立三维局部坐标系,则毛笔中心线上第i个控制节点P<sub>i</sub>的x坐标值和y坐标值计算公式如下:x<sub>i</sub>=(S<sub>9</sub>-S<sub>i</sub>)×cosγ                        (3)y<sub>i</sub>=(S<sub>9</sub>-S<sub>i</sub>)×sinγ                                (4)各个控制节点的z坐标值计算公式如下:z<sub>0</sub>=0<maths num="0002"><![CDATA[<math><mrow><msub><mi>Z</mi><mn>1</mn></msub><mo>=</mo><msub><mi>Z</mi><mn>0</mn></msub><mo>+</mo><msup><mrow><mo>[</mo><msub><mi>P</mi><mn>0</mn></msub><msup><msub><mi>P</mi><mn>1</mn></msub><mn>2</mn></msup><mo>-</mo><msubsup><mi>S</mi><mn>1</mn><mn>2</mn></msubsup><mo>]</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>Z</mi><mn>2</mn></msub><mo>=</mo><msub><mi>Z</mi><mn>1</mn></msub><mo>+</mo><msup><mrow><mo>[</mo><msub><mi>P</mi><mn>1</mn></msub><msup><msub><mi>P</mi><mn>2</mn></msub><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>S</mi><mn>2</mn></msub><mo>-</mo><msub><mi>S</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup></mrow></math>]]></maths>…<maths num="0004"><![CDATA[<math><mrow><msub><mi>Z</mi><mi>i</mi></msub><mo>=</mo><msub><mi>Z</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msup><mrow><mo>[</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><msup><msub><mi>P</mi><mi>i</mi></msub><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>-</mo><msub><mi>S</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup></mrow></math>]]></maths>由此推导出P<sub>i</sub>的z坐标值计算如公式(5)所示:<maths num="0005"><![CDATA[<math><mrow><msub><mi>Z</mi><mi>i</mi></msub><mo>=</mo><msub><mi>Z</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msup><mrow><mo>[</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><msup><msub><mi>P</mi><mi>i</mi></msub><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>S</mi><mi>i</mi></msub><mo>-</mo><msub><mi>S</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>]</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>当毛笔受到压力产生弯曲变形时,椭圆长半轴R<sub>ib</sub>随着力的增加而变大。当压力F为零时,R<sub>ib</sub>=R<sub>i</sub>;当施加于毛笔的压力F最大时,椭圆长半轴R<sub>ib</sub>最大,取此时的椭圆长半轴R<sub>ib</sub>=1.5R<sub>i</sub>,R<sub>ib</sub>的计算公式为:R<sub>ib</sub>=R<sub>i</sub>×(1+p/2)               (6)假定椭圆的面积不变,则椭圆短半轴R<sub>ia</sub>的计算公式为:<maths num="0006"><![CDATA[<math><mrow><msub><mi>R</mi><mi>ia</mi></msub><mo>=</mo><msubsup><mi>R</mi><mi>i</mi><mn>2</mn></msubsup><mo>/</mo><msub><mi>R</mi><mi>ib</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>再根据椭圆长短轴确定毛笔笔道宽度,从而建立毛笔压力与毛笔变形的关系;在虚拟绘制过程中,若要求笔道较宽,加大毛笔的压力,若要求笔道较窄,减小毛笔压力。
地址 116024 辽宁省大连市凌工路2号