Spin-orbit coupling (SOC) reduces the spatial symmetry of ferromagnetic solids. That is, the physical properties of ferromagnetic materials are anisotropic, depending on the magnetization direction. In this thesis, by means of numerical calculations with full-relativistic density functional theory, we studied two kinds of physical properties: surface magnetic anisotropy energy (MAE) and anisotropic thermoelectric power due to Lifshitz transitions. After a short introduction to the full-relativistic density functional theory in Chapter 2, the MAE of ferromagnetic thin films is studied in Chapter 3. For such systems, separation of different contributions, such as bulk magnetocrystalline anisotropy (MCA) energy, shape anisotropy energy, and surface/interface anisotropy energy, is crucial to gain better understanding of experiments. By fitting our calculating results for thick slabs to a phenomenological model, reliable surface MAE could be obtained. Following this idea, we have studied the MAE of Co slabs with different geometries, focusing on the effects of orbital polarization correction (OPC). We found that the surface anisotropy is mainly determined by the geometry. While OPC gives better results of orbital moments, it overestimates the MAE. In the second part of Chapter3, the effects of electric fields on the MAE of L10 ferromagnetic thin films are studied. Using a simple model to simulate the electric field, our calculations are in good agreement with previous experimental results. We predicted that for CoPt, even larger effects exist. Moreover, we found that it is the amount of screening charge that determines the magnetoelectric coupling effects. This gives us some clue about how to achieve electric field control of magnetization direction. In Chapter 4, Lifshitz transitions in L10 FePt caused by a canted magnetic field are studied. We found several Lifshitz transitions in ordered FePt with tiny features in DOS. Using a two-band model, it is demonstrated that at such transitions, the singular behaviour of kinetic properties is due to the interband scattering, and the singularity itself is proportional to the derivative of the singular DOS. For FePt, such singularity will be smeared into anomaly by chemical disorder. Using CPA, we studied the effects of energy level broadening for the critical bands in FePt. We found that for experimentally available FePt thin films, Lifshitz transitions would induce up to a 3% increase of thermopower as the magnetization is rotated from the easy axis to the hard axis. / Spin-Bahn-Kopplung reduziert die Symmetrie ferromagnetischer Festkörper. Das bedeutet, dass die physikalischen Eigenschaften ferromagnetischer Stoffe anisotrop bezüglich der Magnetisierungsrichtung sind. In dieser Dissertation werden mittels numerischer voll-relativistischer Dichtefunktional-Rechnungen zwei Arten physikalischer Eigenschaften untersucht: magnetische Oberflächen-Anisotropieenergie (MAE) und anisotrope Thermokraft durch Lifshitz-Übergänge. Nach einer kurzen Einführung in die relativistische Dichtefunktional-Theorie in Kapitel 2 wird in Kapitel 3 die MAE ferromagnetischer dünner Filme untersucht. In diesen Systemen ist es für ein Verständnis experimenteller Ergebnisse wichtig, verschiedene Beiträge zu separieren: Volumenanteil der magnetokristallinen Anisotropie (MCA), Formanistropie und Oberflächen bzw. Grenzflächenanisotropie. Durch Anpassen berechneter Daten für dicke Schichten an ein phänomenologisches Modell konnten verlässliche Oberflächen Anisotropien erhalten werden. In dieser Weise wurde die MAE von Co- Schichten mit unterschiedlichen Geometrien untersucht, wobei der Einfluss von Orbitalpolarisations-Korrekturen (OPC) im Vordergrund stand. Es wurde gefunden, dass die Oberflächenanisotropie hauptsächlich von der Geometrie bestimmt wird. Während OPC bessere Ergebnisse für die Orbitalmomente liefert, wird die MAE überschätzt. Im zweiten Teil von Kapitel 3 wird der Einfluss elektrischer Felder auf die MAE von dünnen ferromagnetischen Filmen mit L10-Struktur untersucht. Unter Verwendung eines einfachen Modells zur Simulation des elektrischen Feldes liefern die Rechnungen gute Übereinstimmung mit vorliegenden experimentellen Ergebnissen. Es wird vorhergesagt, dass für CoPt ein noch größerer Effekt existiert. Weiterhin wurde gefunden, dass die magnetoelektrische Kopplung von der Größe der Abschirmladung bestimmt wird. Dies ist eine wichtige Einsicht, um die Magnetisierungsrichtung durch ein elektrisches Feld kontrollieren zu können. In Kapitel 4 werden Lifshitz-Übergänge untersucht, die ein gekantetes Magnetfeld hervorruft. Es wurden mehrere Lifshitz-Übergänge in geordnetem FePt gefunden, welche kleine Anomalien in der Zustandsdichte hervorrufen. Mit Hilfe eines Zweiband-Modells wird gezeigt, dass an solchen Übergängen das singuläre Verhalten kinetischer Eigenschaften durch Interband- Streuung verursacht wird und dass die Singularität proportional zur Ableitung der singulären Zustandsdichte ist. In FePt wird durch chemische Unordnung diese Singularität zu einer Anomalie verschmiert. Der Einfluss einer Verbreiterung der Energieniveaus der kritischen Bänder in FePt wurde mittels CPA untersucht. Es wurde gefunden, dass in experimentell verfügbaren dünnen FePt-Filmen Lifshitz-Übergänge bis zu 3% Erhöhung der Thermokraft erzeugen, wenn die Magnetisierung von der leichten in die harte Richtung gedreht wird.

