Der Gegenstand der vorliegenden Arbeit ist die Simulation von Betonstrukturen beliebiger Geometrie unter überwiegender Zugbelastung. Die Modellierung erfolgt auf Makroebene als Kontinuum und zur Lösung des mechanischen Feldproblems wird die Finite-Elemente-Methode verwendet. Ein neues Materialmodell für Beton und eine Erweiterung der Bruchenergetischen Regularisierung werden vorgestellt. Die Arbeit ist in zwei Teile gegliedert. Im ersten Teil wird ein lokales, anisotropes Schädigungsmodell abgeleitet, wobei als Schädigungsvariable ein symmetrischer Tensor zweiter Stufe gewählt wird. Die Verwendung einer Normalenregel im Raum der dissipativen Kräfte zur Bestimmung der Schädigungsevolution und die Definition der Schädigungsgrenzflächen im Raum der dissipativen Kräfte gewährleisten die Gültigkeit der Hauptsätze der Thermodynamik und des Prinzips der maximalen Dissipationsrate. Vorteilhaft ist die Symmetrie der Materialtangente, die sich aus diesem Vorgehen ergibt. Eine Formulierung mit drei entkoppelten Schädigungsgrenzflächen wird vorgeschlagen. Eine wichtige Forderung bei der Ableitung des Materialmodells war die Verwendung einer möglichst geringen Anzahl von Materialparametern, welche darüber hinaus aus wenigen Standardversuchen bestimmbar sein sollten. Das Schädigungsmodell enthält als Materialparameter den Elastizitätsmodul, die Querdehnzahl, die Zugfestigkeit und die auf eine Einheitsfläche bezogene Bruchenergie. Im zweiten Teil der Arbeit stehen Lokalisierung und Regularisierung im Fokus der Betrachtungen. Aufgrund der lokalen Formulierung des Materialmodells tritt bei Finite-Elemente Simulationen eine Netzabhängigkeit der Simulationsergebnisse auf. Um dieser Problematik zu begegnen und netzunabhängige Simulationen zu erreichen, werden Regularisierungstechniken angewendet. In dieser Arbeit wird die Bruchenergetische Regularisierung eingesetzt, die durch die Einführung einer äquivalenten Breite in ein lokal formuliertes Stoffgesetz gekennzeichnet ist. Die spezielle Wahl eines Wertes für die äquivalente Breite beruht auf der Forderung, dass in der Simulation die korrekte Bruchenergie je Einheitsfläche für den Bruchprozess verbraucht wird, d.h. die Energiedissipation der Realität entspricht. In vorliegender Arbeit wird die neue These aufgestellt, dass die Energiedissipation nur für den Fall korrekt abgebildet wird, wenn die im Stoffgesetz enthaltene äquivalente Breite in jedem Belastungsinkrement der Breite des Bereiches entspricht, in dem in der Simulation Energie dissipiert wird. In einer Simulation wird in den Bereichen Energie dissipiert, in denen die Schädigung im aktuellen Belastungsinkrement zunimmt. In vorliegender Arbeit werden die energiedissipierenden Bereiche daher als Pfad der Schädigungsrate bezeichnet. Um Erkenntnisse über die Entwicklung des Pfades der Schädigungsrate über den Belastungsverlauf zu erhalten, wurden umfangreiche Untersuchungen anhand von Simulationen eines beidseitig gekerbten Betonprobekörpers unter kombinierter Zug-Schubbeanspruchung durchgeführt, wobei die gewählten Werte für die äquivalente Breite variiert wurden. Es wurde stets eine Diskretisierung mit linearen Verschiebungselementen verwendet, wobei die Bereiche mit zu erwartender Schädigung feiner und regelmäßig mit Elementen quadratischer Geometrie diskretisiert wurden. Die Ergebnisse der Untersuchungen zeigen, dass die Breite des Pfades der Schädigungsrate abhängig ist von der Schädigung am betrachteten Materialpunkt, dem von Schädigungsrichtung und Elementkante eingeschlossenen Winkel, der Elementgröße und den Materialparametern. Um die geforderte Übereinstimmung von äquivalenter Breite und der Breite des Pfades der Schädigungsrate zu erreichen, werden neue Ansätze für die äquivalente Breite vorgeschlagen, die die erwähnten Einflüsse berücksichtigen. Simulationen unter Verwendung der neuen Ansätze für die äquivalente Breite führen zu einer guten Übereinstimmung von äquivalenter Breite und der Breite des Pfades der Schädigungsrate in der Simulation. Die Ergebnisse der Simulationen, wie z.B. Last-Verformungsbeziehung und Rissverläufe, sind netzunabhängig und stimmen gut mit den experimentellen Beobachtungen überein. Basierend auf den gewonnenen Erkenntnissen wird eine Erweiterung der Bruchenergetischen Regularisierung vorgeschlagen: die Adaptive Bruchenergetische Regularisierung. Im abschließenden Kapitel der Arbeit werden mit der vorgeschlagenen Theorie, dem neuen Schädigungsmodell und der Adaptiven Bruchenergetischen Regularisierung, noch zwei in der Literatur gut dokumentierte Versuche simuliert. Die Simulationsergebnisse entsprechen den experimentellen Beobachtungen. / This doctoral thesis deals with the simulation of predominantly tensile loaded plain concrete structures. Concrete is modeled on the macro level and the Finite Element Method is applied to solve the resulting mechanical field problem. A new material model for concrete based on continuum damage mechanics and an extended regularization technique based on the fracture energy approach are presented. The thesis is subdivided into two parts. In the first part, a local, anisotropic damage model for concrete is derived. This model uses a symmetric second-order tensor as the damage variable, which enables the simulation of orthotropic degradation. The validity of the first and the second law of thermodynamics as well as the validity of the principle of maximum dissipation rate are required. Using a normal rule in the space of the dissipative forces, which are the thermodynamically conjugated variables to the damage variables, and the definition of the loading functions in the space of the dissipative forces guarantee their validity. The suggested formulation contains three decoupled loading functions. A further requirement in the derivation of the model was the minimization of the number of material parameters, which should be determined by a small number of standard experiments. The material parameters of the new damage model are the Young’s modulus, the Poisson’s ratio, the tensile strength and the fracture energy per unit area. The second part of the work focuses on localization and regularization. If a Finite Element simulation is performed using a local material model for concrete, the results of the Finite Element simulation are mesh-dependent. To attain mesh-independent simulations, a regularization technique must be applied. The fracture energy approach, which is characterized by introducing a characteristic length in a locally formulated material model, is used as regularization technique in this work. The choice of a value for the characteristic length is founded by the requirement, that the fracture energy per unit area, which is consumed for the fracture process in the simulation, must be the same as in experiment, i.e. the energy dissipation must be correct. In this dissertation, the new idea is suggested that the correct energy dissipation can be only attained if the characteristic length in the material model coincides in every loading increment with the width of the energy-dissipating zone in the simulation. The energy-dissipating zone in a simulation is formed by the integration points with increasing damage and obtains the name: damage rate path. Detailed investigations based on simulations of a double-edge notched specimen under mixed-mode loading are performed with varying characteristic lengths in order to obtain information concerning the evolution of the damage rate path during a simulation. All simulations were performed using displacement-based elements with four nodes. The range with expected damage was always finer and regularly discretized. The results of the simulations show that the width of the damage rate path depends on the damage at the specific material point, on the angle between damage direction and element edges, on the element size and on the material parameters. Based on these observations, new approaches for the characteristic length are suggested in order to attain the coincidence of the characteristic length with the width of the damage rate path. Simulations by using the new approaches yield a sufficient coincidence of the characteristic length with the width of the damage rate path. The simulations are mesh-independent and the results of the simulation, like load-displacement curves or crack paths, correspond to the experimental results. Based on all new information concerning the regularization technique, an extension of the fracture energy approach is suggested: the adaptive fracture energy approach. The validity and applicability of the suggested theory, the new anisotropic damage model and the adaptive fracture energy approach, are verified in the final chapter of the work with simulations of two additional experiments, which are well documented in the literature. The results of the simulations correspond to the observations in the experiments.

