Cette thèse présente le développement d'un modèle numérique avancé, capable de simuler les interactions entre des vagues de surface de cambrure quelconque et des corps rigides immergés ayant des mouvements de grande amplitude. Fondé sur la théorie potentielle, il propose une résolution couplée de la dynamique vagues/structure par la méthode implicite de Van Daalen (1993), encore appelée méthode du potentiel d'accélération par Tanizawa (1995). La précision du modèle à deux dimensions est testée sur un ensemble d'applications impliquant le mouvement forcé ou libre d'un cylindre horizontal immergé, de section circulaire : diffraction par un cylindre fixe, radiation par un cylindre en mouvement forcé de grande amplitude, absorption des vagues par le cylindre de Bristol. Pour chaque application, les résultats numériques sont comparés à des résultats expérimentaux ou analytiques issus de la théorie linéaire, avec un bon accord en particulier pour les petites amplitudes de mouvement du cylindre et pour les vagues de faibles cambrures. La génération de vagues irrégulières et la prise en compte d'un second corps cylindrique immergé sont ensuite intégrées au modèle, et illustrées sur des applications pratiques avec des systèmes récupérateurs d'énergie des vagues simples. Enfin, le modèle est étendu en trois dimensions avec des premières applications au cas d'une sphère décrivant des mouvements de grande amplitude / This PhD is dedicated to the development of an advanced numerical model for simulating interactions between free surface waves of arbitrary steepness and rigid bodies in high amplitude motions. Based on potential theory, it solves the coupled dynamics of waves and structure with the implicit method by Van Daalen (1993), also named the acceleration potential method by Tanizawa (1995). The precision of this two-dimensional model is tested on a wide range of applications involving the forced motion or free motion of a submerged horizontal cylinder of circular cross-section : diffraction by a fixed cylinder, radiation by a cylinder in specified high amplitude motions, wave absorption by the Bristol cylinder. In each of these applications, numerical results are compared to experimental data or analytical solutions based on the linear wave theory, with a good agreement especially for small amplitude motions of the cylinder and small wave steepnesses. The irregular wave generation by a paddle and the possibility to add an extra circular cylinder are integrated in the model and illustrated on practical applications with simple wave energy converters. The model is finally extended to three dimensions, with preliminary results for a sphere in large amplitude heaving oscillations

La position assise peut être source d'inconfort, particulièrement en avion lors des vols longcourriers. Cet inconfort provient en partie de facteurs mécaniques liés à l'interaction entre le siège et le passager. Disposer de modèles biomécaniques pouvant simuler cette interaction et estimer ces facteurs permettrait d'optimiser le design du siège d'avion lors de sa phase de conception afin d'améliorer son ergonomie et réduire l'inconfort du passager. L'objectif de cette thèse est de développer un outil numérique permettant d'estimer les facteurs mécaniques menant à l'inconfort et la fatigue des passagers. Cet outil combine deux différents types de modèles. Un modèle éléments finis permettant de simuler la déformation des tissus sous-cutanés et un modèle corps rigides permettant d'estimer les efforts musculaires et articulaires. Une méthode de couplage des deux modèles a été développée permettant ainsi de simuler une position à partir de laquelle l'ensemble des facteurs d'inconfort sont estimés. Un modèle éléments finis a été développé après avoir fait une étude de sensibilité sur les différents paramètres de modélisation (maillage, géométrie, lois matériaux). Un modèle corps rigides développé par Anybody a été utilisé pour être couplé avec ce modèle éléments finis. Cette méthode de couplage itératif entre les deux modèles a permis de réaliser un ajustement de la posture initiale dans le siège. Afin de simuler l'ensemble de la population et sa grande diversité morphologique, un modèle surfacique paramétrique a été développé à partir de données 3d expérimentales, ce modèle surfacique permettant ainsi d'obtenir un modèle éléments finis représentant tout type d'anthropométrie. Différents processus de validation ont été effectués à l'aide de données et d'un modèle « sujetspécifique ». Les données de pression externe simulées ont été comparées à des données expérimentales. Une étude expérimentale sous IRM ouvert a permis de mesurer les déformations des différents tissus sous-cutanés afin de les comparer aux données simulées. Un outil numérique est donc aujourd'hui disponible pour simuler l'impact du siège sur les passagers, cependant de futures études devraient se concentrer d'une part sur les modèles en étudiant la variation morphologique interne inter individus, le positionnement dans le siège ainsi que l'influence du temps sur les tissus mous et d'autres part sur la définition de critères d'inconfort et de fatigue (inconfort considéré sur des temps longs représentatifs d'un vol long-courrier) / The sitting position could be a source of discomfort, in particular in a long haul flight. This discomfort comes partially from mechanical factors linked to the interaction between the passenger and the seat. Having biomechanical models which can simulate this interaction and estimate these factors would allow optimizing the seat design in its conception phase to improve it ergonomic quality and reduce the passenger discomfort. The objective of this thesis is to develop a digital tool allowing estimating the mechanical factors leading to discomfort and fatigue of the passenger. This tool assemble two kinds of models, a finite element model allowing to simulate the sub dermal tissue deformation and a multibody model allowing to estimate the muscular and joint forces. A coupling method of the two models have been developed allowing simulating a position from where the all discomfort factors are estimated. A finite element model has been developed after having done a sensitivity analysis on the different model parameters (mesh, geometry, material law). A multibody model developed by Anybody was used to be associated with this finite element model. This iterative coupling method between the two models allowed realizing an adjustment of the initial posture in the seat. To simulate the whole population and is large morphological diversity, a parametric shape model was developed from 3d experimental data, this shape model allowing to obtain a finite element model representing any kind of anthropometry. Different validation processes have been realized with experimental data and subject-specific model. The simulated extern pressures were compared to experimental data. An experimental study done in an open MRI allowed to measure the different subcutaneous tissue to compare it to the simulated data. A digital tool is consequently now available to simulate the impact of the seat on the passenger, however future studies should focus on the one hand on the models studying the internal morphological variations between people, the person positioning in the seat, the time influence on the soft tissue and on the other hand on the discomfort and fatigue criteria

