:py:mod:`femo_alpha.dynamic_rm_shell.plate_sim` =============================================== .. py:module:: femo_alpha.dynamic_rm_shell.plate_sim Module Contents --------------- Classes ~~~~~~~ .. autoapisummary:: femo_alpha.dynamic_rm_shell.plate_sim.BC_list femo_alpha.dynamic_rm_shell.plate_sim.PlateSim Functions ~~~~~~~~~ .. autoapisummary:: femo_alpha.dynamic_rm_shell.plate_sim.extractTipDispDolfinx Attributes ~~~~~~~~~~ .. autoapisummary:: femo_alpha.dynamic_rm_shell.plate_sim.comm .. py:class:: BC_list(W_space, direction=0, custom_bc_func=None) Bases: :py:obj:`object` Setting Dirichlet boundary condition .. !! processed by numpydoc !! .. py:method:: return_affected_dof_tuple() .. py:method:: return_bc_list() .. py:class:: PlateSim(mesh, E, nu, rho, dt, Nsteps, element_wise_thickness=False, custom_bc_func=None, add_self_weight=False, g_factor=None, quad_deg=3, comm=MPI.COMM_WORLD) Bases: :py:obj:`object` Implements St. Venant--Kirchhoff constitutive model for isogeometric shell. .. !! processed by numpydoc !! .. py:method:: SVK_residual() .. py:method:: assembleStrainEnergy(w) .. py:method:: compute_solution_vectors_at_time(disp_arr, i, backwards_timestepping=False) .. py:method:: compute_wdot_vector_at_time_n(w_arr) .. py:method:: constructStrainEnergy(w) .. py:method:: construct_CG2_CG1_interpolation_map() .. py:method:: construct_disp_extraction_mats() .. py:method:: construct_force_to_pressure_map() .. py:method:: construct_nodal_disp_map() .. py:method:: create_mat_from_form(inp_form, apply_bcs=False) .. py:method:: dRdt(svk_res, w_cur, w_old, wdot_old, adjoint=True) .. py:method:: dRdw(svk_res, w_cur, w_old, wdot_old, adjoint=True) This function assembles the derivative matrix dR/dw, where R is the residual vector of all time steps and w contains all displacements and rotations. Inputs: w_cur: Solution vector at current time step w_old: Solution vector at previous time step wdot: Solution time derivative at current time step .. !! processed by numpydoc !! .. py:method:: dRdwdot(svk_res, w_cur, w_old, wdot_old) .. py:method:: eval_fe_basis_all_dolfinx(x, dofmap_adjacencylist, mesh_bbt, basix_element, x_idx, mat_shape) .. py:method:: pnorm_stress(m=1e-06, rho=100, alpha=None, regularization=False) Compute the p-norm of the stress `rho` is the Constraint aggregation factor .. !! processed by numpydoc !! .. py:method:: reset_solution_vectors() .. py:method:: set_solution_vector(func_obj, arr) :staticmethod: .. py:method:: set_solution_vectors_for_adjoint(w, w_old, wdot_old) .. py:method:: set_up_tip_dofs(x_tip, cell_tip) .. py:method:: solve_dynamic_problem(residual, saving_outputs=False, PATH=None, POD_matrix=None, POD_mean=None, timing=False) .. py:method:: update_f(f_array) .. py:method:: update_f_history(f_history_array) .. py:method:: update_materialmodel() .. py:method:: update_nsteps(Nsteps) .. py:method:: update_t(t_array) .. py:method:: update_wdot_vector_at_time(w_vec, w_old_vec, wdot_old_vec) .. py:method:: update_wdot_vector_at_time_reverse(w_vec, w_old_vec, wdot_vec) .. py:method:: von_Mises_stress() .. py:method:: wdot_vector() .. py:function:: extractTipDispDolfinx(w, x_tip=[10.0, 0.0, 0.0], cell_tip=18) .. py:data:: comm