(57)【特許請求の範囲】
【請求項１】
透明電極抵抗と直列抵抗とシャント抵抗とダイオード成分と電流源とを用いて太陽電池セルを表す分布定数等価回路モデルから得た、太陽電池セルの特性を求めるための非線形微分方程式を用いて、太陽電池セルを設計する方法であって、
理想太陽電池における電圧Ｖdと電流ｉiscとの関係を表す曲線を求めるステップと、
前記電流源の電流をｉLとし、前記理想太陽電池の電圧Ｖdと電流ｉiscとの関係を表す曲線において、電圧が０のときの曲線の傾きをｇsc、電圧が開放電圧Ｖocに等しいときの曲線の傾きをｇocとして、該曲線を近似式
【数１】
0005637398.tif 000012
で近似するステップと、
前記非線形微分方程式と、前記近似式とを用いて、前記太陽電池セルの寸法と特性との関係を求めるステップと、
求められた前記太陽電池セルの寸法と特性との関係に基づいて、太陽電池セルの寸法を設計するステップとを含むことを特徴とした太陽電池セルの設計方法。
【請求項２】
前記太陽電池セルの寸法を設計するステップにおいて、前記求められた太陽電池セルの寸法と特性との関係を用いて、前記太陽電池セルに要求される発電電力が得られるように太陽電池セルの寸法を設計する、請求項１に記載の太陽電池セルの設計方法。
【請求項３】
前記太陽電池セルの特性は、電流密度、形状因子又は発電電力である、請求項１又は請求項２に記載の太陽電池セルの設計方法。
【請求項４】
前記太陽電池セルの出力電圧をＶa、前記透明電極抵抗の抵抗値をρs、前記直列抵抗の抵抗値をｒc、前記シャント抵抗の抵抗値をｒshとおき、
【数２】
0005637398.tif 000013
とすると、
前記太陽電池セルの寸法がＬp以下となるように選択される、請求項１から請求項３の何れか一項に記載の太陽電池セルの設計方法。
【請求項５】
前記太陽電池セルの出力電圧をＶa、前記透明電極抵抗の抵抗値をρs、前記直列抵抗の抵抗値をｒc、前記シャント抵抗の抵抗値をｒshとおき、
【数３】
0005637398.tif 000014
とすると、
前記太陽電池セルの分布定数等価回路モデルにおける前記シャント抵抗ｒshに流れる電流が漏れ電流の大部分を占めるとき、前記太陽電池セルの寸法がＬs以下となるように選択される、請求項１から請求項４の何れか一項に記載の太陽電池セルの設計方法。
【請求項６】
透明電極抵抗と直列抵抗とシャント抵抗とダイオード成分と電流源とを用いて太陽電池セルを表す分布定数等価回路モデルから得た、太陽電池セルの特性を求めるための非線形微分方程式を用いて、太陽電池セルの特性を求めるためのシミュレーションプログラムであって、
理想太陽電池の電圧Ｖdと電流ｉiscとの関係を表す曲線のデータを入力させるステップと、
前記電流源の電流をｉL、入力された前記理想太陽電池の電圧Ｖdと電流ｉiscとの関係を表す曲線において、電圧が０のときの曲線の傾きをｇsc、電圧が開放電圧Ｖocに等しいときの曲線の傾きをｇocとして、該曲線を近似式
【数４】
0005637398.tif 000015
で近似するステップと、
前記非線形微分方程式と前記近似式とを用いて、前記太陽電池セルの寸法と特性との関係を求めるステップと、
求められた前記太陽電池セルの寸法と特性との関係を表示装置に表示するステップと、
を実行することを特徴とする太陽電池セルの特性を求めるためのシミュレーションプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、透明電極抵抗と直列抵抗とシャント抵抗とダイオード成分と電流源とを用いて太陽電池セルを表す分布定数等価回路モデルから得た、太陽電池セルの特性を求めるための非線形微分方程式を用いて、太陽電池セルを設計する方法に関する。
【背景技術】
【０００２】
太陽電池セルは、図９に示されているように、大型化に伴って、特性の中で形状因子（以下、単にＦＦとも記載する。）及び短絡電流（以下、単にＪscとも記載する。）が減衰することが知られている（図９では、短絡電流Ｊscを光起電流ＪLで割ることで無次元化している）。そこで、大型の太陽電池パネルを作製する際には、図１０に示されているように、特性の減衰を避けるために補助電極の導入やスクライブによるモジュール化によって大型太陽電池パネルを分割する技術が用いられてきている。なお、図１０では、並列する太陽電池セルの幅方向（図中の左右方向）の分割サイズをＬ2、並列する太陽電池セルの幅方向の有効サイズ（太陽発電が可能な部分の寸法）をＬとしている。
【０００３】
しかしながら、分割によって太陽電池セルのサイズを小さくすれば一つの太陽電池セル当たりの特性は高くなるが、有効サイズＬも小さくなって、モジュール全体としての効率は低下するので、太陽電池セルのサイズとモジュール全体の効率とはトレードオフの関係にある。したがって、開口率すなわちＬ／Ｌ2を高く保ちつつ高いセル特性が得られるように太陽電池セルの寸法を設計することが求められる。このためには、太陽電池セルの寸法とセル特性との関係を解析し、その関係から最適な構造を選択することが必要になる。従来、このような太陽電池セルの寸法とセル特性との関係の解析の際には、非特許文献１〜５に開示されているように、分布定数等価回路モデルを用いることが一般的である。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Alexis de Vos, Solar Cells 12, 311-327(1984).
【非特許文献２】D.V. Rumyantsev and J.A. Rodrguez, Solar Cells 28, 241-252(1990).
【非特許文献３】N. Convers Wyeth, Sol. Stat. Elect., 20(1976)629.
【非特許文献４】G.C. Jain and F.M. Stuber, Adv. Energ. Conv 7(1967)167.
【非特許文献５】M. Burgelman et.al, Sol. Energy Mat. & Sol. Cells 51,(1998)129.CIS
【発明の概要】
【発明が解決しようとする課題】
【０００５】
太陽電池セルの特性を求めるためには、上述した分布定数等価回路モデルから得た微分方程式を解く必要があるが、この微分方程式は非線形であるので、これを解くためには、コンピュータを用いて繰り返し計算を行うことによって数値的に非線形微分方程式を解く数値的手法か、非線形微分方程式を近似する近似手法を用いる必要がある。
【０００６】
しかしながら、数値的手法では、プロトタイプとなる太陽電池セルの寸法構造や特性パラメータを仮定して、それらについて個々に数値計算を行った上で、太陽電池セルの特性曲線を求め、求められた特性曲線から最適な寸法構造を選択するという手順で太陽電池セルの寸法の設計を行う必要がある。したがって、多くの時間と計算資源が必要となるという問題があった。
【０００７】
また、近似手法も、等価回路モデルにおけるダイオード成分のような重要なパラメータを無視してしまうなど、実用性に欠けるという問題があった。
【０００８】
よって、本発明の目的は、従来技術に存する問題を解消して、太陽電池セルの分布定数等価回路モデルのパラメータを全て取り入れながら、簡便に太陽電池セルの寸法と特性との関係を表す特性曲線を求めることを可能にすることにある。
【課題を解決するための手段】
【０００９】
上記目的に鑑み、本発明は、透明電極抵抗と直列抵抗とシャント抵抗とダイオード成分と電流源とを用いて太陽電池セルを表す分布定数等価回路モデルから得た、太陽電池セルの特性を求めるための非線形微分方程式を用いて、太陽電池セルを設計する方法であって、理想太陽電池における電圧Ｖdと電流ｉiscとの関係を表す曲線を求めるステップと、前記電流源の電流をｉLとし、前記理想太陽電池の電圧Ｖdと電流ｉiscとの関係を表す曲線において、電圧が０のときの曲線の傾きをｇsc、電圧が開放電圧Ｖocに等しいときの曲線の傾きをｇocとして、該曲線を近似式
【数１】
0005637398.tif 000002
で近似するステップと、前記非線形微分方程式と前記近似式とを用いて、前記太陽電池セルの寸法と特性との関係を求めるステップと、求められた前記太陽電池セルの寸法と特性との関係に基づいて、太陽電池セルの寸法を設計するステップとを含む太陽電池セルの設計方法を提供する。
【００１０】
上記太陽電池セルの設計方法では、理想太陽電池の電圧と電流との関係を表す曲線を折れ線によって近似することにより、太陽電池セルの分布定数等価回路モデルから得られる、太陽電池セルの特性を求めるための微分方程式を線形化して、数式解として解けるようにし、太陽電池セルの等価回路モデルの特性パラメータや形状設計データから、直接的に太陽電池セルの寸法と特性との関係を表す曲線を求めることができるようにする。
【００１１】
上記太陽電池セルの設計方法では、例えば、前記太陽電池セルの寸法を設計するステップにおいて、前記求められた太陽電池セルの寸法と特性との関係を用いて、前記太陽電池セルに要求される発電電力が得られるように太陽電池セルの寸法を設計することができる。
また、前記太陽電池セルの特性は、例えば電流密度、形状因子又は発電電力とすることができる。
【００１２】
上記太陽電池セルの設計方法では、前記太陽電池セルの出力電圧をＶa、前記透明電極抵抗の抵抗値をρs、前記直列抵抗の抵抗値をｒc、前記シャント抵抗の抵抗値をｒshとおき、
【数２】
0005637398.tif 000003
とすると、
前記太陽電池セルの寸法がＬp以下となるように選択されることが好ましい。
【００１３】
さらに、前記太陽電池セルの出力電圧をＶa、前記透明電極抵抗の抵抗値をρs、前記直列抵抗の抵抗値をｒc、前記シャント抵抗の抵抗値をｒshとおき、
【数３】
0005637398.tif 000004
とすると、
前記太陽電池セルの分布定数回路モデルにおける前記シャント抵抗ｒshに流れる電流が漏れ電流の大部分を占めるとき、前記太陽電池セルの寸法がＬs以下となるように選択されることが好ましい。
【００１４】
また、本発明は、透明電極抵抗と直列抵抗とシャント抵抗とダイオード成分と電流源とを用いて太陽電池セルを表す分布定数等価回路モデルから得た、太陽電池セルの特性を求めるための非線形微分方程式を用いて、太陽電池セルの特性を求めるシミュレーションプログラムであって、理想太陽電池の電圧Ｖdと電流ｉiscとの関係を表す曲線のデータとの入力を求めるステップと、前記電流源の電流をｉL、入力された前記理想太陽電池の電圧Ｖdと電流ｉiscとの関係を表す曲線において、電圧が０のときの曲線の傾きをｇsc、電圧が開放電圧Ｖocに等しいときの曲線の傾きをｇocとして、該曲線を近似式
【数４】
0005637398.tif 000005
で近似するステップと、前記非線形微分方程式と前記近似式とを用いて、前記太陽電池セルの寸法と特性との関係を求めるステップと、求められた前記太陽電池セルの寸法と特性との関係を表示装置に表示するステップとを実行する、太陽電池セルの特性を求めるためのシミュレーションプログラムを提供する。
【発明の効果】
【００１５】
本発明によれば、理想太陽電池セルにおける電圧と電流との関係を表す曲線を折れ線によって近似することにより、分布定数等価回路モデルから得られる非線形微分方程式について、特性パラメータ（電流源、ダイオード成分、並列抵抗、直列抵抗、透明電極抵抗）を全て含みながらも線形化して数式解として解くことを可能にする。したがって、各特性パラメータが太陽電池セルの寸法設計への影響を考慮することを可能にしながら、複雑な計算によることなく、太陽電池セルの分布定数等価回路モデルの特性パラメータや形状設計データから、直接的に太陽電池セルの寸法と特性との関係を表す曲線を求めることができるようになる。
【図面の簡単な説明】
【００１６】
【図１】（ａ）及び（ｂ）は太陽電池セルの分布定数等価回路モデルを示す説明図、（ｃ）は本発明に従って理想太陽電池セルの電流−電圧特性曲線を折れ線で近似した近似曲線を示すグラフである。
【図２】本発明に従って求めた太陽電池セルの電流密度と出力電圧との関係を示すグラフである。
【図３】本発明に従って求めた太陽電池セルの長さと短絡電流及び形状因子との関係を示すグラフである。
【図４】太陽電池セルの長さと短絡電流及び形状因子との関係を表す曲線と共に、太陽電池セルの長さと太陽電池セルの出力電力との関係を表す曲線を示したグラフである。
【図５】（ａ）は補助電極で囲まれた太陽電池セルの全体構成を示す全体構成図であり、（ｂ）は補助電極で囲まれた太陽電池セルの短絡電流と長さＬ及び幅Ｗとの関係を表す、数値的解法で求めた曲線及び本発明による折れ線近似曲線を示すグラフである。
【図６】太陽電池セルの透明電極抵抗が太陽電池セルの特性に与える影響を示すグラフである。
【図７】太陽電池セルの開放電圧が太陽電池セルの特性に与える影響を示すグラフである。
【図８】太陽電池セルのシャント抵抗が太陽電池セルの特性に与える影響を示すグラフである。
【図９】太陽電池セルの電流値及び形状因子と長さとの関係を示すグラフである。
【図１０】太陽電池セルの構成を示す全体構成図である。
【発明を実施するための形態】
【００１７】
以下、図面を参照しながら、本発明による太陽電池セルを設計する方法を説明する。なお、本願における「太陽電池セルを設計する方法」とは、太陽電池モジュールなどで使用される太陽電池セルの寸法（例えば太陽電池モジュールにおいて並列に配置される太陽電池セルの長さ）や形状によってセル特性がどのように変化するかを解析し、解析結果に基づいて、要求される仕様が得られるように太陽電池セルの最適な寸法を設計する方法を意味する。また、以下の説明では、太陽電池セルの長さとは、電流の流れる方向（図１０における左右方向）を意味する。
【００１８】
太陽電池セルは、図１（ａ）に示すように、電流源、ダイオード成分、シャント抵抗（すなわち並列抵抗）、直列抵抗を組み合わせて構成した分布定数等価回路モデルによって表すことができる。詳細には、透明電極と下部電極との間に、電流源、ダイオード成分、シャント抵抗とを並列に配置し且つこれと直列に直列抵抗を配置したものを透明電極抵抗を挟んでｘ方向に並列に複数接続した形態となっている。ここで、シャント抵抗の抵抗値をｒsh、直列抵抗の抵抗値ｒc、透明電極抵抗の抵抗値をρsとし、また、太陽電池セルの長さをＬ、出力電圧をＶaとし、ｘ軸が図中左右方向に延び、左端がｘ＝０、右端がｘ＝Ｌとなるように定義するものとする。
【００１９】
このような分布定数等価回路モデルにおいて、図１（ｂ）に示されるように、位置ｘにおける透明電極上の電位及び電流をそれぞれＶ（ｘ）及びｉ（ｘ）とすると共に、位置ｘにおいて、直列抵抗を流れる電流、シャント抵抗を流れる電流、電流源を流れる電流（光起電流）、ダイオード成分を流れる電流を、それぞれ、ｉc（ｘ）、ｉsh（ｘ）、ｉL、ｉdとすると、局所的な電圧及び電流から、電気素量ｑ、ボルツマン係数ｋB、ダイオード因子ｎ、測定温度Ｔを用いて以下の関係式が得られる。
【数５】
0005637398.tif 000006
【００２０】
さらに、上関係式から以下の微分方程式が導かれる。
【数６】
0005637398.tif 000007
【００２１】
しかしながら、微分方程式（１）は非線形であり、これを数式解として解くことはできない。このため、従来は、コンピュータを用いた数値的解法により上記非線形微分方程式（１）を解いていた。これに対して、本発明による太陽電池セルを設計する方法では、図１（ｃ）に示されているように、ショックレーモデルによる理想太陽電池の電流−電圧特性（図１（ｂ）で破線によって囲まれた部分においてシャント抵抗ｒshを無限大としｉc（ｘ）に相当する電流をｉiscとしたときのｉiscとＶdとの関係）を表す曲線を以下の近似式（２）により近似し、ｉisc＝ｉd＋ｉLと表すことができることと近似式（２）とを利用して、非線形微分方程式（１）を線形化する。
【数７】
0005637398.tif 000008
ここで、ｇscは、理想太陽電池の電流−電圧特性を表す曲線において、電圧が０のときの曲線の傾き、ｇocは、理想太陽電池の電流−電圧特性を表す曲線において、電圧が開放電圧Ｖocに等しいときの曲線の傾きである。
【００２２】
このように近似式（２）を用いて微分方程式（１）を線形化することにより、微分方程式（１）の数式解を得ることができ、太陽電池セルの出力特性を解析的に求めることが可能となる。
【００２３】
図２は、近似式（２）を用いて微分方程式（１）を解くことにより求めた太陽電池セルの電流密度（太陽電池セルの出力電流／太陽電池セルの発電部分の表面積）と出力電圧との関係を表すグラフを示す。近似式（２）を用いると、微分方程式（１）が線形化され、数式解を求めることが可能となるので、太陽電池の長さごとに太陽電池セルの電流密度と出力電圧との関係を表す曲線を容易に求めることができ、太陽電池の寸法の設計が容易になる。また、数式解を求めることができるので、例えば表計算ソフトなどを用いて太陽電池セルの電流密度と出力電圧との関係を表す曲線を求めることができるようになり、コンピュータ資源及び労力の大幅な節約を実現することが可能になる。
【００２４】
図３は、近似式（２）を用いて非線形微分方程式（１）を線形化して解くことにより求めた、短絡電流値（ｊsc／ｉL）及び形状因子（ＦＦ）と太陽電池セルの長さとの関係を表す曲線を示す。図３において、実線で示されているのが、近似式（２）を用いて微分方程式（１）を数式解を求めることにより得られる短絡電流値及び形状因子と太陽電池セルの長さとの関係を表す曲線であり、破線で示されているのが、微分方程式（１）を従来の数値的手法で解くことにより得られる短絡電流値及び形状因子と太陽電池セルの長さとの関係を表す曲線である。図３から、近似式（２）を用いて微分方程式（１）を数式解を求めることにより得られる短絡電流値及び形状因子と太陽電池セルの長さとの関係を表す曲線と、微分方程式（１）を従来の数値的手法で解くことにより得られる短絡電流値及び形状因子と太陽電池セルの長さとの関係を表す曲線とはほぼ一致していることが分かる。なお、正確には、図３では、短絡電流値を理想太陽電池の光起電流ｉLで割ることにより無次元化した値を示している。
【００２５】
また、近似式（２）を用いて微分方程式（１）を解くにあたって、以下のようなパラメータＬp得られる。
【数８】
0005637398.tif 000009
太陽電池セルの長さがこのＬpとなる位置が図３に矢印で示されており、太陽電池セルの長さがＬpを超えると、短絡電流値（実際には短絡電流値ｊscを光起電流ｉLによって割った値）が急激に減衰すると共に、形状因子（ＦＦ）が減衰しなくなり飽和状態になることが分かる。すなわち、太陽電池セルの長さがＬpとなるときに太陽電池セルの最大出力を得ることができると言える。したがって、上記パラメータＬpは太陽電池セルの長さの設計指標として利用することができ、太陽電池セルの長さはＬp以下となるように選択することが好ましい。また、パラメータＬpは長さを所定の値にしたときの太陽電池セルの特性を求めることなく求められるので、Ｌpを設計指標として利用すれば、太陽電池セルの長さの設計が容易になる。
【００２６】
図４は、近似式（２）を用いることにより、短絡電流値（ｊsc／ｉL）及び形状因子（ＦＦ）と太陽電池セルの長さとの関係を表す曲線と共に、太陽電池セルの出力電力（正確には出力電力Ｐを長さＬ＝０の極限における出力電力Ｐ0によって割って無次元化した値）と太陽電池セルの長さとの関係を表す曲線とを示したグラフである。図４においても、図３と同様に、太陽電池セルの長さがＬpとなる位置が矢印で示されている。
【００２７】
図４から、太陽電池セルの出力電力Ｐは、太陽電池セルの長さがＬp以下の場合、直線的に減衰しており、太陽電池セルの長さＬ＝０の極限における形状因子をＦ0とすると、太陽電池セルの出力電力Ｐは、
Ｐ〜Ｐ０（１−１／（４Ｆ0））Ｌ／Ｌp）
と表すことができることが分かる。すなわち、太陽電池セルの長さがＬp以下の場合、太陽電池セルの特性を求めることなく、出力電力を上式から近似して求めることができる。また、Ｌpの前後で曲線の傾きが変わっており、Ｌpは、太陽電池セルの設計の際に、開口率と発電効率とのトレードオフ関係を考えるためにも有用となることが分かる。
【００２８】
また、上記パラメータＬpは、長さ方向だけでなく、幅方向（平面において長さ方向と垂直な方向）を含む２次元の特性解析にも適用することが可能である。図５（ａ）に示されているように、コの字状に補助電極で囲まれた太陽電池セルの寸法を設計する場合において、太陽電池セルの長さをＬ、幅をＷとするとき、Ｌ又はＷ／２をＬp以下になるように設計することにより、電流の急激な減衰を防止することができる。すなわち、図５（ｂ）の黒く塗りつぶされた領域（矢印で折れ線モデルと指された領域）から長さＬ及び幅Ｗの値を選択することにより、発電効率のよい太陽電池セルを設計することができる。また、本発明による太陽電池セルを設計する方法を用いて、補助電極の形状及び寸法に関する設計も行うことが可能になる。
【００２９】
図６〜図８は、近似式（２）を微分方程式（１）に適用して解くことにより得られた数式解から求めた、透明電極抵抗ρs、開放電圧Ｖoc、シャント抵抗ｒshが太陽電池セルの特性に与える影響を示すグラフである。
【００３０】
図６（ａ）及び（ｂ）は、透明電極抵抗ρsの値を変えたときの、太陽電池セルの長さＬと電流密度（図６（ａ））及び形状因子ＦＦ（図６（ｂ））との関係、図６（ｃ）はＬpと透明電極抵抗ρsとの関係を示し、図７（ａ）及び（ｂ）は、開放電圧Ｖocの値を変えたときの、太陽電池セルの長さＬと電流密度（図７（ａ））及び形状因子ＦＦ（図７（ｂ））との関係、図７（ｃ）はＬpと開放電圧Ｖocとの関係を示す。透明電極抵抗ρsや開放電極Ｖocの値に拘わらず、長さの増加に対して電流又は電流密度が急激に減衰する傾向及び長さの増加に対して形状因子ＦＦの減衰が飽和する傾向があり、電流又は電流密度が急激に減衰するポイント及び形状因子ＦＦがほぼ飽和するポイントとＬpとの間に相関がある。したがって、Ｌpを設計指標として使用できることが分かる。
【００３１】
図８（ａ）及び（ｂ）は、シャント抵抗ｒshの値を変えたときの、太陽電池セルの長さＬと電流密度（図８（ａ））及び形状因子ＦＦ（図８（ｂ））との関係、図８（ｃ）はパラメータＬpとシャント抵抗ｒshとの関係を示す。シャント抵抗ｒshの値が大きいとき（例えばｒsh≧５００Ωｃｍのとき）には、長さの増加に対して電流又は電流密度が急激に減衰する傾向及び長さの増加に対して形状因子ＦＦの減衰が飽和する傾向があり、電流又は電流密度が急激に減衰するポイント及び形状因子ＦＦがほぼ飽和するポイントとパラメータＬpとの間に相関がある。したがって、パラメータＬpを設計指標として使用できることが分かる。一方、シャント抵抗ｒshが小さいとき（例えばｒsh＜５００Ωｃｍのとき）には、上記のような傾向があてはまるとは言えない。このようにシャント抵抗ｒshが小さいときには、図８（ｃ）から、以下のような指標Ｌsを導入すれば、Ｌsが設計指標となることが分かる。
【数９】
0005637398.tif 000010
【００３２】
次に、本発明による太陽電池セルの設計方法の手順を説明する。
まず、透明電極抵抗と直列抵抗とシャント抵抗とダイオード成分とを用いた太陽電池セルを表す分布定数等価回路モデルから、太陽電池セルの特性を求めるための非線形微分方程式を求める。非線形微分方程式は例えば上述した式（１）のようになる。
【００３３】
次に、理想太陽電池（図１（ｂ）において破線で囲まれた部分によって構成されると仮定した太陽電池）における電流−電圧特性を求める。理想太陽電池の電流−電圧特性としては、設計しようとする太陽電池セルと同じ材料によって構成された十分に小さな寸法（例えば長さ２ｍｍ）の試験太陽電池セルに様々な電圧を印加したときの出力電流を測定することによって実験的に得た電流−電圧特性を用いることができる。
【００３４】
得られた電流−電圧特性の曲線から、ｉL、Ｖoc、ｒc、ｒsh、ｉ0をさらに定めることができる。ここで、ｉLは分布定数等価回路モデルにおける電流源の電流値である。詳細には、得られた太陽電池セルの電流−電圧曲線において電圧が０のときの電流の値からｉLが、電流が０のときの電圧の値からＶocが定まる。また、ｒsh、ｒcは、それぞれ、得られた試験太陽電池セルの電流−電圧曲線における電圧が０のときの曲線の傾き、電圧が開放電圧Ｖocに等しいときの曲線の傾きに等しいと近似することにより求めることができる。さらに、上述したようにして求めたｉLの値を用いて、得られた試験太陽電池セルの電流−電圧曲線に以下の式（３）をフィッティングすることにより、定数ｉ0を定めることができる。
【数１０】
0005637398.tif 000011
【００３５】
なお、透明電極抵抗ρsは、透明電極の材料の特性から求めればよい。
【００３６】
次に、得られた理想太陽電池の電圧Ｖdと電流ｉiscとの関係を表す曲線において、電圧が０のときの曲線の傾きをｇsc、電圧が開放電圧Ｖocに等しいときの曲線の傾きをｇocとして、該曲線を上述した近似式（２）で近似し、非線形微分方程式（１）と近似式（２）とを用いて、太陽電池セルの寸法と特性との関係、例えば、太陽電池セルの長さと、電流密度、形状因子又は出力電力との関係を求める。そして、求められた関係から、太陽電池セルに要求される発電電力が得られるように太陽電池セルの寸法を設計する。
【００３７】
非線形微分方程式（１）と近似式（２）とから、Ｌpを求め、Ｌp以下となるように、太陽電池セルの寸法を選択してもよい。Ｌp以下となるように太陽電池セルの寸法を選択すれば、発電効率のよい太陽電池セルを設計することができる。
【００３８】
本発明は、透明電極抵抗の抵抗値、理想太陽電池の電流−電圧特性データを入力すると、電池非線形微分方程式（１）と近似式（２）とを用いて、太陽電池セルの特性を求め、求められた太陽電池セルの特性をディスプレイなどの表示装置に表示するシミュレーションプログラムとして利用することもできる。
【図１】
0005637398.tif 000016
【図２】
0005637398.tif 000017
【図３】
0005637398.tif 000018
【図４】
0005637398.tif 000019
【図５】
0005637398.tif 000020
【図６】
0005637398.tif 000021
【図７】
0005637398.tif 000022
【図８】
0005637398.tif 000023
【図９】
0005637398.tif 000024
【図１０】
0005637398.tif 000025