JP7180274B2 - Information processing device, display program and display method - Google Patents

Information processing device, display program and display method Download PDF

Info

Publication number
JP7180274B2
JP7180274B2 JP2018198329A JP2018198329A JP7180274B2 JP 7180274 B2 JP7180274 B2 JP 7180274B2 JP 2018198329 A JP2018198329 A JP 2018198329A JP 2018198329 A JP2018198329 A JP 2018198329A JP 7180274 B2 JP7180274 B2 JP 7180274B2
Authority
JP
Japan
Prior art keywords
shape
cable
probability distribution
linear
flexible object
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.)
Active
Application number
JP2018198329A
Other languages
Japanese (ja)
Other versions
JP2020067694A (en
Inventor
左千夫 小林
大地 島田
秀城 阿部
弘樹 小林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2018198329A priority Critical patent/JP7180274B2/en
Publication of JP2020067694A publication Critical patent/JP2020067694A/en
Application granted granted Critical
Publication of JP7180274B2 publication Critical patent/JP7180274B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本発明は、情報処理装置、表示プログラム及び表示方法に関する。 The present invention relates to an information processing device, a display program, and a display method.

例えば、工場の製造の効率化には、組立作業のロボット化が望まれている。組立作業にはケーブル等の線状柔軟物を扱う工程が含まれているため、ロボット動作の妥当性を短時間で検証するために、柔軟物の変形シミュレーションが用いられる。 For example, in order to improve the efficiency of manufacturing in factories, robotization of assembly work is desired. Since the assembly work includes the process of handling linear flexible objects such as cables, deformation simulation of flexible objects is used to verify the validity of robot motions in a short time.

特開2005-306187号公報JP 2005-306187 A 特開2005-38398号公報JP-A-2005-38398

柔軟物には、長さ、剛性および重量に公差が存在するため、これらの公差を考慮して変形シミュレーションを行うようにすると、柔軟物の存在位置として所定の範囲を示すことができる。 Since flexible objects have tolerances in length, rigidity, and weight, if deformation simulation is performed with these tolerances taken into account, it is possible to indicate a predetermined range as the position of the flexible object.

また、変形シミュレーションにおいて上記公差以外の要素についても考慮することができれば、柔軟物の存在位置を更に精度よくシミュレーションすることができると考えられる。 Also, if it is possible to consider factors other than the tolerances in the deformation simulation, it is possible to simulate the existence position of the flexible object with higher accuracy.

1つの側面では、本発明は、線状柔軟物が存在する可能性のある範囲を精度よく特定し、表示することができる情報処理装置、表示プログラム及び表示方法を提供することを目的とする。 In one aspect, an object of the present invention is to provide an information processing device, a display program, and a display method that can accurately identify and display a range in which a linear flexible object may exist.

一つの態様では、情報処理装置は、線状柔軟物に複数の拘束点を設定する設定部と、設定した前記複数の拘束点のうち前記線状柔軟物上で隣接する2つの拘束点間における前記線状柔軟物の長さ誤差の確率分布を、前記2つの拘束点それぞれの前記線状柔軟物の長さ方向に関する位置ずれの確率分布から求め、求めた前記長さ誤差の確率分布に基づいて、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の尤もらしさを示す指標とを算出する算出部と、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の前記指標とに基づいて、前記線状柔軟物が存在する可能性のある範囲を特定し、表示する表示処理部と、を備えている。 In one aspect, an information processing apparatus includes a setting unit that sets a plurality of constraint points on a linear flexible object; The probability distribution of the length error of the flexible linear object is obtained from the probability distribution of the displacement of the two restraint points in the length direction of the flexible linear object, and based on the obtained probability distribution of the length error. a calculator for calculating a shape that the linear flexible object can take between the two constraint points and an index indicating the likelihood of the shape; a display processing unit that identifies and displays a range in which the linear flexible object may exist based on the shape and the index of the shape.

線状柔軟物が存在する可能性のある範囲を精度よく特定し、表示することができる。 A range in which a linear flexible object may exist can be accurately specified and displayed.

一実施形態に係る情報処理装置の構成を概略的に示す図である。1 is a diagram schematically showing the configuration of an information processing device according to an embodiment; FIG. 情報処理装置の機能ブロック図である。3 is a functional block diagram of an information processing device; FIG. 図3(a)は、ケーブル長さの製造誤差の確率分布を示す図であり、図3(b)は、ケーブル剛性の製造誤差の確率分布を示す図であり、図3(c)は、ケーブル重量の製造誤差の確率分布を示す図である。FIG. 3(a) is a diagram showing the probability distribution of the cable length manufacturing error, FIG. 3(b) is a diagram showing the probability distribution of the cable stiffness manufacturing error, and FIG. FIG. 4 is a diagram showing the probability distribution of manufacturing errors in cable weight; 情報処理装置の処理を示すフローチャートである。4 is a flowchart showing processing of an information processing device; 図5(a)~図5(d)は、2台のロボットがケーブルを電子機器に敷設する手順を示す図である。FIGS. 5(a) to 5(d) are diagrams showing a procedure for two robots to lay a cable to an electronic device. 図6(a)は、把持について説明するための図であり、図6(b)は、固定について説明するための図である。FIG. 6(a) is a diagram for explaining gripping, and FIG. 6(b) is a diagram for explaining fixing. 図4のステップS16の詳細処理を示すフローチャートである。FIG. 5 is a flow chart showing detailed processing of step S16 in FIG. 4 ; FIG. 計算区間を示す図である。It is a figure which shows a calculation area. 図9(a)~図9(c)は、拘束点の長さ誤差の確率分布の計算方法について説明するための図である。FIGS. 9(a) to 9(c) are diagrams for explaining a method of calculating the probability distribution of length errors of constraint points. 図10(a)、図10(b)は、ケーブルの空間的な尤度分布の計算方法を示す図(その1)である。FIGS. 10(a) and 10(b) are diagrams (part 1) showing a method of calculating the spatial likelihood distribution of cables. 図11(a)、図11(b)は、ケーブルの空間的な尤度分布の計算方法を示す図(その2)である。FIGS. 11(a) and 11(b) are diagrams (part 2) showing a method of calculating a spatial likelihood distribution of cables. 図12(a)、図12(b)は、ケーブルの空間的な尤度分布の計算方法を示す図(その3)である。FIGS. 12(a) and 12(b) are diagrams (part 3) showing a method of calculating the spatial likelihood distribution of cables. 図13(a)~図13(c)は、ケーブルが存在する可能性がある範囲を表示する例について説明するための図である。FIGS. 13(a) to 13(c) are diagrams for explaining an example of displaying a range where a cable may exist. 図14(a)、図14(b)は、工法の検討について説明するための図であり、図14(c)は、作業失敗リスクの把握について説明するための図である。FIGS. 14(a) and 14(b) are diagrams for explaining the examination of the construction method, and FIG. 14(c) is a diagram for explaining the understanding of the work failure risk. 図15(a)~図15(c)は、干渉チェックについて説明するための図である。FIGS. 15(a) to 15(c) are diagrams for explaining the interference check.

以下、情報処理装置の一実施形態について、図1~図15に基づいて詳細に説明する。 An embodiment of an information processing apparatus will be described in detail below with reference to FIGS. 1 to 15. FIG.

図1には、一実施形態に係る情報処理装置10のハードウェア構成が示されている。図1の情報処理装置10は、例えば、ロボットが、ケーブル、ワイヤ、糸、ホース、ロープ、チューブ、カテーテルなどの線状柔軟物(ここではケーブルとする)を対象(ここではPC(Personal Computer)などの電子機器とする)に敷設する間における、ケーブルが存在する可能性のある範囲をシミュレーションする装置である。情報処理装置10は、シミュレーションにより得られたケーブルの存在可能性のある範囲を表示する。 FIG. 1 shows the hardware configuration of an information processing device 10 according to one embodiment. In the information processing apparatus 10 of FIG. 1, for example, a robot targets a linear flexible object (here, a cable) such as a cable, wire, thread, hose, rope, tube, catheter (here, a PC (Personal Computer)). It is a device that simulates the range where a cable may exist while laying it in an electronic device such as a cable. The information processing apparatus 10 displays a possible cable existence range obtained by the simulation.

情報処理装置10は、図1に示すように、CPU(Central Processing Unit)90、ROM(Read Only Memory)92、RAM(Random Access Memory)94、記憶部(ここではHDD(Hard Disk Drive))96、ネットワークインタフェース97、表示部93、入力部95、及び可搬型記憶媒体用ドライブ99等を備えている。これら情報処理装置10の構成各部は、バス98に接続されている。表示部93は、液晶ディスプレイ等を含み、入力部95は、キーボードやマウス、タッチパネル等を含む。情報処理装置10では、ROM92あるいはHDD96に格納されているプログラム(表示プログラムを含む)、或いは可搬型記憶媒体用ドライブ99が可搬型記憶媒体91から読み取ったプログラム(表示プログラムを含む)をCPU90が実行することにより、図2に示す各部の機能が実現される。なお、図2の各部の機能は、例えば、ASIC(Application Specific Integrated Circuit)やFPGA(Field Programmable Gate Array)等の集積回路により実現されてもよい。 As shown in FIG. 1, the information processing apparatus 10 includes a CPU (Central Processing Unit) 90, a ROM (Read Only Memory) 92, a RAM (Random Access Memory) 94, a storage unit (here, a HDD (Hard Disk Drive)) 96. , a network interface 97, a display unit 93, an input unit 95, a portable storage medium drive 99, and the like. Each component of the information processing apparatus 10 is connected to the bus 98 . The display unit 93 includes a liquid crystal display and the like, and the input unit 95 includes a keyboard, mouse, touch panel and the like. In the information processing apparatus 10, the CPU 90 executes a program (including a display program) stored in the ROM 92 or HDD 96, or a program (including a display program) read from the portable storage medium 91 by the portable storage medium drive 99. By doing so, the function of each part shown in FIG. 2 is realized. Note that the function of each unit in FIG. 2 may be implemented by an integrated circuit such as an ASIC (Application Specific Integrated Circuit) or an FPGA (Field Programmable Gate Array).

図2には、情報処理装置10の機能ブロック図が示されている。図2に示すように、情報処理装置10のCPU90は、プログラムを実行することにより、入力受信部12、シミュレーション実行部14、処理部16、画像生成部18として機能する。なお、図2には、HDD96等に格納されている各種DB30、32、34、36についても、図示されている。 FIG. 2 shows a functional block diagram of the information processing device 10. As shown in FIG. As shown in FIG. 2, the CPU 90 of the information processing apparatus 10 functions as the input receiving section 12, the simulation executing section 14, the processing section 16, and the image generating section 18 by executing programs. FIG. 2 also shows various DBs 30, 32, 34, and 36 stored in the HDD 96 or the like.

入力受信部12は、ユーザが入力部95を介して入力した情報を受信し、受信した情報を作業動作情報DB30や物性値DB32に格納する。また、入力受信部12は、ユーザが入力したシミュレーション開始指示の情報を受信すると、シミュレーション実行部14にその旨を通知する。ここで、作業動作情報DB30には、各時刻における、ロボットハンドの移動軌跡や姿勢の情報、動作速度、作業の情報(把持や固定などの情報)が格納されている。また、物性値DB32には、ケーブルの長さ、ケーブルの剛性、ケーブルの重量などの情報が格納されている。 The input reception unit 12 receives information input by the user via the input unit 95 and stores the received information in the work movement information DB 30 and the physical property value DB 32 . Further, when the input reception unit 12 receives the information of the simulation start instruction input by the user, the input reception unit 12 notifies the simulation execution unit 14 to that effect. Here, the work motion information DB 30 stores information on the movement trajectory and posture of the robot hand, motion speed, and work information (information on gripping, fixing, etc.) at each time. In addition, the physical property value DB 32 stores information such as cable length, cable rigidity, and cable weight.

シミュレーション実行部14は、入力受信部12からシミュレーション開始指示を受信すると、処理部16に対してシミュレーションの時刻情報(シミュレーション開始から何秒後かを示す情報)を通知することで、処理部16にシミュレーションを実行させる。 When the simulation execution unit 14 receives the simulation start instruction from the input reception unit 12, the simulation execution unit 14 notifies the processing unit 16 of simulation time information (information indicating how many seconds after the start of the simulation). Run the simulation.

処理部16は、シミュレーション実行部14から受信した時刻情報に基づいて、受信した時刻におけるケーブルの存在可能性のある範囲をシミュレーションする。ここで、処理部16は、図2に示すように、設定部としての制約計算部22、算出部としてのケーブル形状計算部24、尤度分布計算部26を有する。 Based on the time information received from the simulation execution unit 14, the processing unit 16 simulates the possible cable existence range at the received time. Here, as shown in FIG. 2 , the processing unit 16 has a constraint calculation unit 22 as a setting unit, a cable shape calculation unit 24 as a calculation unit, and a likelihood distribution calculation unit 26 .

制約計算部22は、シミュレーション実行部14から受信した時刻におけるケーブルの制約を計算する。ケーブルの制約とは、例えば、ケーブルの一部が拘束されたり、拘束が解除されることを意味する。ケーブルの拘束には、ケーブルの一部が電子機器に固定されることや、ケーブルの一部がロボットハンドにより保持されることを含んでいる。 The constraint calculation unit 22 calculates the cable constraint at the time received from the simulation execution unit 14 . Restriction of a cable means, for example, that a part of the cable is restricted or unrestricted. Cable restraint includes a portion of the cable being fixed to the electronic device and a portion of the cable being held by the robot hand.

ケーブル形状計算部24は、ケーブルの拘束点がケーブルの長さ方向に位置ずれすることによる拘束点間のケーブルの長さの変動を考慮して、拘束点間のケーブルがとりうる形状および該形状の尤度(尤もらしさを示す指標)を計算する。具体的には、ケーブル形状計算部24は、制約計算部22が計算したケーブルの制約に基づいて、物性値DB32、製造誤差に基づく確率分布DB34、位置ずれの確率分布DB36に格納されている情報から、ケーブルが取りうる形状(形状サンプル)とその尤度を計算する。ここで、製造誤差に基づく確率分布DB34には、ケーブル長さの製造誤差の確率分布(図3(a))や、ケーブル剛性の製造誤差の確率分布(図3(b))、ケーブル重量の製造誤差の確率分布(図3(c))の各情報が格納されている。また、位置ずれの確率分布DB36には、シミュレーション中のケーブルの拘束点(ロボットハンドに把持されたり、電子機器に固定された点)それぞれのケーブルの長さ方向に関する位置ずれの確率分布(図9(c)参照)が格納される。 The cable shape calculation unit 24 considers the variation in the length of the cable between the restraint points due to the displacement of the restraint points of the cable in the length direction of the cable, and determines the shape that the cable can take between the restraint points and the shape. Calculate the likelihood of Specifically, the cable shape calculation unit 24 calculates information stored in the physical property value DB 32, the probability distribution DB 34 based on manufacturing errors, and the probability distribution DB 36 of positional deviation based on the cable constraints calculated by the constraint calculation unit 22. Then, the shape (shape sample) that the cable can take and its likelihood are calculated. Here, the probability distribution DB 34 based on the manufacturing error contains the probability distribution of the manufacturing error of the cable length (Fig. 3(a)), the probability distribution of the manufacturing error of the cable stiffness (Fig. 3(b)), and the weight of the cable. Each piece of information on the probability distribution of manufacturing errors (FIG. 3(c)) is stored. Further, the probability distribution DB 36 of positional deviation stores the probability distribution of positional deviation (Fig. 9 (c)) is stored.

尤度分布計算部26は、ケーブル形状計算部24が計算したケーブルが取りうる形状(形状サンプル)とその尤度とに基づいて、尤度の空間的な分布を計算する。尤度分布計算部26は、計算した尤度分布をシミュレーション実行部14を介して画像生成部18に送信する。 The likelihood distribution calculator 26 calculates a spatial distribution of likelihoods based on the shapes (shape samples) that the cable can take calculated by the cable shape calculator 24 and the likelihoods thereof. The likelihood distribution calculator 26 transmits the calculated likelihood distribution to the image generator 18 via the simulation execution unit 14 .

画像生成部18は、尤度分布計算部26で計算された尤度の空間的な分布を表す画像を生成する。画像生成部18は、生成した画像を表示部93に表示させる。 The image generator 18 generates an image representing the spatial distribution of likelihoods calculated by the likelihood distribution calculator 26 . The image generation unit 18 causes the display unit 93 to display the generated image.

(情報処理装置10の処理について)
次に、情報処理装置10の処理について、図4のフローチャートに沿って、その他図面を適宜参照しつつ詳細に説明する。
(Processing of information processing device 10)
Next, the processing of the information processing apparatus 10 will be described in detail along the flowchart of FIG. 4 and with appropriate reference to other drawings.

ここで、本実施形態では、図5(a)~図5(d)に示すような手順で、2台のロボットを用いて、ケーブル50を電子機器100に敷設する(組み付ける)ものとする。具体的には、図5(a)に示すように、第1ロボットが、ケーブル50の一方のコネクタをコネクタ接続点Aに接続する。次いで、図5(b)に示すように、第1ロボットがコネクタ接続点Aから離れる方向に移動し、点Bにおいて第2ロボットがケーブル50の一部を電子機器100に固定する。次いで、図5(c)に示すように、第1ロボットがコネクタ接続点Aから離れる方向に更に移動し、点Cにおいて第2ロボットがケーブル50の一部を電子機器100に固定する。そして、図5(d)に示すように、第1ロボットが、ケーブル50の他方のコネクタをコネクタ接続点Dに接続する。情報処理装置10は、図4の処理を実行することで、第1ロボット及び第2ロボットによる図5(a)~図5(d)に示す作業の間におけるケーブル50の存在可能性のある範囲をシミュレーションし、その結果を表示部93に表示する。 Here, in this embodiment, two robots are used to lay (assemble) the cable 50 to the electronic device 100 in the procedure shown in FIGS. 5(a) to 5(d). Specifically, the first robot connects one connector of the cable 50 to the connector connection point A, as shown in FIG. 5(a). Next, as shown in FIG. 5B, the first robot moves away from the connector connecting point A, and the second robot fixes part of the cable 50 to the electronic device 100 at the point B. As shown in FIG. Next, as shown in FIG. 5C, the first robot moves further away from the connector connection point A, and the second robot fixes a portion of the cable 50 to the electronic device 100 at point C. Then, as shown in FIG. Then, the first robot connects the other connector of the cable 50 to the connector connection point D, as shown in FIG. 5(d). The information processing apparatus 10 performs the process of FIG. 4 to determine the possible existence range of the cable 50 during the work shown in FIGS. 5A to 5D by the first robot and the second robot. is simulated, and the result is displayed on the display unit 93 .

図4の処理では、まず、ステップS10において、シミュレーション実行部14が、時刻tを0に初期化する。 In the process of FIG. 4, the simulation execution unit 14 first initializes the time t to 0 in step S10.

次いで、ステップS12では、シミュレーション実行部14が、時刻tをδt進める(t←t+δt)。 Next, in step S12, the simulation execution unit 14 advances the time t by δt (t←t+δt).

次いで、ステップS14では、制約計算部22が、作業動作情報DB30を参照して、制約条件を発生させたり、消滅させたりする。この処理では、制約計算部22は、2台のロボットの動作として、作業動作情報DB30に格納されている以下の作業動作情報を得る。
(作業動作情報1)
作業時のロボットハンド先端の移動軌跡と各位置でのロボットの姿勢の時系列情報。
(作業動作情報2)
時刻δtにおけるケーブルの拘束状態(把持、固定)の変化情報(ON/OFF)。
Next, in step S14, the constraint calculation unit 22 refers to the work motion information DB 30 to generate or eliminate constraint conditions. In this process, the constraint calculation unit 22 obtains the following work motion information stored in the work motion information DB 30 as motions of the two robots.
(Work movement information 1)
Time-series information of the movement trajectory of the tip of the robot hand during work and the posture of the robot at each position.
(Work motion information 2)
Change information (ON/OFF) of the restraint state (holding, fixing) of the cable at time δt.

ここで、拘束状態には、把持と固定がある。把持は、ケーブルをロボットハンドが持っている状態を表現するものであり、図6(a)に示すように、ロボットハンド先端の位置と時刻δtにおいて当該位置と一致しているケーブルの位置とを、以後の時刻でも一致させることを意味する。一方、固定は、組み付けや一時的な固定(テープ等による仮止め)を表現するものである。時刻δtにおいて図6(b)に示すようにケーブルが固定された場合には、固定位置(ロボットで押さえる場合にはロボットハンドの先端位置、ピン留めであればピン位置)にあるケーブル上の点を以後の時刻で固定位置と一致させることを意味する。以降の説明では、上述したような拘束されたケーブル上の点を「拘束点」と呼ぶものとする。 Here, the restrained state includes grasping and fixing. Grasping represents the state in which the robot hand holds the cable. As shown in FIG. , means to match at later times as well. On the other hand, fixation expresses assembly or temporary fixation (temporary fixation with tape or the like). When the cable is fixed at time δt as shown in FIG. 6(b), the point on the cable at the fixed position (the tip position of the robot hand if held by the robot, or the pin position if pinned) to match the fixed position at a later time. In the following description, the points on the cable that are constrained as described above will be referred to as "constraint points."

次いで、ステップS16では、ケーブル形状計算部24が、ケーブルの形状を計算する処理を実行する。このステップS16では、図7のフローチャートに沿った処理を実行する。ここでは、一例として、図5(b)の状態から図5(c)の状態に遷移する間を例にとり、図7の処理(ステップS16の処理)について説明する。なお、図7の処理では、誤差のないケーブルの形状はもちろん、誤差を考慮した場合(誤差がある場合)のケーブルの形状(形状サンプル)についても計算する。 Next, in step S16, the cable shape calculator 24 executes a process of calculating the shape of the cable. In this step S16, the process according to the flowchart of FIG. 7 is executed. Here, as an example, the process of FIG. 7 (the process of step S16) will be described while transitioning from the state of FIG. 5B to the state of FIG. 5C. In the process of FIG. 7, not only the shape of the cable without error but also the shape (shape sample) of the cable when considering the error (when there is an error) is calculated.

図7の処理では、まず、ステップS50において、ケーブル形状計算部24が、両端の拘束点のケーブルの長さ方向に関する位置ずれとケーブル長さに関する製造誤差の確率分布に従い、計算区間における長さ誤差と確率密度を計算する。ここで、図5(b)の状態から図5(c)の状態に遷移する間における「計算区間」とは、図8に示す把持拘束点を挟む前後の区間(計算区間1、2)である。 In the process of FIG. 7, first, in step S50, the cable shape calculation unit 24 calculates the length error in the calculation section according to the probability distribution of the positional deviation of the restraint points at both ends of the cable in the length direction of the cable and the manufacturing error of the cable length. and calculate the probability density. 5(b) to FIG. 5(c) is the interval before and after the grip constraint point shown in FIG. 8 (calculation intervals 1 and 2). be.

(計算区間1について)
計算区間1における長さ誤差と確率密度を求める場合、点Bにおける作業によるケーブルの長さ方向に関する位置ずれの確率分布と、把持拘束点における作業によるケーブルの長さ方向に関する位置ずれの確率分布と、長さに関する製造誤差の確率分布とを用いる。
(Regarding calculation interval 1)
When calculating the length error and the probability density in the calculation interval 1, the probability distribution of the positional deviation in the length direction of the cable due to the work at the point B and the probability distribution of the positional deviation in the length direction of the cable due to the work at the grip constraint point , and the probability distribution of the manufacturing error on the length.

具体的には、計算区間1における長さ誤差T1は、点Bにおける作業によるケーブルの長さ方向に関する位置ずれをTB、把持拘束点における作業によるケーブルの長さ方向に関する位置ずれをTG、長さに関する製造誤差をTMとすると、次式(1)から求められる。なお、L1は、計算区間1の長さ、Lはケーブルの全長を意味する。
1=TB+TG+(L1/L)×TM …(1)
Specifically, the length error T 1 in the calculation interval 1 is defined as T B , the displacement in the length direction of the cable due to work at point B, and T G , the displacement in the length direction of the cable due to work at the gripping constraint point. , and the manufacturing error related to the length is T M , it is obtained from the following equation (1). In addition, L1 means the length of the calculation section 1, and L means the total length of the cable.
T1 = TB + TG +(L1/L)× TM (1)

また、長さ誤差T1の確率密度関数P(T1)は、TB、TG、TMそれぞれの確率分布を表す関数の合成積である。ここで、長さ誤差T1の確率密度関数P(T1)は、誤差の確率分布を連続値の関数として表現したものであり、この確率密度関数が、計算区間1の長さ誤差の確率分布にあたる。なお、長さに関する製造誤差TMの確率分布は、図3(a)に示すものであり、製造誤差に基づく確率分布DB34(図2参照)に格納されている。なお、拘束点のケーブルの長さ方向に関する位置ずれの確率分布は、位置ずれの確率分布DB36(図2参照)に格納されているが、その計算方法については後述する。 Also, the probability density function P(T 1 ) of the length error T 1 is a composite product of functions representing the respective probability distributions of T B , T G , and T M . Here, the probability density function P(T 1 ) of the length error T 1 expresses the probability distribution of the error as a function of continuous values. distribution. Note that the probability distribution of the manufacturing error T M related to the length is shown in FIG. 3(a), and is stored in the probability distribution DB 34 (see FIG. 2) based on the manufacturing error. The probability distribution of the positional displacement of the restraint point in the length direction of the cable is stored in the positional displacement probability distribution DB 36 (see FIG. 2), and the calculation method thereof will be described later.

(計算区間2について)
計算区間2における長さ誤差と確率分布を求める場合、把持拘束点における作業によるケーブルの長さ方向に関する位置ずれの確率分布と、未固定の端点における作業によるケーブルの長さ方向に関する位置ずれの確率分布と、長さに関する製造誤差の確率分布とを用いる。
(Regarding calculation interval 2)
When obtaining the length error and probability distribution in calculation interval 2, the probability distribution of positional displacement in the length direction of the cable due to work at the grip restraint point and the probability of positional displacement in the length direction of the cable due to work at the unfixed end point distribution and the probability distribution of the manufacturing error with respect to length.

具体的には、計算区間2における長さ誤差T2は、把持拘束点における作業によるケーブルの長さ方向に関する位置ずれをTG、未固定の端点における作業によるケーブルの長さ方向に関する位置ずれをTend(ただし、Tend=0)、長さに関する製造誤差をTMとすると、次式(2)から求められる。なお、L2は、計算区間2の長さを意味する。
2=Tend+TG+(L2/L)×TM …(2)
Specifically, the length error T 2 in the calculation interval 2 is the positional deviation in the length direction of the cable due to work at the grip restraint point, and the positional deviation in the length direction of the cable due to work at the unfixed end point. T end (where T end =0), and T M being the manufacturing error related to the length, can be obtained from the following equation (2). Note that L2 means the length of calculation interval 2 .
T 2 =T end +T G +(L2/L)×T M (2)

また、長さ誤差T2の確率密度関数P(T2)は、Tend、TG、TMそれぞれの確率分布を表す関数の合成積であり、この確率密度関数が、計算区間2の長さ誤差の確率分布にあたる。 Moreover, the probability density function P(T 2 ) of the length error T 2 is a composite product of functions representing the probability distributions of T end , T G , and T M , and this probability density function is the length of the calculation interval 2 It corresponds to the probability distribution of the error.

図7に戻り、次のステップS52では、ケーブル形状計算部24が、剛性に関する製造誤差の確率分布に従い、剛性の誤差と確率密度を計算する。例えば、ケーブル形状計算部24は、図3(b)より、計算区間1の剛性の誤差G1とその確率密度PG1を計算する。 Returning to FIG. 7, in the next step S52, the cable shape calculator 24 calculates stiffness errors and probability densities according to the probability distribution of manufacturing errors relating to stiffness. For example, the cable shape calculator 24 calculates the stiffness error G1 of the calculation section 1 and its probability density PG1 from FIG. 3(b).

次いで、ステップS54では、ケーブル形状計算部24が、重量に関する製造誤差の確率分布に従い、計算区間における重量の誤差と確率密度を計算する。例えば、ケーブル形状計算部24は、図3(c)より、計算区間1における重量の誤差M1と、その確率密度PM1を計算する。 Next, in step S54, the cable shape calculation unit 24 calculates the weight error and the probability density in the calculation interval according to the probability distribution of the weight-related manufacturing error. For example, the cable shape calculator 24 calculates the weight error M1 in the calculation section 1 and its probability density PM1 from FIG. 3(c).

次いで、ステップS56では、ケーブル形状計算部24が、長さ、剛性、重量を更新し、各値の確率密度を乗算して尤度を計算する。ステップS50で計算した計算区間1における長さ誤差をL1、この長さ誤差L1の確率密度をPL1とすると、尤度Xは次式(3)より求めることができる。
X=PL1×PG1×PM1 …(3)
Next, in step S56, the cable shape calculator 24 updates the length, stiffness, and weight, multiplies the probability density of each value, and calculates the likelihood. Assuming that the length error in calculation interval 1 calculated in step S50 is L1 and the probability density of this length error L1 is PL1, the likelihood X can be obtained from the following equation (3).
X=PL1×PG1×PM1 (3)

なお、計算区間2についても同様にして尤度を求めることができる。 Note that the likelihood for calculation interval 2 can also be obtained in the same manner.

次いで、ステップS58では、ケーブル形状計算部24が、更新された長さ、剛性、重量で形状サンプルの形状を計算し、尤度を添えて保存する。 Next, in step S58, the cable shape calculator 24 calculates the shape of the shape sample with the updated length, stiffness, and weight, and stores it with the likelihood.

ここで、各区間の形状は、ケーブルをパラメトリック曲線で近似し、拘束点の制約を満たす条件のもとでポテンシャルエネルギを最小化する曲線形状を求めることで計算することができる。なお、形状の計算方法については、「ワイヤ・ハーネスの変形シミュレーションと設計支援システム,情報処理学会論文誌 48(2),909-917,2007-02-15」に記載されている方法などを採用することができる。 Here, the shape of each section can be calculated by approximating the cable with a parametric curve and determining the curve shape that minimizes the potential energy under the condition that the constraints of the constraint points are satisfied. Regarding the calculation method of the shape, the method described in "Deformation simulation of wire harness and design support system, Information Processing Society of Japan 48(2), 909-917, 2007-02-15" is adopted. can do.

次いで、ステップS60では、ケーブル形状計算部24が、回数が上限に達したか否か、すなわち必要数だけ形状サンプルが得られたか否かを判断する。このステップS60の判断が否定された場合には、ステップS50に戻り、上述した処理を繰り返す。一方、ステップS60の判断が肯定された場合には、図7の全処理(ステップS16)を終了し、図4のステップS20に移行する。なお、ステップS20に移行するときには、誤差のないケーブルの形状及びその尤度と、誤差がある場合のケーブルの形状(形状サンプル)及びその尤度が保存された状態となっている。 Next, in step S60, the cable shape calculator 24 determines whether or not the number of times has reached the upper limit, that is, whether or not the necessary number of shape samples have been obtained. If the determination in step S60 is negative, the process returns to step S50 and repeats the above-described processing. On the other hand, if the determination in step S60 is affirmative, the entire process (step S16) in FIG. 7 is terminated, and the process proceeds to step S20 in FIG. Note that when the process proceeds to step S20, the shape of the cable without error and its likelihood, and the shape of the cable with an error (shape sample) and its likelihood are stored.

ここで、ケーブル形状計算部24は、拘束点において固定や把持が発生するたびに、拘束点のケーブルの長さ方向に関する位置ずれの確率分布を計算し、位置ずれの確率分布DB36に保存しておく。以下、図9(a)~図9(c)に基づいて、拘束点のケーブルの長さ方向に関する位置ずれの確率分布の計算方法について説明する。 Here, the cable shape calculation unit 24 calculates the probability distribution of positional deviation in the length direction of the cable at the constraint point each time fixation or gripping occurs at the constraint point, and stores the probability distribution in the positional deviation probability distribution DB 36. back. A method of calculating the probability distribution of the positional deviation of the restraint point in the length direction of the cable will be described below with reference to FIGS. 9(a) to 9(c).

(a)ケーブル形状計算部24は、誤差が無いケーブルにおける拘束点をZとし、拘束点Zに対応する各形状サンプル上の各点を抽出する。図9(a)においては、太実線が誤差が無いケーブルの形状を示しており、その他の線種が誤差がある形状サンプルを示している。また、形状サンプル上において示す点が、拘束点Zに対応する点(始点から同じ位置にある点)を示している。 (a) The cable shape calculator 24 sets Z as a constraint point on a cable with no error, and extracts each point on each shape sample corresponding to the constraint point Z. FIG. In FIG. 9A, the thick solid line indicates the shape of the cable with no error, and the other line types indicate shape samples with errors. Also, the points shown on the shape sample indicate points corresponding to the constraint point Z (points at the same position from the starting point).

(b)次いで、ケーブル形状計算部24は、拘束点Zに対応する各点の位置ずれ量Eを以下の方法で計算する。
(i) 図9(b)に示すように、拘束点Zと各点を結ぶ移動ベクトルvと拘束点Zにおける接線とを特定して内積を求め、その値が0以下であれば位置ずれ量Eの符号を負とし、それ以外であれば正とする。
(ii) 位置ずれ量の絶対値は、移動ベクトルvの長さとする。
(b) Next, the cable shape calculator 24 calculates the positional deviation amount E of each point corresponding to the constraint point Z by the following method.
(i) As shown in FIG. 9(b), determine the inner product of the movement vector v connecting the constraint point Z and each point and the tangent line to the constraint point Z, and if the value is 0 or less, the displacement amount Let the sign of E be negative, otherwise positive.
(ii) The absolute value of the displacement amount is the length of the movement vector v.

(c)次いで、ケーブル形状計算部24は、拘束点Zに対応する点における尤度Xを形状サンプルそれぞれの尤度と同値とする。 (c) Next, the cable shape calculator 24 sets the likelihood X at the point corresponding to the constraint point Z to the same value as the likelihood of each shape sample.

(d)次いで、ケーブル形状計算部24は、拘束点Zに対応する各点について、x軸(横軸)を位置ずれ量E、y軸(縦軸)を尤度としてプロットし、図9(c)に示すような分布関数F(E)を計算する。この分布関数F(E)を適宜スケールした分布を作業により発生する拘束点でのケーブルの長さ方向に関する位置ずれの確率分布とする。 (d) Next, the cable shape calculation unit 24 plots the x-axis (horizontal axis) as the positional deviation amount E and the y-axis (vertical axis) as the likelihood for each point corresponding to the restraint point Z. Calculate the distribution function F(E) as shown in c). The distribution obtained by appropriately scaling this distribution function F(E) is used as the probability distribution of the positional deviation in the length direction of the cable at the restraint point generated by the work.