Im Rahmen dieser Doktorarbeit werden de Haas-van Alphen-Untersuchungen an den nichtmagnetischen Borkarbidsupraleitern LuNi2B2C und YNi2B2C präsentiert. Aus den Quantenoszillationen in der normalleitenden Phase in Kombination mit Bandstrukturrechnungen konnten Informationen über die verzweigte Fermiflächenarchitektur und über die Elektron-Phonon-Kopplung der Borkarbide gewonnen werden. Die Kopplung ist stark anisotrop und fermiflächenabhängig. Dies spricht für einen Mehrbandmechanismus der Supraleitung in der Materialklasse. Zusätzlich konnten de Haas-van-Alphen-Oszillationen mehrerer Fermiflächen unterhalb von Bc2 tief in der Shubnikov-Phase beobachtet werden. Das Verhalten dieser Oszillationen lässt sich nicht mit bisher bekannten Theorien beschreiben. Allerdings weist das Bestehen der Oszillationen weit unterhalb von Bc2 auf ein Bestehen von elektronischen Zuständen in der Shubnikov-Phase hin.

Carbon-based conductors like carbon nanotubes (CNTs) and graphene nanoribbons (GNRs) have many properties, which make them relevant for potential electronic applications. Among them are high conductances and tunable band gap sizes. These properties make CNTs and GNRs useful in many circumstances, e.g. as channel material in transistors or transparent electrodes in solar cells. Plenty of literature can be found on the topic of single linear CNTs/GNRs. Some applications however require a large network of these conductors. In addition, a single conductor has only a small impact on the network conductance, which reduces the need to control the properties of each individual nanotube/-ribbon. This leads to networks being easier to apply. In this work, the conductance of large networks of GNRs is calculated using the quantum-transport formalism (QT). This has not been done before in literature. In order to apply QT to such a large amount of atoms, the recursive Green's function formalism is used. For this the networks are devided into subcells, which are represented by tight-binding matrices. Similar networks are also examined using two different nodal analysis (NA) approaches, where the nanoribbons are treated as ohmic conductors. For NA with one-dimensional conductors, major discrepancies are found in regards to the QT model. However, networks consisting of two-dimensional conductors (NA-2D) have many properties similar to the QT networks. A recipe to approximate the QT results with NA-2D is presented.:1. Introduction 2. Theoretical Principles 2.1 Carbon-based Conductors 2.1.1 Structure and Properties 2.1.2 Networks 2.2 Tight-Binding Model 2.3 Quantum Transport 2.3.1 Introduction 2.3.2 Level Broadening 2.3.3 Current Flow 2.3.4 Transmission 2.4 Nodal Analysis 3. Implementation 3.1 Quantum Tranport 3.1.1 Network Generation 3.1.2 Density-Functional based Tight-Binding Method 3.1.3 Recursive Green's Function Algorithm 3.1.4 Conductance 3.2 Nodal Analysis 3.2.1 One-dimensional Conductors 3.2.2 Two-dimensional Conductors 4. Results 4.1 Quantum Transport 4.1.1 Band Structures and Fermi Energies 4.1.2 Ideal Transmission and Consistency Tests 4.1.3 Percolation 4.1.4 Transmission 4.1.5 Conductance 4.1.6 Power Law Scaling 4.1.7 Size Dependence and Confinement Effects 4.1.8 Calculation Time 4.2 Nodal Analysis 4.2.1 One-dimensional Conductors 4.2.2 Two-dimensional Conductors 4.2.3 Calculation Time 4.3 Approximating QT with NA 4.3.1 Optimal Parameters 4.3.2 Percolation 4.3.3 Conductance 4.3.4 Power Law Scaling 5. Conclusions / Graphenbasierte Leiter wie Kohlenstoff-Nanoröhrchen (engl. 'carbon nanotubes', CNTs) oder Graphen-Nanobänder (engl. 'graphene nanoribbons', GNRs) haben viele Eigenschaften, die sie für potenzielle elektronische Anwendungen interessant machen. Darunter sind hohe Leitfähigkeiten und einstellbare Bandlückengrößen. Dadurch sind CNTs und GNRs in vielen Bereichen nützlich, z.B. als Kanalmaterial in Transistoren oder als transparente Elektroden in Solarzellen. Es gibt viel Literatur über einzelne, lineare CNTs/GNRs. Einige Anwendungen benötigen jedoch ein großes Netzwerk dieser Leiter. Zusätzlich hat ein einzelner Leiter wenig Einfluss auf die Leitfähigkeit des Netzwerks, wodurch die Eigenschaften der einzelnen Nanoröhrchen/-streifen weniger streng kontrolliert werden müssen. Dies führt dazu, dass es einfacher ist Netzwerke zu nutzen. In dieser Arbeit wird die Leitfähigkeit von großen GNR-Netzwerken mittels Quantentransport (QT) berechnet. Dies wurde in der Literatur noch nicht getan. Um QT auf eine so große Menge an Atomen anzuwenden wird der rekursive Greenfunktions-Formalismus benutzt. Dazu werden die Netzwerke in Unterzellen unterteilt, die durch Tight-Binding-Matrizen dargestellt werden. Ähnliche Netzwerke werden auch mit zwei Versionen der Knotenanalyse (engl. 'nodal analysis', NA) untersucht, welche die Nanobänder wie ohmische Leiter behandelt. Die Ergebnisse der NA mit eindimensionalen Leitern weisen deutliche Unterschiede zu den mit QT erzielten Ergebnissen auf. Wenn jedoch zweidimensionale Leiter in NA verwendet werden (NA-2D) gibt es viele parallelen zu den QT Ergebnissen. Zuletzt wird ein Vorgehen präsentiert, mit dem QT Resultate durch NA-2D Rechnungen genähert werden können.:1. Introduction 2. Theoretical Principles 2.1 Carbon-based Conductors 2.1.1 Structure and Properties 2.1.2 Networks 2.2 Tight-Binding Model 2.3 Quantum Transport 2.3.1 Introduction 2.3.2 Level Broadening 2.3.3 Current Flow 2.3.4 Transmission 2.4 Nodal Analysis 3. Implementation 3.1 Quantum Tranport 3.1.1 Network Generation 3.1.2 Density-Functional based Tight-Binding Method 3.1.3 Recursive Green's Function Algorithm 3.1.4 Conductance 3.2 Nodal Analysis 3.2.1 One-dimensional Conductors 3.2.2 Two-dimensional Conductors 4. Results 4.1 Quantum Transport 4.1.1 Band Structures and Fermi Energies 4.1.2 Ideal Transmission and Consistency Tests 4.1.3 Percolation 4.1.4 Transmission 4.1.5 Conductance 4.1.6 Power Law Scaling 4.1.7 Size Dependence and Confinement Effects 4.1.8 Calculation Time 4.2 Nodal Analysis 4.2.1 One-dimensional Conductors 4.2.2 Two-dimensional Conductors 4.2.3 Calculation Time 4.3 Approximating QT with NA 4.3.1 Optimal Parameters 4.3.2 Percolation 4.3.3 Conductance 4.3.4 Power Law Scaling 5. Conclusions

