CN103714205A - Clock spring model simulating rotatable deformation of flexible body - Google Patents
Clock spring model simulating rotatable deformation of flexible body Download PDFInfo
- Publication number
- CN103714205A CN103714205A CN201310713728.4A CN201310713728A CN103714205A CN 103714205 A CN103714205 A CN 103714205A CN 201310713728 A CN201310713728 A CN 201310713728A CN 103714205 A CN103714205 A CN 103714205A
- Authority
- CN
- China
- Prior art keywords
- clockwork spring
- circle
- torsion
- flexible body
- moment
- 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.)
- Granted
Links
Images
Landscapes
- Developing Agents For Electrophotography (AREA)
Abstract
The invention provides a clock spring model simulating rotatable deformation of a flexible body. A plurality of clock springs are sequentially connected in series to form the clock spring model; in an interactive process, signals which are fed back as force touch information which is calculated by the clock spring model and responds to real-time deformation simulation of the flexible body under the action of torque are output, and the sum of the quantities of torsional deformation of all circles of the clock spring model is equivalent to deformation of the surface of the flexible body. The same method is used for calculating the quantities of the torsional deformation of all the circles of the clock spring model, calculation is simple, and the speed of torsional deformation calculation is increased; different types of flexible bodies can be simulated by adjusting the radiuses of rotatable mandrels of the clock springs, the thicknesses of the clock springs and the suspension section widths, and applicability is high; the clock spring model can be applied to the fields of virtual surgery simulation, remote operation robot control, telemedicine and the like.
Description
Technical field
The invention belongs to computer simulation technique field, relate in particular to a kind of clockwork spring model of simulating the rotatable distortion of flexible body.
Background technology
Virtual operation emulation is the important application of virtual reality technology in modern operative training, mainly by visual, the moment of torsion deformation emulating stressed with modeling, human body soft tissue organ of medical data, formed, the true reappearance of surgical scene is provided for user on vision and sense of touch sense organ.Compare with traditional operation training, virtual operation can be simulated various physics and the physiological property of biological tissue, has greatly improved the effect of operative training; The experimental subjects that can recycle that simultaneously can also be different according to need to providing of doctor, greatly reduces the cost of surgeon's operative training.
It is the gordian technique of virtual operation emulation that the human body flexible bodily form becomes.At present, mainly contain the deformation models such as Mass-spring model, finite element model, boundary element model.That spring-particle modeling method has is simple, calculated amount is less, be easy to the advantages such as realization, but has the drawbacks such as topological analysis difficulty.Although finite element modeling method can carry out accurate and quantitative simulation to the deformation of object, relate to a large amount of numerous and diverse calculating, real-time, interactive poor performance.Though it is discrete that boundary element modeling method is carried out the border of modeling method, simplified calculating, aspect stability, there is certain deficiency.More than explanation, all there is the problems such as the comparatively numerous and diverse and simulation accuracy of calculating is not high in these conventional flexible bodies distortion Method of Physical Modeling, has then affected real-time and the validity calculated.Therefore, set up the accurately haptic modeling method of efficient virtual flexible body deformation simulation, become current fictitious force haptic-display system matter of utmost importance urgently to be resolved hurrily.
In view of the above problems, in virtual operation simulation process, the combination of vision and power sense of touch makes virtual emulation truer, credible, has proposed a kind of clockwork spring model of simulating the rotatable deformation simulation of flexible body.
Summary of the invention
Technical matters to be solved by this invention is to overcome the deficiencies in the prior art, and a kind of clockwork spring model of simulating the rotatable distortion of flexible body is provided.
For solving the problems of the technologies described above, the technical solution adopted in the present invention is:
A clockwork spring model of simulating the rotatable distortion of flexible body, comprises the steps:
Step 1, sets up rectangular coordinate system in space, determines the clockwork spring of taking up an official post flexible body surface meaning point place laying, and its process is as follows:
Step 1.1, sets up rectangular coordinate system in space,
Under given moment of torsion M effect, on flexible body surface, a clockwork spring is laid at O place, arbitrfary point, is being that r place arranges the rotatable mandrel points outside of clockwork spring A apart from arbitrfary point O
0, r is the rotatable mandrel radius of clockwork spring, take arbitrfary point O as initial point, ray OA
0place direction is X-axis positive dirction, sets up XYZ space coordinate system;
Step 1.2, arranges clockwork spring by circle,
Wherein, A outside i-1 circle clockwork spring
iX1point, suspension thickness is h
i=h
1+ (i-1) i of d encloses clockwork spring, wherein h
1be the thickness of the 1st circle clockwork spring, d is given constant and d>0, i circle clockwork spring with initial point O (0,0,0) for the center of circle, with
for the circle of radius, i=1,2,3 ... S, S is natural number,
Wide b, the cross sectional moment of inertia of being of every circle clockwork spring
elastic modulus E depends on the material of flexible body, and flexible body material is all identical;
Step 2, determines the moment of torsion that arbitrary circle clockwork spring consumes;
Set active line and the rotatable mandrel points outside of the clockwork spring A of given moment of torsion M
0place circle tangent, and under given moment of torsion M effect, if P circle clockwork spring produces torsional deflection before total in flexible body, P circle clockwork spring is called distortion check loop; And P≤S, namely the number of turns of clockwork spring at least equals P;
According to clockwork spring characteristic, set:
On front P-1 circle clockwork spring, any point is under given moment of torsion M effect, and the maximum twist power that consumes is equal and be F
0; On P circle clockwork spring, any point is under given moment of torsion M effect, the maximum twist power consuming
all equate, and be not more than F
0;
Step 2.1, in P-1 circle clockwork spring, is total to a little the moment of torsion M consuming before determining on i circle
ifor:
Step 2.2, determines that P circle clockwork spring is total to a little the moment of torsion consuming
for:
On P circle clockwork spring, any point is under given moment of torsion M effect, the maximum twist power of consumption
for:
Step 3, calculates on i circle clockwork spring institute a little, at the moment of torsion M of common consumption
iunder effect, the torsion number of turns n of generation
i;
Wherein, l
ibe the effective length of i circle clockwork spring, its expression formula is as follows:
Described a kind of clockwork spring model of simulating the rotatable distortion of flexible body, the time delay time that described front P circle clockwork spring produces torsional deflection total to be needed is not more than 1ms, also meets the requirement that refreshing frequency is not less than 1000Hz.
Described a kind of clockwork spring model of simulating the rotatable distortion of flexible body, in all clockwork springs, the time delay time that every circle produces torsional deflection to be needed forms Geometric Sequence, meets:
t
i=q
i-1t
1
Wherein, t
irepresent the time delay time that i circle clockwork spring produces torsional deflection needs, the common ratio that q is Geometric Sequence, t
1be that the 1st circle clockwork spring produces the time delay time that torsional deflection needs, 1≤i≤P.
Described a kind of clockwork spring model of simulating the rotatable distortion of flexible body, described clockwork spring model, on arbitrary circle clockwork spring, the twisting resistance of any point consumption reaches after maximum twist power, and its next circle starts to produce torsional deflection.
The invention has the beneficial effects as follows: the present invention proposes a kind of clockwork spring model of simulating the rotatable distortion of flexible body, described clockwork spring model is composed in series successively by a plurality of clockwork springs, in reciprocal process, the signal that output feedback is reaction power tactile data of flexible body real-time deformation emulation under torsional interaction of adopting clockwork spring model and calculating, the stack that produces torsional deflection amount sum in this clockwork spring model on all circles is externally equivalent to the distortion on flexible body surface.The torsional deflection amount computing method of the every circle clockwork spring of clockwork spring model of the present invention are identical, calculate simply, accelerated the speed that torsional deflection is calculated; By regulating the rotatable mandrel radius of clockwork spring, the thickness of clockwork spring and suspension cross-sectional width etc., just can simulate dissimilar flexible body, and applicability is wide; Can be applicable to the fields such as virtual surgery emulation, teleoperated vehicle's control, tele-medicine.
Accompanying drawing explanation
Fig. 1 is clockwork spring model schematic diagram.
Fig. 2 is flexible body deformation simulation process flow diagram.
Fig. 3 is clockwork spring model building method process flow diagram.
Fig. 4 is clockwork spring model moment of torsion, the torsion number of turns and time delay time relationship schematic diagram.
Embodiment
A kind of clockwork spring model of simulating the rotatable distortion of flexible body the present invention being proposed below in conjunction with flow process shown in accompanying drawing is elaborated:
Clockwork spring model schematic diagram as shown in Figure 1.A clockwork spring model of simulating the rotatable distortion of flexible body, its concrete steps are as follows:
Step 1, sets up rectangular coordinate system in space, determines the clockwork spring of taking up an official post flexible body surface meaning point place laying, and its process is as follows:
Step 1.1, sets up rectangular coordinate system in space,
Under given moment of torsion M effect, on flexible body surface, a clockwork spring is laid at O place, arbitrfary point, is being that r place arranges the rotatable mandrel points outside of clockwork spring A apart from arbitrfary point O
0, r is the rotatable mandrel radius of clockwork spring, take arbitrfary point O as initial point, ray OA
0place direction is X-axis positive dirction, sets up XYZ space coordinate system;
Step 1.2, arranges clockwork spring by circle,
Wherein, A outside i-1 circle clockwork spring
i-1point, suspension thickness is h
i=h
1+ (i-1) i of d encloses clockwork spring, wherein h
1be the thickness of the 1st circle clockwork spring, d is given constant and d>0, i circle clockwork spring with initial point O (0,0,0) for the center of circle, with
for the circle of radius, i=1,2,3 ... S, S is natural number,
Wide b, the cross sectional moment of inertia of being of every circle clockwork spring
elastic modulus E depends on the material of flexible body, and flexible body material is all identical;
Step 2, determines the moment of torsion that arbitrary circle clockwork spring consumes;
Suppose active line and the rotatable mandrel points outside of the clockwork spring A of given moment of torsion M
0place circle tangent, and under given moment of torsion M effect, if P circle clockwork spring produces torsional deflection before total in flexible body, P circle clockwork spring is called distortion check loop; And P≤S, namely the number of turns of clockwork spring at least equals P;
According to clockwork spring characteristic, set:
On front P-1 circle clockwork spring, any point is under given moment of torsion M effect, and the maximum twist power that consumes is equal and be F
0; On P circle clockwork spring, any point is under given moment of torsion M effect, the maximum twist power consuming
all equate, and be not more than F
0;
Step 2.1, in P-1 circle clockwork spring, is total to a little the moment of torsion M consuming before determining on i circle
ifor:
Step 2.2, determines that P circle clockwork spring is total to a little the moment of torsion consuming
for:
On P circle clockwork spring, any point is under given moment of torsion M effect, the maximum twist power of consumption
for:
Step 3, calculates on i circle clockwork spring institute a little, at the moment of torsion M of common consumption
iunder effect, the torsion number of turns n of generation
i;
Wherein, l
ibe the effective length of i circle clockwork spring, its expression formula is as follows:
Described a kind of clockwork spring model of simulating the rotatable distortion of flexible body, the time delay time that described front P circle clockwork spring produces torsional deflection total to be needed is not more than 1ms, also meets the requirement that refreshing frequency is not less than 1000Hz, 1≤i≤P.
Described a kind of clockwork spring model of simulating the rotatable distortion of flexible body, in all clockwork springs, the time delay time that every circle produces torsional deflection to be needed forms Geometric Sequence, meets:
t
i=q
i-1t
1
Wherein, t
irepresent the time delay time that i circle clockwork spring produces torsional deflection needs, the common ratio that q is Geometric Sequence, t
1be that the 1st circle clockwork spring produces the time delay time that torsional deflection needs.
Described a kind of clockwork spring model of simulating the rotatable distortion of flexible body, described clockwork spring model, on arbitrary circle clockwork spring, the twisting resistance of any point consumption reaches after maximum twist power, and its next circle starts to produce torsional deflection.
Virtual hand and Virtual Cardiac Mode take below as example, enumerate the embodiment of technical solution of the present invention.
In this example, all virtual hand and Virtual Cardiac Mode all directly adopt the OBJ form of deriving from 3DS MAX 2013 softwares, with 1558 particles, 3114 virtual hand and 3910 particles that triangle gridding forms, the Virtual Cardiac Mode of 7814 triangle gridding formations is that example is carried out deformation simulation, and in experimentation, model obtains and revise very convenient; Operating system is Windows 2000, take 3DS MAX 2013, OpenGL shape library as basis, on Microsoft Visual C++2012 Software Development Platform, carries out emulation.
When virtual hand being detected and collide on virtual heart surface any point, in given counterclockwise moment of torsion M=10 * 10
-3under Nm effect, the inner clockwork spring model of filling of regional area that virtual hand and virtual heart are mutual, in reciprocal process, output is fed back to the signal that adopts reaction power tactile data of virtual heart real-time deformation emulation under torsional interaction that clockwork spring model calculates, as shown in Figure 2;
Under given counterclockwise moment of torsion M effect, when virtual hand collides on virtual heart surface arbitrfary point O, on virtual heart surface, a clockwork spring is laid at O place, arbitrfary point, apart from arbitrfary point O, is being r=1 * 10
-3m place arranges the rotatable mandrel points outside of clockwork spring A
0, r is the rotatable mandrel radius of clockwork spring, take arbitrfary point O as initial point, ray OA
0place direction is X-axis positive dirction, sets up XYZ space coordinate system;
Set gradually each circle of clockwork spring, as shown in figures 1 and 3; Wherein, A outside i-1 circle clockwork spring
i-1point, suspension thickness is h
i=h
1+ (i-1) i of the clockwork spring of d encloses, h
1=0.4 * 10
-3m, d=0.1 * 10
-3m, i circle clockwork spring with initial point O (0,0,0) for the center of circle, with
for the circle of radius, i=1,2,3 ... S, S is natural number;
With O (0,0,0) for the center of circle, with
R
1=r+h
1=1 * 10
-3+ 0.4 * 10
-3=1.4 * 10
-3m is the circle of radius, forms the 1st circle clockwork spring;
R
2=r+h
1+ h
2=r+2h
1+ (2-1) d=1 * 10
-3+ 2 * 0.4 * 10
-3+ 1 * 0.1 * 10
-3=1.9 * 10
-3m is the circle of radius, forms the 2nd circle clockwork spring;
R
3=r+h
1+ h
2+ h
3=r+3h
1+ 3d=1 * 10
-3+ 3 * 0.4 * 10
-3+ 3 * 0.1 * 10
-3=2.5 * 10
-3m is the circle of radius, forms the 3rd circle clockwork spring;
R
4=r+h
1+ h
2+ h
3+ h
4=r+4h
1+ 6d=1 * 10
-3+ 4 * 0.4 * 10
-3+ 6 * 0.1 * 10
-3=3.2 * 10
-3m is the circle of radius, forms the 4th circle clockwork spring;
R
5=r+h
1+ h
2+ h
3+ h
4+ h
5=r+5h
1+ 10d=1 * 10
-3+ 5 * 0.4 * 10
-3+ 10 * 0.1 * 10
-3=4 * 10
-3m is the circle of radius, forms the 5th circle clockwork spring;
Suppose that the wide of every circle clockwork spring is b=6 * 10
-3m, elastic modulus E=3.09 * 10
7pa depends on the material of flexible body, all identical;
The pilot process calculating, the equal round off method of last data retain after radix point 3.
Suppose given counterclockwise moment of torsion M=10 * 10
-3the rotatable mandrel points outside of the active line of Nm and clockwork spring A
0place circle tangent, and under given moment of torsion M effect, if P circle clockwork spring produces torsional deflection before total in flexible body, P circle clockwork spring is called distortion check loop;
According to clockwork spring property settings: as shown in Figure 4; On front P-1 circle clockwork spring, any point is under given counterclockwise moment of torsion M effect, and the maximum twist power consuming equates and is F
0=0.989N, on P circle clockwork spring any point under given counterclockwise moment of torsion M effect, the maximum twist power consuming
all equate, and be not more than F
0=0.989N;
The time delay time that front i circle clockwork spring produces torsional deflection total to be needed is not more than 1ms, 1≤i≤P; And in each circle of all clockwork springs, the time delay time formation that every circle produces torsional deflection to be needed produces with the 1st circle clockwork spring the time delay time t that torsional deflection needs
1=10
-5the Geometric Sequence that s is first term, the q=1.2 of take is common ratio; On arbitrary circle clockwork spring, the twisting resistance of any point consumption reaches after maximum twist power, and its next circle starts to produce torsional deflection.
Suppose that haptic feedback refreshing frequency is 1200Hz, the inverse of haptic feedback refreshing frequency
If any point is under given counterclockwise moment of torsion M effect on the 1st circle clockwork spring, the maximum twist power consuming equates also all to arrive F
0, on the 1st circle clockwork spring, be total to a little the moment of torsion M consuming
1for:
M
1=F
0R
1=0.989×1.4×10
-3=1.385×10
-3N·m,
M
1<M,
The effective length l of the 1st circle clockwork spring
1=2 π R
1,
The cross sectional moment of inertia of the 1st circle clockwork spring
The 1st circle clockwork spring on institute a little, at the moment of torsion M of common consumption
1under effect, with the torsion number of turns n of its generation
1between meet:
The 1st circle clockwork spring produces the time delay time T that torsional deflection needs
1=t
1=10
-5s<T;
Therefore, on the 1st circle clockwork spring, be total to a little the moment of torsion M consuming
1<M, and the 1st circle clockwork spring produces the time delay time T that torsional deflection total needs
1=10
-5s<T, meets the requirement of refreshing frequency; Only has the maximum twist power F that reaches consumption when the 1st twisting resistance that encloses any point consumption on clockwork spring
0after, the 2nd circle clockwork spring just starts to produce torsional deflection.
If any point is under given counterclockwise moment of torsion M effect on the 2nd circle clockwork spring, the maximum twist power consuming equates also all to arrive F
0, on the 2nd circle clockwork spring, be total to a little the moment of torsion M consuming
2for:
M
2=F
0R
2=0.989×1.9×10
-3=1.879×10
-3N·m,
On front 2 circle clockwork springs, being total to a little the moment of torsion sum consuming is: M
1+ M
2=1.385 * 10
-3+ 1.879 * 10
-3=3.264 * 10
-3nm<M=10 * 10
-3nm,
The effective length l of the 2nd circle clockwork spring
2=2 π (R
1+ R
2),
The cross sectional moment of inertia of the 2nd circle clockwork spring:
The 2nd circle clockwork spring on institute a little, at the moment of torsion M of common consumption
2under effect, with the torsion number of turns n of its generation
2between meet:
Front 2 circle clockwork springs produce the time delay time T that torsional deflection needs
2=t
1+ t
2=(1+q) t
1=(1+1.2) * 10
-5=2.2 * 10
-5s<T, T is the inverse of haptic feedback refreshing frequency here,
Therefore, on front 2 circle clockwork springs, be total to a little the moment of torsion M consuming
1+ M
2<M, and front 2 circle clockwork springs produce the time delay time T that torsional deflection total needs
2=2.2 * 10
-5s<T, meets the requirement of refreshing frequency; Only has the maximum twist power F that reaches consumption when the 2nd twisting resistance that encloses any point consumption on clockwork spring
0after, the 3rd circle clockwork spring just starts to produce torsional deflection.
If any point is under given counterclockwise moment of torsion M effect on the 3rd circle clockwork spring, the maximum twist power consuming equates also all to arrive F
0, on the 3rd circle clockwork spring, be total to a little the moment of torsion M consuming
3for:
M
3=F
0R
3=0.989×2.5×10
-3=2.473×10
-3N·m,
On front 3 circle clockwork springs, being total to a little the moment of torsion sum consuming is:
M
1+M
2+M
3=1.385×10
-3+1.879×10
-3+2.473×10
-3=5.737×10
-3N·m<M=10×10
-3N·m
The effective length l of the 3rd circle clockwork spring
3=2 π (R
1+ R
2+ R
3),
The cross sectional moment of inertia of the 3rd circle clockwork spring:
The 3rd circle clockwork spring on institute a little, at the moment of torsion M of common consumption
3under effect, with the torsion number of turns n of its generation
3between meet:
Front 3 circle clockwork springs produce the time delay time that torsional deflection needs:
T
3=t
1+t
2+t
3=(1+q+q
2)t
1=(1+1.2+1.2
2)×10
-5=3.64×10
-5s<T,
Here T is the inverse of haptic feedback refreshing frequency,
Therefore, on front 3 circle clockwork springs, be total to a little the moment of torsion M consuming
1+ M
2+ M
3<M, and front 3 circle clockwork springs produce the time delay time T that torsional deflection total needs
3=3.64 * 10
-5s<T, meets the requirement of refreshing frequency; Only has the maximum twist power F that reaches consumption when the 3rd twisting resistance that encloses any point consumption on clockwork spring
0after, the 4th circle clockwork spring just starts to produce torsional deflection.
If any point is under given counterclockwise moment of torsion M effect on the 4th circle clockwork spring, the maximum twist power consuming equates also all to arrive F
0, on the 4th circle clockwork spring, be total to a little the moment of torsion M consuming
4for:
M
4=F
0R
4=0.989×3.2×10
-3=3.165×10
-3N·m,
On front 4 circle clockwork springs, being total to a little the moment of torsion sum consuming is:
M
1+M
2+M
3+M
4=1.385×10
-3+1.879×10
-3+2.473×10
-3+3.165×10
-3,
=8.902×10
-3N·m<M=10×10
-3N·m
The effective length l of the 4th circle clockwork spring
4=2 π (R
1+ R
2+ R
3+ R
4),
The cross sectional moment of inertia of the 4th circle clockwork spring:
The 4th circle clockwork spring on institute a little, at the moment of torsion M of common consumption
4under effect, with the torsion number of turns n of its generation
4between meet:
Front 4 circle clockwork springs produce the time delay time that torsional deflection needs:
T
4=t
1+t
2+t
3+t
4=(1+q+q
2+q
3)t
1=5.368×10
-5s<T,
Here T is the inverse of haptic feedback refreshing frequency,
Therefore, on front 4 circle clockwork springs, be total to a little the moment of torsion M consuming
1+ M
2+ M
3+ M
4<M, and front 4 circle clockwork springs generation torsional deflections amount to T between the time delay needing
4=5.368 * 10
-5s<T, meets the requirement of refreshing frequency; Only has the maximum twist power F that reaches consumption when the 4th twisting resistance that encloses any point consumption on clockwork spring
0after, the 5th circle clockwork spring just starts to produce torsional deflection.
If any point is under given counterclockwise moment of torsion M effect on the 5th circle clockwork spring, the maximum twist power consuming equates also all to arrive F
0, on the 5th circle clockwork spring, be total to a little the moment of torsion M consuming
5for:
M
5=F
0R
5=0.989×4×10
-3=3.956×10
-3N·m,
On front 5 circle clockwork springs, being total to a little the moment of torsion sum consuming is:
M
1+M
2+M
3+M
4+M
5=(1.385+1.879+2.473+3.165+3.956)×10
-3
=12.858×10
-3N·m>M=10×10
-3N·m
Therefore, on front 5 circle clockwork springs, be total to a little the moment of torsion M consuming
1+ M
2+ M
3+ M
4+ M
5>M, is total to a little the moment of torsion sum consuming on front 5 circle clockwork springs and is not less than given counterclockwise moment of torsion M, and the 5th circle clockwork spring, for distortion check loop, does not need to judge whether to meet the requirement of refreshing frequency again;
The 5th circle clockwork spring is out of shape the moment of torsion that is total to a little consumption on check loop
for:
Distortion check loop also on the 5th circle clockwork spring any point under given counterclockwise moment of torsion M effect, the maximum twist power of consumption
for:
According to clockwork spring property settings: on front 4 circle clockwork springs, any point is under given counterclockwise moment of torsion M effect, and the maximum twist power of consumption equates to be also all assumed to F
0=0.989N, distortion check loop also on the 5th circle clockwork spring any point under given counterclockwise moment of torsion M effect, the maximum twist power of consumption
all equate to be also all not more than F
0=0.989N;
The effective length of the 5th circle clockwork spring
The cross sectional moment of inertia of the 5th circle clockwork spring:
The 5th circle clockwork spring is out of shape on check loop institute a little, at the moment of torsion M of common consumption
5under effect, with the torsion number of turns n of its generation
5between meet:
Front 5 circle clockwork springs are in given counterclockwise moment of torsion M=10 * 10
-3under Nm effect, the raw torsion of common property several n that turn-take are:
n=n
1+n
2+n
3+n
4+n
5=1.961+3.211+4.298+5.375+1.274=16.119。
Attention: adopting the clockwork spring model of the rotatable distortion of flexible body to calculate under given counterclockwise torsional interaction, in the process of flexible body real-time deformation emulation, if r, d, h
1, these parameters of b choose excessive, in the modeling method of the clockwork spring model of the rotatable distortion of flexible body, the distortion check loop numerical value of clockwork spring is just few, calculated amount is little, real-time is good, but deformation simulation poor effect; If r, d, h
1, these parameters of b choose too small, in the modeling method of the clockwork spring model of the rotatable distortion of flexible body, the distortion check loop numerical value of clockwork spring is just large, calculated amount is large, real-time is not good, but deformation simulation effect is better; In addition t is being set
1and t
ibetween be related to time, should be taken into account the hardware configuration of program operation computer-chronograph itself, therefore in the process of the whole program of debugging, compromise and select these parameters, constantly repeatedly debug, thereby make deformation effect more true to nature.
For verifying implementation result of the present invention, operator is by the deformation simulation that the handle of PHANTOM OMNI hand controller end touches, perception and control virtual hand are rotated virtual heart, and by the power tactile data Real-time Feedback producing in reciprocal process to operator.Experimental result shows: this model is simply effective, and deformation effect is true to nature, and image is smooth, and power tactile sensation is steady, and operator to the perception of virtual environment and alternately accurately and reliably, can meet the requirement of virtual operation simulation interactive system.
Claims (4)
1. a clockwork spring model of simulating the rotatable distortion of flexible body, is characterized in that, comprises the steps:
Step 1, sets up rectangular coordinate system in space, determines the clockwork spring of taking up an official post flexible body surface meaning point place laying, and its process is as follows:
Step 1.1, sets up rectangular coordinate system in space,
Under given moment of torsion M effect, on flexible body surface, a clockwork spring is laid at O place, arbitrfary point, is being that r place arranges the rotatable mandrel points outside of clockwork spring A apart from arbitrfary point O
0, r is the rotatable mandrel radius of clockwork spring, take arbitrfary point O as initial point, ray OA
0place direction is X-axis positive dirction, sets up XYZ space coordinate system;
Step 1.2, arranges clockwork spring by circle,
Wherein, A outside i-1 circle clockwork spring
i-1point, suspension thickness is h
i=h
1+ (i-1) i of d encloses clockwork spring, wherein h
1be the thickness of the 1st circle clockwork spring, d is given constant and d>0, i circle clockwork spring be with initial point O (0,0,0) for the center of circle, with
for the circle of radius, i=1,2,3 ... S, S is natural number,
Wide b, the cross sectional moment of inertia of being of every circle clockwork spring
elastic modulus E depends on the material of flexible body, and flexible body material is all identical;
Step 2, determines the moment of torsion that arbitrary circle clockwork spring consumes;
Set active line and the rotatable mandrel points outside of the clockwork spring A of given moment of torsion M
0place circle tangent, and under given moment of torsion M effect, if P circle clockwork spring produces torsional deflection before total in flexible body, P circle clockwork spring is called distortion check loop; And P≤S, namely the number of turns of clockwork spring at least equals P;
According to clockwork spring characteristic, set:
On front P-1 circle clockwork spring, any point is under given moment of torsion M effect, and the maximum twist power that consumes is equal and be F
0; On P circle clockwork spring, any point is under given moment of torsion M effect, the maximum twist power consuming
all equate, and be not more than F
0;
Step 2.1, in P-1 circle clockwork spring, is total to a little the moment of torsion M consuming before determining on i circle
ifor:
Step 2.2, determines that P circle clockwork spring is total to a little the moment of torsion consuming
for:
On P circle clockwork spring, any point is under given moment of torsion M effect, the maximum twist power of consumption
for:
Step 3, calculates on i circle clockwork spring institute a little, at the moment of torsion M of common consumption
iunder effect, the torsion number of turns n of generation
i;
Wherein, l
ibe the effective length of i circle clockwork spring, its expression formula is as follows:
2. a kind of clockwork spring model of simulating the rotatable distortion of flexible body according to claim 1, is characterized in that, the time delay time that described front P circle clockwork spring produces torsional deflection total to be needed is not more than 1ms.
3. a kind of clockwork spring model of simulating the rotatable distortion of flexible body according to claim 1, is characterized in that, in all clockwork springs, the time delay time that every circle produces torsional deflection to be needed forms Geometric Sequence, meets:
t
i=q
i-1t
1
Wherein, t
irepresent the time delay time that i circle clockwork spring produces torsional deflection needs, the common ratio that q is Geometric Sequence, t
1be that the 1st circle clockwork spring produces the time delay time that torsional deflection needs, 1≤i≤P.
4. according to a kind of clockwork spring model of simulating the rotatable distortion of flexible body described in claim 1 or 2 or 3, it is characterized in that, described clockwork spring model, on arbitrary circle clockwork spring, the twisting resistance of any point consumption reaches after maximum twist power, and its next circle starts to produce torsional deflection.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310713728.4A CN103714205B (en) | 2013-12-23 | 2013-12-23 | A kind of clockwork spring model simulating the rotatable deformation of flexible body |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310713728.4A CN103714205B (en) | 2013-12-23 | 2013-12-23 | A kind of clockwork spring model simulating the rotatable deformation of flexible body |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103714205A true CN103714205A (en) | 2014-04-09 |
CN103714205B CN103714205B (en) | 2016-10-05 |
Family
ID=50407176
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310713728.4A Expired - Fee Related CN103714205B (en) | 2013-12-23 | 2013-12-23 | A kind of clockwork spring model simulating the rotatable deformation of flexible body |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103714205B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107080586A (en) * | 2017-03-24 | 2017-08-22 | 南京信息工程大学 | It is used for the cross spring model of the real-time deformation of soft tissue in virtual operation |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101964023A (en) * | 2010-09-27 | 2011-02-02 | 南京信息工程大学 | Modeling method of asymmetric plate spring virtual model capable of realizing flexible force tactile representation |
WO2012099709A2 (en) * | 2011-01-20 | 2012-07-26 | Nothwestern University | Further improvements to ankle-foot prosthesis and orthosis capable of automatic adaptation to sloped walking surfaces and methods of use |
-
2013
- 2013-12-23 CN CN201310713728.4A patent/CN103714205B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101964023A (en) * | 2010-09-27 | 2011-02-02 | 南京信息工程大学 | Modeling method of asymmetric plate spring virtual model capable of realizing flexible force tactile representation |
WO2012099709A2 (en) * | 2011-01-20 | 2012-07-26 | Nothwestern University | Further improvements to ankle-foot prosthesis and orthosis capable of automatic adaptation to sloped walking surfaces and methods of use |
Non-Patent Citations (3)
Title |
---|
崔桐等: "一种用于力觉再现的柔性体变形仿真弹簧-质点模型", 《东南大学学报(自然科学版)》, vol. 37, no. 5, 30 September 2007 (2007-09-30), pages 849 - 852 * |
张小瑞等: "用于软组织变形仿真的层状菱形连接模型", 《计算机辅助设计与图形学学报》, vol. 22, no. 3, 31 March 2010 (2010-03-31), pages 473 - 479 * |
张小瑞等: "虚拟物体力/触觉交互算法的研究进展", 《系统仿真学报》, vol. 23, no. 4, 30 April 2011 (2011-04-30), pages 637 - 647 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107080586A (en) * | 2017-03-24 | 2017-08-22 | 南京信息工程大学 | It is used for the cross spring model of the real-time deformation of soft tissue in virtual operation |
CN107080586B (en) * | 2017-03-24 | 2019-05-21 | 南京信息工程大学 | The cross spring model of the real-time deformation of soft tissue is used in virtual operation |
Also Published As
Publication number | Publication date |
---|---|
CN103714205B (en) | 2016-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101587372B (en) | Modeling method for enhanced force tactile of virtual reality human-computer interaction | |
CN101976298B (en) | Modeling method of symmetrical type plate spring virtual model enhancing haptic feedback | |
CN103869983B (en) | A kind of flexible article deformation simulation method mutual for power haptic | |
CN101976299B (en) | Modeling method of symmetrical laminated spring virtual model for flexile force touch reappearance | |
CN103400023A (en) | Soft tissue deformation simulation method | |
CN101964023B (en) | Modeling method of asymmetric plate spring virtual model capable of realizing flexible force tactile representation | |
CN107992672A (en) | A kind of soft tissue deformation modeling method based on virtual spring | |
CN101496028A (en) | Method of simulating deformable object using geometrically motivated model | |
CN110289104B (en) | Soft tissue compression and deformation recovery simulation method | |
CN104063057B (en) | Volute spiral spring simulation system supporting haptic feedback enhancement | |
CN106485000A (en) | Based on virtual spring and spherical harmonic function Force Feedback Model modeling method | |
Zhang et al. | An optimized mass-spring model with shape restoration ability based on volume conservation | |
CN103745058A (en) | Method for simulating tension/deformation on soft tissue epidermis of any shape | |
CN107146288A (en) | The soft tissue model modeling method of pressing deformation in real time is supported in virtual operation | |
CN105955458B (en) | Feedback force generation method and device | |
CN106528993A (en) | Modeling method based on combined spring virtual model composed of disc springs | |
CN107590851B (en) | Soft deformation model construction method based on mixing of mass spring and filler model | |
CN106295085A (en) | The modeling method of the overlapping disk spring dummy model of flexible force tactile sense reproduction | |
CN103714205A (en) | Clock spring model simulating rotatable deformation of flexible body | |
CN107080586A (en) | It is used for the cross spring model of the real-time deformation of soft tissue in virtual operation | |
CN101901565B (en) | Virtual flexible body deformation operation simulation system supporting haptic feedback | |
CN103699753B (en) | A kind of stranded wire helical spring dummy model supporting to strengthen haptic device | |
CN106485028B (en) | Support the flexible body deformation simulation modeling method pressed in real time | |
CN103699741A (en) | Rotationally deformable clockwork spring model of strengthened simulative flexible body | |
CN103678806A (en) | Incense coil spring model simulating flexible body torsional deformation |
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: 20161005 Termination date: 20211223 |