This work aims at dimensional reduction of nonlinear material models in an asymptotically accurate manner. The three-dimensional(3-D) nonlinear material models considered include isotropic, orthotropic and dielectric compressible hyperelastic material models. Hyperelastic materials have potential applications in space-based inflatable structures, pneumatic membranes, replacements for soft biological tissues, prosthetic devices, compliant robots, high-altitude airships and artificial blood pumps, to name a few. Such structures have special engineering properties like high strength-to-mass ratio, low deflated volume and low inflated density. The majority of these applications imply a thin shell form-factor, rendering the problem geometrically nonlinear as well. Despite their superior engineering properties and potential uses, there are no proper analysis tools available to analyze these structures accurately yet efficiently. The development of a unified analytical model for both material and geometric nonlinearities encounters mathematical difficulties in the theory but its results have considerable scope. Therefore, a novel tool is needed to dimensionally reduce these nonlinear material models. In this thesis, Prof. Berdichevsky’s Variational Asymptotic Method(VAM) has been applied rigorously to alleviate the difficulties faced in modeling thin shell structures(made of such nonlinear materials for the first time in the history of VAM) which inherently exhibit geometric small parameters(such as the ratio of thickness to shortest wavelength of the deformation along the shell reference surface) and physical small parameters(such as moderate strains in certain applications). Saint Venant-Kirchhoff and neo-Hookean 3-D strain energy functions are considered for isotropic hyperelastic material modeling. Further, these two material models are augmented with electromechanical coupling term through Maxwell stress tensor for dielectric hyperelastic material modeling. A polyconvex 3-D strain energy function is used for the orthotropic hyperelastic model. Upon the application of VAM, in each of the above cases, the original 3-D nonlinear electroelastic problem splits into a nonlinear one-dimensional (1-D) through-the-thickness analysis and a nonlinear two-dimensional(2-D) shell analysis. This greatly reduces the computational cost compared to a full 3-D analysis. Through-the-thickness analysis provides a 2-D nonlinear constitutive law for the shell equations and a set of recovery relations that expresses the 3-D field variables (displacements, strains and stresses) through thethicknessintermsof2-D shell variables calculated in the shell analysis (2-D). Analytical expressions (asymptotically accurate) are derived for stiffness, strains, stresses and 3-D warping field for all three material types. Consistent with the three types of 2-D nonlinear constitutive laws,2-D shell theories and corresponding finite element programs have been developed. Validation of present theory is carried out with a few standard test cases for isotropic hyperelastic material model. For two additional test cases, 3-Dfinite element analysis results for isotropic hyperelastic material model are provided as further proofs of the simultaneous accuracy and computational efficiency of the current asymptotically-correct dimensionally-reduced approach. Application of the dimensionally-reduced dielectric hyperelastic material model is demonstrated through the actuation of a clamped membrane subjected to an electric field. Finally, the through-the-thickness and shell analysis procedures are outlined for the orthotropic nonlinear material model.