例えば、次式(4)のように、正規分布を定数(k)倍して、関数F(E)とすることができる。 For example, as shown in the following equation (4), the normal distribution can be multiplied by a constant (k) to obtain a function F(E).

Figure 0007180274000001
Figure 0007180274000001

この場合、拘束点Zと対応する点との距離の総和が最小となる分散σと定数kを求めることにより、分布関数F(E)を決定することができる。 In this case, the distribution function F(E) can be determined by obtaining the variance σ and the constant k that minimize the sum of the distances between the constraint point Z and the corresponding points.

そして、F(E)は尤度の関数であるため、F(E)を全区間で積分したときに1となるようにスケールしたP(E)(=F(E)/k)を位置ずれ量Eの確率分布とする。 Since F(E) is a function of likelihood, P(E) (=F(E)/k) scaled to be 1 when F(E) is integrated over the entire interval is used as the displacement Let the probability distribution of quantity E be

なお、上記のようにして計算された拘束点のケーブルの長さ方向に関する位置ずれの確率分布は、位置ずれの確率分布DB36に保存される。一方、拘束状態の制約が解除された場合(OFFになった場合)には、その拘束点の長さ方向に関する位置ずれの確率分布を0に更新し、位置ずれの確率分布DB36に保存する。 It should be noted that the probability distribution of the positional displacement of the restraint point in the length direction of the cable calculated as described above is stored in the positional displacement probability distribution DB 36 . On the other hand, when the constraint in the constraint state is released (turned OFF), the positional deviation probability distribution in the length direction of the constraint point is updated to 0 and stored in the positional deviation probability distribution DB 36 .

