Location: BG Mica models @ e105ad14397f / BG tutorial model solid mechanics 4 Lorenz s-domain.cellml

Author:
Peter H <p.hunter@auckland.ac.nz>
Date:
2017-05-13 17:27:58+12:00
Desc:
update
Permanent Source URI:
https://models.physiomeproject.org/workspace/496/rawfile/e105ad14397fd6fc3de65305e5d042795124869e/BG tutorial model solid mechanics 4 Lorenz s-domain.cellml

<?xml version='1.0'?>
<model name="transfer_function" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
    <!-- transfer function for Lorenz force actuator-->
    <units name="per_s">
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_s2">
        <unit exponent="-2" units="second"/>
    </units>
    <units name="J_per_C">
        <unit units="joule"/>
        <unit exponent="-1" units="coulomb"/>
    </units>
    <units name="J_per_C2">
        <unit units="joule"/>
        <unit exponent="-2" units="coulomb"/>
    </units>
    <units name="Js_per_C2">
        <unit units="joule"/>
        <unit units="second"/>
        <unit exponent="-2" units="coulomb"/>
    </units>
    <units name="Js2_per_C2">
        <unit units="joule"/>
        <unit exponent="2" units="second"/>
        <unit exponent="-2" units="coulomb"/>
    </units>
    <units name="m_per_s">
        <unit units="metre"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="m_per_s2">
        <unit units="metre"/>
        <unit exponent="-2" units="second"/>
    </units>
    <units name="J_per_m">
        <unit units="joule"/>
        <unit exponent="-1" units="metre"/>
    </units>
    <units name="J_per_m2">
        <unit units="joule"/>
        <unit exponent="-2" units="metre"/>
    </units>
    <units name="J2_per_m4">
        <unit exponent="2" units="joule"/>
        <unit exponent="-4" units="metre"/>
    </units>
    <units name="Js_per_m2">
        <unit units="joule"/>
        <unit units="second"/>
        <unit exponent="-2" units="metre"/>
    </units>
    <units name="J2s2_per_m4">
        <unit exponent="2" units="joule"/>
        <unit exponent="2" units="second"/>
        <unit exponent="-4" units="metre"/>
    </units>
    <units name="Js2_per_m2">
        <unit units="joule"/>
        <unit exponent="2" units="second"/>
        <unit exponent="-2" units="metre"/>
    </units>
    <units name="Js_per_C_m">
        <unit units="joule"/>
        <unit units="second"/>
        <unit exponent="-1" units="coulomb"/>
        <unit exponent="-1" units="metre"/>
    </units>
    <component name="main">
        <variable initial_value="0" name="t" units="second"/>
        <variable initial_value="0.1" name="omega" units="per_s"/>
        <variable initial_value="1" name="omega_ref" units="per_s"/>
        <variable name="logOmega" units="dimensionless"/>
        <variable initial_value="0" name="E_1" units="J_per_C2"/>
        <variable initial_value="2000" name="E_2" units="J_per_m2"/>
        <variable initial_value="3.5" name="R_1" units="Js_per_C2"/>
        <variable initial_value="0.4" name="R_2" units="Js_per_m2"/>
        <variable initial_value="0.15" name="L_1" units="Js2_per_C2"/>
        <variable initial_value="0.018" name="L_2" units="Js2_per_m2"/>
        <variable initial_value="10" name="Bl" units="Js_per_C_m"/>
        <variable name="x_1" units="J_per_m2"/>
        <variable name="x_2" units="J2_per_m4"/>
        <variable name="omega_3" units="per_s"/>
        <variable name="logOmega_3" units="dimensionless"/>
        <variable name="G_real" units="Js_per_C2"/>
        <variable name="G_imag" units="Js_per_C2"/>
        <variable name="amplitude" units="Js_per_C2"/>
        <variable initial_value="1" name="amplitude_ref" units="Js_per_C2"/>
        <variable name="phase" units="dimensionless"/>
        <variable initial_value="1" name="phase_ref" units="dimensionless"/>
        <variable name="phase_degrees" units="dimensionless"/>
        <variable name="logAmplitude" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>omega_3</ci>
                <apply>
                    <root/>
                    <apply>
                        <divide/>
                        <ci>E_2</ci>
                        <ci>L_2</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>logOmega_3</ci>
                <apply>
                    <log/>
                    <apply>
                        <divide/>
                        <ci>omega_3</ci>
                        <ci>omega_ref</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>omega</ci>
                </apply>
                <cn cellml:units="per_s2">1</cn>
            </apply>
            <apply>
                <eq/>
                <ci>logOmega</ci>
                <apply>
                    <log/>
                    <apply>
                        <divide/>
                        <ci>omega</ci>
                        <ci>omega_ref</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>x_1</ci>
                <apply>
                    <minus/>
                    <ci>E_2</ci>
                    <apply>
                        <times/>
                        <apply>
                            <power/>
                            <ci>omega</ci>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <ci>L_2</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>x_2</ci>
                <apply>
                    <plus/>
                    <apply>
                        <power/>
                        <ci>x_1</ci>
                        <cn cellml:units="dimensionless">2</cn>
                    </apply>
                    <apply>
                        <power/>
                        <apply>
                            <times/>
                            <ci>omega</ci>
                            <ci>R_2</ci>
                        </apply>
                        <cn cellml:units="dimensionless">2</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>G_real</ci>
                <apply>
                    <plus/>
                    <ci>R_1</ci>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R_2</ci>
                            <apply>
                                <power/>
                                <apply>
                                    <times/>
                                    <ci>omega</ci>
                                    <ci>Bl</ci>
                                </apply>
                                <cn cellml:units="dimensionless">2</cn>
                            </apply>
                        </apply>
                        <ci>x_2</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>G_imag</ci>
                <apply>
                    <times/>
                    <ci>omega</ci>
                    <apply>
                        <plus/>
                        <ci>L_1</ci>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <apply>
                                    <power/>
                                    <ci>Bl</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                </apply>
                                <ci>x_1</ci>
                            </apply>
                            <ci>x_2</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>amplitude</ci>
                <apply>
                    <root/>
                    <apply>
                        <plus/>
                        <apply>
                            <power/>
                            <ci>G_real</ci>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                            <power/>
                            <ci>G_imag</ci>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>phase</ci>
                <apply>
                    <arctan/>
                    <apply>
                        <divide/>
                        <ci>G_imag</ci>
                        <ci>G_real</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>phase_degrees</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>phase</ci>
                        <cn cellml:units="dimensionless">180</cn>
                    </apply>
                    <pi/>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>logAmplitude</ci>
                <apply>
                    <log/>
                    <apply>
                        <divide/>
                        <ci>amplitude</ci>
                        <ci>amplitude_ref</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
</model>