Possibilities of modeling non-linear behavior of concrete within standard room temperatures and increased fire-load values using FEM software ANSYS are studied. Temperature dependences of material models are considered. Fire resistance of reinforced concrete and concrete-steel composite construction is analyzed. Fire loads are defined in accordance with relevant standards. Non-linear structural transient analyses are calculated after temperatures were determined by transient thermal analyses. Results obtained from analyses of simple reinforced concrete structure are compared with approach of isotherm 500 °C method.

The thesis deals with the back analysis of the instrumentated embankment on the subsoil, improved by prefabricated vertical drains. The first part is devoted to the theory of consolidation calculation. Furthermore, the author deals with a parametric study of the analytical calculation of prefabricated drains and the comparsion of the analytical and numerical solution. The last and the largest part of the work is a back analysis of the embankment which was built as a part of the subsoil improvement near the Suez Canal in Egypt by Keller GmbH.

In Stetigfördersystemen mit Gleitketten treten Schwingungen auf, die zu Beeinträchtigungen des Materialflusses führen können. Dazu zählen u. a. Kippen oder Rutschen der Fördergüter oder die Schädigung der Förderkette durch schwellende Belastung. Die vorliegende Arbeit befasst sich mit der Analyse und Simulation der auftretenden Schwingungen, mit dem Ziel, ein Simulationsmodell zu entwickeln, welches die dynamischen Effekte in einer Kunststoff-Gleitkette abbilden kann. Einführend werden Gleitkettenfördersysteme analysiert und die Anregungsursachen, die in dem System zu Schwingungen führen, betrachtet. Die Eigenfrequenz ist eine maßgebliche Größe des Schwingungsverhaltens. Deshalb werden mehrere Ansätze zur Berechnung der Eigenfrequenzen einer Gleitkette aufgestellt und deren Zweckmäßigkeit überprüft. Anschließend erfolgt die Erstellung eines Mehrkörper-Simulationsmodells, mit dem neben den Eigenfrequenzen auch die Amplitudenverläufe der Beschleunigung sowie der Kettenzugkraft berechnet werden können. Weiterhin wird aufgrund des viskoelastischen Materialverhaltens der Kunststoffkette ein geeignetes Materialmodell ermittelt und dessen Parameter über eine Optimierungsrechnung bestimmt. Nach der Validierung des Eigenfrequenz- und des Simulationsmodells an einer Versuchsförderanlage erfolgt eine Parameteranalyse, mit der die Auswirkungen der relevanten Einflussgrößen auf das dynamische Verhalten am Beispiel des Versuchsfördersystems eruiert werden. Abschließend werden Empfehlungen zur Reduzierung von Schwingungen in Gleitkettenfördersystemen gegeben. / In continuous conveyor systems occur vibrations which can lead to impairments of material flow. These includes tilting or slipping of transported material as well as damage to the conveyor chain due to swelling stress. The present thesis is concerned with the analysis and simulation of occurring vibrations with the aim of developing a simulation model which can depict the dynamic effects in a plastic slide chain. Introductorily, conveyor systems with sliding chains are analyzed and causes of excitation are examined, which lead to vibrations in the system. The natural frequency is a significant variable of vibration behavior. Therefore, several approaches for calculating the natural frequencies of a sliding chain are compiled and their suitability is reviewed. A multi-body simulation model is built which can be used to calculate amplitude curves of acceleration and chain traction force besides the natural frequencies. Due to the viscoelastic material behavior of plastic chains, a suitable material model is defined and its parameters are determined via an optimization calculation. After validating the natural frequency and the simulation model at a test conveyor, a parameter analysis is carried out with which effects of relevant influencing variables on the dynamic behavior are determined. Finally, recommendations are given for reducing vibrations in sliding chain conveyor systems.