Les assemblages d'accéléromètres sont utilisés entre autres en biomécanique afin d'estimer le champ des accélérations d'un corps rigide et, de là, les variables cinématiques décrivant les déplacements dudit corps rigide. On remarque toutefois que le progrès de cette technologie a été limité par celui des gyroscopes micro-usinés, lesquels sont plus précis que les assemblages d'accéléromètres dans la plupart des applications. Le but des recherches décrites dans cette thèse est d'améliorer la précision des estimés de vitesse angulaire produits par les assemblages d'accéléromètres. Les développements proposés portent sur deux aspects : une classe d'architectures mécaniques d'accéléromètres est proposée ayant pour but de réduire la sensibilité des capteurs aux accélérations transverses, des algorithmes robustes étant également proposés afin d'estimer la vitesse angulaire à partir de certaines ou de toutes les composantes du champ d'accélération du corps rigide. La classe d'accéléromètres proposée est inspirée des architectures de machines à cinématique parallèle (MCP) en assimilant la masse d'épreuve et le support d'un accéléromètre à l'organe terminal et à la base d'une MCP, respectivement. Une caractéristique commune aux MCP proposées ici est que leurs organes terminaux et leurs bases sont reliées par n+1 chaînes cinématiques simples, n = 1, 2, 3 étant leurs nombres de directions sensibles respectifs. Pour cette raison, on appelle accéléromètres simpliciaux multiaxiaux les capteurs résultant de ces MCP. Une version microusinée de l'accéléromètre simplicial biaxial est conçue, fabriquée et testée. En outre, la théorie permettant l'estimation de la vitesse angulaire d'un corps rigide à partir de mesures d'accélérations est réexaminée. De cette révision, quatre algorithmes permettant l'estimation de la vitesse angulaire à partir de mesures d'accélération centripète sont proposés. Enfin, un algorithme permettant l'estimation de la vitesse angulaire d'un corps rigide à partir de mesures d'accélérations centripète et tangentielle et reposant sur le filtrage de Kalman est proposé.