図4に戻り、次のステップS20では、尤度分布計算部26が、ケーブルの空間的な尤度分布を計算する。以下、ケーブルの空間的な尤度分布の計算方法について、図10~図13に基づいて詳細に説明する。 Returning to FIG. 4, in the next step S20, the likelihood distribution calculator 26 calculates the spatial likelihood distribution of the cable. A method of calculating the spatial likelihood distribution of cables will be described in detail below with reference to FIGS. 10 to 13. FIG.

(ケーブルの空間的な尤度分布の計算方法)
まず、尤度分布計算部26は、図10(a)に示すような形状サンプル、誤差なしの形状の中から、所定の尤度X以上の形状サンプル(誤差なしの形状を含む)を抽出する。そして、尤度分布計算部26は、図10(b)に示すように、抽出した形状サンプルをすべて含む閉曲面S(X)(多面体メッシュ)を計算する。このとき、閉曲面S(X)上の点(黒丸で示す点)は尤度Xとして扱う。
(Calculation method of spatial likelihood distribution of cables)
First, the likelihood distribution calculator 26 extracts shape samples (including shapes without errors) having a predetermined likelihood X or more from shape samples without errors shown in FIG. 10(a). . Then, the likelihood distribution calculator 26 calculates a closed curved surface S(X) (polyhedral mesh) including all the extracted shape samples, as shown in FIG. 10(b). At this time, points on the closed surface S(X) (points indicated by black circles) are treated as likelihood X.

