주요 콘텐츠
결과: 30개
Time-dependent evolution of the continuum undergoing combined elasto-plasto-damage transformation.
interpolate_cbrewer(cbrew_init, interp_method, ncolors)
- % INTERPOLATE_CBREWER - interpolate a colorbrewer map to ncolors levels
visualize_paper1_graphs.m
- shift=10; %for different pictures
platehole_mesh(number_of_refinements)
- current directory stored
entryInWhichRows(A)
- function: entryrows=entryInWhichRows(A)
getEdges_triangles(elements)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
visualize_all.m
- return
refinement_uniform_2D(coordinates,elements,dirichlet,neumann,element_parameters)
- function: [coordinates,elements3,dirichlet]=refinement_uniform(coordinates,elements3,dirichlet)
refinement_uniform(coordinates,elements3,dirichlet,neumann,element_parameters)
- function: [coordinates,elements3,dirichlet]=refinement_uniform(coordinates,elements3,dirichlet)
cbrewer(ctype, cname, ncol, interp_method)
- % CBREWER - This function produces a colorbrewer table (rgb data) for a
evaluate_gv(t)
- Calculating gv
conv_ma2av(ma)
- converts matrix to array of vectors
compute_energy_balance.m
- two sided estimate
dissipation_plastic(Zeta,P1_derivative)
- plastic dissipation not depending on damage
phider
- PHIDER Returns the gradients of the basis functions
bulk_mesh(number_of_refinements,internal_slip_BC)
- Uniforme Verfeinerung
shapeder
- SHAPEDER Returns the gradients of the shape functions with
square_mesh_symmetric(number_of_refinements,slipWidthPercentage,nonsymmetric_BC)
- test2_nonsymmetricBC=-1; %enforces beam 1D with Dirichlet BC on the complete right edge
ring_mesh(number_of_refinements)
- A quarter of a ring geometry and boundary conditions
Computation of Friedrichs and Poincare constants in 2D/3D from generalized eigenvalue problems.
affine_transformations(nodes2coord,elems2nodes)
- AFFINE_TRANSFORMATIONS
basis_P1(p)
- BASIS_P1
show_edges(edge2nodes,coordinates,which)
- function: show_edges(edge,coordinates,whichelements)
phider
- PHIDER Returns the gradients of the basis functions
copy_triu(ama)
- COPY_TRIU
assembly_all_matrices.m
- stiffness matrix assembly for RT0 element
getEdges(elems2nodes)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
mass_matrix_RT0( elems, B_K, B_K_det, signs, material )
- MASS_MATRIX_RT0
compute_constants.m
- eigenvalues computation
inttet(p)
- INTTET
get_edges(elems2nodes)
- function: [element2edges, edge2nodes]=get_edges(elems2nodes)
get_boundary_edges(edges2nodes,bedges2nodes)
- GETBOUNDARYEDGES
inttri(p)
- INTTRI
shapeder
- SHAPEDER Returns the gradients of the shape functions with
basis_RT0(p)
- BASIS_RT0
tetrarefine3(XYZ,Elm,ERef)
- author: Ales Janka, ales.janka@unifr.ch, http://perso.unifr.ch/ales.janka
mass_matrix_P1( elems2nodes, B_K_det, material )
- MASS_MATRIX_P1
signs_edges(elems2nodes)
- SIGNS_EDGES
start_2D.m
- level_CM_zeroEigenvalueMultiplicity=zeros(levels,1);
refinement_uniform_3D(nodes2coord,elems2nodes,dirichlet)
- function: [nodes2coord,elems2nodes,dirichlet]=refinement_uniform_3D(nodes2coord,elems2nodes,dirichlet)
eigs_smallest_positive(K,M)
- two_smallest_eigenvalues_Poincare=eigs(K_P1,M_P1,2,-1e-3);
show_nodal_scalar(nodalValue,nodes2coord,elems2nodes,nodalDisplacement)
- fill3(X(elems2nodes)',Y(elems2nodes)',nodalValue(elems2nodes)',nodalValue(elems2nodes)');
solution_RT0( elems, B_K, B_K_det, signs, u_h, qr )
- SOLUTION_RT0
show_elements(elements,coordinates,which)
- function: show_elements(elements,coordinates,whichelements)
get_midpoints( elems2nodes, nodes2coord )
- MIDPOINTS
mass_matrix_Nedelec0( elems2edges, B_K, B_K_det, signs, material )
- MASS_MATRIX_NEDELEC0
basis_Nedelec0(p)
- BASIS_NEDELEC0
add_paths.m
- This script adds the absolute location of
start_Friedrichs.m
- coarse mesh of the unit square
stiffness_matrix_RT0( elems, B_K_det, signs, material )
- STIFFNESS_MATRIX_RT0
solution_P1( elems2nodes, B_K, u_h, qr)
- SOLUTION_P1
get_boundary_faces(elems2faces,faces2nodes)
- GETBOUNDARYFACES
refinement_uniform(coordinates,elements3,dirichlet)
- function: [coordinates,elements3,dirichlet]=refinement_uniform(coordinates,elements3,dirichlet)
Evaluations of Cournot-Nash and Stackelber-Cournot-Nash equilibria for models with several firms using Gauss method
Contact problem is solved by the Uzawa algorithm and the error of the solution is estimated.
start.m
- SETUP START
energy(v)
- we assume that v is P1 fuction
u_exact(f,phi,x)
- phi is scalar here, not the vector of constraints!
energy_norm(v)
- assume v is a piecewise linear function
compute_exact_error.m
- exact error computation
Animation of a rotating elastic ring subjected to time dependent periodical boundary conditions.
Finite element solution of the double porosity model and its a posteriori error estimate
deleterepeatedrows(matrix)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
flux_optimization.m
- geometry properties (can not be avoided)
entryInWhichRows(A)
- function: entryrows=entryInWhichRows(A)
start.m
- format long
Verfeinerung_Uniform(coordinates,elements3,dirichlet,neumann,nodes2edge)
- New array created, it returns two nodes belonging to a given edge
getEdges(elements)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
Quadrature of the absolute value of a function
trapzAbs(x, y, varargin)
- trapzAbs Trapezoidal numerical integration of absolute value of given
trapzAbs2D(coordinates, fcn, elements, varargin)
- trapzAbs2D Trapezoidal numerical integration of absolute value of
trapzPlus2D(coordinates, fcn, elements, varargin)
- trapzPlus2D Trapezoidal numerical integration of parts of function of
trapzMinus(x, y, varargin)
- trapzMinus Trapezoidal numerical integration of parts of function,
deleterepeatedrows(matrix)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
trapzMinus2D(coordinates, fcn, elements, varargin)
- trapzMinus2D Trapezoidal numerical integration of parts of function of
comparison_1D.m
- This script was created to demonstrate some ways of possible use of
show_elements(elements,coordinates,which)
- function: show_elements(elements,coordinates,whichelements)
example_1D.m
- This is a simple example of using the 'trapzAbs' function. The inputs are
trapzPlusMinus2D(coordinates, fcn, elements, varargin)
- trapzAbs2D Trapezoidal numerical integration of absolute value of
get_edges(elems2nodes)
- function: [element2edges, edge2nodes]=get_edges(elems2nodes)
trapzPlus(x, y, varargin)
- trapzPlus Trapezoidal numerical integration of parts of function,
trapz2D_P1(coordinates,fcn,elements)
- values of P1 functions at midpoints of elements
comparison_2D.m
- This script was created to demonstrate a new function trapzAbs2D.
example_2D.m
- This is a simple example of using the 'trapzAbs2D' function. The inputs
refinement_uniform_2D(nodes2coord,elems2nodes,dirichlet,element_parameters)
- function: [nodes2coord,elems2nodes,dirichlet]=refinement_uniform_2D(nodes2coord,elems2nodes,dirichlet)
An unconstrained minimization of FEM discretized energies.
export_fig
- Exports figures in a publication-quality format
crop_borders
- Crop the borders of an image or stack of images
phider
- Returns the gradients of the basis functions
print2array
- Exports a figure to an image array
Software for a new article - more details later
Detailed implementation of Bogner-Fox-Schmit elements.
shapeder
- Returns the gradients of the shape functions with
print2array
- Exports a figure to an image array
shapefun
- Returns the values of the shape functions.
crop_borders
- Crop the borders of an image or stack of images
export_fig
- Exports figures in a publication-quality format
Finite element solution of a contact problem with a Tresca friction by the semismooth* Newton method.
user_string
- Get/set a user specific string
crop_borders
- Crop the borders of an image or stack of images
pdftops
- Calls a local pdftops executable with the input command
hyperlink
- create a string that is displayable as hyperlink in Matlab console
print2array
- Exports a figure to a bitmap RGB image array
ghostscript
- Calls a local GhostScript executable with the input command
A quasistatic evolution of viscoelastic von Kármán plates.
print2array
- Exports a figure to an image array
export_fig
- Exports figures in a publication-quality format
ghostscript
- Calls a local GhostScript executable with the input command
shapeder
- Returns the gradients of the shape functions with
print2eps
- Prints figures to eps with improved line styles
crop_borders
- Crop the borders of an image or stack of images
eps2pdf
- Convert an eps file to pdf format using ghostscript
screen2jpeg
- Generate a JPEG file of the current figure with
user_string
- Get/set a user specific string
shapefun
- Returns the values of the shape functions.
Implementation of 3D Crouzeix-Raviart mortar finite element
tetrarefine3(XYZ,Elm,ERef)
- author: Ales Janka, ales.janka@unifr.ch, http://perso.unifr.ch/ales.janka
entryInWhichRows(A)
- function: entryrows=entryInWhichRows(A)
deleterepeatedrows(matrix)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
exact_solution(point)
- onesvector=ones(size(x));
rotate_intersection_elements(elements4,coordinates)
- elements4_OnIntersection=elements4(T1_Plane_index,:);
stima_massCR_3D(elements3,coordinates,element2faces)
- Crouzeix-Raviart mass matrix computation
getEdges(elements)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
face_middlepoint(X,F)
- middle_points=zeros(size(F,2),3);
getFaces(elements)
- function: [element2faces, face2nodes]=getFaces(elements)
f(point)
- VolumeForce =ones(size(VolumeForce));
rhs_CR_3D(elements3, coordinates, element2faces, faces)
- GaussPoints2=GaussPoints1;
Calculation of two-phase obstacle problem and its a posteriori error estimate
errorExact(x,u,v,u_derivative)
- computed square H1 seminorm
compute_exact_error.m
- exact error computation
start_demo_Uzawa.m
- SETUP START
u_exact(x)
- middle_branch=find(abs(x)<=r_contact);
start.m
- global x lambda_plus lambda_minus u_left_value u_right_value
energy_norm(v)
- assume v is a piecewise linear function
start_demo_majorant.m
- SETUP START
Algorithm generates edges of a conforming 2D/3D triangulation.
deleterepeatedrows(matrix)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
getEdges(elements)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
entryInWhichRows(A)
- function: entryrows=entryInWhichRows(A)
Calculation of two-phase obstacle problem in 2D and its a posteriori error estimate
refinement_uniform_3D(nodes2coord,elems2nodes,dirichlet)
- function: [nodes2coord,elems2nodes,dirichlet]=refinement_uniform_3D(nodes2coord,elems2nodes,dirichlet)
change_jet.m
- This script help produce a new 'jet'-like colormap based on other RGB reference colors
visualize_solution_mu.m
- colorbar
get_boundary_edges(edges2nodes,bedges2nodes)
- GETBOUNDARYEDGES
show_edges(edge2nodes,coordinates,which)
- function: show_edges(edge,coordinates,whichelements)
basis_RT0(p)
- BASIS_RT0
signs_faces(nodes2coord,elems2faces,faces2nodes,B_K)
- SIGNS_FACES
trapzPlus2D(nodes, elements, varargin)
- trapzPlus2D Trapezoidal numerical integration of parts of function of
shapefun
- SHAPEFUN Returns the values of the shape functions.
interpolate_cbrewer(cbrew_init, interp_method, ncolors)
- % INTERPOLATE_CBREWER - interpolate a colorbrewer map to ncolors levels
trapzMinus2D(nodes, elements, varargin)
- trapzMinus2D Trapezoidal numerical integration of parts of function of
tricont(varargin)
- TRICONT Contours data on a triangular mesh
solution_RT0( elems, B_K, B_K_det, signs, u_h, qr )
- SOLUTION_RT0
intquad( qr, dim )
- INTQUAD
start_majorant.m
- majorant minimization
mass_matrix_Nedelec0( elems2edges, B_K, B_K_det, signs, material )
- MASS_MATRIX_NEDELEC0
shapeder
- SHAPEDER Returns the gradients of the shape functions with
stiffness_matrix_RT0( elems, B_K_det, signs, material )
- STIFFNESS_MATRIX_RT0
energy(v,alpha_plus,alpha_minus,exact_quadrature)
- derivative
mass_matrix_P1( elems2nodes, B_K_det, material )
- MASS_MATRIX_P1
stiffness_matrix_Nedelec0( elems2edges, B_K, B_K_det, signs, material )
- STIFFNESS_MATRIX_NEDELEC0
get_boundary_faces(elems2faces,faces2nodes)
- GETBOUNDARYFACES
basis_P1(p)
- BASIS_P1
print_table.m
- for level=levels_energy_error:levels_level_exact_energy
tridemo.m
- TRIDEMO Demonstration of triangular
numInt_1D
- numInt_1D Trapezoidal numerical integration of parts of function,
numInt_2D
- numInt_2D Trapezoidal numerical integration of parts of function of
Vectorized assembly of FEM matrices and including complete calculations of elastoplasticity benchmarks.
kpde2dmss
- Assembles the mass matrix with P1 finite element
kpde3dstf
- Assembles the stiffness matrix with P1 finite element
kpde2dgphi
- Elements area and gradient of basis functions
kelas3dvmes
- Computes 3D Von Mises effective stress
kelas2dshow
- Visualisation du maillage déformé et/ou des contraintes
stima
- Computes element stiffness matrix for tetraeder.
kelas2dvmes
- Computes 2D Von Mises effective stress
kpde2dstf
- Assembles the stiffness matrix with P1 finite element
Simulation of a mesoscopic model of temperature/external field-driven evolution of magnetization
show_mesh(coordinates,elements)
- colorbar;
entryInWhichRows(A)
- function: entryrows=entryInWhichRows(A)
phider
- PHIDER Returns the gradients of the basis functions
show_elements(elements,coordinates,which)
- function: show_elements(elements,coordinates,whichelements)
prolongate(vNodal,elements,edge2nodes)
- prolongates bilinear function defined on rectangles
shapefun
- SHAPEFUN Returns the values of the shape functions.
refinement_uniform_2D(coordinates,elements,dirichlet,neumann,element_parameters)
- function: [coordinates,elements3,dirichlet]=refinement_uniform(coordinates,elements3,dirichlet)
stifness_matrixP1_2D(elements,coordinates)
- particular part for a given element in a given dimension
getEdges_triangles(elements)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
show_deformed_mesh(U)
- colorbar;
screen2jpeg(filename)
- SCREEN2JPEG Generate a JPEG file of the current figure with
refinement_uniform(coordinates,elements3,dirichlet,neumann,element_parameters)
- function: [coordinates,elements3,dirichlet]=refinement_uniform(coordinates,elements3,dirichlet)
mesh_preparation
- Preparation - extracting Dirichlet nodes
magnet_mesh_old(number_of_refinements)
- Uniforme Verfeinerung
CMRmap(M)
- Create a monochrome-compatible colour map
shapeder
- SHAPEDER Returns the gradients of the shape functions with
Construction of the Lagrange interpolation polynomial from its definition.
lagrangePoly(xp,yp,xx)
- Example (the most important):
Time-dependent 2D elastoplasticity solver for a two-yied model.
mesh_preparation
- Preparation - extracting Dirichlet nodes
evaluate_Phi(U,P1prev,P2prev)
- global yield_type
platehole_mesh(number_of_refinements)
- current directory stored
evaluate_P_on_element(Rlocal,Ulocal,P1prevlocal,P2prevlocal, C, mu, sigmay1, sigmay2, h1, h2, tolerance)
- gamma (epsilon11,epsilon22,epsilon12) on the element
ring_mesh(number_of_refinements)
- A square beam fixed on the left edge in both directions and pulled on the right edge (Neumann conditions)
Lshape_mesh(number_of_refinements)
- Uniforme Verfeinerung
cook_mesh(number_of_refinements)
- current directory stored
beam_2D_symetric_mesh(number_of_refinements)
- A square beam fixed on the left edge in both directions and pulled on the right edge (Neumann conditions)
evaluate_fv(time)
- Calculating fv
refinement_uniform(coordinates,elements3,dirichlet,neumann)
- function: [coordinates,elements3,dirichlet]=refinement_uniform(coordinates,elements3,dirichlet)
evaluate_gv(t)
- Calculating gv
We provide an alternative to the function rotm2eul from Robotics System Toolbox and comparisons.
print2array
- Exports a figure to an image array
export_fig
- Exports figures in a publication-quality format
crop_borders
- Crop the borders of an image or stack of images
Given a figure with (eg. tight) axes, this simple code automatically recalculates and enlarges axes them by a prescribed factor.
Implementation of the semismooth* method.
build_AB
- Construct block matrices A and B from vector a and matrix GRAD
Calculaton of Error identities for variational problems with obstacles
print_figure(filename)
- Saving eps with matlab and then producing pdf and png with system commands
prepare_new_figure(subplotsx,subplotsy)
- parameters for figure and panel size
solution_parabolic_obstacle_1D(f,phi,x)
- phi is scalar here, not the vector of constraints!
solution_classical_obstacle_1D(f,phi,x)
- phi is scalar here, not the vector of constraints!
solution_double_obstacle_1D(x,alphaPlus,alphaMinus)
- An analytical solution is generated for alphaPlus>=2 and alphaMinus>=2 only.
start.m
- gap_primal=-f.*(coincidence_u.*(v-phi)); error_primal_nonlinear=trapz(x,gap_primal);
Evaluation of the HADAMARD’S INEQUALITY IN THE MEAN in 2 space-dimensions.
crop_borders
- Crop the borders of an image or stack of images
print2array
- Exports a figure to an image array
lbmap
- Returns specified Light-Bertlein colormap.
An obstacle problem is solved and the error of the solution is estimated
screen2jpeg(filename)
- SCREEN2JPEG Generate a JPEG file of the current figure with
stifness_matrixP1_2D(elements,coordinates)
- particular part for a given element in a given dimension
shapefun
- SHAPEFUN Returns the values of the shape functions.
compute_majorant(vRectangular, phi_nodal, coordinates, elementsRectangular,
- iteration_display=1;
refinement_uniform_2D(coordinates,elements,dirichlet,neumann,element_parameters)
- function: [coordinates,elements3,dirichlet]=refinement_uniform(coordinates,elements3,dirichlet)
u_exact_ring_spherical_obstacle(xy,f,phi,Lx,Ly,rho)
- figure(1000)
getEdges_rectangles(elements)
- function: [element2edges, edge2nodes]=edge_numbering(elements)
visualize_solution.m
- discrete rectangular solution
phider
- PHIDER Returns the gradients of the basis functions
shapeder
- SHAPEDER Returns the gradients of the shape functions with
start.m
- clear all
prolongate(vNodal,elements,edge2nodes)
- prolongates bilinear function defined on rectangles
mesh_setup.m
- mesh generation
conv_ma2av(ma)
- converts matrix to array of vectors
또한 다음 목록에서 웹사이트를 선택하실 수도 있습니다.
미주
- América Latina (Español)
- Canada (English)
- United States (English)
유럽
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom(English)
아시아 태평양
- Australia (English)
- India (English)
- New Zealand (English)
- 中国
- 日本Japanese (日本語)
- 한국Korean (한국어)