L'objectif de cette thèse est l'étude du comportement cinématique des mécanismes bouclés de corps rigides. Le modèle mathématique d'un tel mécanisme est l'équation de fermeture f (q1,..., qm) = e où q1,...,qm sont des coordonnées articulaires et f est une fonction analytique à valeur dans un groupe de Lie. L'étude des propriétés cinématiques se ramène à celle de l'ensemble des configurations admissibles f-1 ( e ) qui est une sous-variété dans le cas régulier où f est une subimmersion. Par contre, l'étude est beaucoup plus difficile lorsque f possède des singularités. On utilise comme outil fondamental le formalisme de la géométrie différentielle des groupes de Lie pour le groupe des déplacements et la structure de Δ - module de son algèbre de Lie, ceci permet une écriture simple et condensée des équations de la cinématique et facilite leur traitement symbolique. Nous avons montré que l'analyse au deuxième ordre de l'équation de fermeture est suffisante pour les mécanismes 6R paradoxaux. Un algorithme d'évaluation du rang d'un ensemble de champs antisymétriques (équiprojectifs) est développé et est utilisé pour étudier les processus de génération des sous algèbres de Lie. Nous avons proposé également des méthodes de cinématique inverse pour des mécanismes spatiaux, ces méthodes permettent de résoudre l'équation de fermeture indépendamment d'un choix des coordonnées et d'obtenir des conditions nécessaires et suffisantes de résolution : notamment, la méthode simplifie considérablement la procédure de résolution pour les mécanismes 6R spatiaux.

L'objectif de cette thèse est l'étude du comportement dynamique des milieux curvilignes, en grands déplacements. Ce qui introduit une source de non linéarité géométrique qui se manifeste dans le terme d'inertie ainsi que dans le terme de rigidité. Le milieu curviligne considéré est modélisé par une suite continue de sections rigides liées par des milieux élastiques de masse nulle. On n'introduit aucune hypothèse simplificatrice dans la description des efforts intérieurs. Dans le modèle proposé, nous pouvons introduire une loi de comportement élastique non linéaire ce qui rajoute une deuxième source de non linéarité. On utilise ici comme outil fondamental le formalisme de la géométrie différentielle des groupes de Lie, ceci permet une écriture simple et condensée des équations de la dynamique et facilite leur traitement numérique. Les équations sont résolues par un algorithme numérique élaboré dans le même formalisme, ce qui évite l'utilisation "lourde" des paramètres de coordonnées. Enfin, les résultats obtenus sont appliqués à deux exemples concrets : le premier d'origine industrielle concerne le comportement du faisceau de câbles robotiques, le deuxième issu du Génie parasismique traite du comportement dynamique de grands bâtiments.