図11(a)には、閉曲面S(X)に回転楕円体で規定される領域を重ね合わせた状態が示されている。この回転楕円体で既定される領域は、ケーブルが届く最大の領域である。すなわち、計算区間の両端点が固定された一定の長さの曲線が届く範囲が回転楕円体で規定される領域となっている。この回転楕円体の面上の点(白丸で示す点)は尤度0となる。 FIG. 11(a) shows a state in which a region defined by a spheroid is superimposed on the closed curved surface S(X). The area defined by this spheroid is the maximum area that the cable can reach. That is, the range reached by a curve of a certain length with fixed endpoints of the calculation interval is an area defined by the spheroid. Points on the surface of this spheroid (points indicated by white circles) have a likelihood of zero.

尤度分布計算部26は、回転楕円体内部にある任意の点における尤度を、誤差なしの形状、閉曲面S(X)、回転楕円体それぞれの表面の点における尤度を空間的に補間することで計算する。なお、任意の点は、後述する画像生成部18で画像を生成する(可視化する)際に必要となる解像度に対応する数の点を意味する。 The likelihood distribution calculator 26 spatially interpolates the likelihood at any point inside the spheroid, the likelihood at points on the shape without error, the closed surface S(X), and the surface of the spheroid. Calculate by doing Note that the arbitrary points mean the number of points corresponding to the resolution required when an image is generated (visualized) by the image generation unit 18, which will be described later.