Der Gegenstand der vorliegenden Arbeit ist die Simulation von Betonstrukturen beliebiger Geometrie unter überwiegender Zugbelastung. Die Modellierung erfolgt auf Makroebene als Kontinuum und zur Lösung des mechanischen Feldproblems wird die Finite-Elemente-Methode verwendet. Ein neues Materialmodell für Beton und eine Erweiterung der Bruchenergetischen Regularisierung werden vorgestellt. Die Arbeit ist in zwei Teile gegliedert. Im ersten Teil wird ein lokales, anisotropes Schädigungsmodell abgeleitet, wobei als Schädigungsvariable ein symmetrischer Tensor zweiter Stufe gewählt wird. Die Verwendung einer Normalenregel im Raum der dissipativen Kräfte zur Bestimmung der Schädigungsevolution und die Definition der Schädigungsgrenzflächen im Raum der dissipativen Kräfte gewährleisten die Gültigkeit der Hauptsätze der Thermodynamik und des Prinzips der maximalen Dissipationsrate. Vorteilhaft ist die Symmetrie der Materialtangente, die sich aus diesem Vorgehen ergibt. Eine Formulierung mit drei entkoppelten Schädigungsgrenzflächen wird vorgeschlagen. Eine wichtige Forderung bei der Ableitung des Materialmodells war die Verwendung einer möglichst geringen Anzahl von Materialparametern, welche darüber hinaus aus wenigen Standardversuchen bestimmbar sein sollten. Das Schädigungsmodell enthält als Materialparameter den Elastizitätsmodul, die Querdehnzahl, die Zugfestigkeit und die auf eine Einheitsfläche bezogene Bruchenergie. Im zweiten Teil der Arbeit stehen Lokalisierung und Regularisierung im Fokus der Betrachtungen. Aufgrund der lokalen Formulierung des Materialmodells tritt bei Finite-Elemente Simulationen eine Netzabhängigkeit der Simulationsergebnisse auf. Um dieser Problematik zu begegnen und netzunabhängige Simulationen zu erreichen, werden Regularisierungstechniken angewendet. In dieser Arbeit wird die Bruchenergetische Regularisierung eingesetzt, die durch die Einführung einer äquivalenten Breite in ein lokal formuliertes Stoffgesetz gekennzeichnet ist. Die spezielle Wahl eines Wertes für die äquivalente Breite beruht auf der Forderung, dass in der Simulation die korrekte Bruchenergie je Einheitsfläche für den Bruchprozess verbraucht wird, d.h. die Energiedissipation der Realität entspricht. In vorliegender Arbeit wird die neue These aufgestellt, dass die Energiedissipation nur für den Fall korrekt abgebildet wird, wenn die im Stoffgesetz enthaltene äquivalente Breite in jedem Belastungsinkrement der Breite des Bereiches entspricht, in dem in der Simulation Energie dissipiert wird. In einer Simulation wird in den Bereichen Energie dissipiert, in denen die Schädigung im aktuellen Belastungsinkrement zunimmt. In vorliegender Arbeit werden die energiedissipierenden Bereiche daher als Pfad der Schädigungsrate bezeichnet. Um Erkenntnisse über die Entwicklung des Pfades der Schädigungsrate über den Belastungsverlauf zu erhalten, wurden umfangreiche Untersuchungen anhand von Simulationen eines beidseitig gekerbten Betonprobekörpers unter kombinierter Zug-Schubbeanspruchung durchgeführt, wobei die gewählten Werte für die äquivalente Breite variiert wurden. Es wurde stets eine Diskretisierung mit linearen Verschiebungselementen verwendet, wobei die Bereiche mit zu erwartender Schädigung feiner und regelmäßig mit Elementen quadratischer Geometrie diskretisiert wurden. Die Ergebnisse der Untersuchungen zeigen, dass die Breite des Pfades der Schädigungsrate abhängig ist von der Schädigung am betrachteten Materialpunkt, dem von Schädigungsrichtung und Elementkante eingeschlossenen Winkel, der Elementgröße und den Materialparametern. Um die geforderte Übereinstimmung von äquivalenter Breite und der Breite des Pfades der Schädigungsrate zu erreichen, werden neue Ansätze für die äquivalente Breite vorgeschlagen, die die erwähnten Einflüsse berücksichtigen. Simulationen unter Verwendung der neuen Ansätze für die äquivalente Breite führen zu einer guten Übereinstimmung von äquivalenter Breite und der Breite des Pfades der Schädigungsrate in der Simulation. Die Ergebnisse der Simulationen, wie z.B. Last-Verformungsbeziehung und Rissverläufe, sind netzunabhängig und stimmen gut mit den experimentellen Beobachtungen überein. Basierend auf den gewonnenen Erkenntnissen wird eine Erweiterung der Bruchenergetischen Regularisierung vorgeschlagen: die Adaptive Bruchenergetische Regularisierung. Im abschließenden Kapitel der Arbeit werden mit der vorgeschlagenen Theorie, dem neuen Schädigungsmodell und der Adaptiven Bruchenergetischen Regularisierung, noch zwei in der Literatur gut dokumentierte Versuche simuliert. Die Simulationsergebnisse entsprechen den experimentellen Beobachtungen. / This doctoral thesis deals with the simulation of predominantly tensile loaded plain concrete structures. Concrete is modeled on the macro level and the Finite Element Method is applied to solve the resulting mechanical field problem. A new material model for concrete based on continuum damage mechanics and an extended regularization technique based on the fracture energy approach are presented. The thesis is subdivided into two parts. In the first part, a local, anisotropic damage model for concrete is derived. This model uses a symmetric second-order tensor as the damage variable, which enables the simulation of orthotropic degradation. The validity of the first and the second law of thermodynamics as well as the validity of the principle of maximum dissipation rate are required. Using a normal rule in the space of the dissipative forces, which are the thermodynamically conjugated variables to the damage variables, and the definition of the loading functions in the space of the dissipative forces guarantee their validity. The suggested formulation contains three decoupled loading functions. A further requirement in the derivation of the model was the minimization of the number of material parameters, which should be determined by a small number of standard experiments. The material parameters of the new damage model are the Young’s modulus, the Poisson’s ratio, the tensile strength and the fracture energy per unit area. The second part of the work focuses on localization and regularization. If a Finite Element simulation is performed using a local material model for concrete, the results of the Finite Element simulation are mesh-dependent. To attain mesh-independent simulations, a regularization technique must be applied. The fracture energy approach, which is characterized by introducing a characteristic length in a locally formulated material model, is used as regularization technique in this work. The choice of a value for the characteristic length is founded by the requirement, that the fracture energy per unit area, which is consumed for the fracture process in the simulation, must be the same as in experiment, i.e. the energy dissipation must be correct. In this dissertation, the new idea is suggested that the correct energy dissipation can be only attained if the characteristic length in the material model coincides in every loading increment with the width of the energy-dissipating zone in the simulation. The energy-dissipating zone in a simulation is formed by the integration points with increasing damage and obtains the name: damage rate path. Detailed investigations based on simulations of a double-edge notched specimen under mixed-mode loading are performed with varying characteristic lengths in order to obtain information concerning the evolution of the damage rate path during a simulation. All simulations were performed using displacement-based elements with four nodes. The range with expected damage was always finer and regularly discretized. The results of the simulations show that the width of the damage rate path depends on the damage at the specific material point, on the angle between damage direction and element edges, on the element size and on the material parameters. Based on these observations, new approaches for the characteristic length are suggested in order to attain the coincidence of the characteristic length with the width of the damage rate path. Simulations by using the new approaches yield a sufficient coincidence of the characteristic length with the width of the damage rate path. The simulations are mesh-independent and the results of the simulation, like load-displacement curves or crack paths, correspond to the experimental results. Based on all new information concerning the regularization technique, an extension of the fracture energy approach is suggested: the adaptive fracture energy approach. The validity and applicability of the suggested theory, the new anisotropic damage model and the adaptive fracture energy approach, are verified in the final chapter of the work with simulations of two additional experiments, which are well documented in the literature. The results of the simulations correspond to the observations in the experiments.

