(57)【要約】【課題】拡散反射と鏡面反射を伴う光拡散性の強い光反射板を高精度に模擬可能な光強度分布シミュレーション装置、シミュレーション方法、該方法を実行させるプログラムおよび該プログラムを格納した記録媒体を提供する。【解決手段】光反射板模擬部１２２では、鏡面反射成分と拡散反射成分とを分離して模擬する構成としており、それぞれの反射成分を算出するための鏡面／拡散反射成分算出手段１２２ａを備えている。また、鏡面／拡散反射成分算出手段１２２ａは、鏡面反射成分を複数のピークフィッティング関数を用いてフィッティングしたものを用いるように構成されており、これにより鏡面反射成分を高精度に評価することが可能となっている。【選択図】図１
2010020528.tif 000002
【特許請求の範囲】
【請求項１】
光源の一方に光反射板を配置し他方に受光面を配置した構成において、前記光源から出射された光子が前記受光面に前記光反射板で反射されて入射するかあるいは前記光源から直接入射する個数から光強度分布を算出する光反射板の光強度分布シミュレーション装置であって、
前記光子が前記光反射板に入射したときにその反射方向を決定する光反射板模擬部と、
前記光子が前記光源に入射したときにその吸収の判定または反射方向の決定を行う光源模擬部と、
前記光子を発生させ、前記光子が前記光反射板に入射する場合には前記光反射板模擬部を用いて前記光子の反射方向を決定し、前記光源に入射する前記光子を前記光源模擬部を用いて前記光子の吸収または反射方向を決定し、前記受光面に入射する前記光子の個数から前記光強度分布を算出する追跡計算を、前記光子の発生回数が所定数に達するまで行う光子追跡計算部と、を備え、
前記光反射板模擬部は、前記光反射板に入射した前記光子が鏡面反射するか拡散反射するかを所定の確率で決定している
ことを特徴とする光反射板の光強度分布シミュレーション装置。
【請求項２】
前記光反射板模擬部は、前記鏡面反射による光強度分布を２以上の所定の関数を用いて模擬し、前記光強度分布に基づいて反射する光子の反射方向の決定を行う
ことを特徴とする請求項１に記載の光反射板の光強度分布シミュレーション装置。
【請求項３】
前記所定の関数は、ガウス関数、ルジャンドル関数、ピアソン関数のいずれか１つ、あるいは２以上を組み合わせた関数である
ことを特徴とする請求項２に記載の光反射板の光強度分布シミュレーション装置。
【請求項４】
前記鏡面反射するか拡散反射するかの確率は、前記光子の入射角に対する前記鏡面反射の反射率がフレネル則に従うように決定されている
ことを特徴とする請求項１乃至３に記載の光反射板の光強度分布シミュレーション装置。
【請求項５】
前記光反射板は、熱可塑性樹脂発泡体光反射板である
ことを特徴とする請求項１乃至４のいずれか１項に記載の光反射板の光強度分布シミュレーション装置。
【請求項６】
光源の一方に光反射板を配置し他方に受光面を配置した構成において、前記光源から出射された光子が前記受光面に前記光反射板で反射されて入射するかあるいは前記光源から直接入射する個数から光強度分布を算出する光反射板の光強度分布シミュレーション方法であって、
前記光源から光子を出射させる光子発生ステップと、
前記光子が前記光反射板に入射したときにその反射方向を決定する光反射板模擬ステップと、
前記光子が前記光源に入射したときにその吸収の判定または反射方向の決定を行う光源模擬ステップと、
前記光子が前記受光面に入射したときに前記光強度分布を算出する光強度分布算出ステップと、
を有し、
前記光子発生ステップを実行して前記光子を発生させ、前記光反射板に入射する前記光子を前記光反射板模擬ステップを実行して前記光子の反射方向を決定し、前記光源に入射する前記光子を前記光源模擬ステップを実行して前記光子の吸収または反射方向を決定し、前記受光面に入射する前記光子数から前記光強度分布算出ステップを実行する追跡計算を、前記光子の発生回数が所定数に達するまで行い、
前記光反射板模擬ステップでは、前記光反射板に入射した前記光子が鏡面反射するか拡散反射するかを決定する
ことを特徴とする光反射板の光強度分布シミュレーション方法。
【請求項７】
前記光反射板模擬ステップでは、前記鏡面反射による光強度分布を２以上の所定の関数を用いて模擬し、前記光強度分布に基づいて前記反射する光子の反射方向の決定を行う
ことを特徴とする請求項６に記載の光反射板の光強度分布シミュレーション方法。
【請求項８】
前記所定の関数は、ガウス関数、ルジャンドル関数、ピアソン関数のいずれか１つ、あるいは２以上を組み合わせた関数である
ことを特徴とする請求項６に記載の光反射板の光強度分布シミュレーション装置。
【請求項９】
前記鏡面反射するか拡散反射するかの確率は、測定データとフレネル則に基づいて事前にモデル化された鏡面反射率と、拡散反射率との比率に基づいて決定されている
ことを特徴とする請求項６乃至８に記載の光反射板の光強度分布シミュレーション方法。
【請求項１０】
コンピュータにより、請求項６乃至９のいずれか１項に記載の光反射板の光強度分布シミュレーション方法を実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光を反射させる光反射板の光強度分布シミュレーション装置、シミュレーション方法、および該方法を実行させるプログラムに関し、特に熱可塑性樹脂発泡体光反射板の光強度分布シミュレーション技術に関する。
【背景技術】
【０００２】
光源から出射された光が、物体に反射されたり吸収されることにより所定の位置で得られる光強度分布を評価するシミュレーション手法が従来から知られている。特許文献１では、画像読取装置の再反射光分布を評価するシミュレーション装置等が開示されている。
【０００３】
また、ファザード看板、ＬＣＤなどの透過型ディスプレイユニットに用いられるバックライトユニットとして、従来は導光板などを用いるエッジライト方式のものが用いられてきたが、その設計開発には幾何光学をベースとする解析的シミュレーション手法、もしくはモンテカルロ法をベースとするシミュレーション手法が用いられてきた。
【０００４】
中でも、モンテカルロ法を用いたシミュレーション手法は、光強度分布を高精度に評価する方法として知られている。モンテカルロ法を用いたシミュレーション手法では、光源から出射される光を光子の集合として扱い、光子１つずつに対し、その反射・吸収を逐次計算することで光強度分布を評価している。
【０００５】
近年、透過型ディスプレイユニットの大型化にともない、直下型と呼ばれる方式のバックライトが一般化してきている。直下型方式のバックライトにおいては、白色の光拡散性の強い光反射板を光出射面の裏面側に設けるのが一般的である。このようなバックライトユニットの設計開発では、光学系の解析学的な手法による評価は困難であり、モンテカルロ法をベースとしたシミュレーションがより好適な手法となる。白色の光拡散性の強い反射板として、熱可塑性樹脂発泡体光反射板が知られている。
【特許文献１】特開２００８−６５６４７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来のシミュレーション手法では、評価対象を構成する部材の光学特性を、完全拡散面、もしくは反射光をガウス関数等のピークフィッティング関数でフィッティングした反射面等からなるものとして模擬していた。これに対し、熱可塑性樹脂発泡体光反射板のような白色の光拡散性の強い光反射板では、鏡面反射に加えて拡散反射も大きな割合で生じている。そのため、上記従来の模擬方法ではこのような光反射板の光学特性を十分に模擬することができず、熱可塑性樹脂発泡体光反射板を用いたバックライトのシミュレーション結果が実測値から大きく乖離していた。このように、従来のシミュレーション手法は、白色の光拡散性の強い光反射板の評価用として実用に供することはできない、といった問題があった。
【０００７】
そこで、本発明は上記問題を解決するためになされたものであり、拡散反射と鏡面反射を伴う光拡散性の強い光反射板を高精度に模擬可能な光強度分布シミュレーション装置、シミュレーション方法、および該方法を実行させるプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明の光反射板の光強度分布シミュレーション装置の第１の態様は、光源の一方に光反射板を配置し他方に受光面を配置した構成において、前記光源から出射された光子が前記受光面に前記光反射板で反射されて入射するかあるいは前記光源から直接入射する個数から光強度分布を算出する光反射板の光強度分布シミュレーション装置であって、前記光子が前記光反射板に入射したときにその反射方向を決定する光反射板模擬部と、前記光子が前記光源に入射したときにその吸収の判定または反射方向の決定を行う光源模擬部と、前記光子を発生させ、前記光子が前記光反射板に入射する場合には前記光反射板模擬部を用いて前記光子の反射方向を決定し、前記光源に入射する前記光子を前記光源模擬部を用いて前記光子の吸収または反射方向を決定し、前記受光面に入射する前記光子の個数から前記光強度分布を算出する追跡計算を、前記光子の発生回数が所定数に達するまで行う光子追跡計算部と、を備え、前記光反射板模擬部は、前記光反射板に入射した前記光子が鏡面反射するか拡散反射するかを所定の確率で決定していることを特徴とする。
【０００９】
本発明の光反射板の光強度分布シミュレーション装置の他の態様は、前記光反射板模擬部は、前記鏡面反射による光強度分布を２以上の所定の関数を用いて模擬し、前記光強度分布に基づいて反射する光子の反射方向の決定を行うことを特徴とする。
【００１０】
本発明の光反射板の光強度分布シミュレーション装置の他の態様は、前記所定の関数は、ガウス関数、ルジャンドル関数、ピアソン関数のいずれか１つ、あるいは２以上を組み合わせた関数であることを特徴とする。
【００１１】
本発明の光反射板の光強度分布シミュレーション装置の他の態様は、前記鏡面反射するか拡散反射するかの確率は、前記光子の入射角に対する前記鏡面反射の反射率がフレネル則に従うように決定されていることを特徴とする。
【００１２】
本発明の光反射板の光強度分布シミュレーション装置の他の態様は、前記光反射板は、熱可塑性樹脂発泡体光反射板であることを特徴とする。
【００１３】
本発明の光反射板の光強度分布シミュレーション方法の第１の態様は、光源の一方に光反射板を配置し他方に受光面を配置した構成において、前記光源から出射された光子が前記受光面に前記光反射板で反射され入射するかあるいは前記光源から直接入射する個数から光強度分布を算出する光反射板の光強度分布シミュレーション方法であって、前記光源から光子を出射させる光子発生ステップと、前記光子が前記光反射板に入射したときにその反射方向を決定する光反射板模擬ステップと、前記光子が前記光源に入射したときにその吸収の判定または反射方向の決定を行う光源模擬ステップと、前記光子が前記受光面に入射したときに前記光強度分布を算出する光強度分布算出ステップと、を有し、前記光子発生ステップを実行して前記光子を発生させ、前記光反射板に入射する前記光子を前記光反射板模擬ステップを実行して前記光子の反射方向を決定し、前記光源に入射する前記光子を前記光源模擬ステップを実行して前記光子の吸収または反射方向を決定し、前記受光面に入射する前記光子数から前記光強度分布算出ステップを実行する追跡計算を、前記光子の発生回数が所定数に達するまで行い、前記光反射板模擬ステップでは、前記光反射板に入射した前記光子が鏡面反射するか拡散反射するかを決定することを特徴とする。
【００１４】
本発明の光反射板の光強度分布シミュレーション方法の他の態様は、前記光反射板模擬ステップでは、前記鏡面反射による光強度分布を２以上の所定の関数を用いて模擬し、前記光強度分布に基づいて前記反射する光子の反射方向の決定を行うことを特徴とする。
【００１５】
本発明の光反射板の光強度分布シミュレーション方法の他の態様は、前記所定の関数は、ガウス関数、ルジャンドル関数、ピアソン関数のいずれか１つ、あるいは２以上を組み合わせた関数であることを特徴とする。
【００１６】
本発明の光反射板の光強度分布シミュレーション方法の他の態様は、前記鏡面反射するか拡散反射するかの確率は、測定データとフレネル則に基づいて事前にモデル化された鏡面反射率と、拡散反射率との比率に基づいて決定されていることを特徴とする。
【００１７】
本発明の他の態様は、コンピュータにより、上記の光反射板の光強度分布シミュレーション方法を実行させるプログラムである。
【発明の効果】
【００１８】
本発明によれば、拡散反射と鏡面反射を伴う光拡散性の強い光反射板を高精度に模擬可能な光強度分布シミュレーション装置、シミュレーション方法、および該方法を実行させるプログラムを提供することができる。
【発明を実施するための最良の形態】
【００１９】
本発明の好ましい実施の形態における光反射板の光強度分布シミュレーション装置、シミュレーション方法、および該方法を実行させるプログラムについて、図面を参照して詳細に説明する。同一機能を有する各構成部については、図示及び説明簡略化のため、同一符号を付して示す。
【００２０】
以下では、拡散反射と鏡面反射を伴う光拡散性の強い光反射板の一例として、熱可塑性樹脂発泡体光反射板を用いたときの実施形態について説明する。熱可塑性樹脂発泡体光反射板は、微細発泡を有して白色の光拡散性の強い光反射板であり、図２に示すような構造を有している。図２（ａ）は、熱可塑性樹脂発泡体光反射板（以下では単に光反射板という）１０の平面図を、図２（ｂ）は光反射板１０の側面図を、図２（ｃ）は光反射板１０の断面の一部を拡大した部分断面図を、それぞれ示している。
【００２１】
図２（ｃ）では、光反射板１０の上部に光源２０を１つだけ配置した例を示している。光反射板１０は、最表面（以下では第１表面とする）１１から所定厚さの未発泡層１３と、それより内側の発泡層１４とで構成されており、未発泡層１３と発泡層１４との間に界面（以下では第２表面とする）１２が形成されている。第２表面１２は、第１表面１１と略平行に形成されているものとする。
【００２２】
光源２０から光反射板１０に入射するときの入射角は、第１表面１１の法線と入射光とのなす角度αである。また、光反射板１０の法線方向を０度とし水平方向を９０度とする緯度θ（図２（ｂ））と、光反射板１０の水平面上の方位角である経度φ（図２（ａ））とを定義して以下で用いるものとする。
【００２３】
図２のように構成された光反射板１０では、入射光の一部が第１表面１１または第２表面１２で鏡面反射されるとともに、残りの入射光は発泡層１４まで達して拡散反射される。鏡面反射される入射光と拡散反射される入射光との割合は、光反射板１０への入射角αによって異なっており、一例として入射角が６０°のときには、鏡面反射が約１７％となり拡散反射が約８３％となる。
【００２４】
本発明の第１の実施の形態に係る光強度分布シミュレーション装置およびシミュレーション方法を、図１に示すブロック図を用いて説明する。図１は、本実施形態の光強度分布シミュレーション装置１００の構成を示すブロック図である。本実施形態の光強度分布シミュレーション装置およびシミュレーション方法は、上記の熱可塑性樹脂発泡体光反射板１０のような拡散反射と鏡面反射とを伴う光拡散性の強い光反射板を対象に、これを用いたときの光強度分布を高精度に評価することが可能となっている。
【００２５】
以下では、図３に示す光源２０および光反射板１０を評価対象の一例として、本実施形態の光強度分布シミュレーション装置１００およびシミュレーション方法を説明する。図３では、光反射板１０の一辺に略等しい長さの線状の光源２０が、光反射板１０の上部に複数本平行に配列されている。光反射板１０は、図２に示した熱可塑性樹脂発泡体光反射板と同じ構造のものとしている。また、光強度分布を算出する位置を、光反射板１０の上部の所定高さの位置に設定しており、以下ではこれを受光面３０としている。なお、ここでは簡単化のために評価対象として光反射板１０および光源２０のみを表示しているが、その他の部材として、たとえばアクリル板、光学フィルム、プリズムシート、偏光フィルム等を含めることができる。
【００２６】
本実施形態の光強度分布シミュレーション装置１００は、入出力装置１１０、演算処理装置１２０、および記憶装置１３０を備えている。入出力装置１１０は、シミュレーション計算を実行させるための各種条件等を入力するとともに、シミュレーション計算の結果を出力する手段を有している。入力する各種条件として、シミュレーション対象の光学部材である光反射板１０、光源２０等の配置データや光学特性データ等がある。また、光強度分布を評価する受光面３０の位置もここで設定することができる。入力された各種条件は、演算処理装置１２０を介して記憶装置１３０に記憶させ、必要に応じて記憶装置１３０から読み込んで用いている。
【００２７】
演算処理装置１２０は、光子追跡計算部１２１、光反射板模擬部１２２、および光源模擬部１２３を有している。さらに別の光学部材を有するときは、それを模擬する模擬部を追加することができる。たとえば、別の光学部材としてアクリル板、光学フィルム、プリズムシート、および偏光フィルムを有するときは、それぞれを模擬するアクリル板模擬部、光学フィルム模擬部、プリズムシート模擬部、および偏光フィルム模擬部を追加する。光子追跡計算部１２１は、たとえばモンテカルロ法を用いて任意の光源２０から光子を出射させ、これが光反射板１０や別の光源２０に入射して反射したり吸収される軌跡を追跡計算する。
【００２８】
光反射板模擬部１２２は、入出力装置１１０から入力された光反射板１０の配置データや光学特性データ等に基づいて、光反射板１０に入射した光子の軌跡を決定するのに用いられる。光源模擬部１２３は、入出力装置１１０から入力された光源２０の配置データや光学特性データ等に基づいて、光源２０からの光子の出射や入射してきた光子の軌跡等を決定するのに用いられる。
【００２９】
光反射板模擬部１２２および光源模擬部１２３を用いて光子追跡計算部１２１でシミュレーション計算された光強度分布等の結果は、記憶装置１３０に保存されるとともに入出力装置１１０に適宜表示させることができる。
【００３０】
本実施形態の光反射板模擬部１２２は、拡散反射と鏡面反射を伴う光拡散性の強い光反射板１０を高精度に模擬できるように構成されている。従来は、光反射板による反射光の光強度を所定の関数を用いて近似するように構成されていた。すなわち、光強度の測定データを所定の関数でフィッティングし（以下では従来フィッティングと呼ぶ）、これをシミュレーション計算に用いていた。ここで、光強度の測定データは、例えば１つの光源から光反射板に光を照射して反射される光強度を測定したデータである。なお、フィッティングに用いられる関数は、光強度のピークをフィッティングするのに好適な関数（ピークフィッティング関数）であり、たとえばガウス関数、ルジャンドル関数、ピアソン関数等を用いることができる。
【００３１】
しかしながら、このような従来フィッティングでは、光反射板による反射が経度に依存しない拡散反射成分を有しているときには、これを適切に模擬することができなかった。そのため、測定データと従来フィッティングで算出した光強度とに大きな乖離が見られた。
【００３２】
これに対し、本実施形態の光反射板模擬部１２２では、鏡面反射成分と拡散反射成分とを分離して模擬する構成としており、それぞれの反射成分の確率を算出するための鏡面／拡散反射成分算出手段１２２ａを備えている。また、鏡面／拡散反射成分算出手段１２２ａは、鏡面反射成分を複数の所定の関数を用いてフィッティングしたものを用いるように構成されており、これにより鏡面反射成分を高精度に評価することが可能となっている。所定の関数として、たとえばガウス関数、ルジャンドル関数、ピアソン関数等のピークフィッティング関数を用いることができる。あるいは、これらのピークフィッティング関数を２以上組み合わせて用いるようにしてもよい。
【００３３】
予め測定されたデータを上記のように高精度にフィッティングして用いる本実施形態の鏡面／拡散反射成分算出手段１２２ａによれば、鏡面反射成分および拡散反射成分のそれぞれの割合を高精度に算出してシミュレーション計算に用いるようにすることができる。鏡面反射成分および拡散反射成分のそれぞれの割合は光の入射角によって異なることから、０〜９０°の全ての入射角で光強度を測定し、鏡面／拡散反射成分算出手段１２２ａではこの測定データをフィッティングしたピークフィッティング関数を用いるのが好ましい。しかしながら、入射角全体で測定データを取得することは難しく、特に高い入射角の測定データが得られないことが多い。
【００３４】
そこで、本実施形態では、鏡面反射がフレネル則に従うことを利用して、測定された所定の入射角までの光強度データをもとに、測定できない入射角に対する反射成分をフレネル則に基づいて補間して用いるようにしている。すなわち、測定可能な入射角に対しては、鏡面反射成分および拡散反射成分のそれぞれの割合を測定データから求め、測定できない入射角に対しては、鏡面反射成分の割合をフレネル則を用いて算出している。一例として、０〜６０°の入射角の反射成分割合が測定データから得られる場合には、６０〜９０°間の鏡面反射成分割合をフレネル則に従って算出している。このようにして得られた０〜９０°の光強度のデータを、複数のピークフィッティング関数を用いてフィッティングしている。
【００３５】
入射角αで反射板に入射した入射光は、図４に示すように、反射面において反射角αで鏡面反射する成分と屈折角βで透過する成分とに分けられる。このとき、フレネル則は次式で表される。
2010020528.tif 000003
2010020528.tif 000004
2010020528.tif 000005
2010020528.tif 000006
ここで、ｒｐ、ｒｓがそれぞれｐ偏光、ｓ偏光の反射率、ｔｐ、ｔｓがそれぞれｐ偏光、ｓ偏光の透過率を表し、ｎ１、ｎ２が図４に示す第１層２１、第２層２２のそれぞれの屈折率を表す。但し、ｐ偏光およびｓ偏光は、それぞれ入射面内方向で入射方向に垂直方向の偏光、および入射面に垂直方向の偏光を表している。
【００３６】
熱可塑性樹脂発泡体光反射板１０では、上記の透過率で第１表面１１を透過した光は、第２表面１２で再びフレネル則に基づいて鏡面反射するか、あるいは第２表面１２を透過して発泡層１４に達する。発泡層１４に達した光は、ここで拡散反射される。熱可塑性樹脂発泡体光反射板１０の各層における入射角に対する反射率の変化の一例を図５に示す。同図では、第１表面１１および第２表面１２におけるそれぞれの鏡面反射率をＲ１ｍ、Ｒ２ｍで示している。また、第１表面１１での反射分を除いた第２表面での実質的な鏡面反射率をＲ２で示し、第１表面１１での鏡面反射率Ｒ１ｍと第２表面での実質的な鏡面反射率Ｒ２とを加算した全鏡面反射率をＲｍで示している。
【００３７】
図５では、さらに発泡層１４における拡散反射の反射率をＲｄで示している。同図より、入射角が５０°を超えたあたりから９０°に近づくにつれて、鏡面反射の反射率Ｒｍが大きくなることがわかる。それに伴って、拡散反射の反射率Ｒｄが急激に低下している。高緯度の光強度は、大きな入射角の光が鏡面反射した影響を強く受けることから、高緯度の光強度の測定データが得られない場合には、フレネル則に基づく反射率で高緯度の光強度を補間するのが好ましい。
【００３８】
一方、発泡層１４における拡散反射成分については、理論ランバート拡散〔完全拡散〕を模擬することにより、拡散反射による光強度Ｉｄを次式で比較的良好に近似することができる。
Ｉｄ＝ｃｏｓ（θ）
ここで、θは緯度を表す。
【００３９】
上記説明の本実施形態の光強度分布シミュレーション装置１００で行われるシミュレーション計算の方法を、図６に示す流れ図を用いてさらに詳細に説明する。本実施形態のシミュレーション方法は、図６に示す流れ図に基づいて作成されたプログラムを図１に示した演算処理装置１２０に装荷して実行される。このプログラムは、ＦＤ等の記録媒体に保存することができる。
【００４０】
本実施形態のシミュレーション方法では、モンテカルロ法を用いて受光面３０における光強度分布を算出している。モンテカルロ法では、一つ一つの光子の軌跡をランダムに決定しているが、多数の光子全体で見たときは、所定の光学特性に従うように決定している。以下の説明では、このような光子の軌跡の決定方法を「所定の条件に従ってランダムに決定する」と表現する。
【００４１】
まずステップＳ１において、シミュレーション計算に必要な初期値設定を行う。すなわち、シミュレーション評価の対象となる光反射板１０や光源２０等の寸法、配置、光学特性等を設定する。光反射板１０の光学特性の一つとして、鏡面反射成分のフィッティング式がある。また、別の光学部材としてアクリル板、光学フィルム、プリズムシート、偏光フィルム等を有するときは、それぞれの寸法、配置、光学特性等を設定する。初期値設定では、光強度分布を評価する受光面３０の位置設定も行う。さらに、モンテカルロ法の初期値設定として、光子の数ｎｐを０に初期化する。
【００４２】
ステップＳ２では、光子の数ｎｐに１を加算して光子数を積算する。続くステップＳ３では、１つの光子（以下ではＰ１とする）について、これを出射させる光源２０の位置および出射方向を所定の条件に従ってランダムに決定する。
【００４３】
ステップＳ４では、ステップＳ３で発生させた光子Ｐ１の入射先を、光子Ｐ１の出射位置（光源２０の位置）及び出射方向から算出する。その結果、光子Ｐ１が受光面３０に入射するかそれ以外の光反射板１０等に入射するかをステップＳ５で判定し、受光面３０に入射せずに反射板１０またはその他部材に入射する場合には、次のステップＳ６で光子Ｐ１の追跡計算を行う。
【００４４】
ステップＳ６の光子追跡計算では、光子Ｐ１が光反射板１０に入射したと判定されると、光反射板模擬部１２２を用いて光子Ｐ１が反射される方向を上記説明のピークフィッティング関数で模擬した鏡面反射と拡散反射の確率に従ってランダムに決定する。また、光子Ｐ１が光源２０に入射したと判定されると、光源模擬部１２３を用いて光子Ｐ１が吸収されるか反射されるかを決定し、反射される場合にはその方向を鏡面反射の条件に従ってランダムに決定する。さらに別の光学部材を有するときには、光子Ｐ１が入射した部材に対応する模擬部を用いて吸収または反射方向を決定する。
【００４５】
ステップＳ７において、ステップＳ６の光子追跡計算の結果光子Ｐ１がいずれかの方向に反射されたと判定される場合には、再びステップＳ４に戻ってステップＳ６までの処理を繰り返し行い、光子Ｐ１の軌跡を評価していく。またステップＳ７において、光子Ｐ１がいずれかで吸収されると判定されると、次のステップＳ９に進む。
一方、ステップＳ５で光子Ｐ１が受光面３０に入射したと判定されると、ステップＳ８において受光面３０での光強度分布を算出する。
【００４６】
ステップＳ８までの処理により、光子Ｐ１がいずれかの部材に吸収されると、ステップＳ９において、ステップＳ３でそれまでに発生させた光子数ｎｐが所定の数Ｎｍａｘに達したか否かを判定し、Ｎｍａｘに達していない場合にはステップＳ２に戻って光子数ｎｐを１だけ加算し、次のステップＳ３で別の光子を発生させてシミュレーション計算を継続する。一方、光子数ｎｐが所定の数Ｎｍａｘに達すると、シミュレーション計算を終了する。光子数Ｎｍａｘの光子を発生させたときにステップＳ８で算出される受光面３０での光強度分布が、シミュレーション結果として記憶装置１３０に記憶されるとともに、入出力装置１１０に出力される。
【００４７】
シミュレーション結果の一例を図７に示す。同図は、横軸を光反射板１０に沿った位置、縦軸を光反射板１０による反射光の輝度としたときのシミュレーション結果の一例を示しており、隣接する２つの光源２０近傍の輝度を示している。同図では、本実施形態のシミュレーション装置１００によるシミュレーション結果を符号４１で示しており、これに加えて実測データを符号４２、従来のシミュレーション結果を符号４３で示している。従来のシミュレーション結果は、光反射板１０による反射光の光強度を、拡散反射成分と鏡面反射成分とに分離せずに１つのガウス関数でフィッティングしたものを用いたときのシミュレーション結果である。
【００４８】
図７に示すように、従来フィッティングを用いた従来のシミュレーション結果４３では、輝度が実測データ４２から大きく乖離している位置が多数見られるのに対し、本実施形態のシミュレーション装置１００によるシミュレーション結果４１は、実測データ４２と全体的によい一致を示している。特に、光源と光源との間における輝度については、本実施形態のシミュレーション結果４１が従来のシミュレーション結果４３に比べて大幅に改善されており、拡散反射による光強度への寄与が適切にシミュレーションされていることが示されている。
【００４９】
本発明の第２の実施の形態に係る光強度分布シミュレーション装置およびシミュレーション方法について、以下に説明する。本実施形態の光強度分布シミュレーション装置では、光反射板模擬部１２２において光反射板１０による反射光を拡散反射成分と鏡面反射成分とに分離し、鏡面反射成分のピークを１つのピークフィッティング関数を用いてフィッティングしたものを用いている。本実施形態でも、ピークフィッティング関数として、たとえばガウス関数、ルジャンドル関数、ピアソン関数等を用いることができる。
【００５０】
上記の第１の実施形態では、光反射板模擬部１２２が鏡面反射成分を複数のピークフィッティング関数を用いてフィッティングしたものを用いるように構成されていたが、本実施形態では１つのピークフィッティング関数だけで鏡面反射成分をフィッティングするように構成されている。この場合でも、光反射板１０による反射光を拡散反射成分と鏡面反射成分とに分離して評価するようにしていることから、従来フィッティングを用いた場合に比べて光強度分布を高精度に評価することが可能となっている。とくに、光源と光源との間の位置における輝度については、本実施形態によるシミュレーション結果が、従来のシミュレーション結果に比べて大幅に改善されている。
【００５１】
上記説明のように本発明によれば、拡散反射と鏡面反射を伴う光拡散性の強い光反射板を高精度に模擬可能な光強度分布シミュレーション装置、シミュレーション方法、該方法を実行させるプログラムおよび該プログラムを格納した記録媒体を提供することができる。上記では、シミュレーション対象として熱可塑性樹脂発泡体光反射板を用いて説明したが、これに限らず、鏡面反射に加えて拡散反射を伴う光反射板を対象にシミュレーションを行うのに好適なシミュレーション装置およびシミュレーション方法を提供することができる。
【００５２】
なお、本実施の形態における記述は、本発明に係る反射板の光強度分布シミュレーション装置、シミュレーション方法、該方法を実行させるプログラムおよび該プログラムを格納した記録媒体の一例を示すものであり、これに限定されるものではない。本実施の形態における光強度分布シミュレーション装置等の細部構成及び詳細な動作等に関しては、本発明の趣旨を逸脱しない範囲で適宜変更可能である。
【図面の簡単な説明】
【００５３】
【図１】本発明の第１の実施形態に係る光反射板の光強度分布シミュレーション装置の構成を示すブロック図である。
【図２】熱可塑性樹脂発泡体光反射板の構造の一例を示す図である。
【図３】シミュレーション対象の光源および光反射板の配置構成を示す図である。
【図４】フレネル則の説明に用いる反射面での反射成分、透過成分を示す図である。
【図５】入射角に対する鏡面反射の反射率及び拡散反射の反射率の変化の一例を示す図である。
【図６】第１の実施形態の光強度分布シミュレーション方法の処理を示す流れ図である。
【図７】シミュレーション結果の一例として、輝度の分布を示す図である。
【符号の説明】
【００５４】
１０      熱可塑性樹脂発泡体光反射板
１１      第１表面
１２      第２表面
１３      未発泡層
１４      発泡層
２０      光源
２１      第１層
２２      第２層
３０      受光面
１００    光強度分布シミュレーション装置
１１０    入出力装置
１２０    演算処理装置
１２１    光子追跡計算部
１２２    光反射板模擬部
１２２ａ  鏡面／拡散反射成分算出手段
１２３    光源模擬部
１３０    記憶装置
【図１】
2010020528.tif 000007
【図２】
2010020528.tif 000008
【図３】
2010020528.tif 000009
【図４】
2010020528.tif 000010
【図５】
2010020528.tif 000011
【図６】
2010020528.tif 000012
【図７】
2010020528.tif 000013