Le docking moléculaire est une méthode permettant de prédire l'orientation d'une molécule donnée relativement à une autre lorsque celles-ci forment un complexe. Le premier algorithme de docking moléculaire a vu jour en 1990 afin de trouver de nouveaux candidats face à la protéase du VIH-1. Depuis, l'utilisation de protocoles de docking est devenue une pratique standard dans le domaine de la conception de nouveaux médicaments. Typiquement, un protocole de docking comporte plusieurs phases. Il requiert l'échantillonnage exhaustif du site d'interaction où les éléments impliqués sont considérées rigides. Des algorithmes de clustering sont utilisés afin de regrouper les candidats à l'appariement similaires. Des méthodes d'affinage sont appliquées pour prendre en compte la flexibilité au sein complexe moléculaire et afin d'éliminer de possibles artefacts de docking. Enfin, des algorithmes d'évaluation sont utilisés pour sélectionner les meilleurs candidats pour le docking. Cette thèse présente de nouveaux algorithmes de protocoles de docking qui facilitent la prédiction des structures de complexes protéinaires, une des cibles les plus importantes parmi les cibles visées par les méthodes de conception de médicaments. Une première contribution concerne l‘algorithme Docktrina qui permet de prédire les conformations de trimères protéinaires triangulaires. Celui-ci prend en entrée des prédictions de contacts paire-à-paire à partir d'hypothèse de corps rigides. Ensuite toutes les combinaisons possibles de paires de monomères sont évalués à l'aide d'un test de distance RMSD efficace. Cette méthode à la fois rapide et efficace améliore l'état de l'art sur les protéines trimères. Deuxièmement, nous présentons RigidRMSD une librairie C++ qui évalue en temps constant les distances RMSD entre conformations moléculaires correspondant à des transformations rigides. Cette librairie est en pratique utile lors du clustering de positions de docking, conduisant à des temps de calcul améliorés d'un facteur dix, comparé aux temps de calcul des algorithmes standards. Une troisième contribution concerne KSENIA, une fonction d'évaluation à base de connaissance pour l'étude des interactions protéine-protéine. Le problème de la reconstruction de fonction d'évaluation est alors formulé et résolu comme un problème d'optimisation convexe. Quatrièmement, CARBON, un nouvel algorithme pour l'affinage des candidats au docking basés sur des modèles corps-rigides est proposé. Le problème d'optimisation de corps-rigides est vu comme le calcul de trajectoires quasi-statiques de corps rigides influencés par la fonction énergie. CARBON fonctionne aussi bien avec un champ de force classique qu'avec une fonction d'évaluation à base de connaissance. CARBON est aussi utile pour l'affinage de complexes moléculaires qui comportent des clashes stériques modérés à importants. Finalement, une nouvelle méthode permet d'estimer les capacités de prédiction des fonctions d'évaluation. Celle-ci permet d‘évaluer de façon rigoureuse la performance de la fonction d'évaluation concernée sur des benchmarks de complexes moléculaires. La méthode manipule la distribution des scores attribués et non pas directement les scores de conformations particulières, ce qui la rend avantageuse au regard des critères standard basés sur le score le plus élevé. Les méthodes décrites au sein de la thèse sont testées et validées sur différents benchmarks protéines-protéines. Les algorithmes implémentés ont été utilisés avec succès pour la compétition CAPRI concernant la prédiction de complexes protéine-protéine. La méthodologie développée peut facilement être adaptée pour de la reconnaissance d'autres types d'interactions moléculaires impliquant par exemple des ligands, de l'ARN… Les implémentations en C++ des différents algorithmes présentés seront mises à disposition comme SAMSON Elements de la plateforme logicielle SAMSON sur http://www.samson-connect.net ou sur http://nano-d.inrialpes.fr/software. / Molecular docking is a method that predicts orientation of one molecule with respect to another one when forming a complex. The first computational method of molecular docking was applied to find new candidates against HIV-1 protease in 1990. Since then, using of docking pipelines has become a standard practice in drug discovery. Typically, a docking protocol comprises different phases. The exhaustive sampling of the binding site upon rigid-body approximation of the docking subunits is required. Clustering algorithms are used to group similar binding candidates. Refinement methods are applied to take into account flexibility of the molecular complex and to eliminate possible docking artefacts. Finally, scoring algorithms are employed to select the best binding candidates. The current thesis presents novel algorithms of docking protocols that facilitate structure prediction of protein complexes, which belong to one of the most important target classes in the structure-based drug design. First, DockTrina - a new algorithm to predict conformations of triangular protein trimers (i.e. trimers with pair-wise contacts between all three pairs of proteins) is presented. The method takes as input pair-wise contact predictions from a rigid-body docking program. It then scans and scores all possible combinations of pairs of monomers using a very fast root mean square deviation (RMSD) test. Being fast and efficient, DockTrina outperforms state-of-the-art computational methods dedicated to predict structure of protein oligomers on the collected benchmark of protein trimers. Second, RigidRMSD - a C++ library that in constant time computes RMSDs between molecular poses corresponding to rigid-body transformations is presented. The library is practically useful for clustering docking poses, resulting in ten times speed up compared to standard RMSD-based clustering algorithms. Third, KSENIA - a novel knowledge-based scoring function for protein-protein interactions is developed. The problem of scoring function reconstruction is formulated and solved as a convex optimization problem. As a result, KSENIA is a smooth function and, thus, is suitable for the gradient-base refinement of molecular structures. Remarkably, it is shown that native interfaces of protein complexes provide sufficient information to reconstruct a well-discriminative scoring function. Fourth, CARBON - a new algorithm for the rigid-body refinement of docking candidates is proposed. The rigid-body optimization problem is viewed as the calculation of quasi-static trajectories of rigid bodies influenced by the energy function. To circumvent the typical problem of incorrect stepsizes for rotation and translation movements of molecular complexes, the concept of controlled advancement is introduced. CARBON works well both in combination with a classical force-field and a knowledge-based scoring function. CARBON is also suitable for refinement of molecular complexes with moderate and large steric clashes between its subunits. Finally, a novel method to evaluate prediction capability of scoring functions is introduced. It allows to rigorously assess the performance of the scoring function of interest on benchmarks of molecular complexes. The method manipulates with the score distributions rather than with scores of particular conformations, which makes it advantageous compared to the standard hit-rate criteria. The methods described in the thesis are tested and validated on various protein-protein benchmarks. The implemented algorithms are successfully used in the CAPRI contest for structure prediction of protein-protein complexes. The developed methodology can be easily adapted to the recognition of other types of molecular interactions, involving ligands, polysaccharides, RNAs, etc. The C++ versions of the presented algorithms will be made available as SAMSON Elements for the SAMSON software platform at http://www.samson-connect.net or at http://nano-d.inrialpes.fr/software.

