FR 2811766 A1 20020118 FR 0009366 A 20000713 La présente invention concerne un procédé de caractérisation de la compacité d'objets tels que des cigarettes ou des filtres ainsi que du pouvoir de remplissage d'un volume par des produits tels que du tabac. D'une manière générale, on sait que les méthodes actuelles de caractérisation de la compacité des cigarettes, des filtres et du pouvoir de remplissage reposent sur un même principe physique: le produit à caractériser est placé entre deux mâchoires qui le soumettent à une force capable de provoquer un écrasement. L'amplitude de la déformation résultante est alors mesurée après un certain temps et est directement reliée aux propriétés caractéristiques communément appelée compacité et pouvoir de remplissage. Les méthodes actuelles assimilent donc ces propriétés caractéristiques à une amplitude d'écrasement. La forme des mâchoires, la valeur de la force appliquée et le temps au bout duquel est mesurée l'amplitude d'écrasement, sont fortement variables d'un fournisseur d'appareil à un autre, mais le principe reste le même. La caractérisation de la compacité ou du pouvoir de remplissage par une amplitude d'écrasement constitue un moyen simple mais est loin de fournir 2 2811766 une description détaillée du comportement du produit lorsqu'il est soumis à une force. En effet, l'écrasement dépend des déformations plastique et élastique, et de la relaxation des forces à l'intérieur du produit. Or, l'amplitude d'écrasement après un temps donné ne fournit aucune information sur ces trois phénomènes physiques. Les méthodes actuelles sont ainsi peu discriminantes puisque peu descriptives des propriétés mécaniques du produit. L'invention a donc plus particulièrement pour but de supprimer ces inconvénients. Elle propose à cet effet un procédé de caractérisation reposant sur une modélisation physique des phénomènes et l'ajustement de ce modèle à des données expérimentales de manière à caractériser le produit en ce qui concerne la compacité ou le pouvoir de remplissage. D'une façon plus précise, ce procédé comprend les étapes suivantes * La modélisation physique des phénomènes rentrant dans le processus d'écrasement, cette modélisation prenant en compte la déformation élastique, la déformation plastique et la relaxation. * L'ajustement du modèle obtenu lors de la phase de modélisation à des données expérimentales, et * La caractérisation du produit à partir du modèle ajusté grâce à des paramètres autres que la valeur d'amplitude d'écrasement, relatifs à la dissipation d'énergie lors de la déformation, la relaxation des forces ou même l'accumulation de l'énergie lors de la déformation. Avantageusement, la susdite phase d'ajustement pourra comprendre: ò L'application au produit d'une force connue d'écrasement, à l'aide d'un dispositif comprenant des mâchoires, un moyen d'appliquer une force sur ces mâchoires et un système d'enregistrement du déplacement relatif desdites mâchoires, ò Une phase d'acquisition des données dans laquelle on enregistre le déplacement au cours du temps, * Une phase de caractérisation physique par ajustement du modèle aux données expérimentales, à l'aide d'un outil de multirégression, * Une phase de description du produit à l'aide d'au moins trois paramètres caractéristiques de la déformation plastique et élastique et du phénomène de relaxation. On constate que ce procédé qui fournit des informations relatives à la déformation élastique, à la déformation plastique et à la relaxation lors du processus d'écrasement est plus discriminante que les méthodes utilisées jusqu'ici. De plus, selon ce procédé, l'information est extraite dès les premières secondes d'écrasement et est donc obtenue beaucoup plus rapidement qu'elle ne l'est avec les procédés actuels basés sur des mesures d'amplitude d'écrasement. Le procédé selon l'invention sera mieux compris à la lecture de la description ci-après donnée à titre d'exemple non limitatif, avec référence aux dessins annexés dans lesquels: La figure 1 est un schéma d'un circuit électrique constituant un modèle du comportement d'une cigarette, ou d'un échantillon de tabac, soumise à une force d'écrasement; La figure 2 est un diagramme montrant les courbes de modélisation de déplacements en fonction du temps de la mâchoire d'écrasement de la cigarette, ou de l'échantillon de tabac, sur laquelle on applique la force d'écrasement; La figure 3 montre une courbe d'écrasement en fonction du temps illustrant l'ajustement du modèle aux mesures d'écrasement enregistrées expérimentalement. Comme précédemment mentionné, le procédé selon l'invention comprend la modélisation du comportement de l'objet à étudier soumis à une force d'écrasement. A cet effet, la demanderesse a recherché les analogies qui peuvent exister entre les phénomènes physiques mis en oeuvre dans le processus d'écrasement et les phénomènes électriques dont les lois sont mieux connues et qui sont donc plus facilement exploitables. A ce titre, elle s'est basée sur les constatations suivantes * la déformation plastique, qui induit une perte d'énergie fournie, peut être représentée par une résistance électrique RI qui dissipe de l'énergie lorsqu'un courant la traverse, * La déformation élastique, qui représente un stockage d'énergie peut être représentée par une capacité C qui se charge ou se décharge lorsqu'elle est soumise à des variations de tension, * La relaxation qui constitue une perte d'énergie accumulée peut être représentée par une résistance R2 branchée en parallèle de la capacité induisant un courant de fuite. Le tableau I ci-après fait apparaître les analogies entre les phénomènes mécaniques et les phénomènes électriques précédemment évoqués. Tableau I Phénomènes mécaniques Phénomènes électriques Force F(t) Tension u(t) Vitesse v(t) Courant i(t) Energie fournie: ft F(t).v(t).dt Energie fournie ft u(t).i(t).dt Dissipation d'énergie Résistance R: ft R.I2 (t). dt Accumulation d'énergie Capacité 1 x C x u2 (t) Le modèle électrique équivalent à un objet tel qu'une cigarette, un filtre ou un échantillon de tabac, soumise à une force d'écrasement est représenté sur la figure 1. Ce modèle se présente sous la forme d'un circuit électrique monté entre deux bornes sur lesquelles est appliquée une tension u(t) et comprenant en série, d'une part, une résistance RI et, d'autre part, une capacité C sur laquelle est montée en parallèle une résistance R2. Le calcul de la réponse d'une cigarette ou du tabac à une force F(t) revient à calculer la réponse du circuit électrique à une tension u(t). Dans le cas général, ce calcul peut être effectué à l'aide de la transformation de Laplace. En supposant que la force d'amplitude F est appliquée instantanément, c'est à dire que la tension u(t) soit un échelon de tension, le courant i(t) circulant dans la résistance RI est de la forme: 6 2811766 jiU R2 R,+R, i(t) _ 1-eR xR2XC soit, par analogie, pour la vitesse de déplacement: Pl- +1P2 l v(t) =- P2 x 1- e PlxP2xp P3 Pl Pl +P2 et pour déplacement: d(t) = ft v(t).dt It i(t).dt o o d(t)= FP2'P3 x 1-e Pl P2 P3, + F xt (P. +p2)2 Pl +P2 avec pi: paramètre caractéristique de la dissipation d'énergie lors de la déformation P2: paramètre caractéristique de la relaxation des forces p3' paramètre caractéristique de l'accumulation d'énergie lors de la déformation Le modèle électrique proposé ci-dessus permet de calculer le déplacement quelle que soit la force appliquée. En voici un autre exemple: si nous supposons que la force F(t) n'est pas appliquée instantanément, mais qu'elle augmente linéairement avant d'atteindre sa valeur maximale F au bout d'un temps to, on obtient une relation plus complexe: 2 PIP2Xt2 axpxp' XP xa 2 a xP2 XP3 d(t) = a X Pl X p2 X p x e p p1 + x t + (Pl +P2 t)2.(p, +P2) (Pl +P2) 3 2 ( Pl+P2>"t-t0) 2 axp xp32X Px eP P2P3 2-1 +ax p2xp3 (t-to)xH(t-to) (pl +p2)3 2.(p +P2) (Pl +P2) Avec a: Vitesse d'augmentation de la tension entre t--=0 et t=--to (a=F/to) H(t): Fonction d'Heaviside ou d'échelon La modélisation du déplacement est tracée pour les deux cas de figure précédents, pour F=355, pl=4, P2=400, p3=0.4 et to=3s, sur la figure 2 dans laquelle la courbe CI correspond à l'application d'un échelon de force tandis que la courbe C2 est la réponse à une application progressive de la force. On remarque sur cette figure 2 qu'après une déformation majeure, le phénomène de relaxation induit un déplacement linéaire. Le modèle étant établi, il s'agit maintenant de l'exploiter pour caractériser le produit. Le dispositif expérimental employé reste le même que ceux déjà utilisés par les méthodes actuelles de mesure de compacité ou de pouvoir de remplissage. Par contre, avec la méthode proposée, il est indispensable de connaître ou d'enregistrer la force appliquée et le déplacement tout au long du processus d'écrasement. Connaissant la force appliquée au cours du temps, on simule alors à l'aide du modèle, le déplacement théorique, que l'on ajuste au déplacement mesuré. Il faut toutefois préciser qu'une loi d'application de force simple, comme un échelon de force par exemple, simplifie considérablement la modélisation du déplacement et donc l'étape de l'ajustement. L'exploitation des données se fait donc en ajustant le modèle théorique aux données expérimentales. L'ajustement du modèle théorique aux données expérimentales est réalisé à l'aide d'un outil informatique de multirégression non linéaire. Cet outil peut être basé sur différents algorithmes de minimisation tel que celui de Newton Raphson, de la plus grande pente ou celui de Levenberg-Marquardt. Dans ce dernier cas, il est possible d'exprimer l'incertitude sur les trois paramètres Pl, P2 et p3 du modèle. Un exemple est donné sur la figure 3, dans le cas d'une application d'un échelon de force et pour lequel les mesures enregistrées sont indiquées sur le tableau II ciaprès. Tableau II Temps Ecrasement 0,00 0,00 1,00 55,20 2,00 82,20 3,00 98,70 4,00 113,90 ,00 125,70 6,00 131,60 7,00 139,80 8,00 145,70 9,00 151,60 ,00 154,40 Sur la figure 3, la courbe C3 correspond au modèle tandis que les points P correspondent aux valeurs obtenues expérimentalement. Les valeurs des paramètres du modèle permettant d'ajuster les données expérimentales sont, à 95%, les suivantes pl = 4.4 +0.2 P2 = 48+ 5 P3 = 0.40 + 0.01 avec F=300. Ces trois paramètres caractérisent le comportement du produit lorsqu'il est soumis à l'écrasement. pl est un paramètre relatif à la déformation plastique. Plus il est grand et plus la déformation plastique est importante. P2 est un paramètre relatif à la relaxation. Plus il est grand et plus ce phénomène est faible. Enfin, P3 est un paramètre relatif à la déformation élastique. Plus il est grand et plus la déformation élastique est importante. Ces trois paramètres, obtenus grâce à ce nouveau procédé de mesure, permettent de caractériser beaucoup plus complètement le produit que ne le fait une valeur d'amplitude d'écrasement après un temps donné. De plus, on peut remarquer que le nouveau procédé permet de calculer a posteriori le résultat que donneraient les méthodes actuelles. La réciproque est fausse. Ceci démontre bien le caractère plus descriptif du procédé selon l'invention. On peut préciser que l'incertitude sur les paramètres peut être réduite en augmentant le nombre de mesures. Revendications 1. Procédé de caractérisation de la compacité d'objets tels que des cigarettes ou des filtres ainsi que du pouvoir de remplissage d'un volume par un produit tel que du tabac, caractérisé en ce qu'il comprend les étapes suivantes * La modélisation physique des phénomènes rentrant dans le processus d'écrasement de l'objet ou du produit, cette modélisation prenant en compte la déformation élastique, la déformation plastique et la relaxation. a L'ajustement du modèle obtenu lors de la phase de modélisation à des données expérimentales, et * La caractérisation de l'objet ou du produit à partir du modèle ajusté grâce à des paramètres autres que la valeur d'amplitude d'écrasement, relatifs à la dissipation d'énergie lors de la déformation, la relaxation des forces ou même l'accumulation de l'énergie lors de la déformation. 2. Procédé selon la revendication 2, caractérisé en ce que le modèle utilisé est celui d'un circuit électrique monté entre deux bornes sur lesquelles est appliquée une tension u(t) et comprenant en série, d'une part, une résistance RI et, d'autre part, une capacité sur laquelle est montée en parallèle une résistance R2. 3. Procédé selon l'une des revendications précédentes, caractérisé en ce que les susdites phases d'ajustement et de caractérisation comprennent: * L'application au produit d'une force connue d'écrasement, à l'aide d'un dispositif comprenant des mâchoires, un moyen d'appliquer une force sur 1l ces mâchoires et un système d'enregistrement du déplacement relatif desdites mâchoires, * Une phase d'acquisition des données dans laquelle on enregistre le déplacement au cours du temps, * Une phase de caractérisation physique par ajustement du modèle aux données expérimentales, à l'aide d'un outil de multirégression, * Une phase de description du produit à l'aide d'au moins trois paramètres caractéristiques de la déformation plastique et élastique et du phénomène de relaxation. 4. Procédé selon les revendications précédentes, caractérisé en ce que, dans le cas de l'application d'un échelon de force d'amplitude F, le susdit modèle est de la forme: _plp2 X f.2 -- xI d(t) = F.p2p3 x 1 -e p. p2p3 + X t (Pl + P2)2 Pl +p2 avec Pi: paramètre caractéristique de la dissipation d'énergie lors de la déformation P2: paramètre caractéristique de la relaxation des forces p3: paramètre caractéristique de l'accumulation d'énergie lors de la déformation d(t): déplacement en fonction du temps F: valeur maximale de la force appliquée F(t) 5. Procédé selon la revendication 3, caractérisé en ce que l'ajustement du susdit modèle aux données expérimentales est réalisé à l'aide d'un outil informatique de multirégression non linéaire basé sur un algorithme de multirégression tel que celui de Newton Raphson, de la plus grande pente ou celui de Levenberg-Marquardt permettant d'exprimer l'incertitude sur les trois paramètres Pl, P2, p3 du modèle.