CN103345773B - 一种基于力反馈技术的毛笔建模方法 - Google Patents
一种基于力反馈技术的毛笔建模方法 Download PDFInfo
- Publication number
- CN103345773B CN103345773B CN201310248722.4A CN201310248722A CN103345773B CN 103345773 B CN103345773 B CN 103345773B CN 201310248722 A CN201310248722 A CN 201310248722A CN 103345773 B CN103345773 B CN 103345773B
- Authority
- CN
- China
- Prior art keywords
- writing brush
- brush
- writing
- pressure
- spring
- 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.)
- Expired - Fee Related
Links
Abstract
一种基于力反馈技术的毛笔建模方法,属于三维虚拟绘制技术领域,是在三维虚拟绘制中引入力反馈技术,毛笔几何模型将毛笔定义为毛笔中心折线以及过各节点上的不同半径的圆,毛笔笔尖为一个点,当毛笔受力发生变形时,各节点上的圆发生变形,形成椭圆;毛笔力学模型为弹簧振子模型:在毛笔笔头安装一个垂直于纸面的弹簧,弹簧始终沿纸面法向变形。通过力反馈设备向毛笔施压弹簧发生变形产生反馈力,由反馈力计算毛笔变形后的中心折线及变形后各椭圆的长短轴,由此确定毛笔笔道宽度,从而建立毛笔压力与毛笔变形的关系。本发明揭示了三维虚拟绘制过程中力对毛笔变形的影响机理,满足了基于力反馈技术的虚拟绘制要求。
Description
技术领域
本发明涉及三维虚拟绘制技术领域,特别是涉及三维虚拟绘制中的一种基于力反馈技术的毛笔建模方法。
背景技术
毛笔建模是三维虚拟绘制的基础,在国外,Strassman在文献“S.Strassmann.HairyBrushes[C].SIGGRAPH’86Proceedings.USA:ACMPress,1986:225-232.”中提出了一个虚拟画笔模型,通过鼠标确定控制点,插值形成笔迹路径,键盘输入压力值控制笔迹宽度,但未引入力反馈设备。2001年BillBaxter等在文献“BillBaxter,VincentScheib,MingC.Lin,DineshManocha.DAB:InteractiveHapticPaintingwith3DVirtualBrushes[C]//ProceedingsoftheACMSIGGRAPHConferenceonComputerGraphics,2001:461-468”中采用一种弹簧振子的半隐式积分算法和细分表面进行毛笔建模,应用该模型不能准确的建立画笔变形与力的关系,参见文献WilliamV.BaxterandMingC.Lin.AVersatileInteractive3DBrushModel[C]//Proceedingsofthe12thPacificConferenceonComputerGraphicsandApplications,Seoul,SouthKorea,2004:319-328。2002年Jeng-shengYeh等在文献“Jeng-shengYeh,Ting-yuLien,MingOuhyoung.OntheEffectsofHapticDisplayinBrushandInkSimulationforChinesePaintingandCalligraphy[C]//Proceedingsofthe10thPacificConferenceonComputerGraphicsandApplication.2002:439-441”中改进了BillBaxter的画笔建模方法,采用可弯曲的弹簧构建毛笔骨架,建立了一个八方向笔毛的毛笔模型,但论文未给出力和毛笔变形的关系。2003年宓晓峰等在文献“宓晓峰,唐敏,林建贞,董金祥.基于经验的虚拟毛笔模型[J].计算机研究与发展,2003,40(8):1244-1251”中提出了一种基于经验的虚拟毛笔模型,但未考虑力和毛笔变形的关系。2008年孙美君在文献“孙美君,孙济洲,王征,丁兆伟.3D实用毛笔模型的物理仿真[J].天津大学学报,2008,41(3):293-299”中构建了一个实用的3D毛笔模型,提出势能-阀值理论控制毛笔变形。
综上所述,现有的毛笔建模方法不能满足基于力反馈技术的三维虚拟绘制要求。
发明内容
本发明提出了一种基于力反馈技术的毛笔建模方法,毛笔模型分为毛笔几何模型和毛笔力学模型。毛笔几何模型将毛笔定义为毛笔中心折线以及过各节点上的不同半径的圆,毛笔笔尖为一个点,当毛笔受力发生变形时,各节点上的圆发生变形,形成椭圆;毛笔力学模型为弹簧振子模型:在毛笔笔头安装一个垂直于纸面的弹簧,弹簧始终沿纸面法向变形。通过力反馈设备向毛笔施压,弹簧发生变形,从而产生反馈力,由反馈力计算毛笔变形后的中心折线及变形后各椭圆的长短轴,由此确定毛笔笔道宽度,从而建立毛笔压力与毛笔变形的关系。在虚拟绘制过程中,若要求笔道较宽,可以加大毛笔的压力,若笔道要求较窄,可减小毛笔压力。该毛笔模型可用于三维虚拟绘制过程中。
一种基于力反馈技术的毛笔建模方法,包括以下部分:
(1)建立基于力反馈技术的毛笔几何模型
如附图1所示,首先将毛笔定义为毛笔中心折线以及过各节点上的不同半径Ri的圆,毛笔笔尖为一个点,当毛笔受力发生变形时,各节点上的圆发生变形,形成椭圆。毛笔笔头根部节点上的圆变形较小,该圆在虚拟绘制过程中形状不发生变化。
设毛笔模型中心折线的长度为L,在中心折线上从毛笔笔尖到笔头根部依次设立控制节点P0,P1,P2,…,P9,每相邻两个控制节点之间由线段相连形成折线,由于毛笔尖部相较于毛笔笔头根部更容易产生变形,因此设定毛笔笔尖到笔头根部每相邻两个控制节点之间的线段长度逐渐增大,呈等差数列分布,公差其中P0P1是控制节点P0和P1之间的线段长度,n为中心折线上线段的个数,n=9。
(2)建立毛笔力学模型
毛笔力学模型为弹簧振子模型:在毛笔笔头安装一个垂直于纸面的弹簧,弹簧始终沿纸面法向变形,如附图2所示,在控制节点P9和P9在纸面上的投影点处设置虚拟弹簧,毛笔未受力作用时,与P0重合;当毛笔受力作用时,弹簧受力向下运动,压力F的计算公式为:
其中λ为压力因子,单位为N/mm;X为毛笔受压接触纸面后向下移动的距离,即弹簧变形量,单位为mm;H为毛笔硬度,与F成正比,变化范围为:0<H≤1。W为毛笔湿度,即毛笔内部剩余水墨量与毛笔所蘸含的总水墨量的比值,变化范围为:0≤W≤1,F随着湿度的减小而增大。
(3)确定力与毛笔变形的关系
压力使毛笔产生弯曲变形,如附图3所示,Si为毛笔笔头中心线P0到Pi线段在平面上的投影,一般情况下,压力越大,Si越大,当施加于笔上的压力为0时,笔尖P0与纸面刚刚接触,Si为0;当压力最大时,毛笔笔头的中心线全部投影在纸面上。设p为压力影响因子,p值为压力与最大输出力的比值,将其转换为角度α,其所用的力反馈设备能够提供的最大输出力为7.9牛顿,因此p=F/7.9,变化范围为:0≤p≤1。
Si=(P0P1+P1P2+...+Pi-1Pi)×cosα(2)
其中Pi-1Pi是毛笔中心线控制节点Pi-1到Pi之间的线段长度。
设毛笔沿运动方向移动微小距离,附图3中从A运动到B,由x轴正方向逆时针旋转到向量方向时所经过的角度为γ角,即由x轴正方向逆时针旋转到运动方向的反方向时,所经过的角度为γ角;其取值范围为0≤γ<2π,建立如附图3所示的三维局部坐标系,则毛笔中心线上第i个控制节点Pi的x坐标值和y坐标值计算公式如下:
xi=(S9-Si)×cosγ(3)
yi=(S9-Si)×sinγ(4)
各个控制节点的z坐标值计算公式如下:
z0=0
…
由此可以推导出Pi的z坐标值计算如公式(5)所示:
当毛笔受到压力产生弯曲变形时,椭圆长半轴Rib随着力的增加而变大。当压力F为零时,Rib=Ri;当施加于毛笔的压力F最大时,椭圆长半轴Rib最大,经实际观察,取此时的椭圆长半轴Rib=1.5Ri,Rib的计算公式为:
Rib=Ri×(1+p/2)(6)
本发明假定椭圆的面积不变,则椭圆短半轴Ria的计算公式为:
再根据椭圆长短轴确定毛笔笔道宽度,从而建立毛笔压力与毛笔变形的关系。在虚拟绘制过程中,若要求笔道较宽,可以加大毛笔的压力,若要求笔道较窄,可减小毛笔压力。
本发明在三维虚拟绘制中引入了力反馈技术,提出了一种基于力反馈技术的毛笔建模方法,建立了毛笔压力与毛笔变形的关系。在虚拟绘制过程中,若要求笔道较宽,可以加大毛笔的压力,若要求笔道较窄,可减小毛笔压力。本发明揭示了三维虚拟绘制过程中力对毛笔变形的影响机理,满足了基于力反馈技术的虚拟绘制要求。
附图说明
图1是基于力反馈技术的毛笔几何模型示意图。
图2是毛笔弹簧振子力学模型示意图。
图3是毛笔变形后中心折线控制节点计算示意图。
图4是基于力反馈技术的毛笔建模流程图。
具体实施方式
本具体实施例中,使用的力反馈设备为美国SensAbleTechnologies公司的PhantomDesktop,最大输出力Fmax=7.9牛顿,以虚拟绘制过程中三次采样时间内P3控制节点位置变化与相应椭圆变化为例说明力和毛笔变形的关系,具体过程如下:
1、建立毛笔几何模型,取毛笔的硬度系数H=0.8,毛笔中心线长度为27mm,设定P0P1=2mm,根据公差公式求得d=0.25mm,从而可以得到每相邻两个控制节点之间的线段P0P1,P1P2,…,P8P9长度,本实施例中取R3=3mm。
2、计算压力
本实例中取压力因子λ=0.25N/mm,第一次采样时毛笔湿度W1为0.8,毛笔下压位移X1为6.75mm,根据公式(1)求得第一次采样时力反馈设备施加于毛笔的压力F1=0.75N;第二次采样时毛笔湿度W2为0.5,毛笔下压位移X2为12.5mm,同理求得F2=1.67N;第三次采样时毛笔湿度W3为0.2,毛笔下压位移X3为19.25mm,同理求得F3=3.21N。
3、计算毛笔变形
主要计算三次采样时间内P3控制节点位置坐标和相应的椭圆长半轴和短半轴变化。本实施例中假定三次采样过程中毛笔沿与x轴夹角为π/4的直线运动,因此设γ=π/4;第一次采样时根据压力求得力影响因子p=0.095,根据公式(2)求得S3=1.008,根据公式(3)、(4)、(5)求得P3的坐标为(2.1358,2.1358,6.674),根据公式(6)和(7)求得椭圆长半轴R3b=3.1425mm,短半轴R3a=2.864mm;同理可以求得第二次采样时P3的坐标为(4.6757,4.6757,6.46),椭圆长半轴R3b=3.3165mm,短半轴R3a=2.7137mm;第三次采样时P3的坐标为(8.534,8.534,5.422),椭圆长半轴R3b=3.609mm,短半轴R3a=2.494mm。
由此可以看出毛笔变形后毛笔笔道的宽度随着毛笔压力的变大而变宽。
Claims (1)
1.一种基于力反馈技术的毛笔建模方法,其特征包括以下步骤:
(1)建立基于力反馈技术的毛笔几何模型
将毛笔定义为毛笔中心折线以及过各节点上的不同半径Ri的圆,毛笔笔尖为一个点,当毛笔受力发生变形时,各节点上的圆发生变形,形成椭圆;假设该圆在虚拟绘制过程中形状不发生变化;
设毛笔模型中心折线的长度为L,在中心折线上从毛笔笔尖到笔头根部依次设立控制节点P0,P1,P2,…,P9,每相邻两个控制节点之间由线段相连形成折线,由于毛笔尖部相较于毛笔笔头根部更容易产生变形,因此设定毛笔笔尖到笔头根部每相邻两个控制节点之间的线段长度逐渐增大,呈等差数列分布,公差其中P0P1是控制节点P0和P1之间的线段长度,n为中心折线上线段的个数,n=9;
(2)建立毛笔力学模型,毛笔力学模型为弹簧振子模型
在毛笔笔头安装一个垂直于纸面的弹簧,弹簧始终沿纸面法向变形,在控制节点P9和P9在纸面上的投影点P'9处设置虚拟弹簧,毛笔未受力作用时,P'9与P0重合;当毛笔受力作用时,弹簧受力向下运动,压力F的计算公式为:
其中λ为压力因子,单位为N/mm;X为毛笔受压接触纸面后向下移动的距离,即弹簧变形量,单位为mm;H为毛笔硬度,与F成正比,变化范围为:0<H≤1;W为毛笔湿度,即毛笔内部剩余水墨量与毛笔所蘸含的总水墨量的比值,变化范围为:0≤W≤1,F随着湿度的减小而增大;
(3)确定力与毛笔变形的关系
压力使毛笔产生弯曲变形,Si为毛笔笔头中心线P0到Pi线段在平面xP'9y上的投影,压力越大,Si越大,当施加于笔上的压力为0时,笔尖P0与纸面刚刚接触,Si为0;当压力最大时,毛笔笔头的中心线全部投影在纸面上;设p为压力影响因子,p值为压力与最大输出力的比值,将其转换为角度α,力反馈设备能够提供的最大输出力为7.9牛顿,因此p=F/7.9,变化范围为:0≤p≤1;
Si=(P0P1+P1P2+...+Pi-1Pi)×cosα(2)
其中Pi-1Pi是毛笔中心线控制节点Pi-1到Pi之间的线段长度;
设毛笔沿运动方向移动微小距离,由x轴正方向逆时针旋转到运动方向的反方向时,所经过的角度为γ角,其取值范围为0≤γ<2π,建立三维局部坐标系,则毛笔中心线上第i个控制节点Pi的x坐标值和y坐标值计算公式如下:
xi=(S9-Si)×cosγ(3)
yi=(S9-Si)×sinγ(4)
各个控制节点的z坐标值计算公式如下:
z0=0
…
由此推导出Pi的z坐标值计算如公式(5)所示:
当毛笔受到压力产生弯曲变形时,椭圆长半轴Rib随着力的增加而变大;当压力F为零时,Rib=Ri;当施加于毛笔的压力F最大时,椭圆长半轴Rib最大,取此时的椭圆长半轴Rib=1.5Ri,Rib的计算公式为:
Rib=Ri×(1+p/2)(6)
假定椭圆的面积不变,则椭圆短半轴Ria的计算公式为:
再根据椭圆长短轴确定毛笔笔道宽度,从而建立毛笔压力与毛笔变形的关系;在虚拟绘制过程中,若要求笔道较宽,加大毛笔的压力,若要求笔道较窄,减小毛笔压力。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310248722.4A CN103345773B (zh) | 2013-06-20 | 2013-06-20 | 一种基于力反馈技术的毛笔建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310248722.4A CN103345773B (zh) | 2013-06-20 | 2013-06-20 | 一种基于力反馈技术的毛笔建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103345773A CN103345773A (zh) | 2013-10-09 |
CN103345773B true CN103345773B (zh) | 2015-11-18 |
Family
ID=49280566
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310248722.4A Expired - Fee Related CN103345773B (zh) | 2013-06-20 | 2013-06-20 | 一种基于力反馈技术的毛笔建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103345773B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105404481A (zh) * | 2014-06-23 | 2016-03-16 | 联想(北京)有限公司 | 一种图像信息处理方法及电子设备 |
CN108805968B (zh) * | 2018-05-31 | 2020-07-10 | 西安电子科技大学 | 基于视频的虚拟毛笔建模方法 |
TWI730697B (zh) * | 2020-03-27 | 2021-06-11 | 胡厚飛 | 毛筆檢測機台及毛筆檢測方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11224080A (ja) * | 1998-02-05 | 1999-08-17 | Adk Fuji System Kk | 毛筆文字描画方法 |
CN1773516A (zh) * | 2005-07-08 | 2006-05-17 | 天津大学 | 基于毛笔和宣纸的水墨传输模型与仿真算法 |
CN1870052A (zh) * | 2005-07-08 | 2006-11-29 | 天津大学 | 基于压感笔等装置的三维毛笔模型与仿真方法 |
CN102831292A (zh) * | 2012-07-09 | 2012-12-19 | 大连理工大学 | 一种基于力反馈技术的虚拟水墨扩散范围计算方法 |
-
2013
- 2013-06-20 CN CN201310248722.4A patent/CN103345773B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11224080A (ja) * | 1998-02-05 | 1999-08-17 | Adk Fuji System Kk | 毛筆文字描画方法 |
CN1773516A (zh) * | 2005-07-08 | 2006-05-17 | 天津大学 | 基于毛笔和宣纸的水墨传输模型与仿真算法 |
CN1870052A (zh) * | 2005-07-08 | 2006-11-29 | 天津大学 | 基于压感笔等装置的三维毛笔模型与仿真方法 |
CN102831292A (zh) * | 2012-07-09 | 2012-12-19 | 大连理工大学 | 一种基于力反馈技术的虚拟水墨扩散范围计算方法 |
Non-Patent Citations (3)
Title |
---|
3D实用毛笔模型的物理仿真;孙美君 等;《天津大学学报》;20080331;第41卷(第3期);293-299 * |
DAB: Interactive Haptic Painting with 3D Virtual Brushes;Bill Baxter et al;《ACM Siggraph Conference on Computer Graphics》;20011231;461-468 * |
On the Effects of Haptic Display in Brush and Ink Simulation for Chinese Painting and Calligraphy;Jeng-sheng Yeh et al;《The 10th Pacific Conference on Computer Graphics and Applications》;20021231;439-441 * |
Also Published As
Publication number | Publication date |
---|---|
CN103345773A (zh) | 2013-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102520849B (zh) | 电子毛笔书写方法及系统 | |
CN102436348B (zh) | 一种基于力信息的毛笔书法表达方法 | |
Bixler et al. | Analysis of a swimmer's hand and arm in steady flow conditions using computational fluid dynamics | |
CN103869983B (zh) | 一种用于力触觉人机交互的柔性物体变形仿真方法 | |
CN103559732B (zh) | 一种生成毛笔笔迹的方法 | |
CN106600534A (zh) | 一种基于多尺度小波支持向量机优化的克里金空间插值方法 | |
Fu et al. | Analytical modeling and numerical simulation for three-roll bending forming of sheet metal | |
CN105159459A (zh) | 一种可用于移动终端的虚拟物体三维形状触觉再现方法 | |
CN103345773B (zh) | 一种基于力反馈技术的毛笔建模方法 | |
CN104834380A (zh) | 一种应用于移动终端的柔性物体的触觉建模与表达方法 | |
CN103412978A (zh) | 一种基于工件三维扫描的拉弯工艺模型修正方法 | |
CN103729564A (zh) | 一种基于粒子图像测速技术的压力场计算方法和装置 | |
CN102496176B (zh) | 一种基于力信息的钢笔书法表达方法 | |
CN106384384A (zh) | 一种三维产品模型的形状优化算法 | |
CN103366056A (zh) | 基于有限元分析的机械零件圆柱度误差的测量方法 | |
CN102521856B (zh) | 电子毛笔书写笔形变形方法 | |
CN101635063A (zh) | 一种椭圆截面弹性柱体自由扭转/变形的模拟方法 | |
CN109902418A (zh) | 一种计算eacld中心刚体-悬臂梁模型动力学响应的仿真方法 | |
CN103942397A (zh) | 基于幂函数的修形齿轮数字化建模方法 | |
Zhang et al. | Fuzzy-adaptive control method for off-road vehicle guidance system | |
CN104182569A (zh) | 产品装配误差累积过程的三维可视化模拟方法 | |
CN103543694B (zh) | 一种基于嵌入式智能终端的数控车削加工仿真方法 | |
CN103824322A (zh) | 利用透视关系交互创建三维曲线的方法及系统 | |
CN102332173B (zh) | 一种基于向量代数的简明软笔效果生成方法 | |
CN102708293B (zh) | 一种电极模型和头模型配准的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151118 Termination date: 20190620 |