Die Dissertation liefert einen Beitrag zur Modellierung und optimalen Ansteuerung von vernetzten hybriden Energiesystemen. Die Arbeit beschreibt die Entwicklung einer modellprädiktiven Regelung (MPC) für konkrete Energiesysteme. Dafür wird eine Betrachtung zu berücksichtigender wirtschaftlicher und technischer Rahmenbedingungen vorgenommen, die zur Formulierung notwendiger Nebenbedingungen für die MPC genutzt wird. Für den Umgang mit dem ansteigenden Rechenbedarf der MPC bei steigender Systemzahl wird ein alternativer Ansatz auf Basis eines auktionsbasierten Algorithmus vorgestellt. Die Modellierung der Energiesysteme wird ausgehend von einer bestehenden Laboranlage vorgenommen. Die Erprobung der vorgestellten Ansätze erfolgt in einer Simulationsumgebung, die die Untersuchung verschiedener Szenarien erlaubt. Im Rahmen der Simulationsszenarien mit unterschiedlicher Systemzahl und Zusammensetzung der Energie-systeme wird eine Sensibilitätsanalyse der vorgestellten MPC vorgenommen. Die Interpretation der Ergebnisse erfolgt auf Basis numerischer und empirischer Bewertungskriterien.

In der vorliegenden Arbeit wird die Entwicklung eines kinetischen Energy Harvesters vorgestellt, der auf Grundlage des elektrostatischen Wandlerprinzips aus Vibrationen elektrische Energie generiert. Für die Umsetzung wurde eine Siliziummikrostruktur entworfen, die für Arbeitsfrequenzen unter 100 Hz ausgelegt ist. Die Zahnstruktur der verwendeten Elektroden ermöglicht Spaltabstände im Submikrometerbereich und folglich große Kapazitätsänderungen, die durch die Elektrodengeometrie zusätzlich mit einer höheren Frequenz als die mechanische Bewegung stattfinden. Vergleichsweise große Leistungsausbeuten und geringe Quellimpedanzen sind dadurch erreichbar. Die geometrischen Parameter der Elektroden wurden unter Berücksichtigung der auftretenden Fertigungstoleranzen und Wechselwirkungen zueinander optimiert. Für die Ausnutzung einer ausreichend großen Inertialmasse wurde ein feinwerktechnisch hergestellter Hebelmechanismus an die Mikrostruktur angekoppelt. Über diesen wird zusätzlich ein neuer Ansatz zur Abstimmung der Eigenfrequenz des Harvesters umgesetzt. Experimentelle Untersuchungen zeigten Ausgangsleistungen im einstelligen Mikrowattbereich bei Anregungen im Zehntel m/s²-Bereich. Durch fortschreitende Optimierungen der Fertigungstechnologie sind noch deutliche Leistungssteigerungen um mindestens zwei Größenordnungen möglich. Weiterhin wird ein Energiemanagementsystem vorgestellt, welches die effiziente Übertragung der Energie auf den Verbraucher ermöglicht. / In this work the development of a kinetic energy harvester using the electrostatic conversion principle is presented. The silicon microstructure is designed to work in frequency ranges below 100Hz. Its toothed electrode structure enables gap distances in the sub micrometer range and consequently high changes of capacitance. Additionally, due to the electrode geometry the frequency of the capacitance changes is higher then the frequency of the mechanical movement. Thus high power outputs and low source impedances can be reached. The electrodes geometric parameters were optimized considering manufacturing tolerances and interactions of the parameters. To reach a sufficient inertial mass, a lever mechanism manufactured by precision engineering was connected to the microstructure. This mechanism also allows the implementation of a new method of frequency tuning. In experimental tests power outputs in the single digit microwatt range under excitations of 0.3 m/s² were reached. In accordance of further optimizations of the manufacturing technology significantly higher outputs, by at least two orders of magnitude, are possible,. Furthermore an energy management system is presented, that allows the efficient transfer of the electrical energy to the consumer.