Le docking moléculaire est une méthode permettant de prédire l'orientation d'une molécule donnée relativement à une autre lorsque celles-ci forment un complexe. Le premier algorithme de docking moléculaire a vu jour en 1990 afin de trouver de nouveaux candidats face à la protéase du VIH-1. Depuis, l'utilisation de protocoles de docking est devenue une pratique standard dans le domaine de la conception de nouveaux médicaments. Typiquement, un protocole de docking comporte plusieurs phases. Il requiert l'échantillonnage exhaustif du site d'interaction où les éléments impliqués sont considérées rigides. Des algorithmes de clustering sont utilisés afin de regrouper les candidats à l'appariement similaires. Des méthodes d'affinage sont appliquées pour prendre en compte la flexibilité au sein complexe moléculaire et afin d'éliminer de possibles artefacts de docking. Enfin, des algorithmes d'évaluation sont utilisés pour sélectionner les meilleurs candidats pour le docking. Cette thèse présente de nouveaux algorithmes de protocoles de docking qui facilitent la prédiction des structures de complexes protéinaires, une des cibles les plus importantes parmi les cibles visées par les méthodes de conception de médicaments. Une première contribution concerne l‘algorithme Docktrina qui permet de prédire les conformations de trimères protéinaires triangulaires. Celui-ci prend en entrée des prédictions de contacts paire-à-paire à partir d'hypothèse de corps rigides. Ensuite toutes les combinaisons possibles de paires de monomères sont évalués à l'aide d'un test de distance RMSD efficace. Cette méthode à la fois rapide et efficace améliore l'état de l'art sur les protéines trimères. Deuxièmement, nous présentons RigidRMSD une librairie C++ qui évalue en temps constant les distances RMSD entre conformations moléculaires correspondant à des transformations rigides. Cette librairie est en pratique utile lors du clustering de positions de docking, conduisant à des temps de calcul améliorés d'un facteur dix, comparé aux temps de calcul des algorithmes standards. Une troisième contribution concerne KSENIA, une fonction d'évaluation à base de connaissance pour l'étude des interactions protéine-protéine. Le problème de la reconstruction de fonction d'évaluation est alors formulé et résolu comme un problème d'optimisation convexe. Quatrièmement, CARBON, un nouvel algorithme pour l'affinage des candidats au docking basés sur des modèles corps-rigides est proposé. Le problème d'optimisation de corps-rigides est vu comme le calcul de trajectoires quasi-statiques de corps rigides influencés par la fonction énergie. CARBON fonctionne aussi bien avec un champ de force classique qu'avec une fonction d'évaluation à base de connaissance. CARBON est aussi utile pour l'affinage de complexes moléculaires qui comportent des clashes stériques modérés à importants. Finalement, une nouvelle méthode permet d'estimer les capacités de prédiction des fonctions d'évaluation. Celle-ci permet d‘évaluer de façon rigoureuse la performance de la fonction d'évaluation concernée sur des benchmarks de complexes moléculaires. La méthode manipule la distribution des scores attribués et non pas directement les scores de conformations particulières, ce qui la rend avantageuse au regard des critères standard basés sur le score le plus élevé. Les méthodes décrites au sein de la thèse sont testées et validées sur différents benchmarks protéines-protéines. Les algorithmes implémentés ont été utilisés avec succès pour la compétition CAPRI concernant la prédiction de complexes protéine-protéine. La méthodologie développée peut facilement être adaptée pour de la reconnaissance d'autres types d'interactions moléculaires impliquant par exemple des ligands, de l'ARN… Les implémentations en C++ des différents algorithmes présentés seront mises à disposition comme SAMSON Elements de la plateforme logicielle SAMSON sur http://www.samson-connect.net ou sur http://nano-d.inrialpes.fr/software. / Molecular docking is a method that predicts orientation of one molecule with respect to another one when forming a complex. The first computational method of molecular docking was applied to find new candidates against HIV-1 protease in 1990. Since then, using of docking pipelines has become a standard practice in drug discovery. Typically, a docking protocol comprises different phases. The exhaustive sampling of the binding site upon rigid-body approximation of the docking subunits is required. Clustering algorithms are used to group similar binding candidates. Refinement methods are applied to take into account flexibility of the molecular complex and to eliminate possible docking artefacts. Finally, scoring algorithms are employed to select the best binding candidates. The current thesis presents novel algorithms of docking protocols that facilitate structure prediction of protein complexes, which belong to one of the most important target classes in the structure-based drug design. First, DockTrina - a new algorithm to predict conformations of triangular protein trimers (i.e. trimers with pair-wise contacts between all three pairs of proteins) is presented. The method takes as input pair-wise contact predictions from a rigid-body docking program. It then scans and scores all possible combinations of pairs of monomers using a very fast root mean square deviation (RMSD) test. Being fast and efficient, DockTrina outperforms state-of-the-art computational methods dedicated to predict structure of protein oligomers on the collected benchmark of protein trimers. Second, RigidRMSD - a C++ library that in constant time computes RMSDs between molecular poses corresponding to rigid-body transformations is presented. The library is practically useful for clustering docking poses, resulting in ten times speed up compared to standard RMSD-based clustering algorithms. Third, KSENIA - a novel knowledge-based scoring function for protein-protein interactions is developed. The problem of scoring function reconstruction is formulated and solved as a convex optimization problem. As a result, KSENIA is a smooth function and, thus, is suitable for the gradient-base refinement of molecular structures. Remarkably, it is shown that native interfaces of protein complexes provide sufficient information to reconstruct a well-discriminative scoring function. Fourth, CARBON - a new algorithm for the rigid-body refinement of docking candidates is proposed. The rigid-body optimization problem is viewed as the calculation of quasi-static trajectories of rigid bodies influenced by the energy function. To circumvent the typical problem of incorrect stepsizes for rotation and translation movements of molecular complexes, the concept of controlled advancement is introduced. CARBON works well both in combination with a classical force-field and a knowledge-based scoring function. CARBON is also suitable for refinement of molecular complexes with moderate and large steric clashes between its subunits. Finally, a novel method to evaluate prediction capability of scoring functions is introduced. It allows to rigorously assess the performance of the scoring function of interest on benchmarks of molecular complexes. The method manipulates with the score distributions rather than with scores of particular conformations, which makes it advantageous compared to the standard hit-rate criteria. The methods described in the thesis are tested and validated on various protein-protein benchmarks. The implemented algorithms are successfully used in the CAPRI contest for structure prediction of protein-protein complexes. The developed methodology can be easily adapted to the recognition of other types of molecular interactions, involving ligands, polysaccharides, RNAs, etc. The C++ versions of the presented algorithms will be made available as SAMSON Elements for the SAMSON software platform at http://www.samson-connect.net or at http://nano-d.inrialpes.fr/software.

