JPWO2010021287A1 - 基板の反り予測方法、基板の反り予測システムおよび基板の反り予測プログラム - Google Patents
基板の反り予測方法、基板の反り予測システムおよび基板の反り予測プログラム Download PDFInfo
- Publication number
- JPWO2010021287A1 JPWO2010021287A1 JP2010525675A JP2010525675A JPWO2010021287A1 JP WO2010021287 A1 JPWO2010021287 A1 JP WO2010021287A1 JP 2010525675 A JP2010525675 A JP 2010525675A JP 2010525675 A JP2010525675 A JP 2010525675A JP WO2010021287 A1 JPWO2010021287 A1 JP WO2010021287A1
- Authority
- JP
- Japan
- Prior art keywords
- warpage
- layer
- elastic modulus
- printed wiring
- wiring board
- 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.)
- Pending
Links
Images
Classifications
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05K—PRINTED CIRCUITS; CASINGS OR CONSTRUCTIONAL DETAILS OF ELECTRIC APPARATUS; MANUFACTURE OF ASSEMBLAGES OF ELECTRICAL COMPONENTS
- H05K3/00—Apparatus or processes for manufacturing printed circuits
- H05K3/0005—Apparatus or processes for manufacturing printed circuits for designing circuits by computer
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/18—Chip packaging
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05K—PRINTED CIRCUITS; CASINGS OR CONSTRUCTIONAL DETAILS OF ELECTRIC APPARATUS; MANUFACTURE OF ASSEMBLAGES OF ELECTRICAL COMPONENTS
- H05K2201/00—Indexing scheme relating to printed circuits covered by H05K1/00
- H05K2201/06—Thermal details
- H05K2201/068—Thermal details wherein the coefficient of thermal expansion is important
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05K—PRINTED CIRCUITS; CASINGS OR CONSTRUCTIONAL DETAILS OF ELECTRIC APPARATUS; MANUFACTURE OF ASSEMBLAGES OF ELECTRICAL COMPONENTS
- H05K2201/00—Indexing scheme relating to printed circuits covered by H05K1/00
- H05K2201/09—Shape and layout
- H05K2201/09009—Substrate related
- H05K2201/09136—Means for correcting warpage
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
基板全体の反り予測精度を向上させるとともに、配線パターンに起因して発生する局所変形も高精度に予測する基板の反り予測システムを提供する。多層のプリント配線基板の配線パターン、各層の弾性率と線膨張係数と厚さ、温度プロファイルデータを少なくとも含む解析データを取り込み、前記取り込んだ解析データを基に、前記多層のプリント配線基板の前記各層ごとにメッシングを行って構造体モデルを作成し、作成した前記構造体モデルから、前記プリント配線基板の配線層と樹脂層のペアであって前記配線層と前記樹脂層との接合面での要素の節点が一致せず不整合となっている不整合ペア層を抽出し、抽出した前記不整合ペア層について、導体領域と樹脂領域とを含む配線層の等価弾性率を算出し、算出した前記等価弾性率を基に、各接合面の反り変化を表す反り変化データを算出し、前記構造体モデルと前記反り変化データと温度変化を基に、前記多層のプリント配線基板の反りを算出する。
Description
本発明は、BGA(Ball Grid Array)やCSP(Chip Size Package)等のLSI(Large Scale Integration)パッケージのモールドの際に発生する反りをはじめ、樹脂を使用した電子部品全般の反りを高精度に予測し、かつ、これらの反りを有した電子部品をはんだ材料を用いてプリント配線基板に実装するリフロー工程において、プリント配線基板および電子部品に発生する反りに起因した接続信頼性の低下を防止するための最適条件を求める基板の反り予測方法、基板の反り予測システムおよび基板の反り予測プログラムに関する。特に、本発明は、軽薄短小化の傾向が進み、曲げ剛性が低下して反りの発生が顕著となってきた携帯用電子機器に使用するプリント配線基板のリフロー時における反りを高精度に予測し、設計段階で反り防止対策を実施するための基板の反り予測方法、基板の反り予測システムおよび基板の反り予測プログラムに関する。
近年、電子機器の軽薄短小化が急速に進み、実装技術についてはより高密度かつ高信頼性が強く求められている。また、高機能化のためLSIパッケージは多ピン化し、それに伴ってサイズは大きくかつ薄くなる傾向にあるため、LSIパッケージをリフロー工程ではんだ付けするプリント配線基板側も薄くする必要があり、曲げ剛性の低下に伴って、反りが増大する傾向にある。特に、多層のプリント配線基板においては、樹脂層とこれに挟まれる導体(主に銅)の配線パターン層の熱膨張係数が大きく異なり、このミスマッチに起因した基板全体の反りが問題となる。その上、配線パターンは必ずしも層内で均一に存在していないので、パターンのある部分とない部分とのミスマッチも加わり、局所的な反りも発生するという問題もある。このような原因で発生した反りは、リフロー時におけるはんだとの接続不良を起こす可能性があり、接続信頼性を大きく低下させる原因となっている。さらに、反ったままのプリント配線基板を薄い筐体に組み込むことは難しく、歩留り低下の主要な原因ともなっている。
これらの問題を解決するため、有限要素法(FEM)を用いたシミュレーションが関連する企業、大学、研究機関などで行われるようになってきた。しかし、何層もの異なる配線パターンを全て正確にメッシング(要素切り)すると要素数が莫大となって、大規模な情報処理装置が必要となり、設計段階でフィードバックしながら使用するには程遠いものとなる。このため、特許文献1のように配線パターンをその密度に応じた大きさの領域に分割し、これらの領域の等価弾性率を算出したり、特許文献2のように配線パターンの中で周期性のある領域を抽出し、この領域に均質化処理を加えたりして、配線パターンを考慮しつつ、要素数を低減させる工夫が提案されている。また、非特許文献1では、配線CAD(Computer Aided Design)からのガーバーデータをメッシュデータに変換する際、ガーバーデータを幾つかのグループに分け、各グループ内で導体と樹脂との含有率から、そのグループを導体とするか、樹脂とするかを決めて、要素数を低減させる対策が行われている。また、FEMを使用しない方法も研究されており、非特許文献2の多層ばり理論を拡張させて、非特許文献3のように多層基板の反りを短時間で予測する方法も発表されている。
伊東ら,"プリント配線板の反り解析技術の研究",Mate2006論文集,pp.461−466,2006.
尾田ら,"多層ばり理論によるプリント基板の応力・変形の評価",日本機械学会論文集, 59巻563号,pp.203−208,1993.
平田,"多層基板の反り応力算出ツールの開発および応力極小値発生原因の究明",エレクトロニクス実装学術講演大会論文集,pp.93−94,2005.
有限要素法では通常、接合面の要素は一致(整合)させておく必要がある。このため、配線パターンとの接合面の節点は一致させなければならず、配線層以外の全ての層も同じ細かさの要素で作成する必要がでてくる。FEMには通常は上記のような制約が存在し、プリント配線基板では異なる配線パターンを有する配線層が何層も存在するので、全体としての要素数が膨大になっていた。このため、配線パターンは細かいまま、全体モデルの要素数は低減させるという、相反する要求を満足できる手段はこれまで存在しなかった。
例えば、特許文献1で示されている方法では配線パターンの密度に応じた等価弾性率を使用しているため解析精度が低いという課題があった。また、特許文献2では均質化法を用いているので、特許文献1のように配線パターンの密度に応じた等価弾性率を使用する場合よりも精度は高くなるが、周期性のある配線パターン部にしか使用できないという課題があった。さらに、非特許文献1では、配線CADからのガーバーデータをメッシュデータに変換する際、導体と樹脂との含有率から、配線パターンを導体とするか、樹脂とするかを決めるため、この方法でも解析精度が低下するという課題があった。なお、ここでいう解析精度とは、プリント配線基板の反りの最大と最小の差と、配線パターンに起因する局所的な反りや応力の予測も含んでおり、何らかの方法で配線パターンを簡易化する限り、この両者の予測精度は大幅に低下してしまう。さらに、非特許文献2を発展させた非特許文献3のような有限要素法そのものを使用しない提案もあり、この方法は確かにメッシングを行う必要がなく、計算時間も有限要素法よりも非常に早くなる。しかし、この方法で扱えるのは2次元の多層基板なので、配線パターンを扱うことはできないという課題もあった。
例えば、特許文献1で示されている方法では配線パターンの密度に応じた等価弾性率を使用しているため解析精度が低いという課題があった。また、特許文献2では均質化法を用いているので、特許文献1のように配線パターンの密度に応じた等価弾性率を使用する場合よりも精度は高くなるが、周期性のある配線パターン部にしか使用できないという課題があった。さらに、非特許文献1では、配線CADからのガーバーデータをメッシュデータに変換する際、導体と樹脂との含有率から、配線パターンを導体とするか、樹脂とするかを決めるため、この方法でも解析精度が低下するという課題があった。なお、ここでいう解析精度とは、プリント配線基板の反りの最大と最小の差と、配線パターンに起因する局所的な反りや応力の予測も含んでおり、何らかの方法で配線パターンを簡易化する限り、この両者の予測精度は大幅に低下してしまう。さらに、非特許文献2を発展させた非特許文献3のような有限要素法そのものを使用しない提案もあり、この方法は確かにメッシングを行う必要がなく、計算時間も有限要素法よりも非常に早くなる。しかし、この方法で扱えるのは2次元の多層基板なので、配線パターンを扱うことはできないという課題もあった。
本発明は、上述した課題を解決するためになされたものであり、基板全体の反り予測精度を向上させるとともに、配線パターンに起因して発生する局所変形も高精度に予測する技術を実現し、これらの予測結果を基にプリント配線基板の反り低減方法の究明に設計段階から役立てて、電子機器のさらなる薄型化を可能とする基板の反り予測方法、基板の反り予測システムおよび基板の反り予測プログラムを提供することを目的とする。
上述した課題を解決するために、本発明に係る基板の反り予測方法では、多層のプリント配線基板の配線パターン、各層の弾性率と線膨張係数と厚さ、温度プロファイルデータを少なくとも含む解析データをコンピュータに取り込み、前記コンピュータで、前記取り込んだ解析データを基に、前記多層のプリント配線基板の前記各層ごとにメッシングを行って構造体モデルを作成し、前記コンピュータで、作成した前記構造体モデルから、前記プリント配線基板の配線層と樹脂層のペアであって前記配線層と前記樹脂層との接合面での要素の節点が一致せず不整合となっている不整合ペア層を抽出し、前記コンピュータで、抽出した前記不整合ペア層について、導体領域と樹脂領域とを含む配線層の等価弾性率を算出し、前記コンピュータで、算出した前記等価弾性率を基に、各接合面の反り変化を表す反り変化データを算出し、前記コンピュータで、前記構造体モデルと前記反り変化データと温度変化を基に、前記多層のプリント配線基板の反りを算出する。
また、本発明に係る基板の反り予測システムは、多層のプリント配線基板の配線パターン、各層の弾性率と線膨張係数と厚さ、温度プロファイルデータを少なくとも含む解析データを取り込む解析データ取込部と、前記解析データ取込部が取り込んだ前記解析データを基に、前記多層のプリント配線基板の前記各層ごとにメッシングを行って構造体モデルを作成する構造体モデル作成部と、前記構造体モデル作成部が作成した前記構造体モデルから、前記プリント配線基板の配線層と樹脂層のペアであって前記配線層と前記樹脂層との接合面での要素の節点が一致せず不整合となっている不整合ペア層を抽出する不整合ペア層抽出部と、前記不整合ペア層抽出部が抽出した前記不整合ペア層について、導体領域と樹脂領域とを含む配線層の等価弾性率を算出する等価弾性率算出部と、算出した前記等価弾性率を基に、各接合面の反り変化を表す反り変化データを算出するペア層反り予測部と、前記構造体モデルと前記反り変化データと温度変化を基に、前記多層のプリント配線基板の反りを算出する演算部と、を備えている。
また、本発明に係る基板の反り予測プログラムは、多層のプリント配線基板の配線パターン、各層の弾性率と線膨張係数と厚さ、温度プロファイルデータを少なくとも含む解析データを取り込む手段、前記取り込んだ解析データを基に、前記多層のプリント配線基板の前記各層ごとにメッシングを行って構造体モデルを作成する手段、作成した前記構造体モデルから、前記プリント配線基板の配線層と樹脂層のペアであって前記配線層と前記樹脂層との接合面での要素の節点が一致せず不整合となっている不整合ペア層を抽出する手段、抽出した前記不整合ペア層について、導体領域と樹脂領域とを含む配線層の等価弾性率を求める手段、算出した前記等価弾性率を基に、各接合面の反り変化を表す反り変化データを算出する手段、前記構造体モデルと前記反り変化データと温度変化を基に、前記多層のプリント配線基板の反りを算出する手段、としてコンピュータを機能させる。
本発明では、配線パターンの要素は細かくしたまま、パターンがないその他の樹脂層は比較的粗い要素でモデル化し、要素の節点がずれたままの接合面の反りは反りの理論式の演算結果で置換えるようにし、モデル全体はFEMで解析することができる。したがって、全体の要素数を低減できる。また、各接合面は反りの理論式で変化させるため、解析精度の低下を防ぐことができる。これまでは、配線パターンが原因で局所的な変形が発生すると、その近傍に実装した電子部品の接続部に高い応力が発生して接続信頼性が低下するという課題があった。これに対して、本発明では局所的な変形が高精度に予測できるので、電子部品の配置を変更することで、電子部品の接続信頼性の低下を未然に防ぐことができる。また、大規模な情報処理装置などを使用せず、PC(Personal Computer)レベルのマシンを用いて、プリント配線基板の高精度な反り予測とこの予測結果を基にした反りの効果的な低減対策とを設計上流段階で行うことが可能となる。
また、全体の要素数を低減させることができるため、解析時間が短縮できる。
また、全体の要素数を低減させることができるため、解析時間が短縮できる。
以下、図面を参照して本発明の実施形態を説明する。
図1を参照すると、本発明の第1の実施形態に係る基板の反り予測システムは、入力装置1と、データ処理装置2と、出力装置4と、記憶装置3とから構成される。データ処理装置2は、解析データ取込部21と、構造体モデル作成部22と、不整合ペア層抽出部23と、等価弾性率算出部24と、ペア層反り予測部25と、演算部26とを有する。
図1を参照すると、本発明の第1の実施形態に係る基板の反り予測システムは、入力装置1と、データ処理装置2と、出力装置4と、記憶装置3とから構成される。データ処理装置2は、解析データ取込部21と、構造体モデル作成部22と、不整合ペア層抽出部23と、等価弾性率算出部24と、ペア層反り予測部25と、演算部26とを有する。
次に各部の結びつき、役割、機能について説明する。多層のプリント配線基板の少なくとも、配線パターン、各層の弾性率と線膨張係数と厚さ、温度プロファイルを、データライブラリから入力装置1を介して入力する。または、解析者が、入力装置1を介して、これらのデータを入力する。これにより、解析に必要なデータを、解析データ取込部21が取り込む。次に、構造体モデル作成部22が各層ごとにメッシングを行って構造体モデルを作成する。この際、配線パターンは複雑な形状を有しているため、要素は細かくなり要素数は増加するが、配線パターンの特徴をできるだけ維持できるようにモデル化する。また、殆どが樹脂のみの樹脂層については、配線パターン層のメッシングには捉われずに比較的粗いメッシングを行う。さらに、不整合ペア層抽出部23は、配線パターン層と樹脂層のペアを全て抽出する。等価弾性率算出部24は、各配線パターン層の導体領域と樹脂領域の比率から等価弾性率を算出する。
次に、ペア層反り予測部25は、算出した等価弾性率を反りの理論式に代入し、全体モデルに温度が印加された場合、各接合面の反り変形を理論式の演算結果に置き換えられるようにユーザーサブルーチンを作成してFEMソフトウェアにインストールする。次に、演算部26は、全体モデルと反りの理論式に同一の温度変化を印加し、FEMソフトウェアで反りを算出する。
次に、図1の基板の反り予測システムについて、さらに詳細に説明する。入力装置1は、例えば、キーボードとマウスなどの入力デバイスが用いられる。データ処理装置2は、CPU(中央処理装置)やメモリ等を用いて構成される。記憶装置3は、コンピュータ読み取り可能なデータを記憶するものであり、例えば、磁気ディスク記憶媒体が適用される。出力装置4は、例えば、ディスプレイ装置である。
このように、基板の反り予測システムは、入力装置1と、データ処理装置2と、出力装置4と、記憶装置3とから構成されている。また、データ処理装置2は、解析データ取込部21と、構造体モデル作成部22と、不整合ペア層抽出部23と、等価弾性率算出部24と、ペア層反り予測部25と、演算部26とを含んでいる。
ここで、本発明の各実施形態において特に重要となる等価弾性率算出式およびバイメタル反り算出式、ならびに、両者を連立させる新しい方法に関して説明する。
図4は、不整合ペア層抽出部23で抽出したペア層の模式図を示している。抽出されたペア層は、配線層31と樹脂層32のペア層であり、接合面33での要素の節点は、ずれて一致していない。さらに、配線層31は、導体領域34と樹脂領域35との混合層となっている。配線層31は、このように複数の材料から成っており、この配線層31の等価弾性率は、複合則を用いて
1/E=Vr/Er+Vc/Ec (1)
と表される。ここでEは、等価弾性率、VrとVcは、それぞれ樹脂と導体の体積比率、ErとEcは、それぞれ樹脂と導体の弾性率である。また、バイメタルの反りの理論式は、チモシェンコの理論に従えば、
h=a1+a2, m=a1/a2, n=E1/E2 とおくと、以下に示す式(2)及び式(3)で表される。
図4は、不整合ペア層抽出部23で抽出したペア層の模式図を示している。抽出されたペア層は、配線層31と樹脂層32のペア層であり、接合面33での要素の節点は、ずれて一致していない。さらに、配線層31は、導体領域34と樹脂領域35との混合層となっている。配線層31は、このように複数の材料から成っており、この配線層31の等価弾性率は、複合則を用いて
1/E=Vr/Er+Vc/Ec (1)
と表される。ここでEは、等価弾性率、VrとVcは、それぞれ樹脂と導体の体積比率、ErとEcは、それぞれ樹脂と導体の弾性率である。また、バイメタルの反りの理論式は、チモシェンコの理論に従えば、
h=a1+a2, m=a1/a2, n=E1/E2 とおくと、以下に示す式(2)及び式(3)で表される。
ここで、E1、E2はバイメタルのそれぞれの層の弾性率、a1、a2はそれぞれの層の厚さ、α1、α2はそれぞれの層の線膨張係数、Lはバイメタルの長さ、T0は初期温度、T1は最終温度、ρはプリント配線基板の曲率、Wはプリント配線基板のたわみであり、θはプリント配線基板のたわみを円弧と想定した場合の前記円弧に対する中心角の半分である。半径ρの円弧を考え、弦の長さをLとする。また、中心から弦の中央を通る垂直線を引くと、弦の中央から円弧の垂直線との交点までの長さがいわゆる正矢、すなわち、たわみWとなる。中心角の半分をθとすると、W=ρ−sqrt(ρ2−L2/4)(ただし、sqrtは平方根)となり、これを解くと、W=ρ(1−cosθ)≒ρθ2/2となる。さらに、L/2≒ρθであるから、式(3)を得る。
上記の式(1)〜(3)は既知の式であるが、これらをプリント配線基板の抽出したペア層の変形予測に使用する。すなわち、等価弾性率を求める式(1)を用い、導体領域34と樹脂領域35との混合層となっている配線層31の等価弾性率Eを求める。次に、この等価弾性率Eをバイメタル式(2)における配線パターン層(ここでは2層目)の弾性率E2に代入する。すなわち、
n=E1/E (4)
で表される式(4)を式(2)に代入する。そして、この関係式を、配線層31とこの配線層31とペアとなっている樹脂層32との接合面33の変化として用い、温度変化を全体モデルと式(2)に同時に印加する。これにより、精度低下の原因となる接合部分は理論式で高精度に変形し、その変形を利用しながらFEMでプリント配線基板全体の反りを解析することができる。したがって、配線パターンの形状を残しながら、樹脂層は粗い要素でモデル化することが可能となり、全体としての要素数の低減が実現できる。また、解析精度の低下原因となる不整合な接合面は理論式に沿って変形させることができ、高精度な反り解析が可能となる。
上記の式(1)〜(3)は既知の式であるが、これらをプリント配線基板の抽出したペア層の変形予測に使用する。すなわち、等価弾性率を求める式(1)を用い、導体領域34と樹脂領域35との混合層となっている配線層31の等価弾性率Eを求める。次に、この等価弾性率Eをバイメタル式(2)における配線パターン層(ここでは2層目)の弾性率E2に代入する。すなわち、
n=E1/E (4)
で表される式(4)を式(2)に代入する。そして、この関係式を、配線層31とこの配線層31とペアとなっている樹脂層32との接合面33の変化として用い、温度変化を全体モデルと式(2)に同時に印加する。これにより、精度低下の原因となる接合部分は理論式で高精度に変形し、その変形を利用しながらFEMでプリント配線基板全体の反りを解析することができる。したがって、配線パターンの形状を残しながら、樹脂層は粗い要素でモデル化することが可能となり、全体としての要素数の低減が実現できる。また、解析精度の低下原因となる不整合な接合面は理論式に沿って変形させることができ、高精度な反り解析が可能となる。
次に、データ処理装置2の構成について、さらに説明する。
解析データ取込部21は、入力装置1から入力されるプリント配線基板の少なくとも配線パターンデータと材料特性値、具体的には各層の弾性率、線膨張係数、さらに、形状データとして各層の厚さ、モデルに印加する温度プロファイルなどのデータを取り込む。
構造体モデル作成部22は、形状データを基にメッシングを行う。ここでは、配線パターンに起因するプリント配線基板の局所的な反り及び変形も高精度に予測することを目指し、配線層はパターンの特徴が十分残る程度に要素を細かくし、一方、それ程細かくする必要のない樹脂層は粗く要素分割して、後述する方法で接合させる。そのために、構造体モデル作成部22は、各層と接合された相手側(すなわち、配線層に着目した場合は樹脂層、樹脂層に着目した場合は配線層)の要素分割状態に捉われず、配線層は細かく、樹脂層は粗く、を基本方針としてメッシングする。
解析データ取込部21は、入力装置1から入力されるプリント配線基板の少なくとも配線パターンデータと材料特性値、具体的には各層の弾性率、線膨張係数、さらに、形状データとして各層の厚さ、モデルに印加する温度プロファイルなどのデータを取り込む。
構造体モデル作成部22は、形状データを基にメッシングを行う。ここでは、配線パターンに起因するプリント配線基板の局所的な反り及び変形も高精度に予測することを目指し、配線層はパターンの特徴が十分残る程度に要素を細かくし、一方、それ程細かくする必要のない樹脂層は粗く要素分割して、後述する方法で接合させる。そのために、構造体モデル作成部22は、各層と接合された相手側(すなわち、配線層に着目した場合は樹脂層、樹脂層に着目した場合は配線層)の要素分割状態に捉われず、配線層は細かく、樹脂層は粗く、を基本方針としてメッシングする。
メッシングが終了したプリント配線基板モデルの各接合面は不整合となっているので、不整合ペア層抽出部23は、この不整合境界を有するペア(2層の組)を全て抽出する。等価弾性率算出部24は、各ペア層の配線層側において、導体領域と樹脂領域の弾性率と体積比率を式(1)に代入して各配線層の等価弾性率を求める。ペア層反り予測部25は、等価弾性率算出部24で求めた配線層側の等価弾性率を、バイメタル式(2)の配線層側の弾性率として代入する。具体的には、式(4)を式(2)に代入する。演算前のこの時点では、印加する温度のみが未定数として残り、他のパラメータは全て決定することができる。この式(2)をバイメタルの反りの式(3)と連成させ、ペア層の接合界面の反りの変化をこの理論式の演算結果に置換えるユーザーサブルーチンを作成して、FEMソフトウェアにインストールする。
演算部26は、メッシングしたプリント配線基板の全体モデルに温度変化を与えるとともに、サブルーチンとして組み込んだバイメタルの理論式(2)にも同一温度を同一時間ステップで与え、この結果を式(3)に入力して接合面の温度に起因した反り変化を算出する。FEMソフトウェアは、この接合面の変化を利用して全体の反りと局所的変形を算出する。
以上のように、解析データ取込部21は、基板の配線パターン、各層の弾性率と線膨張係数と厚さ、温度プロファイルなどのデータを入力装置1から取り込む。構造体モデル作成部22は、取り込んだ解析データを基に構造体モデルを作成する。不整合ペア層抽出部23は、作成したモデルから多層プリント配線基板の各接合ペア層を抽出する。等価弾性率算出部24は、抽出した接合ペア層の内、配線と樹脂を含む複合材料から成る配線層側の等価弾性率を算出する。ペア層反り予測部25は、算出した等価弾性率をバイメタルの反り理論式に代入して、接合面の反り変化をバイメタル式の演算結果で置き換える。演算部26は、バイメタル式とFEMモデル全体とに同一温度変化(温度プロファイル)を印加してFEMで多層プリント配線基板の反りと局所変形を算出する。
上記各部が行う処理を順次実施することで、解析精度の低下要因となる各接合面は精度の高い反り理論式で変形させながら、モデル全体はFEM解析することが可能となり、解析精度を保ちつつ、要素数の削減も実現できる。このため、設計上流段階で、PCレベルの解析マシンを用いて、設計段階のプリント配線基板のリフローにおける反り予測を高精度に行うことが可能となる。さらに、この予測結果から、材料特性、層厚さなどを調整することで、反りの低減が可能となる。
さらに、解析モデルに電子部品モデルを追加すれば、配線パターンに起因する局所変形のために接続部の応力・ひずみが高くなる部品を予測でき、配置の変更を行って接続信頼性の低下を未然に防止できるという特徴がある。
さらに、解析モデルに電子部品モデルを追加すれば、配線パターンに起因する局所変形のために接続部の応力・ひずみが高くなる部品を予測でき、配置の変更を行って接続信頼性の低下を未然に防止できるという特徴がある。
次に、本発明の第2の実施形態について図面を参照して詳細に説明する。
図2は、本発明の第2の実施形態に係る基板の反り予測システムの構成を示すブロック図である。図2を参照すると、本発明の第2の実施形態に係る反り予測システムは、本発明の第1の実施形態に係る反り予測システムと同様に、入力装置1、データ処理装置2(図1のデータ処理装置2と同様の構成)、記憶装置3、出力装置4を備える。また、図2において、等価弾性率算出プログラム・バイメタル反り算出プログラム5は、データ処理装置2に読み込まれたデータの処理、および、データ処理装置2を構成する各部の動作を制御し、データ処理装置2を構成する各部における処理結果を記憶装置3に記憶する。
以上のように、データ処理装置2は、等価弾性率算出プログラム・バイメタル反り算出プログラム5の制御に従って、第1の実施形態におけるデータ処理装置2が行う処理と同一の処理を実行する。
図2は、本発明の第2の実施形態に係る基板の反り予測システムの構成を示すブロック図である。図2を参照すると、本発明の第2の実施形態に係る反り予測システムは、本発明の第1の実施形態に係る反り予測システムと同様に、入力装置1、データ処理装置2(図1のデータ処理装置2と同様の構成)、記憶装置3、出力装置4を備える。また、図2において、等価弾性率算出プログラム・バイメタル反り算出プログラム5は、データ処理装置2に読み込まれたデータの処理、および、データ処理装置2を構成する各部の動作を制御し、データ処理装置2を構成する各部における処理結果を記憶装置3に記憶する。
以上のように、データ処理装置2は、等価弾性率算出プログラム・バイメタル反り算出プログラム5の制御に従って、第1の実施形態におけるデータ処理装置2が行う処理と同一の処理を実行する。
次に、データ処理装置2の動作(特に、不整合ペア層抽出部23、等価弾性率算出部24、ペア層反り予測部25、及び、演算部26)について、図3のフローチャートを用いて説明する。
ペア層の番号を示すパラメータiを0に初期化した後(ステップA1)、ペア層番号に1を加算する工程を経て(ステップA2)、加算後のペア層番号が総ペア層数Nを越えたか否か判断する(ステップA3)。加算後のペア層番号が総ペア層数Nを越えたならば処理をステップA8に移行し、越えていない場合は処理をステップA4へ移行する。
ステップA4では、現在のペア層番号に対応する不整合ペア層を抽出し、ステップA5において、導体領域と樹脂領域とを含む配線層側について式(1)を用いて等価弾性率を算出する。
さらに、ステップA6では、バイメタル式(2)に式(4)の関係式を代入する。すなわち、ステップA5で求めた等価弾性率をバイメタル式(2)の配線層側の弾性率として代入する。ステップA7では、印加温度のみが未定な式(2)と、これを用いて反りを算出する式(3)に従って現在のペア層の接合面が変化するようユーザーサブルーチンを作成して、FEMソフトウェアにインストールする。
ペア層の番号を示すパラメータiを0に初期化した後(ステップA1)、ペア層番号に1を加算する工程を経て(ステップA2)、加算後のペア層番号が総ペア層数Nを越えたか否か判断する(ステップA3)。加算後のペア層番号が総ペア層数Nを越えたならば処理をステップA8に移行し、越えていない場合は処理をステップA4へ移行する。
ステップA4では、現在のペア層番号に対応する不整合ペア層を抽出し、ステップA5において、導体領域と樹脂領域とを含む配線層側について式(1)を用いて等価弾性率を算出する。
さらに、ステップA6では、バイメタル式(2)に式(4)の関係式を代入する。すなわち、ステップA5で求めた等価弾性率をバイメタル式(2)の配線層側の弾性率として代入する。ステップA7では、印加温度のみが未定な式(2)と、これを用いて反りを算出する式(3)に従って現在のペア層の接合面が変化するようユーザーサブルーチンを作成して、FEMソフトウェアにインストールする。
次に、ステップA8では、バイメタル式(2)、(3)とFEMモデル全体に同一温度変化を印加する。ステップA9では、FEM解析を用いて反りを演算する。この際、解析精度の低下要因となる各接合面は、精度の高い反り理論式で変形させることができ、モデル全体はFEM解析することになる。
このように、本実施形態における特徴は、ステップA5で導体領域と樹脂領域とを含む配線層の等価弾性率を理論式(1)で算出し、ステップA6において、バイメタル式(2)にステップA5で求めた等価弾性率を代入して各ペア層の反り変化を理論式で算出し、ステップA7で各接合面の反り変化を前記理論式で算出した結果に置換え、ステップA8でバイメタル式(2)とFEMモデル全体に同一温度変化を印加し、それにより、理論式で算出した接合面の反りの結果を用いながらFEM解析を行えることにある。
弾性率などの材料特性がクリープなどの特殊な挙動を示す場合、ユーザーサブルーチンでプログラム化してFEMソフトウェアに組み込んで解析することはこれまで一般に行われていた。しかし、複数の材料から成る層の等価弾性率を算出し、この等価弾性率の値を用いた理論式から得られる演算結果を接合面の反り変化と見做しながらFEMで全体の反り解析を行うという例は今まで見られていない。なお、もしステップA5〜A8がなければ、各層は接合すべき境界面が不整合のままなので、面内方向の伸縮だけが発生し、反りを求めることはできなくなる。
また、何らかの方法で接合面を接合できたとしても、図5に示す関連する技術に係る配線層と樹脂層から成るバイメタルの変形後の模式図のように、FEM解析では、要素が粗い樹脂層の反りの変化に要素が細かい配線層の変化が従属してしまう。このため、高精度な解析は期待できない。
なお、他の実施形態として、バイメタルの式の代りに、スプライン曲線などの曲線、曲面方程式を適用しても同様な効果を得ることができる。
なお、他の実施形態として、バイメタルの式の代りに、スプライン曲線などの曲線、曲面方程式を適用しても同様な効果を得ることができる。
上述した実施形態は、携帯電子機器を短期間に開発していく必要のあるプリント基板設計部門において、LSIパッケージをリフローを用いてプリント配線基板に実装する際の基板の反り、局所変形を高精度に予測し、反りや局所変形の低減対策を検討するといった用途に適用できる。
なお、上記で説明した各処理ステップを実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより、基板の反り予測に係る上述した種々の処理を行ってもよい。
ここでいうコンピュータシステムとは、OS(オペレーティングシステム)や周辺機器等のハードウェアを含むものであってもよい。また、コンピュータシステムは、WWW(ワールドワイドウェブ)システムを利用している場合であれば、ホームページ提供環境(あるいは表示環境)も含むものとする。
コンピュータ読み取り可能な記録媒体とは、フレキシブルディスク、光磁気ディスク、ROM(Read Only Memory)、フラッシュメモリ等の書き込み可能な不揮発性メモリ、CD(コンパクトディスク)−ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。また、コンピュータ読み取り可能な記録媒体とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ(例えばDRAM(Dynamic Random Access Memory))のように、一定時間プログラムを保持しているものも含むものとする。
上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介してあるいは伝送媒体中の伝送波により、他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する伝送媒体は、インターネット等のネットワーク(通信網)や電話回線等の通信回線(通信線)のように情報を伝送する機能を有する媒体のことをいう。上記プログラムは、前述した機能の一部を実現するためのものであっても良い。また、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル(差分プログラム)であっても良い。
以上、実施形態を参照して本発明を説明したが、本発明は上述した実施形態に限定されるものではない。本発明の構成や詳細には、本発明の範囲内で当業者が理解し得る様々な変更を加えることができる。
なお、本発明において、複合則を用いて、配線層の等価弾性率を導体領域と樹脂領域の比率から算出するようにしても良い。
また、有限要素法(FEM)を用いた多層基板のモデル化において、材料特性が異なる層間の境界面要素の節点位置を整合させずにメッシングし、外力が加えられた時の不整合な2面の形状変化を理論式の演算結果と一致させながら他の部分はFEMで演算し、複数材料から成る層については、等価弾性率を用いるようにしても良い。
また、ペア層の番号を示すパラメータを初期化する第1の工程と、ペア層番号に1を加算する第2の工程と、加算後のペア層番号が総ペア層数を越えたか否か判断し、総ペア層数を越えたら後述する第8の工程に処理を引渡し、越えていない場合は続く第4の工程へ処理を引き渡す第3の工程と、現在のペア層番号に対応する不整合ペア層を抽出する第4の工程と、導体領域と樹脂領域とを含む配線層の等価弾性率を求める第5の工程と、接合面の反りを算出する反りの理論式に前記等価弾性率を代入する第6の工程と、現在のペア層の接合面の反り変形を前記反りの理論式の演算結果に置き換える第7の工程と、前記反りの理論式とFEMモデル全体に同一温度変化を印加する第8の工程と、FEM解析を用いて反りを演算する第9の工程とを有するものであっても良い。
また、ペア層の番号を示すパラメータを初期化する第1の工程と、ペア層番号に1を加算する第2の工程と、加算後のペア層番号が総ペア層数を越えたか否か判断し、総ペア層数を越えたら後述する第8の工程に処理を引渡し、越えていない場合は続く第4の工程へ処理を引き渡す第3の工程と、現在のペア層番号に対応する不整合ペア層を抽出する第4の工程と、導体領域と樹脂領域とを含む配線層の等価弾性率を求める第5の工程と、接合面の反りを算出する反りの理論式に前記等価弾性率を代入する第6の工程と、現在のペア層の接合面の反り変形を前記反りの理論式の演算結果に置き換える第7の工程と、前記反りの理論式とFEMモデル全体に同一温度変化を印加する第8の工程と、FEM解析を用いて反りを演算する第9の工程とを有するものであっても良い。
また、前記第5の工程において、ペア層の内、導体領域と樹脂領域とを含む配線層の等価弾性率を求めるため、導体領域と樹脂領域の比率から配線層の等価弾性率を算出する複合則を用いるようにしても良い。
また、前記第6の工程において、前記第5の工程で算出した配線層の等価弾性率を反りの理論式に代入し、かつ、理論式としてバイメタルの反りの式を用いるようにしても良い。
また、前記第7の工程において、各ペア層に温度が印加された場合、それぞれのペア層の接合面の反りの変化をバイタルの反りの理論式から求める演算結果に置き換えるようにしても良い。
また、前記第8の工程において、前記第7の工程で各接合面の反り変化を算出するのに使用するバイメタルの反りの理論式と、プリント配線基板全体モデルの両者に対し、時刻歴で変化する温度プロファイルを同時に印加するようにしても良い。
また、前記第9の工程において、各接合面の反りは反りの理論式に従い変化させつつ、同時にプリント配線基板の反りをFEM解析するようにしても良い。
また、前記第6の工程において、接合面の反り変形を高精度に予測する方法として、前記第5の工程を行わず、バイメタルの反り理論式に代えて、スプライン曲線などの曲線、曲面方程式を用いるようにしても良い。
また、入力装置から送られる少なくとも多層のプリント配線基板の配線パターン、各層の弾性率と線膨張係数と厚さ、温度プロファイルデータを取り込む解析データ取込部と、各層ごとにメッシングを行う構造体モデル作成部と、配線層と樹脂層のペアを全て抽出する不整合ペア層抽出部と、各配線層の導体領域と樹脂領域の比率から等価弾性率を算出する等価弾性率算出部と、算出した等価弾性率を反りの理論式に代入し、各接合面の反り変化を理論式の演算結果に置き換えるペア層反り予測部と、全体モデルと反りの理論式に同一の温度変化を印加し、反りを算出する演算部とを備えるようにしても良い。
また、入力装置から渡される形状モデルデータを処理するデータ処理装置と、前記データ処理した結果を表示する出力装置と、データ処理結果を記録する記憶装置とを備えるようにしても良い。
この出願は、2008年8月18日に出願された日本出願特願2008−209809号を基礎とする優先権を主張し、その開示の全てをここに取り込む。
本発明では、配線パターンの要素は細かくしたまま、全体の要素数を低減できる。また、各接合面は反りの理論式で変化させることで解析精度の低下を防ぐことができる。また、局所的な変形が高精度に予測できるので、電子部品の配置を変更することで、電子部品の接続信頼性の低下を未然に防ぐことができる。また、大規模な情報処理装置などを使用せず、PCレベルのマシンを用いて、プリント配線基板の高精度な反り予測とこの予測結果を基にした反りの効果的な低減対策とを設計上流段階で行うことが可能となる。また、全体の要素数を低減させることができるため、解析時間が短縮できる。
1 入力装置
2 データ処理装置
3 記憶装置
4 出力装置
5 等価弾性率算出プログラム・バイメタル反り算出プログラム
21 解析データ取込部
22 構造体モデル作成部
23 不整合ペア層抽出部
24 等価弾性率算出部
25 ペア層反り予測部
26 演算部
31 配線層
32 樹脂層
33 接合面
34 導体領域
35 樹脂領域
E 導体と樹脂との混合層の等価弾性率
Vc 導体の体積比率
Vr 樹脂の体積比率
Ec 導体の弾性率
Er 樹脂の弾性率
a1、a2 バイメタルの各層の厚さ
α1、α2 バイメタルの各層の線膨張係数
L バイメタルの長さ
T0 初期温度
T1 最終温度
ρ プリント配線基板の曲率
W プリント配線基板のたわみ
A1 ペア層の番号を示すパラメータを初期化する工程
A2 ペア層番号に1を加算する工程
A3 加算後のペア層番号に対して、総ペア層数を越えたか否か判断し、総ペア層数を越えたら工程A8に処理を引渡し、越えていない場合は次工程A4へ処理を引き渡す工程
A4 現在のペア層番号に対応する不整合ペア層を抽出する工程
A5 配線層側の等価弾性率を算出する工程
A6 バイメタル式に等価弾性率を代入する工程
A7 接合面の反り変化をバイメタル式の演算結果に置換えする工程
A8 バイメタル式とFEMモデル全体に同一温度変化を印加する工程
A9 FEM解析を用いて反りを演算する工程
2 データ処理装置
3 記憶装置
4 出力装置
5 等価弾性率算出プログラム・バイメタル反り算出プログラム
21 解析データ取込部
22 構造体モデル作成部
23 不整合ペア層抽出部
24 等価弾性率算出部
25 ペア層反り予測部
26 演算部
31 配線層
32 樹脂層
33 接合面
34 導体領域
35 樹脂領域
E 導体と樹脂との混合層の等価弾性率
Vc 導体の体積比率
Vr 樹脂の体積比率
Ec 導体の弾性率
Er 樹脂の弾性率
a1、a2 バイメタルの各層の厚さ
α1、α2 バイメタルの各層の線膨張係数
L バイメタルの長さ
T0 初期温度
T1 最終温度
ρ プリント配線基板の曲率
W プリント配線基板のたわみ
A1 ペア層の番号を示すパラメータを初期化する工程
A2 ペア層番号に1を加算する工程
A3 加算後のペア層番号に対して、総ペア層数を越えたか否か判断し、総ペア層数を越えたら工程A8に処理を引渡し、越えていない場合は次工程A4へ処理を引き渡す工程
A4 現在のペア層番号に対応する不整合ペア層を抽出する工程
A5 配線層側の等価弾性率を算出する工程
A6 バイメタル式に等価弾性率を代入する工程
A7 接合面の反り変化をバイメタル式の演算結果に置換えする工程
A8 バイメタル式とFEMモデル全体に同一温度変化を印加する工程
A9 FEM解析を用いて反りを演算する工程
Claims (9)
- 多層のプリント配線基板の配線パターン、各層の弾性率と線膨張係数と厚さ、温度プロファイルデータを少なくとも含む解析データをコンピュータに取り込み、
前記コンピュータで、前記取り込んだ解析データを基に、前記多層のプリント配線基板の前記各層ごとにメッシングを行って構造体モデルを作成し、
前記コンピュータで、作成した前記構造体モデルから、前記プリント配線基板の配線層と樹脂層のペアであって前記配線層と前記樹脂層との接合面での要素の節点が一致せず不整合となっている不整合ペア層を抽出し、
前記コンピュータで、抽出した前記不整合ペア層について、導体領域と樹脂領域とを含む配線層の等価弾性率を算出し、
前記コンピュータで、算出した前記等価弾性率を基に、各接合面の反り変化を表す反り変化データを算出し、
前記コンピュータで、前記構造体モデルと前記反り変化データと温度変化を基に、前記多層のプリント配線基板の反りを算出する基板の反り予測方法。 - 複合則を用いて、前記配線層の前記等価弾性率を前記導体領域と前記樹脂領域の比率から算出する請求項1記載の基板の反り予測方法。
- 前記配線層を前記樹脂層に比べて細かくメッシングする請求項1または請求項2に記載の基板の反り予測方法。
- 前記等価弾性率を前記配線層の弾性率として用いた反りの理論式と前記解析データとに基づいた演算結果を前記反り変化データとして算出する請求項1〜3のいずれかに記載の基板の反り予測方法。
- 前記反りの論理式は、チモシェンコの理論に従ったバイメタルの理論式である請求項4に記載の基板の反り予測方法。
- 前記反りの論理式は、曲線方程式または曲面方程式である請求項4に記載の基板の反り予測方法。
- 前記構造体モデルと前記反りの理論式に同一の温度変化を同時に適用し、前記反り変化データを用いた有限要素法で前記反りを算出する請求項4〜6のいずれかに記載の基板の反り予測方法。
- 多層のプリント配線基板の配線パターン、各層の弾性率と線膨張係数と厚さ、温度プロファイルデータを少なくとも含む解析データを取り込む解析データ取込部と、
前記解析データ取込部が取り込んだ前記解析データを基に、前記多層のプリント配線基板の前記各層ごとにメッシングを行って構造体モデルを作成する構造体モデル作成部と、
前記構造体モデル作成部が作成した前記構造体モデルから、前記プリント配線基板の配線層と樹脂層のペアであって前記配線層と前記樹脂層との接合面での要素の節点が一致せず不整合となっている不整合ペア層を抽出する不整合ペア層抽出部と、
前記不整合ペア層抽出部が抽出した前記不整合ペア層について、導体領域と樹脂領域とを含む配線層の等価弾性率を算出する等価弾性率算出部と、
算出した前記等価弾性率を基に、各接合面の反り変化を表す反り変化データを算出するペア層反り予測部と、
前記構造体モデルと前記反り変化データと温度変化を基に、前記多層のプリント配線基板の反りを算出する演算部と、
を備えた基板の反り予測システム。 - 多層のプリント配線基板の配線パターン、各層の弾性率と線膨張係数と厚さ、温度プロファイルデータを少なくとも含む解析データを取り込む手段、
前記取り込んだ解析データを基に、前記多層のプリント配線基板の前記各層ごとにメッシングを行って構造体モデルを作成する手段、
作成した前記構造体モデルから、前記プリント配線基板の配線層と樹脂層のペアであって前記配線層と前記樹脂層との接合面での要素の節点が一致せず不整合となっている不整合ペア層を抽出する手段、
抽出した前記不整合ペア層について、導体領域と樹脂領域とを含む配線層の等価弾性率を求める手段、
算出した前記等価弾性率を基に、各接合面の反り変化を表す反り変化データを算出する手段、
前記構造体モデルと前記反り変化データと温度変化を基に、前記多層のプリント配線基板の反りを算出する手段、
としてコンピュータを機能させるための基板の反り予測プログラム。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008209809 | 2008-08-18 | ||
JP2008209809 | 2008-08-18 | ||
PCT/JP2009/064263 WO2010021287A1 (ja) | 2008-08-18 | 2009-08-12 | 基板の反り予測方法、基板の反り予測システムおよび基板の反り予測プログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JPWO2010021287A1 true JPWO2010021287A1 (ja) | 2012-01-26 |
Family
ID=41707165
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010525675A Pending JPWO2010021287A1 (ja) | 2008-08-18 | 2009-08-12 | 基板の反り予測方法、基板の反り予測システムおよび基板の反り予測プログラム |
Country Status (2)
Country | Link |
---|---|
JP (1) | JPWO2010021287A1 (ja) |
WO (1) | WO2010021287A1 (ja) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5637820B2 (ja) * | 2010-11-08 | 2014-12-10 | 富士通株式会社 | 実装基板解析装置,実装基板解析方法および実装基板解析プログラム |
US9247636B2 (en) | 2013-03-12 | 2016-01-26 | International Business Machines Corporation | Area array device connection structures with complimentary warp characteristics |
CN105512360B (zh) * | 2015-11-24 | 2018-11-23 | 广州兴森快捷电路科技有限公司 | 预测封装基板蚀刻后翘曲的方法 |
KR20220015244A (ko) | 2020-07-30 | 2022-02-08 | 삼성전자주식회사 | 인터포저의 워피지 예측 방법 및 이를 사용한 반도체 패키지의 제조 방법 |
CN113068326B (zh) * | 2021-03-29 | 2022-09-30 | 北京小米移动软件有限公司 | 一种焊接质量处理方法及装置、电路板 |
CN115392142B (zh) * | 2022-10-31 | 2023-04-18 | 深圳市城市交通规划设计研究中心股份有限公司 | 沿海环境简支梁弹性模量预测方法、电子设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004013437A (ja) * | 2002-06-05 | 2004-01-15 | Toshiba Corp | 基板の反り解析方法及びそのシステム、基板の反り解析プログラム |
JP2006053747A (ja) * | 2004-08-11 | 2006-02-23 | Fujitsu Ltd | 基板変形予測装置、基板変形予測方法および基板変形予測プログラム |
JP2007080111A (ja) * | 2005-09-16 | 2007-03-29 | Matsushita Electric Ind Co Ltd | 板状体解析方法 |
JP2007199961A (ja) * | 2006-01-25 | 2007-08-09 | Nec Corp | 有限要素法解析モデルの解析方法、解析システム、及び解析プログラム |
-
2009
- 2009-08-12 JP JP2010525675A patent/JPWO2010021287A1/ja active Pending
- 2009-08-12 WO PCT/JP2009/064263 patent/WO2010021287A1/ja active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004013437A (ja) * | 2002-06-05 | 2004-01-15 | Toshiba Corp | 基板の反り解析方法及びそのシステム、基板の反り解析プログラム |
JP2006053747A (ja) * | 2004-08-11 | 2006-02-23 | Fujitsu Ltd | 基板変形予測装置、基板変形予測方法および基板変形予測プログラム |
JP2007080111A (ja) * | 2005-09-16 | 2007-03-29 | Matsushita Electric Ind Co Ltd | 板状体解析方法 |
JP2007199961A (ja) * | 2006-01-25 | 2007-08-09 | Nec Corp | 有限要素法解析モデルの解析方法、解析システム、及び解析プログラム |
Non-Patent Citations (1)
Title |
---|
JPN6013030540; 尾田十八 他: '多層ばり理論によるプリント基板の応力・変形の評価' 日本機械学会論文集 第59巻,第563号, 19930725, p.203-208 * |
Also Published As
Publication number | Publication date |
---|---|
WO2010021287A1 (ja) | 2010-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Darveaux | Effect of simulation methodology on solder joint crack growth correlation and fatigue life prediction | |
KR101142881B1 (ko) | 해석 장치, 해석 방법 및 해석 프로그램을 기록한 컴퓨터로 판독가능한 기록 매체 | |
WO2010021287A1 (ja) | 基板の反り予測方法、基板の反り予測システムおよび基板の反り予測プログラム | |
WO2005088483A1 (ja) | 電子パッケージ信頼性予測装置および電子パッケージ信頼性予測プログラム | |
Syed | Updated life prediction models for solder joints with removal of modeling assumptions and effect of constitutive equations | |
JP4774810B2 (ja) | 基板反り解析方法そのシステム、プログラム及び記録媒体 | |
JP2007199961A (ja) | 有限要素法解析モデルの解析方法、解析システム、及び解析プログラム | |
Vandevelde et al. | Parameterized modeling of thermomechanical reliability for CSP assemblies | |
JP2011039741A (ja) | 基板反り解析方法、基板反り解析システム、および基板反り解析プログラム | |
JPWO2009044699A1 (ja) | 3次元基板反り解析システム、3次元基板反り解析装置、3次元基板反り解析方法およびプログラム。 | |
Ren et al. | Shell-based simplified electronic package model development and its application for reliability analysis | |
JPWO2007086120A1 (ja) | 情報処理装置、シミュレーション方法、情報処理プログラム | |
JP5024636B2 (ja) | 基板又は電子部品の反り解析方法、基板又は電子部品の反り解析システム及び基板又は電子部品の反り解析プログラム | |
Upreti et al. | Stochastic Finite Element Thermal Analysis of a Ball Grid Array Package | |
JP2007027152A (ja) | プリント基板の設計方法 | |
JP5071081B2 (ja) | 半導体装置に対する同時動作信号ノイズ見積り方法における同時動作信号ノイズ基礎特性取得方法、及びプログラム | |
US20110010153A1 (en) | Model analysis system, model analysis method, and model analysis program | |
JP2006278803A (ja) | 基板反り解析方法およびそのシステム、基板の反り解析プログラム。 | |
Ren et al. | Application of ABAQUS/Explicit submodeling technique in drop simulation of system assembly | |
JP2011159870A (ja) | 設計条件算出装置、基板の反りと応力の最適化方法、及び基板の反りと応力の最適化プログラム | |
JP5040363B2 (ja) | 解析装置、解析方法及び解析プログラム | |
JP2010033493A (ja) | レイアウト配線混雑予測装置およびその方法、並びにプログラム | |
JP4876777B2 (ja) | 有限要素法解析モデルの要素分割システム、要素分割装置及び要素分割方法並びにそのプログラム | |
Hossain et al. | Design optimization and reliability of PWB level electronic package | |
Davuluri et al. | Thermomechanical durability of high i/o bga packages |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20120713 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130625 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20131022 |