In einer Einleitung werden verschiedene Möglichkeiten der Kopplung unterschiedlicher physikalischer Domänen diskutiert. Begriffe wie Kopplung auf Systemebene über Terminals und auf Feldebene über Lastvektor oder Matrixkopplung werden zugeordnet. Wie diese Kopplungsmöglichkeiten sich in der ANSYS Simulationsumgebung wiederfinden, wird aufgezeigt. Am Beispiel eines Akuators wird erläutert, welche physikalischen Domänen gekoppelt betrachtet werden müssen, um die hier vorhandenen temperaturabhängigen Materialeigenschaften zu berücksichtigen. In einem Beispiel zu einer Steckverbindung wird aufgezeigt, wie eine vom Kontaktdruck abhängige Leitfähigkeit berücksichtigt wird und Ergebnisse aus dieser Simulation werden diskutiert.

In this thesis we investigate the formation of patterns in a simple activator-inhibitor model supplemented with an inhibitory nonlocal coupling term. This model exhibits a wave instability for slow inhibitor diffusion, while, for fast inhibitor diffusion, a Turing instability is found. For moderate values of the inhibitor diffusion these two instabilities occur simultaneously at a codimension-2 wave-Turing instability. We perform a weakly nonlinear analysis of the model in the neighbourhood of this codimension-2 instability. The resulting amplitude equations consist in a set of coupled Ginzburg-Landau equations. These equations predict that the model exhibits bistability between travelling waves and Turing patterns. We present a study of interfaces separating wave and Turing patterns arising from the codimension-2 instability. We study theoretically and numerically the dynamics of such interfaces in the framework of the amplitude equations and compare these results with numerical simulations of the model near and far away from the codimension-2 instability. Near the instability, the dynamics of interfaces separating small amplitude Turing patterns and travelling waves is well described by the amplitude equations, while, far from the codimension-2 instability, we observe a locking of the interface velocities. This locking mechanism is imposed by the absence of defects near the interfaces and is responsible for the formation of drifting pattern domains, i.e. moving localised patches of travelling waves embedded in a Turing pattern background and vice versa.

