毛髪の断面形状を測定する形状測定ステップと、前記毛髪の一端を固定し、前記断面形状の所定の方向に沿った前記毛髪の荷重と変位との関係を測定する試験ステップと、前記毛髪の荷重と変位との関係から、前記毛髪の柔軟性を評価する評価ステップとを、有することを特徴とする毛髪の柔軟性評価方法。前記形状測定ステップは、前記断面形状を楕円近似したときの長軸と短軸とを求め、前記試験ステップは、少なくとも前記長軸および前記短軸のいずれか一方に沿った前記毛髪の荷重と変位との関係を測定することを、特徴とする請求項1記載の毛髪の柔軟性評価方法。前記試験ステップは、前記毛髪の固定位置から所定の距離lの位置での前記毛髪の荷重Pと変位との関係を測定し、前記評価ステップは、変位とPl3/3IIは荷重Pの負荷方向に対する前記毛髪の断面二次モーメントとの関係から、前記毛髪の柔軟性を評価することを、特徴とする請求項1または2記載の毛髪の柔軟性評価方法。前記毛髪は複数の層から成る断面構造を有し、前記評価ステップは、横軸を変位、縦軸をPl3/3Iとして測定値をプロットしたときのグラフの傾きSEを求め、前記傾きSEから前記毛髪の各層の弾性率および/または厚みを評価することを、特徴とする請求項3記載の毛髪の柔軟性評価方法。前記形状測定ステップは、非接触で前記毛髪の断面形状を測定し、前記試験ステップは、前記毛髪の固定位置から1~3mm離れた位置で、150μN以下の荷重範囲での前記毛髪の荷重と変位との関係を測定することを、特徴とする請求項1、2、3または4記載の毛髪の柔軟性評価方法。非特許文献7に記載のように、毛髪の曲げ試験では、柔軟性を評価するときに、毛髪の断面形状を考慮することが重要であることは指摘されているが、実際に毛髪の断面形状を考慮して柔軟性を高精度に評価したものは存在しない。また、曲げ試験により直接測定される弾性率と毛髪の内部構造との関連性については、これまで具体的に考慮されたことはなかった。本発明は、このような課題に着目してなされたもので、毛髪の断面形状を考慮して柔軟性を高精度に評価することができ、さらに弾性率と毛髪の内部構造とを関連付けて評価を行うことができる毛髪の柔軟性評価方法を提供することを目的としている。上記目的を達成するために、本発明に係る毛髪の柔軟性評価方法は、毛髪の断面形状を測定する形状測定ステップと、前記毛髪の一端を固定し、前記断面形状の所定の方向に沿った前記毛髪の荷重と変位との関係を測定する試験ステップと、前記毛髪の荷重と変位との関係から、前記毛髪の柔軟性を評価する評価ステップとを、有することを特徴とする。本発明に係る毛髪の柔軟性評価方法は、測定した毛髪の断面形状に基づいて、その断面形状の所定の方向に沿った毛髪の荷重と変位との関係を測定するため、毛髪の断面形状を考慮して柔軟性を評価することができる。毛髪の断面形状に対して、常に所定の方向に沿った測定を行うことにより、測定間の幾何学的条件を揃えることができ、統一した基準での毛髪の柔軟性評価を行うことができる。これにより、毛髪処理等に対する毛髪の柔軟性のわずかな変化でも評価することができ、評価精度を高めることができる。本発明に係る毛髪の柔軟性評価方法で、毛髪の荷重と変位との関係は、毛髪の曲げ試験により測定することができる。試験ステップで、毛髪の一端を固定する方法は、いかなる方法であってもよく、例えば、硬化性の樹脂を使用する方法であってもよい。また、断面形状の所定の方向に沿って毛髪の荷重と変位との関係を測定するときには、測定精度を高めるために、顕微鏡で毛髪の形状と方向とを観察しながら行うのが好ましい。なお、ここでの毛髪とは、頭髪だけでなく、眉毛等も含むものである。本発明に係る毛髪の柔軟性評価方法で、前記形状測定ステップは、例えば、前記断面形状を楕円近似したときの長軸と短軸とを求め、前記試験ステップは、少なくとも前記長軸および前記短軸のいずれか一方に沿った前記毛髪の荷重と変位との関係を測定することが好ましい。毛髪の断面形状は、人種や年齢、性別、毛髪の痛み具合などにより、さまざまに異なっていると考えられ、適切な断面形状を仮定して寸法の測定と試験とを実施すると共に、断面形状を考慮して断面二次モーメントを求めて最終的な柔軟性の評価を行うことが肝要である。本発明に係る毛髪の柔軟性評価方法で、前記試験ステップは、前記毛髪の固定位置から所定の距離lの位置での前記毛髪の荷重Pと変位との関係を測定し、前記評価ステップは、変位とPl3/3IIは荷重Pの負荷方向に対する前記毛髪の断面二次モーメントとの関係から、前記毛髪の柔軟性を評価することが好ましい。特に、本発明に係る毛髪の柔軟性評価方法で、前記毛髪は複数の層から成る断面構造を有し、前記評価ステップは、横軸を変位、縦軸をPl3/3Iとして測定値をプロットしたときのグラフの傾きSEを求め、前記傾きSEから前記毛髪の各層の弾性率および/または厚みを評価してもよい。この場合、本発明に係る毛髪の柔軟性評価方法は、以下の測定原理により毛髪の柔軟性を評価することができる。図3に示すように、一端が固定された、片持ちはり構造に対する毛髪の曲げ試験を考える。曲げ試験に使用する測定装置は、本発明者等により開発された特許文献1に記載の荷重測定装置を使用することができる。ここで、毛髪は、断面形状が楕円形状を成し、3層から成る断面構造を有するものとする。毛髪の変形が微小であり、荷重点の移動がないものすれば、以下の式が成立する。ここで、ρはxにおける曲率半径、Mは曲げモーメント、Aは曲げ剛性、E、Iはそれぞれ縦弾性率、断面二次モーメントである。下添え字1、2、3はそれぞれ3層構造の中心部、中間層、最外層を表す。長軸2a、短軸2bの楕円断面に対する短軸方向負荷に対する断面二次モーメントIは、である。なお、長軸方向負荷に対する断面二次モーメントIは、aとbとを入れ替えて得られる。微分方程式1に式2を代入し、これを境界条件dy/dx|x=0=y|x=0のもとで解くと、次式が得られる。ここで、荷重Pと変位yとの関係が負荷方向に依存することに注意が必要である。lは毛髪の固定端から荷重点までの距離である。得られる毛髪に対する荷重-変位線図の傾きより、Aを得ることができる。しかしながら、式2に示すように、Aには3つのEE1、E2、E3に加えて、3つのII1、I2、I3が含まれており、Aの物理的意味を理解するのは困難である。仮に、6つの寸法各層の長軸および短軸を断面観察等により知り得たとしても、単一試験より3つのEを分離決定することは困難である。そこで、直接得られるAの力学的意味を考えることとする。最外層の長軸、および短軸をそれぞれ2a、2bとする。また、中間層および中心部の長軸、および短軸の長さがそれぞれ2aと2bのα倍、β倍定数であるとする。この場合、各層の断面二次モーメントは、I3=I1-α4、I2=Iα4-β4、I1=Iβ4となる。また、各層の縦弾性率を規準となるEの値、ESを用いて、E3=ξES、E2=ψES、E1=ζESと表す。ここでξ、ψ、ζは定数である。この場合、Aは次式で表される。式5において、A/Iに着目すれば、これは負荷方向によらない毛髪構造の弾性率を示すことがわかる。そこでこれをSEstructuralelastisityと定義する。このように定義したSEは、Pl3/3I-変位関係の傾きに一致し、実験により求めることができる。SEは、各層のEと厚みとに関係する。外形寸法の相違は考慮されているため、SEは個々の毛髪サンプルを比較するのにも適している。曲げ試験時には毛髪外表面に高い応力場が形成されるが、測定されるSEには内部層の特性も反映されることに注意が必要である。このように、新たな評価基準SEを用い、弾性率と毛髪の内部構造とを関連付けて、毛髪の柔軟性を評価することができる。本発明に係る毛髪の柔軟性評価方法で、前記形状測定ステップは、非接触で前記毛髪の断面形状を測定し、前記試験ステップは、前記毛髪の固定位置から1~3mm離れた位置で、150μN以下の荷重範囲での前記毛髪の荷重と変位との関係を測定することが好ましい。この場合、試験試料の毛髪を損傷することなく、断面形状を測定することができ、試験ステップでの測定精度の低下を防ぐことができる。また、試験ステップで、毛髪が弾性変形する極微小区間での測定を行うことができる。断面形状の測定は、非接触であればいかなる方法で行ってもよく、例えば、レーザーを用いて行うことができる。本発明によれば、毛髪の断面形状を考慮して柔軟性を評価することができ、さらに弾性率と毛髪の内部構造とを関連付けて評価を行うことができる毛髪の柔軟性評価方法を提供することができる。以下、図面に基づき、本発明の実施の形態について説明する。図1は、本発明の実施の形態の毛髪の柔軟性評価方法で使用する荷重測定装置を示している。荷重測定装置は、本発明者等により開発された特許文献1に記載の荷重測定装置であって、毛髪の曲げ試験を実施可能である。図1に示すように、荷重測定装置10は、支持部材11と梁部材12と荷重部13と変位センサ14とを有している。支持部材11は、3方向に位置を調整可能なXYZステージ21と、XYZステージ21の上に載置された、水平2方向に微小移動可能なピエゾステージ22と、ピエゾステージ22の上に載置された、L型の側面形状を成す長さ調整ジグ23と、長さ調整ジグ23の上に取り付けられたセンサホルダ24とを有している。梁部材12は、一端が連結された二重はり構造カレンチバーから成り、他端が長さ調整ジグ23の下端部に取り付けられている。荷重部13は、針状のチップから成り、梁部材12の一端に固定されている。梁部材12は、長さ方向に沿ってスライドして、固定位置と荷重部13との距離を変更可能に、支持部材11に支持されている。変位センサ14は、静電容量センサから成り、センサホルダ24の先端に取り付けられている。変位センサ14は、梁部材12の一端部の表面との間に隙間をあけて、その表面に対向するよう取り付けられ、その表面との距離の変化を測定可能になっている。これにより、変位センサ14は、荷重部13の突出方向の変位を測定可能になっている。以下に、本発明の実施の形態の毛髪の柔軟性評価方法について、実際の毛髪の試験例に基づいて説明する。試験試料Sampleとして、同一人物の3本の毛髪からそれぞれ5mm程度の3試料合計9試料を採取した。シャープペンシルのペン先に硬化性樹脂をポンプ注入した後、毛髪試料を挿入することにより、毛髪試料の一端を固定した。シャープペンシルを固定し、顕微鏡で観察しながら、一端が固定された毛髪試料を36度刻みで回転させて、断面形状を測定した。測定した断面形状を楕円近似し、各毛髪試料の長軸2aと短軸2bとを求めた。その結果を、表1に示す。次に、図1に示す荷重測定装置10を使用して、各毛髪試料の曲げ試験を行った。荷重部13により、各毛髪試料の固定端からの距離lが2mmおよび2.5mmの位置に荷重を作用させ、梁部材12のたわみを変位センサ14で測定した。荷重範囲を150μN以下とし、毛髪が弾性変形する範囲で曲げ試験を行った。各毛髪試料の変形は、y=δS-δCで与えられる。ここで、δSはピエゾステージ22の変位であり、δCは梁部材12のたわみである。毛髪試料に作用する力は、P=kδCで与えられる。ここで、kはばね定数で、試験時のkは、26.8μN/μmであった。曲げ試験は、各毛髪試料に対して、短軸および長軸に平行な2方向で、それぞれ2箇所固定端からの距離l=2mmおよび2.5mmの位置への負荷で行った。負荷をかける際には、顕微鏡で毛髪の形状と方向とを観察しながら行った。曲げ試験で得られた荷重Pと変位yとの関係の一例を、図2aに示す。得られた荷重Pと変位yとの関係から、横軸を変位y、縦軸をPl3/3Iとして測定値をプロットしたときのグラフを描く。図2aに示したものについて描いたグラフを、図2bに示す。図2bに示すように、異なる負荷方向により得られたデータが良く一致しているのが確認できる。全ての毛髪試料について、Pl3/3I-変位y曲線の傾きSEを求め、表1に示す。ここで、l=2mmまたは2.5mmであり、Iは荷重Pの負荷方向に対する毛髪の断面二次モーメントであり、測定された毛髪の断面形状に基づいて、3式から算出することができる。図2aおよびbから、異なる負荷方向に対する荷重-変位曲線の差異は、毛髪試料の断面二次モーメントによるものであり、断面二次モーメントを考慮した新たな柔軟性の評価基準SEを使用すると、毛髪の柔軟性が負荷方向によらないことが確認された。6式に示すように、SEは毛髪の内部構造を反映しているため、SEを用いることにより、弾性率と毛髪の内部構造とを関連付けて、毛髪の柔軟性を評価することができる。また、SEは、外形寸法の相違も考慮されているため、個々の毛髪サンプルを比較するのにも適している。表1に示すように、毛髪試料A1~A3のSEの値はそれぞれ、5.12±0.18、5.49±0.22、5.16±0.18GPaであり、誤差5未満の高精度なSE測定を行うことができた。全毛髪試料のSEは、5.26±0.26GPaであり、値のバラつきは5未満であった。また、毛髪試料A2のSEは、他のSEに比べて少し大きいことが確認できた。このことから、毛髪処理等によるSEの僅かな変化を検出するには、同一の毛髪を使用する必要があるといえる。ここで、毛髪の形状がSEに及ぼす影響について検討する。6式中のα=0.9、β=0.1、ξ=1、ψ=0.25、ζ=0.4の場合を考える。この場合、6式から、SE=0.50ESとなる。ここで、αの値のみ0.8と変化させると、SE=0.69ESとなり、毛髪の各層の幾何学的形状変化により、SEの値が大きく変動することがわかる。このことから、SEを用いることにより、毛髪の各層の弾性率や厚みを評価することができることがわかる。本発明の実施の形態の毛髪の柔軟性評価方法を使用して、毛髪処理による毛髪の柔軟性の変化について評価を行った。評価には、3つの毛髪試料A1-2、A2-1、A2-2を使用し、各毛髪試料に毛髪処理を施した。毛髪処理として、WT30秒の水処理、ST30秒のシャンプー処理と30秒の水処理、CT30秒のコンディショナー処理と30秒の水処理を行った。各毛髪処理が終わってから40分経過した後に、各毛髪試料についてSEを測定した。測定したSEを、処理前の初期値で規格化した値これをRとするを、表2に示す。表2に示すように、WT後のRはおよそ0.9であり、WTにより毛髪は10程度柔らかくなっている。また、ST後のRもおよそ0.9であり、STによってRはそれほど変わっていない。CT後のRはおよそ0.75である。CT後のSEの変化の内、10はWTによりもたらされたものである。従って、15がCTによる効果である。なお、SEの値は、少なくとも18日経過後には、元に戻った。また、WTを3回繰返したところ、Rは0.9から変化しなかった。ここで、毛髪の内部構造との関係について検討する。各毛髪処理により毛髪の各層に幾何学的変化がないと仮定する。すなわち、SEの変化が各層の弾性率の変化によるものとする。また、簡単のため、二層構造β=0を考える。この場合、毛髪の最外層および中間層に対する変化モデルは、それぞれとなる。ここで、ξ0およびψ0は、毛髪処理前のξおよびψの初期値である。α=0.9、ξ0=1、ψ0=0.25と仮定する。この場合、7式および8式から、R=0.9、0.75に対して、ξ/ξ0=0.85、0.63、ψ/ψ0=0.69、0.23と求まる。このことから、実際に観察されるSEの変化以上の変化が各層に生じていることがわかる。