.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/EM/em_rlc_isopotential.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_EM_em_rlc_isopotential.py: RLC circuit by isopotential =========================== This example shows how to use the corresponding ``contype`` in the ``connect_isopotential()`` method to define an RLC circuit as an inlet boundary condition. The executable file for LS-DYNA is ``ls-dyna_smp_d_R13.1_138-g8429c8a10f_winx64_ifort190.exe``. .. GENERATED FROM PYTHON SOURCE LINES 10-114 .. rst-class:: sphx-glr-script-out .. code-block:: none '/server/output' | .. code-block:: Python import os import sys from em_set_data import rlc_rogoseg from ansys.dyna.core.pre import examples, launch_dynapre from ansys.dyna.core.pre.dynaem import ( DynaEM, EMType, Isopotential, Isopotential_ConnType, NodeSet, RogoCoil, SegmentSet, SolidFormulation, SolidPart, ) from ansys.dyna.core.pre.dynamaterial import EMMATTYPE, MatRigid from ansys.dyna.core.pre.misc import check_valid_ip hostname = "localhost" if len(sys.argv) > 1 and check_valid_ip(sys.argv[1]): hostname = sys.argv[1] solution = launch_dynapre(ip=hostname) fns = [] path = examples.em_rlc_isopotential + os.sep fns.append(path + "em_rlc_isopotential.k") solution.open_files(fns) solution.set_termination(termination_time=0.01) solution.create_database_binary(dt=1e-4) emobj = DynaEM() solution.add(emobj) emobj.set_timestep(tssfac=1, timestep_size_for_mass_scaled=1e-4) emobj.analysis.set_timestep(timestep=1e-4) emobj.analysis.set_em_solver(type=EMType.RESISTIVE_HEATING) matrigid = MatRigid( mass_density=7000, young_modulus=2e11, center_of_mass_constraint=1, translational_constraint=7, rotational_constraint=7, ) matrigid.set_em_permeability_equal(material_type=EMMATTYPE.CONDUCTOR, initial_conductivity=1e4) part1 = SolidPart(1) part1.set_material(matrigid) part1.set_element_formulation(SolidFormulation.CONSTANT_STRESS_SOLID_ELEMENT) emobj.parts.add(part1) nset1 = NodeSet( [ 429, 433, 437, 441, 445, 449, 453, 457, 461, 465, 469, 473, 477, 481, 485, 489, 493, 497, 501, 505, 509, 513, 517, 521, 525, ] ) nset2 = NodeSet( [26, 31, 36, 41, 46, 51, 56, 61, 66, 71, 76, 81, 86, 91, 96, 101, 106, 111, 116, 121, 126, 131, 136, 141, 146] ) isopos_conn1 = Isopotential(nset1) isopos_conn2 = Isopotential(nset2) emobj.connect_isopotential( contype=Isopotential_ConnType.RLC_CIRCUIT, isopotential1=isopos_conn1, value=5e-4, inductance=7.8e-5, capacity=0.0363, initial_voltage=5000, ) emobj.connect_isopotential(contype=Isopotential_ConnType.VOLTAGE_SOURCE, isopotential1=isopos_conn2) emobj.add(RogoCoil(SegmentSet(rlc_rogoseg))) emobj.create_em_output(mats=2, matf=2, sols=2, solf=2) solution.save_file() .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 0.033 seconds) .. _sphx_glr_download_examples_EM_em_rlc_isopotential.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: em_rlc_isopotential.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: em_rlc_isopotential.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: em_rlc_isopotential.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_