In the first two chapters we presented the basics of density functional theory and semiempirical LSD+U approximation, which was implemented in the full-potential local-orbital (FPLO) minimal-basis calculation scheme. In the third chapter we tested the implemented version of LSDA+U on 3d transitional metal monoxides. Essential improvement of the spectroscopic properties was obtained. A simple model describing the value and direction of the magnetic moment of a transition metal ion was presented. The model visualizes the interplay of the spin-orbit coupling and crystal field splitting. In the fourth chapter we calculated the electronic spectrum of the single Zn impurity in CuO2 plane considered as a vacancy in Cu 3d states. The analytic solution for the states of different symmetry was obtained. Depending on the strength of perturbation induced by the impurity on the neighboring Cu ions, the states are either resonant or localized. The critical values of the perturbation were computed. In the fifth chapter we presented the calculations for three novel vanadates: MgVO3, Sb2O2VO3 and VOMoO4. The tight-binding parameters and the exchange integrals were computed. The magnesium and antimony vanadates appeared to be spin-½ one-dimensional systems, the latter having much stronger one-dimensional character and being probably the best realization of inorganic spin-Peierls system. The molybdenum vanadate was found to be two-dimensional spin-½ system. The Mo 4d orbitals play an important role in the electronic transfer.

Die Arbeit befasst sich mit der Elektronenspinresonanz (ESR) stark korrelierter Elektronensysteme mit ferromagnetischen Wechselwirkungen. Es wurden dafür Messungen an den Kondogitter-Systemen CeRuPO und CeOsPO, der Dotierungsreihe CeFeAs1-xPxO, den niederdimensionalen frustrierten Quadratgittern AA’VO(PO4)2 sowie in dem schwach ferromagnetischen Metall Nb1-yFe2+y durchgeführt. Alle Verbindungen zeigen entweder eine ferromagnetische Ordnung oder befinden sich in der Nähe einer ferromagnetischen Instabilität, die die Eigenschaften des stark korrelierten Systems beeinflusst.:Abkürzungsverzeichnis xi Einleitung 1 1 Theoretische Modelle 5 1.1 Auswirkungen des Kristallfeldes 1.2 Phasenübergänge, Landau-Ginzburg Theorie und Quantenphasenübergänge 1.3 Physik Niederdimensionaler Spin-Systeme 1.3.1 Das zweidimensionale XY-Modell 1.3.2 Das J1-J2-Modell auf dem Quadratgitter 1.4 4f-Magnetismus, Kondoeffekt und Kondogitter 1.4.1 Vom freien Elektronengas zur Landau-Fermi-Flüssigkeit 1.4.2 Kondo-Effekt, RKKY-Wechselwirkung und Schwere Fermionen Systeme 1.4.3 Cer: 4f-Magnetismus 1.5 Schwach magnetische Metalle 1.5.1 Pauli-Suszeptibilität und spontan aufgespaltene Bänder 1.5.2 Spinfluktuationen in itineranten Magneten und marginale Fermiflüssigkeit 1.6 Dichtefunktionaltheorie 2 Grundlagen der Elektronenspinresonanz 31 2.1 Prinzip der magnetischen Resonanz und ESR Parameter 2.1.1 Linienform 2.1.2 Intensität 2.1.3 g-Faktor 2.1.4 Linienbreite und Spin-Relaxation 2.2 Experimenteller Aufbau 2.3 ESR in niederdimensionalen Systemen 2.4 ESR an Metallen und Kondogitter Systemen 2.4.1 Auswirkungen der Wechselwirkungen zwischen lokalen Momenten und Leitungselektronen auf die ESR 2.4.2 ESR in Kondogittersystemen 2.5 Leitungselektronenspinresonanz 2.5.1 Linienform in der LESR 2.5.2 Elliot-Yafet Theorie und ihre Verallgemeinerung 2.6 ESR an magnetischen Phasenübergängen 3 ESR in CeTPO (T=Ru, Os) und CeFeAs1-xPxO 59 3.1 Die Eigenschaften von CeTPO (T=Ru, Os) & CeFeAs1-xPxO 3.1.1 Das ferromagnetische Kondogitter CeRuPO 3.1.2 Das antiferromagnetische Kondogitter CeOsPO 3.1.3 Die Dotierungsreihe CeFeAs1-xPxO 3.2 CeRuPO und CeOsPO: Die Bedeutung FM Fluktuationen für die Beobachtbarkeit der ESR in Kondo-Gitter Systemen 3.3 ESR an CeRuPO Einkristallen 3.3.1 Anisotropie und Temperaturverhalten des g-Faktors 3.3.2 Anisotropie und Temperaturverhalten der Linienbreite 3.4 Entwicklung der ESR in der Dotierungsreihe CeFeAs1-xPxO 3.4.1 CeFeAs0.7P0.3O - Supraleitung und Ferromagnetismus 3.5 Zusammenfassung ESR in CeTPO und CeFeAs1-xPxO 4 ESR in den frustrierten Quadratgittern AA’VO(PO4)2 97 4.1 Struktur und magnetische Eigenschaften AA’VO(PO4)2 4.2 ESR an polykristallinen SrZnVO(PO4)2 & BaCdVO(PO4)2 4.2.1 SrZnVO(PO4)2 4.2.2 BaCdVO(PO4)2 4.3 ESR an Pb2VO(PO4)2-Einkristallen 4.3.1 Anisotropie der Resonanz 4.3.2 Temperaturabhängigkeit 4.4 Zusammenfassung der ESR in den frustrierten Quadratgittern 7 4.4.1 Analyse der Linienbreitendivergenz am Phasenübergang 5 Leitungselektronenspinresonanz in Nb1-yFe2+y 125 5.1 Das magnetische Phasendigramm von Nb1-yFe2+y 5.2 LESR in polykristallinen Nb1-yFe2+y 5.2.1 Proben mit ferromagnetischer Ordnung 5.2.2 Proben mit SDW-Ordnung 5.2.3 Zusammenfassung polykristalline Proben 5.3 Bandstruktur und verallgemeinerte Elliot-Yafet-Theorie in NbFe2 5.4 LESR in einkristallinen Proben von Nb1-yFe2+y 5.4.1 Fe-reicher Einkristall: y=0.016 5.4.2 Stöchiometrischer Einkristall: NbFe2 5.4.3 Nb-reicher Einkristall: y=-0.008 5.5 Zusammenfassung LESR in Nb1-yFe2+y 6 Zusammenfassung und Ausblick Anhang A.1 Das zweite Moment der Dipol-Dipol-Wechselwirkung A.2 Untersuchte CeFeAs1-xPxO-Proben Literaturverzeichnis