ここで、任意の点における尤度の補間計算について詳細に説明する。本実施形態では、回転楕円体、閉曲面の各頂点の位置に閉曲面の尤度を持ったスカラー場を考える。この頂点群を3次元Delaunay三角形分割すると、回転楕円体の内部が複数の4面体の集合で表現される。この場合、任意の点pが図11(a)において白丸〇で示すような回転楕円体の表面(又は外側)にある場合には尤度は0となる。一方、任意の点pが回転楕円体の内部にある場合には、図11(b)に示すように、点pを内部に含む4面体を探索し、点pと4面体の各頂点を結ぶことで4つの4面体に分割する。そして、分割された各4面体の体積をV0、V1、V2、V3、元の4面体の体積をVとし、次式(5)から、点pの尤度Xpを求める。 The interpolation calculation of the likelihood at arbitrary points will now be described in detail. In this embodiment, a scalar field having the likelihood of a closed surface at each vertex position of a spheroid and a closed surface is considered. When this vertex group is subjected to three-dimensional Delaunay triangulation, the inside of the spheroid is represented by a set of tetrahedrons. In this case, the likelihood is 0 when an arbitrary point p is on the surface (or outside) of the spheroid indicated by the white circle ◯ in FIG. 11(a). On the other hand, if an arbitrary point p is inside the spheroid, as shown in FIG. divide it into four tetrahedrons. Then, let the volumes of the divided tetrahedrons be V 0 , V 1 , V 2 , and V 3 and the volume of the original tetrahedron be V, and obtain the likelihood X p of the point p from the following equation (5).

Figure 0007180274000002
Figure 0007180274000002

なお、尤度分布計算部26は、閉曲面S(X)を2つ以上計算することができる。例えば、尤度分布計算部26は、図12(a)に示すように、尤度がXa以上の形状サンプルをすべて含む閉曲面(空間)S(Xa)と、尤度がXb(>Xa)以上の形状サンプルをすべて含む閉曲面S(Xb)を計算する。この場合、閉曲面S(Xb)は閉曲面S(Xa)の内部に含まれている。このように複数の閉曲面を計算した場合にも、図12(b)に示すような回転楕円体上の尤度(=0)と複数の閉曲面それぞれの尤度に基づいて上式(5)から任意の点の尤度を求めるようにすればよい。 The likelihood distribution calculator 26 can calculate two or more closed surfaces S(X). For example, as shown in FIG. 12(a), the likelihood distribution calculator 26 calculates a closed surface (space) S(Xa) containing all shape samples with a likelihood of Xa or more and a surface (space) S(Xa) with a likelihood of Xb (>Xa). A closed curved surface S(Xb) including all the above shape samples is calculated. In this case, the closed curved surface S(Xb) is included inside the closed curved surface S(Xa). Even when a plurality of closed surfaces are calculated in this way, the above equation (5 ) to find the likelihood of an arbitrary point.

ここで、上記のようなケーブルの空間的な尤度の計算を行う際には、尤度分布計算部26は、画像生成に必要な位置の尤度を全て計算する。 Here, when calculating the spatial likelihood of the cable as described above, the likelihood distribution calculator 26 calculates all the likelihoods of positions necessary for image generation.

例えば、図13(a)に示すように複数断面でのコンター図としてケーブルの存在可能性のある範囲を表現する場合には、図13(b)に示すように各断面を格子状に分割した際にできる四角形それぞれの中心の位置の尤度を計算する。また、図13(c)に示すように、尤度が高いほど、透過率を低くした形状として表現する場合には、可視化領域を直方体で表し、その直方体を直交メッシュとして小さい長方形に分割した上で、小さい長方形の中心の位置の尤度を計算する。 For example, as shown in FIG. 13(a), when expressing a range in which a cable may exist as a contour diagram of multiple cross sections, each cross section is divided into a grid as shown in FIG. 13(b). Calculate the likelihood of the location of the center of each possible quadrangle. In addition, as shown in FIG. 13C, when representing a shape with a lower transmittance as the likelihood is higher, the visualization region is represented by a rectangular parallelepiped, and the rectangular parallelepiped is divided into small rectangles as an orthogonal mesh. , compute the likelihood of the location of the center of the small rectangle.

なお、尤度分布計算部26は、上述した処理により得られた情報を、シミュレーション実行部14を介して画像生成部18に送信する。 The likelihood distribution calculator 26 transmits the information obtained by the processing described above to the image generator 18 via the simulation execution unit 14 .

図4に戻り、次のステップS22では、画像生成部18が、表示用の画像を作成して表示部93上に表示する。この場合、画像生成部18は、尤度分布計算部26から受信した情報に基づいて、誤差なしのときのケーブルの形状に対して尤度分布を色付けして重畳した図13(a)や図13(c)に示すような画像を作成して表示部93上に表示する。 Returning to FIG. 4, in the next step S22, the image generation unit 18 creates an image for display and displays it on the display unit 93. FIG. In this case, based on the information received from the likelihood distribution calculation unit 26, the image generation unit 18 colors the likelihood distribution with respect to the shape of the cable when there is no error, and superimposes it in FIG. 13(c) is created and displayed on the display unit 93. FIG.

次いで、ステップS24では、シミュレーション実行部14が、終了か否かを判断する。すなわち、時刻tが予め定められている所定時刻(シミュレーション終了時刻)に到達したか否かを判断する。このステップS24の判断が否定された場合には、ステップS12に戻り、ステップS12~S24を実行することで、先ほどの時刻よりもδt秒後についてのシミュレーションを実行する。このようにステップS12~S24を繰り返し実行することにより、ケーブルが存在する可能性のある範囲を表示部93に表示することができる。一方、ステップS24の判断が肯定された場合には、図4の全処理が終了する。 Next, in step S24, the simulation execution unit 14 determines whether or not the simulation is finished. That is, it is determined whether or not the time t has reached a predetermined time (simulation end time). If the determination in step S24 is negative, the process returns to step S12, and steps S12 to S24 are executed to execute the simulation for δt seconds after the previous time. By repeatedly executing steps S12 to S24 in this manner, the range in which the cable may exist can be displayed on the display unit 93. FIG. On the other hand, if the determination in step S24 is affirmative, all the processing in FIG. 4 ends.

(利用シーン)
ここで、画像生成部18が表示部93に表示する画像の利用シーンについて説明する。
(Use scene)
Here, the use scene of the image displayed on the display unit 93 by the image generation unit 18 will be described.

(a)工法の検討
図14(a)、図14(b)は、同一の作業を行う2つの手順それぞれを採用した場合について、ケーブルの存在可能性のある範囲を可視化することで、作業の成功率の高さを見積もることができる。図14(a)の例では、1台のロボットハンドでケーブル50を敷設している。この場合、図14(a)の画像を見たユーザは、ケーブル50のたるみが大きくなり、ケーブル50の位置のばらつきも大きくなるためケーブル50を目標位置に適切に組み付けられない可能性があると判断することができる。一方、図14(b)の例では、2つのロボットハンドを用い、一方のロボットハンドがケーブル50を押さえて、他方のロボットハンドがケーブル50を張った状態で敷設している。この場合、図14(b)の画像を見たユーザは、ケーブル50の位置のばらつきが小さくなるため、目標位置に適切に組み付けることができると判断することができる。
(a) Examination of construction methods Figures 14(a) and 14(b) show how the work can be performed by visualizing the possible range of existence of cables when two procedures for performing the same work are adopted. High success rate can be estimated. In the example of FIG. 14(a), the cable 50 is laid by one robot hand. In this case, the user who saw the image of FIG. 14( a ) thinks that the slack of the cable 50 increases and the variation in the position of the cable 50 increases, so that the cable 50 may not be properly assembled at the target position. can judge. On the other hand, in the example of FIG. 14B, two robot hands are used, one of which holds down the cable 50 and the other of which stretches the cable 50 for laying. In this case, the user who sees the image in FIG. 14(b) can judge that the cable 50 can be appropriately assembled at the target position because the variation in the position of the cable 50 is small.

