梁 宇,李 慧
(北京服装学院 信息工程学院,北京 100029)
随着人们生活水平的不断提高,服装的压力舒适性也开始备受关注。关于服装压力舒适性的研究主要分为真实服装压力和虚拟服装压力2类, 目前关于真实和虚拟的服装压力研究,主要集中在文胸、瑜伽服、束缚带、袜子等贴身衣物上,研究的人体部位主要集中在胸、腰、臀、腿等,虽然有一定的研究成果,但在服装行业,服装设计依然主要依赖于设计师的实践经验和大量的试错。通过仿真研究建立可靠的物理模型,可以对服装压力进行仿真预测。推进服装压力仿真研究,可以减少服装设计中不必要的试错,很大程度避免一些危害人体健康的试验。
纵观前人成果,服装压力舒适性的研究涉及力学、人体工学、纺织材料学等学科领域。关于单方面的生理研究在日本和一些欧美国已经有比较丰硕的成果,真实服装压力研究主要包括服装压力分布研究、服装压力舒适阈值研究及服装压力舒适性的影响因素研究。真实服装压力研究是虚拟服装压力研究的基础。
1991年,MAKAB[1-2]对人体穿着文胸后多个部位服装压力进行测量及分析,找出8个舒适性不佳的部位。Shi mizu等[3]对文胸的压力分布进行了研究,比较在静态和动态下人体8个部位的压力。研究发现,在运动和静止状态下,人体最大受压点略有不同,并且同一个部位的压力值在运动状态下要大于该部位在静止状态的压力值[4]。李世霞[5]研究发现文胸对人体压力较大的部位为钢托、肩部、侧缝和后侧等。
由于人体曲率的改变会影响到各部位的压力值[6-7],不同体型、不同运动状态下人体的各部位受压情况也存在着较大差异。Denton[8]主要对令人体感觉舒适性不佳的服装压力进行研究,研究发现令人感觉舒适度不佳的范围为5.88~9.8 k Pa,与人体皮表毛细血管压力平均值7.85 k Pa相吻合。当服装压力超过皮表毛细血管压力临界值时,血液被迫流向人体穿戴压力较小的身体部位,一般为腿部等较低部位,从而造成这些部位肿胀[9]。刘红[10]主要对女性穿着运动背心时胸部的服装压力进行测量与分析研究,得到胸部的压力舒适范围为0.95~1.36 k Pa,背部为0.59~0.9 k Pa,肩部为1.135~1.64 k Pa,侧缝处为0.48~0.73 k Pa,腹部为0.29~0.43 k Pa。岑司竹[11]结合血流量和运动心率方面的理论知识,对女性穿着文胸时胸部的服装压力进行测量与分析,研究发现侧颈点服装舒适压力范围为4.15~5.67 k Pa,后背肩胛骨提肌点为3.23~4.48 k Pa,BP点为2.56~3.97 k Pa,乳房侧下点为2.02~3.33 k Pa,乳房下缘点为1.69~2.74k Pa,侧中点为1.63~2.18 k Pa,下胸围上侧缝到脊椎1/4点为3.98~4.79 k Pa,下胸围上脊椎与侧缝的中点为4.04~4.91 k Pa。张文斌等[12]发现在穿戴不同款式服装的时候,压力舒适性范围存在不同,在对部分服装款式的压力舒适性范围开展研究后得到泳衣压力舒适范围为0.98~1.96 k Pa,紧身胸衣舒适范围为2.94~4.90 k Pa,医用长袜舒适范围为2.94~5.88k Pa,紧身服舒适范围为小于1.96 k Pa。王永荣等[13]选取女性上身8处测试点,利用束带对舒适服装压力阈值进行研究,肩部、胸部、腹部3个部位测得的压力舒适阈限分别为1.7、0.98、1.34 k Pa,小于以往学者测得的阈限。图1为近年来学者研究着衣后不同人体部位舒适服装压力平均值。
图1 不同人体部位舒适服装压力平均值
有限元法[14-17](Finite Ele ment Method,简称FEM),起源于航空工程中的矩阵分析,此算法的本质是一种求解偏微分方程的算法。其中,偏微分方程通常被用于求解时空相关问题的物理过程。绝大多数情况下,这些偏微分方程的解可能无法求出。不过,可以通过对所需求解域进行离散化构造近似的方程,来获取与偏微分方程类似的数值模型方程,并通过数值方法进行求解。由此得到的解便称作近似解,与偏微分方程真实解相对应,有限元法的主要工作便是求解这些近似解。
网格划分越稠密时,求得近似解越贴近真实解。有限单元算法的核心内容是通过离散化分解,将复杂的求解域模型剖分成有限个细小单元拼合而成的结构。计算机通过对每个单元的力学特性进行分析,将剖分好的单元在有限元算法的规则下,拼接为一个整体,通过计算得到所需求解模型的力学特性。有限元算法的主要研究对象是单元,在进行网格划分时最需要注意的一个问题是相邻单元的参数要求具有一定的连续性。
根据前人研究,织物的多种性能对服装压力都会产生一定影响[18-20],但对人体下胸围受压的应力应变分析中,需将复杂的实际情况抽象为相对简单的力学模型,否则问题将会变得极其复杂难以处理,采用有限元法计算出相对准确可靠的结果,步骤如下。
步骤1:合理的网格剖分,对该模型待求解区域进行网格划分,将整个区域分割成若干离散的单元。对于模型在各方向梯度下的离散化公式如式(1)和(2)所示。
步骤2:单元分析,将每个连接在一起的单元节点作为插值点,进行分片插值,建立线性插值函数。设人体应力为S,应变为∈,模型各方向应力求解偏微分方程为:
模型各方向应变求解偏微分方程为:
式中:u是人体受应力后的位移。
步骤3:求解近似变分方程,输出近似解求解结果。
21世纪初,便有学者将有限元理论的仿真研究应用在服装压力舒适性领域中,取得了一定的成果,近几年由于仿真技术的发展,对虚拟服装压力的仿真研究也更进一步。
Yeung K W 等[21]结合弹性力学理论和三维建模技术,建立了人体及服装动态平衡三维力学模型,其中人体模型由硬骨头、皮肤和软组织3层构成,人体在动态情况下服装压力分布可由该数学模型模拟出来,Yeung K W 通过该模型得到了穿着紧身裤时从脚踝到腰部的压力分布。李毅等[22]基于织物力学特性分析人与服装动态接触过程中的力学分布,研发出服装机械力学系统来进行服装压力预测。覃蕊等[23]运用有限元方法来模拟男短袜袜口对小腿部施压后压力与位移的关系。在研究中,人体被视为弹性体,腿部和袜口之间的接触被视为弹性接触。朱珊[24]通过有限元法对人体腰部所受服装压力做了分析研究,在ANSYS 软件中对人体腰部的应力分析进行了仿真研究。Zhang M 等[25]通过位移分布评价压力舒适度,利用有限元法和曲线拟合探索小腿横截面上压力与位移的多重关系。
2019年Dan R 等[26]以连裤袜腰部为研究对象,构建人体与连裤袜的接触为弹力接触,用ANSYS 模拟了弹力连裤袜样品在穿衣过程中腰部相应的位移分布。2021年Dan R 等[27]通过有限元模拟曲线拟合,深入研究人体穿着弹性连裤袜在运动状态下,压力/位移比与角度之间的函数关系。Shichen Zhang 等[28]在ANSYS中通过有限元模型的全因子分析,评价了不同设计特征对文胸肩带的影响。得出拉伸对文胸肩带运动位移影响最大,对杨氏模量影响最小的结论。
基于不同应用领域,目前常见的有限元软件有COMSOL Multiphysics、ANSYS、L USAS、Al gor、MSC Nastran、ABAQUS、Hyper mesh、L MS-Sam tech、Fe map/NX Nastran和FEPG 等软件。ANSYS是目前最流行的有限元软件[29],其最大特征是可以根据具体需求选择是否使用有限元法,但在几何建模方面存在不足。它通过代码对几何进行建模,这使得所建模型不够准确,并且在求解所用公式的修改方面极为不便。而上述ANSYS的不足之处,在COMSOL中都可以得到很好的解决,在COMSOL 中仅使用有限元法,它本身具有开发几何研究的界面。COMSOL 多物理仿真模拟也为研究复杂问题提供了极强的解决方案。在实际情况中,通过物理过程描述的现象并不是一个简单的在一个物理场下的问题。例如:考虑影响服装压的环境因素,当人体处在温度变化的条件下,在求解时便需要进行固体力学模型和传热模型耦合。图2为在COMSOL中模拟人体下胸围侧中点受压后,该点周围的压力分布,图3 为根据仿真所得数据在MATLAB中绘制人体下胸围侧中点加载后应力分布三维图。
图2 在一点上加载后模型的应力分布
图3 人体最大受压点加载后应力分布
近几年关于真实和虚拟服装压力的研究虽然有了一定的成果,仍存在着不足。
(1)真实服装压力测量方面。现有研究中对关于真实的压力测量研究主要还是静态压力,由于测试压力设备的局限,对运动状态下服装压力的研究还没有普遍开展。
(2)服装压力仿真方面。早期关于虚拟压力的仿真研究普遍存在着模拟的运动过程单一、预测的压力分布误差较大等问题。近年来,由于三维扫描技术的发展,在建模精度方面有了一定的提高,但是大部分学者多是在ANSYS、ABAQUS及CLO 中进行一些虚拟服装压力的研究,目前为止很少有人通过COMSOL对服装压力进行仿真研究,对多物理场耦合仿真也很少涉足。由于人体体型复杂,人体的内部生理环境是一个极为复杂的生态系统,所以现有的仿真研究大都将人体看作刚体来建模,仿真结果也会与真实服装压力存在较大误差。
如今服装行业竞争越来越激烈,服装是否有“黑科技”、是否有足够的专业理论支撑也是相当一部分消费者在购买服装时重点考虑的因素,推进仿真研究,为服装产品提供可靠的理论依据,可以大大减少产品研发中人力资源成本。在以后的研究中还需运用科学手段获取运动状态下的服装压力,提高建模的精确性、准确性和完整性,对人体与服装的关系进行深入研究。