‘Coupled process’ implies that one process affects the initiation and progress of the others and vice versa. The deformation and damage behaviors of rock under loading process change the fluid flow field within it, and lead to altering in permeable characteristics; on the other side inner fluid flow leads to altering in pore pressure and effective stress of rock matrix and flow by influencing stress strain behavior of rock. Therefore, responses of rock to natural or man-made perturbations cannot be predicted with confidence by considering each process independently. As far as hydro-mechanical behavior of rock is concerned, the researchers have always been making efforts to develop the model which can represent the permeable characteristics as well as stress-strain behaviors during the entire damage process. A brittle low porous granite was chosen as the study object in this thesis, the aim is to establish a corresponding constitutive law including the relation between permeability evolution and mechanical deformation as well as the rock failure behavior under hydro-mechanical coupled conditions based on own hydro-mechanical coupled lab tests. The main research works of this thesis are as follows: 1. The fluid flow and mechanical theoretical models have been reviewed and the theoretical methods to solve hydro-mechanical coupled problems of porous medium such as flow equations, elasto-plastic constitutive law, and Biot coupled control equations have been summarized. 2. A series of laboratory tests have been conducted on the granite from Erzgebirge–Vogtland region within the Saxothuringian segment of Central Europe, including: permeability measurements, ultrasonic wave speed measurements, Brazilian tests, uniaxial and triaxial compression tests. A hydro-mechanical coupled testing system has been designed and used to conduct drained, undrained triaxial compression tests and permeability evolution measurements during complete loading process. A set of physical and mechanical parameters were obtained. 3. Based on analyzing the complete stress-strain curves obtained from triaxial compression tests and Hoek-Brown failure criterion, a modified elemental elasto-plastic constitutive law was developed which can represent strength degradation and volume dilation considering the influence of confining pressure. 4. The mechanism of HM-coupled behavior according to the Biot theory of elastic porous medium is summarized. A trilinear evolution rule for Biot’s coefficient based on the laboratory observations was deduced to eliminate the error in predicting rock strength caused by constant Biot’s coefficient. 5. The permeability evolution of low porous rock during the failure process was described based on literature data and own measurements, a general rule for the permeability evolution was developed for the laboratory scale, a strong linear relation between permeability and volumetrical strain was observed and a linear function was extracted to predict permeability evolution during loading process based on own measurements. 6. By combining modified constitutive law, the trilinear Biot’s coefficient evolution model and the linear relationship between permeability and volumetrical strain, a fully hydro-mechanical coupled numerical simulation scheme was developed and implemented in FLAC3D. A series of numerical simulations of triaxial compression test considering the hydro-mechanical coupling were performed with FLAC3D. And a good agreement was found between the numerical simulation results and the laboratory measurements under 20 MPa confining pressure and 10 MPa fluid pressure, the feasibility of this fully hydro-mechanical coupled model was proven.