(b)作業失敗リスクの把握
図14(c)に示すようにケーブル50を敷設する場合、ケーブル50には一定量の余長が必要である。すなわち、ケーブルの存在可能性を示す範囲が狭いと、ケーブルの長さが不足するおそれがある。したがって、図14(c)のような画像を見たユーザは、ケーブル50の存在可能性を示す範囲が狭いため、長さ不足により作業が失敗する可能性があると判断することができる。
(b) Ascertaining work failure risk When laying the cable 50 as shown in FIG. That is, if the range indicating the possibility of existence of a cable is narrow, the length of the cable may be insufficient. Therefore, a user who sees the image as shown in FIG. 14(c) can determine that there is a possibility that the operation will fail due to the short length of the cable 50 because the range indicating the possibility of existence of the cable 50 is narrow.

(c)干渉チェック
図15(a)に示すように、ケーブルの存在可能性が高い範囲に周辺の部品52が存在している場合には、ユーザは、ケーブル50と部品52とが干渉する可能性があると判断することができる。
(c) Interference check As shown in FIG. 15(a), when there is a peripheral component 52 in the range where the possibility of cable existence is high, the user is aware that the cable 50 and the component 52 may interfere with each other. can be determined to be viable.

また、図15(b)に示すようにケーブルを溝54に収める必要がある場合には、ユーザは、ケーブル50の存在可能性が高い範囲と溝54との関係をチェックする。この場合、ユーザは、ケーブル50の存在可能性が高い範囲が溝54から多くはみ出ている場合に、ケーブル50が溝に収まらない可能性が高いと判断することができる。 15(b), when the cable needs to be accommodated in the groove 54, the user checks the relationship between the range where the cable 50 is likely to exist and the groove 54. FIG. In this case, the user can determine that there is a high possibility that the cable 50 will not fit in the groove when the range where the cable 50 is likely to exist protrudes largely from the groove 54 .

また、図15(c)に示すように、ケーブル50をフック56内に収める必要がある場合も、ユーザは、図15(b)と同様に判断すればよい。 Also, as shown in FIG. 15(c), when the cable 50 needs to be housed in the hook 56, the user can make a decision in the same manner as in FIG. 15(b).

これまでの説明から明らかなように、本実施形態では、尤度分布計算部26と、画像生成部18とにより、拘束点間のケーブルがとりうる形状と該形状の尤度とに基づいて、ケーブルが存在する可能性のある範囲を特定し、表示する表示処理部としての機能が実現されている。 As is clear from the description so far, in this embodiment, the likelihood distribution calculator 26 and the image generator 18 operate based on the shape that the cable between the restraint points can take and the likelihood of the shape. A function as a display processing unit that identifies and displays a range where a cable may exist is realized.

以上、詳細に説明したように、本実施形態によると、制約計算部22は、時刻tにおけるケーブルの制約を発生又は消滅させることでケーブルの拘束点を設定する(S14)。また、ケーブル形状計算部24は、拘束点のケーブルの長さ方向に関する位置ずれの確率分布(図9(c))から、ケーブル上の2つの拘束点間におけるケーブルの長さ誤差の確率分布を求める。さらに、ケーブル形状計算部24は、求めた長さ誤差の確率分布に基づいて、2つの拘束点間においてケーブルがとりうる形状と該形状の尤度とを算出する(S16)。そして、尤度分布計算部26は、ケーブルがとりうる形状と該形状の尤度とに基づいて、ケーブルが存在する可能性のある範囲を特定し、画像生成部18は、尤度分布計算部26が特定した範囲の画像を表示部93上に表示する。ここで、ロボットがケーブルを扱う場合に、ケーブルのたるみやロボットがケーブルを把持した際の滑りなどにより、ロボットが把持すべきケーブルの位置がずれたり、電子機器等に固定すべきケーブルの位置がずれたりする。これに対し、本実施形態では、上述のように拘束点の位置ずれの確率分布から求まるケーブル上の2つの拘束点間におけるケーブルの長さ誤差の確率分布を用いてケーブルの存在する可能性のある範囲を特定して表示する。これにより、ロボットがケーブルを扱う際のケーブルの位置ずれを考慮して、ケーブルの存在する可能性のある範囲を精度よく特定し、表示することが可能となる。 As described above in detail, according to the present embodiment, the constraint calculator 22 sets the cable constraint point by generating or eliminating the cable constraint at time t (S14). In addition, the cable shape calculator 24 calculates the probability distribution of the cable length error between two constraint points on the cable from the probability distribution of the positional deviation of the constraint points in the length direction of the cable (Fig. 9(c)). Ask. Further, the cable shape calculator 24 calculates the shape that the cable can take between the two restraint points and the likelihood of that shape based on the calculated probability distribution of the length error (S16). Then, the likelihood distribution calculation unit 26 identifies a range in which the cable may exist based on the shape that the cable can take and the likelihood of the shape, and the image generation unit 18 uses the likelihood distribution calculation unit The image of the range specified by 26 is displayed on the display unit 93 . Here, when the robot handles the cable, the position of the cable to be gripped by the robot shifts due to slack in the cable or slipping when the robot grips the cable, and the position of the cable to be fixed to the electronic device etc. deviate. On the other hand, in the present embodiment, the probability distribution of the cable length error between two constraint points on the cable, which is obtained from the probability distribution of the displacement of the constraint points as described above, is used to determine the possibility of the existence of the cable. Identify and display a certain range. As a result, it is possible to accurately identify and display the range in which the cable may exist, taking into consideration the positional deviation of the cable when the robot handles the cable.

また、本実施形態では、ケーブル形状計算部24は、2つの拘束点間におけるケーブルの長さ誤差の確率分布を、2つの拘束点それぞれのケーブルの長さ方向に関する位置ずれの確率分布とケーブルの長さの製造誤差とに基づいて求める(上式(1)参照)。これにより、ケーブルの長さの製造誤差を考慮して、ケーブルの存在する可能性のある範囲を精度よく特定し、表示することができる。 Further, in this embodiment, the cable shape calculation unit 24 calculates the probability distribution of the length error of the cable between the two constraint points, the probability distribution of the positional deviation in the length direction of the cable at each of the two constraint points, and the probability distribution of the cable. It is determined based on the manufacturing error of the length (see the above formula (1)). As a result, it is possible to accurately specify and display the range in which the cable may exist, taking into consideration manufacturing errors in the length of the cable.

また、本実施形態では、ケーブル形状計算部24は、拘束点間における長さ誤差の確率分布と、ケーブルの重量の製造誤差と、ケーブルの剛性の製造誤差と、に基づいて、ケーブルがとりうる形状と尤度を算出する。これにより、拘束点間におけるケーブルの長さ誤差の確率分布と、ケーブルの重量や剛性の製造誤差を考慮して、ケーブルの存在する可能性のある範囲を精度よく特定し、表示することができる。 In addition, in this embodiment, the cable shape calculator 24 calculates the possible length of the cable based on the probability distribution of the length error between the restraint points, the manufacturing error of the weight of the cable, and the manufacturing error of the stiffness of the cable. Compute shape and likelihood. As a result, considering the probability distribution of cable length error between restraint points and the manufacturing error of cable weight and stiffness, it is possible to accurately identify and display the range where the cable may exist. .

また、本実施形態では、画像生成部18は、ケーブルが存在する可能性のある範囲のうち存在可能性の高い箇所と低い箇所を異なる態様(例えば透過率を異ならせるなど)で表示する。これにより、ユーザは、ケーブルが存在する可能性の高低に基づいて、工法の検討、作業失敗リスクの把握、干渉チェックなどを行うことができる。なお、画像生成部18は、ケーブルが存在する可能性の高い箇所と低い箇所を色分けや透過率の違いにより表現する場合について説明したが、これに限られるものではない。例えば、画像生成部18は、模様の違い、明暗の違いなどにより、ケーブルが存在する可能性の高い箇所と低い箇所を区別可能に表示することとしてもよい。 Further, in the present embodiment, the image generation unit 18 displays, in a range in which a cable may exist, areas with a high possibility of existence and areas with a low possibility of existence in different modes (for example, different transmittances). As a result, the user can examine the construction method, grasp the risk of work failure, and check for interference based on the possibility that the cable exists. Although the image generating unit 18 has described the case where the locations where the cable is likely to exist and the locations where the cable is unlikely to exist are represented by different colors or different transmittances, the present invention is not limited to this. For example, the image generation unit 18 may display a location with a high possibility of a cable being present and a location with a low possibility of the presence of a cable so as to be distinguishable from each other by means of differences in patterns, differences in brightness, and the like.