Le sujet de la thèse est l'étude et la construction de méthodes numériques géométriques pour les équations différentielles, qui préservent des propriétés géométriques du flot exact, notamment la symétrie, la symplecticité des systèmes hamiltoniens, la conservation d'intégrales premières, la structure de Poisson, etc.<br />Dans la première partie, on introduit une nouvelle approche de construction d'intégrateurs numériques géométriques d'ordre élevé en s'inspirant de la théorie des équations différentielles modifiées. Le cas des méthodes développables en B-séries est spécifiquement analysé et on introduit une nouvelle loi de composition sur les B-séries. L'efficacité de cette approche est illustrée par la construction d'un nouvel intégrateur géométrique d'ordre élevé pour les équations du mouvement d'un corps rigide. On obtient également une méthode numérique précise pour le calcul de points conjugués pour les géodésiques du corps rigide.<br />Dans la seconde partie, on étudie dans quelle mesure les excellentes performances des méthodes symplectiques, pour l'intégration à long terme en astronomie et en dynamique moléculaire, persistent pour les problèmes de contrôle optimal. On discute également l'extension de la théorie des équations modifiées aux problèmes de contrôle optimal.<br />Dans le même esprit que les équations modifiées, on considère dans la dernière partie des méthodes de pas fractionnaire (splitting) pour les systèmes hamiltoniens perturbés, utilisant des potentiels modifiés. On termine par la construction de méthodes de splitting d'ordre élevé avec temps complexes pour les équations aux dérivées partielles paraboliques, notamment les problèmes de réaction-diffusion en chimie.