Indiana University-Purdue University Indianapolis (IUPUI) / In this thesis, an explicit time integration code which integrates multibody dynamics and the discrete element method is used for modelling the excavation and moving operation of cohesive soft soil (such as mud and snow) by bulldozers. A soft cohesive soil material model (that includes normal and tangential inter-particle force models) is used that can account for soil compressibility, plasticity, fracture, friction, viscosity and gain in cohesive strength due to compression. In addition, a time relaxation sub-model for the soil plastic deformation and cohesive strength is added in order to account for loss in soil cohesive strength and reduced bulk density due to tension or removal of the compression. This is essential in earth moving applications since the soil that is dug typically becomes loose soil that has lower shear strength and lower bulk density (larger volume) than compacted soil. If the model does not account for loss of soil shear strength then the dug soil pile in front of the blade of a bulldozer will have an artificially high shear strength. A penalty technique is used to impose joint and normal contact constraints. An asperity-based friction model is used to model contact and joint friction. A Cartesian Eulerian grid contact search algorithm is used to allow fast contact detection between particles. A recursive bounding box contact search algorithm is used to allow fast contact detection between the particles and polygonal contact surfaces. A multibody dynamics bulldozer model is created which includes the chassis/body, C-frame, blade, wheels and hydraulic actuators. The components are modelled as rigid bodies and are connected using revolute and prismatic joints. Rotary actuators along with PD (Proportional-Derivative) controllers are used to drive the wheels. Linear actuators along with PD controllers are used to drive the hydraulic actuators. Polygonal contact surfaces are defined for the tires and blade to model the interaction between the soil and the bulldozer. Simulations of a bulldozer performing typical shallow digging operations in a cohesive soil are presented. The simulation of a rear wheel drive bulldozer shows that, it has a limited digging capacity compared to the 4-wheel drive bulldozer. The effect of the relaxation parameter can be easily observed from the variation in the Bulldozer's velocity. The higher the relaxation parameter, the higher is the bulldozer's velocity while it is crossing over the soil patch. For the low penetration depth run the bulldozer takes less time compared to high penetration depth. Also higher magnitudes of torques at front and rear wheels can be observed in case of high penetration depth. The model is used to predict the wheel torque, wheel speed, vehicle speed and actuator forces during shallow digging operations on three types of soils and at two blade penetration depths. The model presented can be used to predict the motion, loads and required actuators forces and to improve the design of the various bulldozer components such as the blade, tires, engine and hydraulic actuators.

Diploma thesis is focused on the description of the parameters of nonlinear material models of concrete, which are implemented in a computational system LS-DYNA, interacting with performance of nonlinear test calculations in system LS-DYNA on selected problems, which are formed mainly by simulations of tests of mechanical and physical properties of concrete in uniaxial compressive and tensile on cylinders with applying different boundary conditions and by simulation of bending slab, with subsequent comparison of some results of test calculations with results of the experiment. The thesis includes creation of appropriate geometric models of selected problems, meshing of these geometric models, description of parameters and application of nonlinear material models of concrete on selected problems, application of loads and boundary conditions on selected problems and performance of nonlinear calculations in a computational system LS-DYNA. Evaluation of results is made on the basis of stress-strain diagrams and load-displacement diagrams based on nonlinear calculations taking into account strain rate effects and on the basis of hysteresis curves based on nonlinear calculations in case of application of cyclic loading on selected problems. Verification of nonlinear material models of concrete is made on the basis of comparison of some results of test calculations with results obtained from the experiment.