<?xml version="1.0" encoding="UTF-8"?>
<!-- Siconos is a program dedicated to modeling, simulation and control-->
<!-- of non smooth dynamical systems.-->
<!---->
<!-- Copyright 2020 INRIA.-->
<!---->
<!-- Licensed under the Apache License, Version 2.0 (the "License");-->
<!-- you may not use this file except in compliance with the License.-->
<!-- You may obtain a copy of the License at-->
<!---->
<!-- http://www.apache.org/licenses/LICENSE-2.0-->
<!---->
<!-- Unless required by applicable law or agreed to in writing, software-->
<!-- distributed under the License is distributed on an "AS IS" BASIS,-->
<!-- WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.-->
<!-- See the License for the specific language governing permissions and-->
<!-- limitations under the License.-->
<!--	-->

<!-- SICONOS Model XML file -->

<SiconosModel>
	<Title>Dry Friction Oscillator</Title>
	<Author>F.Perignon</Author>
	<Description>Oscillator with Coulomb friction </Description>
	<Date>date</Date>
	<Time>
		<t0>0.0</t0>
		<T>40.0</T>
	</Time>
	<NSDS bvp='false'>

		<!-- DSs defined in the problem -->
		<DS_Definition>
			<LagrangianLinearTIDS Name='Oscillator'>
				<q0 vectorSize='3'>
					0.1 0.0 0.0
				</q0>
				<Velocity0  vectorSize='3'>
					0.0 0.0 10.0
				</Velocity0>
				<FExt vectorPlugin="BeadPlugin:FExt"/>
				<Mass matrixRowSize='3' matrixColSize='3'>
					<row>1.0 0.0 0.0</row>
					<row>0.0 1.0 0.0</row>
					<row>0.0 0.0 1.0</row>
				</Mass>
				<K matrixRowSize='3' matrixColSize='3'>
					<row>0.0 0.0 0.0</row>
					<row>0.0 0.0 0.0</row>
					<row>0.0 0.0 0.0</row>
				</K>
				<C matrixRowSize='3' matrixColSize='3'>
					<row>0.0 0.0 0.0</row>
					<row>0.0 0.0 0.0</row>
					<row>0.0 0.0 0.0</row>
				</C>
			</LagrangianLinearTIDS>
		</DS_Definition>
		<!-- Interactions defined in the problem -->
		<Interaction_Definition>
			<Interaction Name='1'>
				<!-- List of couple of DS concerned by this interaction -->
				<DS_Concerned>Oscillator</DS_Concerned>
				<!-- Relation of this interaction -->
				<Interaction_Content>
					<LagrangianLinearTIR>
						<C matrixRowSize='2' matrixColSize='3'>
							<row> 0.0 1.0 0.0 </row>
							<row> 1.0 0.0 0.0</row>
						</C>
					</LagrangianLinearTIR>
					<!-- NS Law of this interaction	-->
					<NewtonImpactFrictionNSL size='2'>
						<en>0.0</en>
						<mu>0.08</mu>
					</NewtonImpactFrictionNSL>
				</Interaction_Content>
			</Interaction>
		</Interaction_Definition>
	</NSDS>
	<!-- Simulation to use in order to solve the problem -->
	<Simulation type='TimeStepping'>
		<TimeDiscretisation>
			<h>0.005</h>
		</TimeDiscretisation>
		<!-- One Step Integrators -->
		<OneStepIntegrator_Definition>
			<MoreauJeanOSI>
				<DS_Concerned>Oscillator</DS_Concerned>
				<Theta all="0.5"></Theta>
			</MoreauJeanOSI>
		</OneStepIntegrator_Definition>
		<!-- One Step NS Problem of th Simulation-->
		<OneStepNSProblems_List>
			<FrictionContact Type='2'>
			</FrictionContact>
		</OneStepNSProblems_List>
	</Simulation>
</SiconosModel>
