[go: up one dir, main page]

JP7282051B2 - 演算装置、演算方法およびプログラム - Google Patents

演算装置、演算方法およびプログラム Download PDF

Info

Publication number
JP7282051B2
JP7282051B2 JP2020045337A JP2020045337A JP7282051B2 JP 7282051 B2 JP7282051 B2 JP 7282051B2 JP 2020045337 A JP2020045337 A JP 2020045337A JP 2020045337 A JP2020045337 A JP 2020045337A JP 7282051 B2 JP7282051 B2 JP 7282051B2
Authority
JP
Japan
Prior art keywords
article
resin
stiffness
tomographic images
unit
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
JP2020045337A
Other languages
English (en)
Other versions
JP2021146514A (ja
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.)
Mitsubishi Heavy Industries Ltd
Original Assignee
Mitsubishi Heavy Industries 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 Mitsubishi Heavy Industries Ltd filed Critical Mitsubishi Heavy Industries Ltd
Priority to JP2020045337A priority Critical patent/JP7282051B2/ja
Priority to US17/801,411 priority patent/US12367572B2/en
Priority to PCT/JP2021/005665 priority patent/WO2021186988A1/ja
Publication of JP2021146514A publication Critical patent/JP2021146514A/ja
Application granted granted Critical
Publication of JP7282051B2 publication Critical patent/JP7282051B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • G06T7/0008Industrial image inspection checking presence/absence
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B29WORKING OF PLASTICS; WORKING OF SUBSTANCES IN A PLASTIC STATE IN GENERAL
    • B29CSHAPING OR JOINING OF PLASTICS; SHAPING OF MATERIAL IN A PLASTIC STATE, NOT OTHERWISE PROVIDED FOR; AFTER-TREATMENT OF THE SHAPED PRODUCTS, e.g. REPAIRING
    • B29C64/00Additive manufacturing, i.e. manufacturing of three-dimensional [3D] objects by additive deposition, additive agglomeration or additive layering, e.g. by 3D printing, stereolithography or selective laser sintering
    • B29C64/10Processes of additive manufacturing
    • B29C64/106Processes of additive manufacturing using only liquids or viscous materials, e.g. depositing a continuous bead of viscous material
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B29WORKING OF PLASTICS; WORKING OF SUBSTANCES IN A PLASTIC STATE IN GENERAL
    • B29CSHAPING OR JOINING OF PLASTICS; SHAPING OF MATERIAL IN A PLASTIC STATE, NOT OTHERWISE PROVIDED FOR; AFTER-TREATMENT OF THE SHAPED PRODUCTS, e.g. REPAIRING
    • B29C64/00Additive manufacturing, i.e. manufacturing of three-dimensional [3D] objects by additive deposition, additive agglomeration or additive layering, e.g. by 3D printing, stereolithography or selective laser sintering
    • B29C64/30Auxiliary operations or equipment
    • B29C64/386Data acquisition or data processing for additive manufacturing
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B33ADDITIVE MANUFACTURING TECHNOLOGY
    • B33YADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
    • B33Y50/00Data acquisition or data processing for additive manufacturing
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B33ADDITIVE MANUFACTURING TECHNOLOGY
    • B33YADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
    • B33Y70/00Materials specially adapted for additive manufacturing
    • B33Y70/10Composites of different types of material, e.g. mixtures of ceramics and polymers or mixtures of metals and biomaterials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G01N23/046Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/12Geometric CAD characterised by design entry means specially adapted for CAD, e.g. graphical user interfaces [GUI] specially adapted for CAD
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/22Moulding
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/24Sheet material
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/26Composites
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30144Printing quality

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Chemical & Material Sciences (AREA)
  • Materials Engineering (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Manufacturing & Machinery (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Health & Medical Sciences (AREA)
  • Optics & Photonics (AREA)
  • Mechanical Engineering (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Composite Materials (AREA)
  • Human Computer Interaction (AREA)
  • Architecture (AREA)
  • Structural Engineering (AREA)
  • Civil Engineering (AREA)
  • Ceramic Engineering (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pulmonology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Description

本開示は、演算装置、演算方法およびプログラムに関する。
材料の強度解析の手法として、有限要素法(Finite Element Method、FEM)が知られている。
特許文献1には、複合材の微視的観察画像を複数の領域に分割し、複合材の有限要素モデルを生成する技術が開示されている。
特許文献2には、タイヤの形状を再現した初期解析モデルに対して、タイヤの残留応力を取得して釣り合い計算を実行して修正解析モデルを作成することにより、構造体の形状及び内部の応力を精度よく再現する技術が開示されている。
特開平9-180002号公報 特許第4581539号公報
強化繊維と樹脂を含む繊維強化プラスチックが知られている。繊維強化プラスチックの強度は、繊維の方向や密度によって異なる。
一方で、近年、繊維強化プラスチックを形成する3Dプリンタ等が検討されており、繊維強化プラスチックで構成される物品は複雑な内部構造を有する場合がある。この場合、物品の内部の剛性は一様でない。そのため、有限要素法により適切に強度評価を行う場合、複雑な内部形状に基づいて、物品を切り分けた複数の小領域それぞれについて剛性に係るパラメータを求めなければならない。
本開示の目的は、上述した課題を解決する演算装置、演算方法およびプログラムを提供することにある。
本開示に係る演算装置は、強化繊維と樹脂を含む物品を平行に切り分けたときの複数の断面における強化繊維の形状および樹脂の形状を表す複数の断層像を取得する断層像取得部と、取得した複数の断層像に基づいて、物品を分割した複数の小領域のそれぞれに係る剛性に係るパラメータを算出する演算部と、を備える。
本開示に係る演算方法は、強化繊維と樹脂を含む物品を平行に切り分けたときの複数の断面における強化繊維の形状および樹脂の形状を表す複数の断層像を取得するステップと、取得した複数の断層像に基づいて、物品を分割した複数の小領域のそれぞれに係る剛性に係るパラメータを算出するステップと、を有する。
本開示に係るプログラムは、コンピュータを、強化繊維と樹脂を含む物品を平行に切り分けたときの複数の断面における強化繊維の形状および樹脂の形状を表す複数の断層像を取得するステップと、取得した複数の断層像に基づいて、物品を分割した複数の小領域のそれぞれに係る剛性に係るパラメータを算出するステップと、として実行する。
上記態様のうち少なくとも1つの態様によれば、物品の内部構造を反映させて算出された剛性パラメータを用いて、強度評価を行うためのデータを生成することができる。
一実施形態に係る演算システムの構成を示す図である。 一実施形態に係るCADデータの一例を示す図である。 一実施形態に係る断層像の一例を示す図である。 一実施形態に係る演算装置を示す概略ブロック図である。 一実施形態に係る計算格子データの一例を示す図である。 一実施形態に係る演算システムの動作を示すフローチャートである。 一実施形態に係る演算システムの構成を示す図である。 一実施形態に係る演算装置の構成を示す概略ブロック図である。 少なくとも1つの実施形態に係るコンピュータの構成を示す概略ブロック図である。
〈第1の実施形態〉
《演算システムの構成》
以下、図面を参照しながら実施形態について詳しく説明する。
図1は、第1の実施形態に係る演算システム10を示す図である。
演算システム10は、3Dプリンタ100が形成した物品300の強度評価を行うシステムである。演算システム10は、3Dプリンタ100と、演算装置200と、を備える。3Dプリンタ100は、演算装置200と有線又は無線で接続されている。物品300は、強化繊維と樹脂とからなる繊維強化プラスチック材で構成される。
《3Dプリンタの構成》
以下、3Dプリンタの構成について説明する。
3Dプリンタ100は、物品300の三次元形状を示すCAD(Computer-Aided Design)データ310の入力を受け付け、当該CADデータに基づいて物品300を形成する。3Dプリンタ100は、強化繊維を出力するヘッダと樹脂を出力するヘッダをそれぞれ備え、ヘッダが出力した強化繊維および樹脂の層を積層することで、物品300を形成する。
3Dプリンタ100は、CADデータに基づいて、物品300を形成する各層を形成するためのヘッダの経路を計算し、当該経路に基づいてヘッダを駆動させる。3Dプリンタ100は、各層におけるヘッダの経路を表す複数の画像を、利用者による参照のために出力する機能を有する。
強化繊維とは物品300の作製に用いられ、物品300の強度を向上させる材料である。強化繊維の例としては、炭素繊維、ガラス繊維等が挙げられる。第1の実施形態における強化繊維は炭素繊維である。また、樹脂とは物品300の作製に用いられ、物品300を形成させるための材料である。樹脂の例としては、ナイロン樹脂、エポキシ樹脂、ウレタン樹脂、ポリエーテルエーテルケトン(Polyetheretherketone,PEEK)、ポリエーテルケトンケトン(Polyetherketoneketone,PEKK)等が挙げられる。
3Dプリンタ100の詳細な動作の一例については、以下に説明する。
まず、3Dプリンタ100は、物品300を示すCADデータを受け入れる。図2は、第1の実施形態に係るCADデータ310の一例を示す図である。CADデータ310は物品300の形状を三次元で示す画像である。
3Dプリンタ100は、CADデータ310を受け入れると、当該CADデータが示す三次元形状を3Dプリンタ100の積層厚さごとに切り分けることで、物品300の断面形状を特定する。積層厚さは、ユーザによって設定されても良いし、3Dプリンタ100に予め設定されていても良い。
3Dプリンタ100は、積層厚さ毎に切り分けた各層の断面形状に基づいて、各層における強化繊維を出力するヘッダのパスを決定する。強化繊維のヘッダのパスは、予め定められた制約条件に基づいて決定される。制約条件は、例えば、許容曲げ角度や密度などが挙げられる。
次に、3Dプリンタ100は、各層の断面形状と強化繊維のヘッダのパスとに基づいて、樹脂を出力するヘッダのパスを決定する。樹脂を出力するヘッダのパスは、断面形状のうち強化繊維のヘッダが通らない部分を充填するように決定される。3Dプリンタ100は、各層における強化繊維のヘッダのパスおよび樹脂のヘッダのパスを示す複数の断層像320を生成する。図3は、第1の実施形態に係る断層像320の一例を示す図である。断層像320は、強化繊維のヘッダのパス322と、樹脂のヘッダのパス321とを異なる色で表す画像である。すなわち、断層像320は、強化繊維と樹脂とを異なる値の画素で示すラスタ画像である。図3においては、強化繊維のヘッダのパス322を白の点線で示し、樹脂のヘッダのパス321を白の実線で示す。
図3に示す例においては、3Dプリンタ100は、CADデータ310が示す三次元形状を、積層厚さの間隔で5つに切り分けて断層像320A、断層像320B、断層像320C、断層像D、断層像Eを生成する。
《演算装置の構成》
以下、演算装置200の構成について説明する。
図4は、演算装置200の構成を示す概略ブロック図である。演算装置200は、3Dプリンタ100が出力する複数の断層像320と、物品300の三次元形状を複数の小領域に切り分けた計算格子とに基づいて、物品300の強度解析を行う。演算装置200は、取得部210と、断層像特定部220と、変換部230と、材料情報特定部240と、演算部250と、モデル生成部260と、解析部270と、を備える。
取得部210は、3Dプリンタ100から複数の断層像320を取得する。取得部210は断層像取得部の一例である。例えば、取得部210は図3に示すような断層像320Aと、断層像320Bと、断層像320Cと、断層像320Dと、断層像320Eと、を3Dプリンタ100から取得する。
また、取得部210は演算システム10のユーザから計算格子データ330を取得する。計算格子データ330とは、物品300が切り分けられた複数のSolid要素331を示す情報である。図5は計算格子データ330の一例を示す図である。図5に示す計算格子データ330には、Solid要素331が高さ方向に3個、奥行き方向に6個、幅方向に22個が並ぶ。すなわち図5に示す例においては、計算格子は、物品300の三次元形状を合計396個のSolid要素331に分割したものである。なお、Solid要素331は、物品300の三次元形状を分割した小領域の一例である。
断層像特定部220は、取得部210が取得した断層像320と、取得部210が取得した計算格子データ330のSolid要素331とを位置合わせして、Solid要素331に対応する断層像320の部分を特定する。すなわち、断層像特定部220は、複数の断層像320の座標系を計算格子データ330の座標系に変換する。3Dプリンタ100の積層厚さは予め分かっているため、断層像特定部220は、複数の断層像320と積層厚さとに基づいて物品300の三次元形状を計算することができる。そのため、断層像特定部220は、断層像320から生成した三次元形状と計算格子データ330が示す三次元形状との誤差が最小となるような姿勢を特定することで、断層像320と計算格子データ330とを位置合わせすることができる。例えば、断層像特定部220は、図3に示すような断層像320と、図5に示すような計算格子データのSolid要素331Aとを位置合わせして、Solid要素331Aに対応する断層像320の部分を特定する。これにより、断層像特定部220は、物品300におけるSolid要素331Aの位置に対応する断層像320Eの部分と、断層像320Dの部分を特定する。
図5に示すように、計算格子データ330は高さ方向に3個に分けられている。他方、図3に示すように、断層像320は高さ方向に5個に分けられている。このように、通常は、計算格子データ330が高さ方向に分けられた分割数に比べて、断層像320が高さ方向に分けられた枚数が多い。そのため、1つのSolid要素331Aに対応する断層像320は1つではなく複数である。
変換部230は、ラスタ画像である各断層像320から、強化繊維をパス要素で表すベクタ画像を生成する。例えば、変換部230は、断層像320の強化繊維を示す画素からなる図形を特定し、当該図形の細線化処理を行うことで、ベクタ画像を生成しても良い。
材料情報特定部240は、断層像特定部220により特定された各Solid要素331に対応する断層像320の部分画像に基づいて、当該Solid要素331における強化繊維と樹脂との面積比率を特定する。例えば、材料情報特定部240は、各部分画像における強化繊維に係る画素の数と樹脂に係る画素の数の比率を算出する。また、材料情報特定部240は、変換部230により変換されたベクタ画像に基づいて、各小領域における強化繊維の方向を特定する。
演算部250は、材料情報特定部240により特定された強化繊維と樹脂との面積比率及び強化繊維の方向に基づいて、Solid要素331の剛性に係るパラメータを算出する。剛性に係るパラメータの例としては、X軸方向、Y軸方向、Z軸方向などの各方向の縦弾性係数、せん断弾性率、ポアソン比、等のパラメータが挙げられる。上記縦弾性係数は、方向によりその値が異なる。演算部250が、剛性に係るパラメータを算出する動作を以下に説明する。
演算部250は、以下の動作により剛性に係るパラメータを算出する。
演算部250は、材料情報特定部240が特定した強化繊維と樹脂との面積比率を受け入れる。その後、演算部250は、炭素繊維強化プラスチック(Carbon Fiber Reinforced Plastic、CFRP)の剛性と樹脂の剛性を上記面積比率で均質化して、特性を算出する。また、演算部250は、強化繊維の角度を検出し、当該角度と剛性の関係式を用いて、算出された上記特定を角度に応じて変換する。ここで、上記繊維の角度とは、予め設定された基準角度からの角度を示す。これにより、演算部250は剛性に係るパラメータを算出できる。
また、演算部250は、1つのSolid要素331の内部に含まれる全ての画像からそれぞれ計算した上記剛性に係るパラメータを均質化する。例えば、演算部250は、等価剛性マトリックスの各種公式を用いて、剛性パラメータを均質化する計算を行う。
モデル生成部260は、取得部210が取得した計算格子データ330に、演算部250が算出した剛性に係るパラメータを格納することで、有限要素法に用いるSolidモデルを生成する。
解析部270は、演算部250により算出されたSolid要素331の剛性に係るパラメータに基づいて、有限要素法による解析を行う。すなわち、解析部270は、モデル生成部260が生成したSolidモデルに基づいて有限要素法による解析で強度評価を行う。
《演算システムの動作》
以下、演算システム10の動作について説明する。
図6は、演算システム10の動作を示すフローチャートである。
取得部210は3Dプリンタ100から複数の断層像320を取得し、演算システム10のユーザから計算格子データ330を取得する(ステップS1)。計算格子データ330は、予めユーザによってCAD情報310が示す三次元形状をメッシュに切り分けることで生成される。
変換部230は、ラスタ画像である各断層像320から、強化繊維をパス要素で表すベクタ画像を生成する(ステップS2)。
断層像特定部220は、ステップS2で取得した断層像320と、計算格子データ330のSolid要素331とを位置合わせして、各Solid要素331に対応する断層像320の部分画像を特定する(ステップS3)。
演算装置200は、計算格子データ330に含まれる複数のSolid要素331を1つずつ選択し、各小領域について以下のステップS5からステップS7の処理を実行する(ステップS4)。
材料情報特定部240は、ステップS4で選択されたSolid要素331に対応する断層像320の部分画像に基づいて、当該Solid要素331における強化繊維と樹脂との面積比率を特定する(ステップS5)。
材料情報特定部240は、ステップS2で生成されたベクタ画像のうち、ステップS3で選択されたSolid要素331に対応する部分に基づいて、ステップS4で選択されたSolid要素331における強化繊維の方向を特定する(ステップS6)。
演算部250は、ステップS5で特定された強化繊維と樹脂との面積比率、及びステップS6で特定された強化繊維の方向に基づいて、ステップS4で選択されたSolid要素331の剛性に係るパラメータを算出する(ステップS7)。すなわち、演算部250は、剛性に係るパラメータを算出して、等価剛性マトリックスの各種公式に基づいて上記剛性に係るパラメータを均質化する。
計算格子データ330に含まれる複数のSolid要素331のすべてについて剛性に係るパラメータが算出されると、モデル生成部260は、ステップS1で取得した計算格子データに、ステップS7で算出した剛性に係るパラメータを格納することで、有限要素法に用いるSolidモデルを生成する(ステップS8)。
解析部270は、ステップS8で生成したSolidモデルに基づいて、有限要素法による解析を行う。すなわち、解析部270は、物品300の全てのSolid要素331の剛性に係るパラメータに基づいて、有限要素法による解析で強度評価を行う(ステップS9)。
上記の動作により、演算システム10は断層像320及び計算格子データ330を用いて剛性に係るパラメータを算出して有限要素法により、物品300の強度評価を行う。演算システム10のユーザは、物品の内部構造を反映させた有限要素法による強度評価を行うことができる。
《作用・効果》
本開示に係る演算装置200は、強化繊維と樹脂を含む物品を平行に切り分けたときの複数の断面における強化繊維の形状および樹脂の形状を表す複数の断層像320を取得する断層像取得部と、取得した複数の断層像320に基づいて、物品を分割した複数の小領域のそれぞれに係る剛性に係るパラメータを算出する演算部250と、を備える。
演算装置200は、断層像320及び計算格子データ330を用いて剛性に係るパラメータを算出する。これにより、演算装置200のユーザは、物品300の内部構造を反映させて算出された剛性パラメータを用いて強度評価を行うことができる。
また、演算装置200は、取得した複数の断層像320に基づいて、各小領域における強化繊維と樹脂との面積比率と、強化繊維の方向とを特定する材料情報特定部240と、を備え、演算部250は、特定された面積比率及び方向に基づいて、複数の小領域のそれぞれの剛性に係るパラメータを算出する。
演算装置200は強化繊維と樹脂との面積比率と、強化繊維の方向を特定して、剛性に係るパラメータを算出する。これにより、演算装置200のユーザは、物品300の内部構造を反映させた強度評価を行うことができる。
また、演算装置200は、複数の断層像320は、それぞれ強化繊維と樹脂とを異なる値の画素で示すラスタ画像によって表され、複数の断層像320のそれぞれから、強化繊維の形状を表すパスを有するベクタ画像を生成する変換部230を備え、材料情報特定部240は、ベクタ画像に基づいて強化繊維の方向を特定する。
演算装置200は、ベクタ画像に基づいて強化繊維の方向を特定して剛性に係るパラメータを算出する。これにより、演算装置200のユーザは、物品300の内部構造を反映させた強度評価を行うことができる。
また、演算装置200の複数の断層像320は、強化繊維と樹脂とを出力可能な3Dプリンタ100が物品300を形成するときのヘッダの経路を表す。
演算装置200は、3Dプリンタ100が生成した断層像320を用いて剛性に係るパラメータを算出する。これにより、演算装置200のユーザは、物品300の内部構造を反映させた強度評価を行うことができる。
また、演算装置200は、算出された剛性に係るパラメータに基づいて、物品の強度解析を行う解析部270と、を備える。
演算装置200のユーザは、物品の内部構造を反映させて算出された剛性パラメータに基づいて強度評価を行うことができる。
〈第2の実施形態〉
以下、第2の実施形態に係る演算システム10について説明する。
第2の実施形態に係る演算システム10は、3Dプリンタ100と、演算装置200に加えて、X線CT装置400を備える。X線CT装置400は、演算装置200と有線又は無線で接続される。図7は、第2の実施形態に係る演算システム10の構成を示す図である。
X線CT装置400は、物品300にX線を照射して、ボクセルデータで表されるX線CT(Computed Tomography)像を生成する。つまり、ユーザは、3Dプリンタ100が形成した物品300をX線CT装置400に投入し、X線CT装置400にX線CT像を生成させる。ユーザは、X線CT装置400によるX線CT処理の完了後に、3Dプリンタ100が生成した断層像320と、X線CT装置400が生成したX線CT像とを、演算装置200に入力する。
図8は、第2の実施形態に係る演算装置200の構成を示す概略ブロック図である。第2の実施形態に係る演算装置200は、第1の実施形態に係る演算装置200の構成に加えて、欠陥判定部211と、格子編集部212と、を備える。
第2の実施形態に係る取得部210は、断層像取得部及び欠陥データ取得部の一例である。取得部210は、3Dプリンタ100から断層像320を取得し、演算システム10のユーザから計算格子データ330を取得する。また、取得部210は、X線CT装置400から、X線CT装置400が生成した物品300のX線CT像を取得する。取得部210が取得するX線CT像は、物品300の欠陥箇所を示す欠陥データの一例である。
欠陥判定部211は計算格子データ330の各Solid要素331が、欠陥箇所に該当するか否かを判定する。欠陥判定部211の詳細な動作は以下に説明する。
欠陥判定部211は、取得部210が取得したX線CT像と、取得部210が取得した計算格子データ330のSolid要素331とを位置合わせして、Solid要素331に対応するX線CT像の部分を特定する。すなわち、欠陥判定部211は、X線CT像の座標系を計算格子データ330の座標系に変換する。その後、欠陥判定部211は、X線CT像から生成した三次元形状と計算格子データ330が示す三次元形状との誤差が最小となるような姿勢を特定することで、X線CT像と計算格子データ330とを位置合わせすることができる。
Solid要素331に対応するX線CT像の部分を特定した後、欠陥判定部211は、各Solid要素331に対応するX線CT像に含まれたボクセルデータを、予め設定された閾値情報に照らし合わせて、各Solid要素331が、欠陥箇所に該当するか否かを判定する。欠陥の例としては、空隙、ウィークボンド、ひび割れ等が挙げられる。その後、欠陥判定部211は、各Solid要素331と、当該Solid要素331が欠陥箇所に該当するか否かと、を関連付けた情報である欠陥データを生成する。
格子編集部212は、欠陥判定部211が生成した欠陥データに基づいて、取得部210が取得した計算格子データ330のSolid要素331のうち、欠陥箇所に該当するSolid要素331を削除する。
演算システム10は、物品300の欠陥箇所に該当するSolid要素331を削除して、剛性に係るパラメータを生成して強度評価を行う。このため、演算システム10のユーザは、物品300の欠陥情報を含む内部構造を反映させた剛性パラメータを算出して、強度評価を行うことができる。
演算システム10は欠陥箇所に該当するSolid要素331を削除する代わりに、当該Solid要素331の剛性に係るパラメータを編集させることにより、欠陥箇所を反映させた強度評価を行っても良い。
《作用・効果》
本開示に係る演算装置200は、物品300の欠陥箇所を示す欠陥データを取得する欠陥取得部を備え、演算部250は、複数の断層像320及び欠陥データに基づいて、複数の小領域それぞれの剛性に係るパラメータを算出する。
演算装置200は、物品300の欠陥データに基づいて剛性に係るパラメータを生成して強度評価を行う。このため、演算システム10のユーザは、物品300の欠陥データを含む内部構造を反映させた剛性パラメータを算出して、強度評価を行うことができる。
〈他の実施形態〉
以上、図面を参照して一実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、様々な設計変更等をすることが可能である。
上記の実施形態における計算格子データ330は、Solid要素331により構成されるが、これに代えてShell要素により構成されるものであっても良い。
図9は、少なくとも1つの実施形態に係るコンピュータの構成を示す概略ブロック図である。
コンピュータ1100は、プロセッサ1110、メインメモリ1120、ストレージ1130、インタフェース1140を備える。
上述の演算装置200は、コンピュータ1100に実装される。そして、上述した各処理部の動作は、プログラムの形式でストレージ1130に記憶されている。プロセッサ1110は、プログラムをストレージ1130から読み出してメインメモリ1120に展開し、当該プログラムに従って上記処理を実行する。また、プロセッサ1110は、プログラムに従って、上述した各記憶部に対応する記憶領域をメインメモリ1120に確保する。
プログラムは、コンピュータ1100に発揮させる機能の一部を実現するためのものであってもよい。例えば、プログラムは、ストレージ1130に既に記憶されている他のプログラムとの組み合わせ、または他の装置に実装された他のプログラムとの組み合わせによって機能を発揮させるものであってもよい。なお、他の実施形態においては、コンピュータ1100は、上記構成に加えて、または上記構成に代えてPLD(Programmable Logic Device)などのカスタムLSI(Large Scale Integrated Circuit)を備えてもよい。PLDの例としては、PAL(Programmable Array Logic)、GAL(Generic Array Logic)、CPLD(Complex Programmable Logic Device)、FPGA(Field Programmable Gate Array)が挙げられる。この場合、プロセッサ1110によって実現される機能の一部または全部が当該集積回路によって実現されてよい。
ストレージ1130の例としては、磁気ディスク、光磁気ディスク、半導体メモリ等が挙げられる。ストレージ1130は、コンピュータ1100のバスに直接接続された内部メディアであってもよいし、インタフェース1140または通信回線を介してコンピュータに接続される外部メディアであってもよい。また、このプログラムが通信回線によってコンピュータ1100に配信される場合、配信を受けたコンピュータ1100が当該プログラムをメインメモリ1120に展開し、上記処理を実行してもよい。少なくとも1つの実施形態において、ストレージ1130は、一時的でない有形の記憶媒体である。
また、当該プログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、当該プログラムは、前述した機能をストレージ1130に既に記憶されている他のプログラムとの組み合わせで実現するもの、いわゆる差分ファイル(差分プログラム)であってもよい。
〈付記〉
各実施形態に記載の演算装置200は、例えば以下のように把握される。
(1)本開示に係る演算装置200は、強化繊維と樹脂を含む物品を平行に切り分けたときの複数の断面における強化繊維の形状および樹脂の形状を表す複数の断層像320を取得する断層像取得部と、取得した複数の断層像320に基づいて、物品を分割した複数の小領域のそれぞれに係る剛性に係るパラメータを算出する演算部250と、を備える。
演算装置200は、断層像320及び計算格子データ330を用いて剛性に係るパラメータを算出する。これにより、演算装置200のユーザは、物品300の内部構造を反映させて算出された剛性パラメータを用いて強度評価を行うことができる。
(2)また、演算装置200は、取得した複数の断層像320に基づいて、各小領域における強化繊維と樹脂との面積比率と、強化繊維の方向とを特定する材料情報特定部240と、を備え、演算部250は、特定された面積比率及び方向に基づいて、複数の小領域のそれぞれの剛性に係るパラメータを算出する。
演算装置200は強化繊維と樹脂との面積比率と、強化繊維の方向を特定して、剛性に係るパラメータを算出する。これにより、演算装置200のユーザは、物品300の内部構造を反映させた強度評価を行うことができる。
(3)また、演算装置200は、複数の断層像320は、それぞれ強化繊維と樹脂とを異なる値の画素で示すラスタ画像によって表され、複数の断層像320のそれぞれから、強化繊維の形状を表すパスを有するベクタ画像を生成する変換部230を備え、材料情報特定部240は、ベクタ画像に基づいて強化繊維の方向を特定する。
演算装置200は、ベクタ画像に基づいて強化繊維の方向を特定して剛性に係るパラメータを算出する。これにより、演算装置200のユーザは、物品300の内部構造を反映させた強度評価を行うことができる。
(4)また、演算装置200の複数の断層像320は、強化繊維と樹脂とを出力可能な3Dプリンタ100が物品300を形成するときのヘッダの経路を表す。
演算装置200は、3Dプリンタ100が生成した断層像320を用いて剛性に係るパラメータを算出する。これにより、演算装置200のユーザは、物品300の内部構造を反映させた強度評価を行うことができる。
(5)また、演算装置200は、算出された剛性に係るパラメータに基づいて、物品の強度解析を行う解析部270と、を備える。
演算装置200のユーザは、物品の内部構造を反映させて算出された剛性パラメータに基づいて強度評価を行うことができる。
(6)本開示に係る演算装置200は、物品300の欠陥箇所を示す欠陥データを取得する欠陥取得部を備え、演算部250は、複数の断層像320及び欠陥データに基づいて、複数の小領域それぞれの剛性に係るパラメータを算出する。
演算装置200は、物品300の欠陥データに基づいて剛性に係るパラメータを生成して強度評価を行う。このため、演算システム10のユーザは、物品300の欠陥データを含む内部構造を反映させた剛性パラメータを算出して、強度評価を行うことができる。
(7)本開示に係る演算方法は、強化繊維と樹脂を含む物品を平行に切り分けたときの複数の断面における強化繊維の形状および樹脂の形状を表す複数の断層像320を取得するステップと、取得した複数の断層像320に基づいて、物品を分割した複数の小領域のそれぞれに係る剛性に係るパラメータを算出するステップと、を有する。
演算方法を用いるユーザは、断層像320及び計算格子データ330を用いて剛性に係るパラメータを算出することができる。これにより、演算方法のユーザは、物品300の内部構造を反映させて算出された剛性パラメータを用いて強度評価を行うことができる。
(8)本開示に係るプログラムは、強化繊維と樹脂を含む物品を平行に切り分けたときの複数の断面における強化繊維の形状および樹脂の形状を表す複数の断層像320を取得するステップと、取得した複数の断層像320に基づいて、物品を分割した複数の小領域のそれぞれに係る剛性に係るパラメータを算出するステップと、として実行する。
プログラムのユーザは、断層像320及び計算格子データ330を用いて剛性に係るパラメータを算出することができる。これにより、プログラムのユーザは、物品300の内部構造を反映させて算出された剛性パラメータを用いて強度評価を行うことができる。
強化繊維とは物品300の作製に用いられ、物品300の強度を向上させる材料である。強化繊維の例としては、炭素繊維、ガラス繊維等が挙げられる。第1の実施形態における強化繊維は炭素繊維である。また、樹脂とは物品300の作製に用いられ、物品300を形成させるための材料である。樹脂の例としては、ナイロン樹脂、エポキシ樹脂、ウレタン樹脂、ポリエーテルエーテルケトン、ポリエーテルケトンケトン等が挙げられる。
計算格子データ330とは、物品300が切り分けられた複数のSolid要素331を示す情報又は、物品300が切り分けられた複数のShell要素を示す情報である。
100 3Dプリンタ
200 演算装置
210 取得部
211 欠陥判定部
212 格子編集部
220 断層像特定部
230 変換部
240 材料情報特定部
250 演算部
260 モデル生成部
270 解析部
300 物品
310 CADデータ
320 断層像
321 樹脂のヘッダのパス
322 強化繊維のヘッダのパス
330 計算格子データ
331 Solid要素
400 X線CT装置
1100 コンピュータ
1110 プロセッサ
1120 メインメモリ
1130 ストレージ
1140 インタフェース

Claims (8)

  1. 強化繊維と樹脂を含む物品を平行に切り分けたときの複数の断面における前記強化繊維の形状および前記樹脂の形状を表す複数の断層像を取得する断層像取得部と、
    取得した前記複数の断層像に基づいて、前記物品を分割した複数の小領域のそれぞれに係る剛性に係るパラメータを算出する演算部と、
    を備える演算装置。
  2. 取得した前記複数の断層像に基づいて、各小領域における前記強化繊維と前記樹脂との面積比率と、前記強化繊維の方向とを特定する材料情報特定部と、を備え、
    前記演算部は、特定された前記面積比率及び前記方向に基づいて、前記複数の小領域のそれぞれの前記剛性に係るパラメータを算出する
    請求項1に記載の演算装置。
  3. 前記複数の断層像は、それぞれ前記強化繊維と前記樹脂とを異なる値の画素で示すラスタ画像によって表され、
    前記複数の断層像のそれぞれから、前記強化繊維の形状を表すパスを有するベクタ画像を生成する変換部を備え、
    前記材料情報特定部は、前記ベクタ画像に基づいて前記強化繊維の方向を特定する
    請求項2に記載の演算装置。
  4. 前記複数の断層像は、前記強化繊維と前記樹脂とを出力可能な3Dプリンタが前記物品を形成するときのヘッダの経路を表す
    請求項1から請求項3の何れか1項に記載の演算装置。
  5. 前記物品の欠陥箇所を示す欠陥データを取得する欠陥取得部を備え、
    前記演算部は、前記複数の断層像及び前記欠陥データに基づいて、前記複数の小領域それぞれの剛性に係るパラメータを算出する
    請求項1に記載の演算装置。
  6. 算出された前記剛性に係るパラメータに基づいて、前記物品の強度解析を行う解析部と、
    を備える請求項1から請求項5の何れか1項に記載の演算装置。
  7. 強化繊維と樹脂を含む物品を平行に切り分けたときの複数の断面における前記強化繊維の形状および前記樹脂の形状を表す複数の断層像を取得するステップと、
    取得した前記複数の断層像に基づいて、前記物品を分割した複数の小領域のそれぞれに係る剛性に係るパラメータを算出するステップと、
    を有する演算方法。
  8. コンピュータを、
    強化繊維と樹脂を含む物品を平行に切り分けたときの複数の断面における前記強化繊維の形状および前記樹脂の形状を表す複数の断層像を取得するステップと、
    取得した前記複数の断層像に基づいて、前記物品を分割した複数の小領域のそれぞれに係る剛性に係るパラメータを算出するステップと、
    として実行するプログラム。
JP2020045337A 2020-03-16 2020-03-16 演算装置、演算方法およびプログラム Active JP7282051B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2020045337A JP7282051B2 (ja) 2020-03-16 2020-03-16 演算装置、演算方法およびプログラム
US17/801,411 US12367572B2 (en) 2020-03-16 2021-02-16 Computing device and computing method for evaluating strength of article formed by 3D printer
PCT/JP2021/005665 WO2021186988A1 (ja) 2020-03-16 2021-02-16 演算装置、演算方法およびプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2020045337A JP7282051B2 (ja) 2020-03-16 2020-03-16 演算装置、演算方法およびプログラム

Publications (2)

Publication Number Publication Date
JP2021146514A JP2021146514A (ja) 2021-09-27
JP7282051B2 true JP7282051B2 (ja) 2023-05-26

Family

ID=77770831

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020045337A Active JP7282051B2 (ja) 2020-03-16 2020-03-16 演算装置、演算方法およびプログラム

Country Status (3)

Country Link
US (1) US12367572B2 (ja)
JP (1) JP7282051B2 (ja)
WO (1) WO2021186988A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117434060B (zh) * 2023-12-20 2024-04-09 深圳市智能派科技有限公司 光固化3d打印用光敏树脂的测试系统、方法及控制设备

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09180002A (ja) 1995-12-25 1997-07-11 Matsushita Electric Works Ltd 複合材料の微視的構造の有限要素モデル作成方法
JPH10124705A (ja) * 1996-10-17 1998-05-15 Nagano Pref Gov 長尺積層材料の設計支援システム
JP4519996B2 (ja) * 2000-07-04 2010-08-04 美津濃株式会社 繊維強化プラスチック製パイプの設計支援方法
JP2003094559A (ja) * 2001-09-21 2003-04-03 Toray Ind Inc ハイブリッド材の設計方法
JP4581539B2 (ja) 2004-07-29 2010-11-17 横浜ゴム株式会社 構造体の解析モデル作成装置及び構造体の解析モデル作成用コンピュータプログラム
JP5583929B2 (ja) * 2009-06-17 2014-09-03 ポリプラスチックス株式会社 配向状態予測方法及び変形挙動解析方法
GB201003065D0 (en) * 2010-02-23 2010-04-07 Simpleware Ltd Image processing method and method of three-dimensional printing incorporating the same
US9916684B2 (en) * 2016-07-27 2018-03-13 Adobe Systems Incorporated Generating portable three-dimensional print-preview renderings of three-dimensional models
US9978173B2 (en) * 2016-07-27 2018-05-22 Adobe Systems Incorporated Generating views of three-dimensional models illustrating defects
WO2018031594A1 (en) * 2016-08-09 2018-02-15 Arevo, Inc. Systems and methods for structurally analyzing and printing parts
US20180093419A1 (en) * 2016-09-30 2018-04-05 Velo3D, Inc. Three-dimensional objects and their formation
US10885407B2 (en) * 2016-11-23 2021-01-05 Simbionix Ltd. Method and system for three-dimensional print oriented image segmentation
EP3643475A4 (en) * 2017-06-20 2021-07-07 Toray Engineering Co., Ltd. METHOD FOR PREDICTING THE STRENGTH OF A STRUCTURE, METHOD FOR MANUFACTURING A STRUCTURE, METHOD TO SUPPORT THE STRUCTURE OF A STRUCTURE AND PROGRAM
JP2019093701A (ja) * 2017-11-25 2019-06-20 正義 平井 装飾体、装飾体照明設備、装飾体製造装置及び装飾体製造方法
CN109383029A (zh) * 2017-08-04 2019-02-26 三纬国际立体列印科技股份有限公司 立体打印设备以及立体打印方法
CN108381916B (zh) * 2018-02-06 2020-10-27 西安交通大学 一种非接触识别缺损形貌的复合3d打印系统及方法
JP2020017097A (ja) * 2018-07-26 2020-01-30 株式会社Subaru 設計支援システム、設計支援方法及び設計支援プログラム

Also Published As

Publication number Publication date
WO2021186988A1 (ja) 2021-09-23
US20230082432A1 (en) 2023-03-16
US12367572B2 (en) 2025-07-22
JP2021146514A (ja) 2021-09-27

Similar Documents

Publication Publication Date Title
US9697639B2 (en) Three-dimensional model data generation device, method and program
CN110001067B (zh) 一种连续纤维增强复合材料3d打印路径规划方法
US10529126B2 (en) Generating slice data from a voxel representation
CN106373168A (zh) 一种基于医疗图像的分割与三维重建方法、3d打印系统
JPH11108859A5 (ja)
CN100421128C (zh) 用于对断层图像数据分段的方法和图像处理系统
JP2013507679A (ja) 三次元物体モデルの3dプリントが可能な方法及びシステム
Sun et al. Reasoning Boolean operation based modeling for heterogeneous objects
US10914690B2 (en) Method for characterising a part
Vatani et al. An enhanced slicing algorithm using nearest distance analysis for layer manufacturing
JP7282051B2 (ja) 演算装置、演算方法およびプログラム
CN107240103B (zh) 一种基于图像分割的数字体积相关算法中边界处理方法
Kang et al. Accuracy assessment of image-based surface meshing for volumetric computed tomography images in the craniofacial region
Fourrier et al. Automated conformal mesh generation chain for woven composites based on CT-scan images with low contrasts
Warkhedkar et al. Material-solid modeling of human body: A heterogeneous B-spline based approach
Clerici et al. Anisotropic adapted meshes for image segmentation: Application to three-dimensional medical data
EP4266032A1 (en) Method and apparatus for analyzing computer tomography data
CN106097457A (zh) 基于stl缺损模型的切割修整方法
Denos et al. SIMULATIONS OF MECHANICAL PROPERTIES FOR DISCONTINUOUS PREPREG COMPOSITES-075
Ficzere et al. Simplified CAD model of human metacarpal for implantation
Thakur et al. Optimization of Direct Slicing Process Using DICOMS for Additive Manufacturing
Remes et al. CT-scan based mechanical finite element analysis of inter-filament voids in fused deposition modelling
Timpano Investigation of Advanced Material Through Digital Volume Correlation Measurement Techniques
Vänskä Voxel-based modelling in additive manufacturing
Hu et al. Developing and optimizing a finite element model of phalange using CT images

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220916

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: 20230425

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230516

R150 Certificate of patent or registration of utility model

Ref document number: 7282051

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150