なお、上記実施形態では、ケーブル形状計算部24が、拘束点間における長さ誤差の確率分布と、ケーブルの重量及び剛性の製造誤差と、に基づいて、ケーブルがとりうる形状と尤度を算出することとしたが、これに限られるものではない。例えば、ケーブルがとりうる形状と尤度の算出において、ケーブルの重量の製造誤差及び剛性の製造誤差の少なくとも一方を用いないこととしてもよい。 In the above embodiment, the cable shape calculation unit 24 calculates the shape and likelihood that the cable can take based on the probability distribution of the length error between the restraint points and the manufacturing error of the weight and stiffness of the cable. However, it is not limited to this. For example, in calculating the shape and likelihood that the cable can take, at least one of the cable weight manufacturing error and stiffness manufacturing error may not be used.

なお、上記の処理機能は、コンピュータによって実現することができる。その場合、処理装置が有すべき機能の処理内容を記述したプログラムが提供される。そのプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。処理内容を記述したプログラムは、コンピュータで読み取り可能な記憶媒体(ただし、搬送波は除く)に記録しておくことができる。 Note that the processing functions described above can be realized by a computer. In that case, a program is provided that describes the processing contents of the functions that the processing device should have. By executing the program on a computer, the above processing functions are realized on the computer. A program describing the processing content can be recorded in a computer-readable storage medium (excluding carrier waves).

プログラムを流通させる場合には、例えば、そのプログラムが記録されたDVD(Digital Versatile Disc)、CD-ROM(Compact Disc Read Only Memory)などの可搬型記憶媒体の形態で販売される。また、プログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することもできる。 When a program is distributed, it is sold in the form of a portable storage medium such as a DVD (Digital Versatile Disc) or a CD-ROM (Compact Disc Read Only Memory) on which the program is recorded. It is also possible to store the program in the storage device of the server computer and transfer the program from the server computer to another computer via the network.

プログラムを実行するコンピュータは、例えば、可搬型記憶媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、自己の記憶装置に格納する。そして、コンピュータは、自己の記憶装置からプログラムを読み取り、プログラムに従った処理を実行する。なお、コンピュータは、可搬型記憶媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することもできる。また、コンピュータは、サーバコンピュータからプログラムが転送されるごとに、逐次、受け取ったプログラムに従った処理を実行することもできる。 A computer that executes a program stores, for example, a program recorded on a portable storage medium or a program transferred from a server computer in its own storage device. The computer then reads the program from its own storage device and executes processing according to the program. The computer can also read the program directly from the portable storage medium and execute processing according to the program. In addition, the computer can also execute processing in accordance with the received program each time the program is transferred from the server computer.

上述した実施形態は本発明の好適な実施の例である。但し、これに限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変形実施可能である。 The embodiments described above are examples of preferred implementations of the present invention. However, the present invention is not limited to this, and various modifications can be made without departing from the spirit of the present invention.

なお、以上の実施形態の説明に関して、更に以下の付記を開示する。
(付記1) 線状柔軟物に拘束点を設定する設定部と、
前記拘束点の前記線状柔軟物の長さ方向に関する位置ずれの確率分布から、前記線状柔軟物上の2つの拘束点間における前記線状柔軟物の長さ誤差の確率分布を求め、求めた前記長さ誤差の確率分布に基づいて、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の尤もらしさを示す指標とを算出する算出部と、
前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の前記指標とに基づいて、前記線状柔軟物が存在する可能性のある範囲を特定し、表示する表示処理部と、
を備える情報処理装置。
(付記2) 前記算出部は、前記長さ誤差の確率分布を、前記2つの拘束点それぞれの前記線状柔軟物の長さ方向に関する位置ずれの確率分布と、前記線状柔軟物の長さの製造誤差と、に基づいて求める、ことを特徴とする付記1に記載の情報処理装置。
(付記3) 前記算出部は、求めた前記長さ誤差の確率分布と、前記線状柔軟物の重量の製造誤差及び剛性の製造誤差の少なくともいずれかと、に基づいて、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の尤度を算出する、ことを特徴とする付記1又は2に記載の情報処理装置。
(付記4) 前記表示処理部は、前記範囲のうち前記線状柔軟物の存在可能性の高い箇所と低い箇所を異なる態様で表示する、ことを特徴とする付記1~3のいずれかに記載の情報処理装置。
(付記5) 線状柔軟物に拘束点を設定し、
前記拘束点の前記線状柔軟物の長さ方向に関する位置ずれの確率分布から、前記線状柔軟物上の2つの拘束点間における前記線状柔軟物の長さ誤差の確率分布を求め、求めた前記長さ誤差の確率分布に基づいて、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の尤もらしさを示す指標とを算出し、
前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の前記指標とに基づいて、前記線状柔軟物が存在する可能性のある範囲を特定し、表示する、
処理をコンピュータに実行させるための表示プログラム。
(付記6) 前記算出する処理では、前記長さ誤差の確率分布を、前記2つの拘束点それぞれの前記線状柔軟物の長さ方向に関する位置ずれの確率分布と、前記線状柔軟物の長さの製造誤差と、に基づいて求める、ことを特徴とする付記5に記載の表示プログラム。
(付記7) 前記算出する処理では、求めた前記長さ誤差の確率分布と、前記線状柔軟物の重量の製造誤差及び剛性の製造誤差の少なくともいずれかと、に基づいて、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の尤度を算出する、ことを特徴とする付記5又は6に記載の表示プログラム。
(付記8) 前記表示する処理では、前記範囲のうち前記線状柔軟物の存在可能性の高い箇所と低い箇所を異なる態様で表示する、ことを特徴とする付記5~7のいずれかに記載の表示プログラム。
(付記9) 線状柔軟物に拘束点を設定し、
前記拘束点の前記線状柔軟物の長さ方向に関する位置ずれの確率分布から、前記線状柔軟物上の2つの拘束点間における前記線状柔軟物の長さ誤差の確率分布を求め、求めた前記長さ誤差の確率分布に基づいて、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の尤もらしさを示す指標とを算出し、
前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の前記指標とに基づいて、前記線状柔軟物が存在する可能性のある範囲を特定し、表示する、
処理をコンピュータが実行することを特徴とする表示方法。
In addition, the following additional remarks will be disclosed with respect to the above description of the embodiment.
(Appendix 1) A setting unit that sets constraint points on a linear flexible object;
The probability distribution of the length error of the linear flexible object between two restraint points on the linear flexible object is obtained from the probability distribution of the positional deviation of the constraint points in the length direction of the linear flexible object. a calculation unit that calculates a shape that the linear flexible object can take between the two restraint points and an index that indicates the likelihood of the shape, based on the probability distribution of the length error;
a display processing unit that identifies and displays a range in which the flexible linear object may exist based on the shape that the flexible linear object can take between the two constraint points and the index of the shape; ,
Information processing device.
(Supplementary note 2) The calculation unit calculates the probability distribution of the length error by combining the probability distribution of the positional deviation of the two restraint points in the length direction of the flexible linear object with the length of the flexible linear object. The information processing apparatus according to Supplementary Note 1, wherein the information processing apparatus is obtained based on a manufacturing error of .
(Supplementary Note 3) The calculation unit calculates the distance between the two restraint points based on the obtained probability distribution of the length error and at least one of the weight manufacturing error and the rigidity manufacturing error of the linear flexible object. 3. The information processing apparatus according to appendix 1 or 2, wherein the shape that the linear flexible object can take and the likelihood of the shape are calculated.
(Supplementary Note 4) The display processing unit according to any one of Supplementary Notes 1 to 3, wherein the display processing unit displays a location with a high possibility of existence of the linear flexible object and a location with a low possibility of existence in the range in different manners. information processing equipment.
(Appendix 5) Set a constraint point on the linear flexible object,
The probability distribution of the length error of the linear flexible object between two restraint points on the linear flexible object is obtained from the probability distribution of the positional deviation of the constraint points in the length direction of the linear flexible object. Based on the probability distribution of the length error, calculating a shape that the linear flexible object can take between the two restraint points and an index indicating the likelihood of the shape,
identifying and displaying a range in which the linear flexible object may exist based on the shape that the linear flexible object can take between the two constraint points and the index of the shape;
A display program for executing processing on a computer.
(Supplementary Note 6) In the calculation process, the probability distribution of the length error is determined by combining the probability distribution of the positional deviation of the two constraint points in the length direction of the flexible linear object with the length of the flexible linear object. The display program according to appendix 5, wherein the display program is obtained based on the manufacturing error of the thickness.
(Appendix 7) In the calculating process, the two constraint points are determined based on the obtained probability distribution of the length error and at least one of the weight manufacturing error and the stiffness manufacturing error of the linear flexible object. 7. The display program according to appendix 5 or 6, wherein the shape that the linear flexible object can take in between and the likelihood of the shape are calculated.
(Additional remark 8) According to any one of additional remarks 5 to 7, characterized in that, in the displaying process, a location with a high possibility of existence of the linear flexible object and a location with a low possibility of existence of the linear flexible object are displayed in different manners. display program.
(Appendix 9) Set a constraint point on the linear flexible object,
The probability distribution of the length error of the linear flexible object between two restraint points on the linear flexible object is obtained from the probability distribution of the positional deviation of the constraint points in the length direction of the linear flexible object. Based on the probability distribution of the length error, calculating a shape that the linear flexible object can take between the two restraint points and an index indicating the likelihood of the shape,
identifying and displaying a range in which the linear flexible object may exist based on the shape that the linear flexible object can take between the two constraint points and the index of the shape;
A display method characterized in that processing is executed by a computer.

10 情報処理装置
18 画像生成部(表示処理部の一部)
22 制約計算部(設定部)
24 ケーブル形状計算部(算出部)
26 尤度分布計算部(表示処理部の一部)
10 information processing device 18 image generation unit (part of display processing unit)
22 Constraint calculation part (setting part)
24 cable shape calculator (calculator)
26 Likelihood distribution calculator (part of display processor)

Claims (6)

線状柔軟物に複数の拘束点を設定する設定部と、
設定した前記複数の拘束点のうち前記線状柔軟物上で隣接する2つの拘束点間における前記線状柔軟物の長さ誤差の確率分布を、前記2つの拘束点それぞれの前記線状柔軟物の長さ方向に関する位置ずれの確率分布から求め、求めた前記長さ誤差の確率分布に基づいて、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の尤もらしさを示す指標とを算出する算出部と、
前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の前記指標とに基づいて、前記線状柔軟物が存在する可能性のある範囲を特定し、表示する表示処理部と、
を備える情報処理装置。
a setting unit for setting a plurality of constraint points on a linear flexible object;
The probability distribution of the length error of the linear flexible object between two constraint points adjacent to each other on the linear flexible object among the plurality of set constraint points is calculated as follows: Obtained from the probability distribution of positional deviation in the length direction , and based on the obtained probability distribution of the length error, showing the shape that the linear flexible object can take between the two restraint points and the likelihood of the shape a calculation unit that calculates the index;
a display processing unit that identifies and displays a range in which the flexible linear object may exist based on the shape that the flexible linear object can take between the two constraint points and the index of the shape; ,
Information processing device.
前記算出部は、前記長さ誤差の確率分布を、前記2つの拘束点それぞれの前記線状柔軟物の長さ方向に関する位置ずれの確率分布と、前記線状柔軟物の長さの製造誤差と、に基づいて求める、ことを特徴とする請求項1に記載の情報処理装置。 The calculation unit calculates the probability distribution of the length error by combining the probability distribution of the positional deviation of the two restraint points in the length direction of the flexible linear object and the manufacturing error of the length of the flexible linear object. 2. The information processing apparatus according to claim 1, wherein the calculation is performed based on . 前記算出部は、求めた前記長さ誤差の確率分布と、前記線状柔軟物の重量の製造誤差及び剛性の製造誤差の少なくともいずれかと、に基づいて、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の尤度を算出する、ことを特徴とする請求項1又は2に記載の情報処理装置。 Based on the obtained probability distribution of the length error and at least one of the weight manufacturing error and the rigidity manufacturing error of the linear flexible article, the calculation unit calculates the length of the linear flexible article between the two restraint points. 3. The information processing apparatus according to claim 1, wherein the shape that the flexible object can take and the likelihood of the shape are calculated. 前記表示処理部は、前記範囲のうち前記線状柔軟物の存在可能性の高い箇所と低い箇所を異なる態様で表示する、ことを特徴とする請求項1~3のいずれか一項に記載の情報処理装置。 4. The display processing unit according to any one of claims 1 to 3, wherein the display processing unit displays a location with a high possibility of existence of the linear flexible object and a location with a low possibility of existence of the linear flexible object in different modes. Information processing equipment. 線状柔軟物に複数の拘束点を設定し、
設定した前記複数の拘束点のうち前記線状柔軟物上で隣接する2つの拘束点間における前記線状柔軟物の長さ誤差の確率分布を、前記2つの拘束点それぞれの前記線状柔軟物の長さ方向に関する位置ずれの確率分布から求め、求めた前記長さ誤差の確率分布に基づいて、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の尤もらしさを示す指標とを算出し、
前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の前記指標とに基づいて、前記線状柔軟物が存在する可能性のある範囲を特定し、表示する、
処理をコンピュータに実行させるための表示プログラム。
Set multiple constraint points on the linear flexible object,
The probability distribution of the length error of the linear flexible object between two constraint points adjacent to each other on the linear flexible object among the plurality of set constraint points is calculated as follows: Obtained from the probability distribution of positional deviation in the length direction , and based on the obtained probability distribution of the length error, showing the shape that the linear flexible object can take between the two restraint points and the likelihood of the shape Calculate the index and
identifying and displaying a range in which the linear flexible object may exist based on the shape that the linear flexible object can take between the two constraint points and the index of the shape;
A display program for executing processing on a computer.
線状柔軟物に複数の拘束点を設定し、
設定した前記複数の拘束点のうち前記線状柔軟物上で隣接する2つの拘束点間における前記線状柔軟物の長さ誤差の確率分布を、前記2つの拘束点それぞれの前記線状柔軟物の長さ方向に関する位置ずれの確率分布から求め、求めた前記長さ誤差の確率分布に基づいて、前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の尤もらしさを示す指標とを算出し、
前記2つの拘束点間において前記線状柔軟物がとりうる形状と該形状の前記指標とに基づいて、前記線状柔軟物が存在する可能性のある範囲を特定し、表示する、
処理をコンピュータが実行することを特徴とする表示方法。
Set multiple constraint points on the linear flexible object,
The probability distribution of the length error of the linear flexible object between two constraint points adjacent to each other on the linear flexible object among the plurality of set constraint points is calculated as follows: Obtained from the probability distribution of positional deviation in the length direction , and based on the obtained probability distribution of the length error, showing the shape that the linear flexible object can take between the two restraint points and the likelihood of the shape Calculate the index and
identifying and displaying a range in which the linear flexible object may exist based on the shape that the linear flexible object can take between the two constraint points and the index of the shape;
A display method characterized in that processing is executed by a computer.
JP2018198329A 2018-10-22 2018-10-22 Information processing device, display program and display method Active JP7180274B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018198329A JP7180274B2 (en) 2018-10-22 2018-10-22 Information processing device, display program and display method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018198329A JP7180274B2 (en) 2018-10-22 2018-10-22 Information processing device, display program and display method

Publications (2)

Publication Number Publication Date
JP2020067694A JP2020067694A (en) 2020-04-30
JP7180274B2 true JP7180274B2 (en) 2022-11-30

Family

ID=70390326

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018198329A Active JP7180274B2 (en) 2018-10-22 2018-10-22 Information processing device, display program and display method

Country Status (1)

Country Link
JP (1) JP7180274B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114722552B (en) * 2022-06-09 2022-11-08 深圳荣耀智能机器有限公司 Cable length verification method and electronic equipment

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004127925A (en) 2002-09-09 2004-04-22 Yazaki Corp Variation prediction method for wire harness, its device, and its program
JP2005038398A (en) 2003-07-01 2005-02-10 Kanto Auto Works Ltd Deformation analysis device for soft object
JP2005306187A (en) 2004-04-21 2005-11-04 Kanto Auto Works Ltd Assembling simulation device of linear flexible material
JP2010238115A (en) 2009-03-31 2010-10-21 Furukawa Electric Co Ltd:The Method for generation of analysis model initial shape, and program

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004127925A (en) 2002-09-09 2004-04-22 Yazaki Corp Variation prediction method for wire harness, its device, and its program
JP2005038398A (en) 2003-07-01 2005-02-10 Kanto Auto Works Ltd Deformation analysis device for soft object
JP2005306187A (en) 2004-04-21 2005-11-04 Kanto Auto Works Ltd Assembling simulation device of linear flexible material
JP2010238115A (en) 2009-03-31 2010-10-21 Furukawa Electric Co Ltd:The Method for generation of analysis model initial shape, and program

Also Published As

Publication number Publication date
JP2020067694A (en) 2020-04-30

Similar Documents

Publication Publication Date Title
US10394979B2 (en) Method and device for elastic object deformation modeling
RU2017136560A (en) SYSTEM AND METHOD FOR INTERACTIVE HELP WITH COGNITIVE MEANS IN SOLVING TASKS
US11292123B2 (en) Continuum robot, modification method of kinematic model of continuum robot, and control method of continuum robot
US8988435B1 (en) Deforming a skin representation using muscle geometries
US20070231779A1 (en) Systems and Methods for Simulation of Organ Dynamics
JP6771995B2 (en) Methods and systems for obtaining vision measurements for digital human models
RU2308764C2 (en) Method for moving a virtual jointed object in virtual space with prevention of collisions of jointed object with elements of environment
US20100305908A1 (en) Harness verification apparatus, harness verification method and storage medium
JP7180274B2 (en) Information processing device, display program and display method
JP4866896B2 (en) Weather forecasting device, weather forecasting method, and weather forecasting program
Spillmann et al. Inextensible elastic rods with torsional friction based on Lagrange multipliers
Hauenstein et al. An a posteriori certification algorithm for Newton homotopies
JP6718098B2 (en) Position estimation apparatus and method
Li et al. A novel FEM-based numerical solver for interactive catheter simulation in virtual catheterization
JP2004507833A (en) How to provide bitwise constraints for test generation
CN105095555B (en) It is a kind of based on particle image velocimetry method and device of the velocity field without scattered smoothing processing
JP6778631B2 (en) Building design information correction support device, building design information correction support method, and program
CN106650235B (en) Method and device for calculating length of jumper cable
CN113034603B (en) Method and device for determining calibration parameters
US8731880B2 (en) Invertible contact model
Endo et al. An application of a digital hand to ergonomic assessment of handheld information appliances
JP4973784B2 (en) Design support program, design support method, and design support apparatus
US11514211B2 (en) Method and system for performing a simulation of a retraction cable motion
JP2013101107A (en) Method for specifying measurement axis having error of three-axis sensor
CN110353852B (en) Method for measuring position relation between two teeth based on computer

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210610

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220524

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220705

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220829

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20221018

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221031

R150 Certificate of patent or registration of utility model

Ref document number: 7180274

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150